Analysis of Germline Variants in CDH1, IGFBP3, MMP1, MMP3, STK15 and VEGF in Familial and Sporadic Renal Cell Carcinoma

BACKGROUND:The investigation of rare familial forms of kidney cancer has provided important insights into the biology of sporadic renal cell carcinoma (RCC). In particular, the identification of the von Hippel Lindau (VHL) familial cancer syndrome gene (VHL) provided the basis for the discovery that VHL is somatically inactivated in most sporadic clear cell RCC. Many cases of familial RCC do not have mutations in known RCC susceptibility genes and there is evidence that genetic modifiers may influence the risk of RCC in VHL disease patients. Hence we hypothesised that low-penetrance functional genetic variants in pathways related to the VHL protein (pVHL) function might (a) modify the phenotypic expression of VHL disease and/or (b) predispose to sporadic RCC. METHODOLOGY/PRINCIPAL FINDINGS:We tested this hypothesis for functional polymorphisms in CDH1 (rs16260), IGFBP3 (rs2854744), MMP1 (rs1799750), MMP3 (rs679620), STK15 (rs2273535) and VEGF (rs1570360). We observed that variants of MMP1 and MMP3 were significant modifiers of RCC risk (and risks of retinal angioma and cerebellar haemangioblastoma) in VHL disease patients. In addition, higher frequencies of the MMP1 rs1799750 2G allele (p = 0.017, OR 1.49, 95%CI 1.06-2.08) and the MMP1/MMP3 rs1799750/rs679620 2G/G haplotype (OR 1.45, 95%CI 1.01-2.10) were detected in sporadic RCC patients than in controls (n = 295). CONCLUSIONS/SIGNIFICANCE:These findings (a) represent the first example of genetic modifiers of RCC risk in VHL disease, (b) replicate a previous report of an association between MMP1/MMP3 variants and sporadic RCC and (c) further implicate MMP1/MMP3-related pathways in the pathogenesis of familial and sporadic RCC

Using Classical Population Genetics Tools with Heterochroneous Data: Time Matters!

BACKGROUND:New polymorphism datasets from heterochroneous data have arisen thanks to recent advances in experimental and microbial molecular evolution, and the sequencing of ancient DNA (aDNA). However, classical tools for population genetics analyses do not take into account heterochrony between subsets, despite potential bias on neutrality and population structure tests. Here, we characterize the extent of such possible biases using serial coalescent simulations. METHODOLOGY/PRINCIPAL FINDINGS:We first use a coalescent framework to generate datasets assuming no or different levels of heterochrony and contrast most classical population genetic statistics. We show that even weak levels of heterochrony ( approximately 10% of the average depth of a standard population tree) affect the distribution of polymorphism substantially, leading to overestimate the level of polymorphism theta, to star like trees, with an excess of rare mutations and a deficit of linkage disequilibrium, which are the hallmark of e.g. population expansion (possibly after a drastic bottleneck). Substantial departures of the tests are detected in the opposite direction for more heterochroneous and equilibrated datasets, with balanced trees mimicking in particular population contraction, balancing selection, and population differentiation. We therefore introduce simple corrections to classical estimators of polymorphism and of the genetic distance between populations, in order to remove heterochrony-driven bias. Finally, we show that these effects do occur on real aDNA datasets, taking advantage of the currently available sequence data for Cave Bears (Ursus spelaeus), for which large mtDNA haplotypes have been reported over a substantial time period (22-130 thousand years ago (KYA)). CONCLUSIONS/SIGNIFICANCE:Considering serial sampling changed the conclusion of several tests, indicating that neglecting heterochrony could provide significant support for false past history of populations and inappropriate conservation decisions. We therefore argue for systematically considering heterochroneous models when analyzing heterochroneous samples covering a large time scale

Reduced Plasmodium vivax Erythrocyte Infection in PNG Duffy-Negative Heterozygotes

BACKGROUND: Erythrocyte Duffy blood group negativity reaches fixation in African populations where Plasmodium vivax (Pv) is uncommon. While it is known that Duffy-negative individuals are highly resistant to Pv erythrocyte infection, little is known regarding Pv susceptibility among heterozygous carriers of a Duffy-negative allele (+/−). Our limited knowledge of the selective advantages or disadvantages associated with this genotype constrains our understanding of the effect that interventions against Pv may have on the health of people living in malaria-endemic regions. METHODS AND FINDINGS: We conducted cross-sectional malaria prevalence surveys in Papua New Guinea (PNG), where we have previously identified a new Duffy-negative allele among individuals living in a region endemic for all four human malaria parasite species. We evaluated infection status by conventional blood smear light microscopy and semi-quantitative PCR-based strategies. Analysis of a longitudinal cohort constructed from our surveys showed that Duffy heterozygous (+/−) individuals were protected from Pv erythrocyte infection compared to those homozygous for wild-type alleles (+/+) (log-rank tests: LM, p = 0.049; PCR, p = 0.065). Evaluation of Pv parasitemia, determined by semi-quantitative PCR-based methods, was significantly lower in Duffy +/− vs. +/+ individuals (Mann-Whitney U: p = 0.023). Overall, we observed no association between susceptibility to P. falciparum erythrocyte infection and Duffy genotype. CONCLUSIONS: Our findings provide the first evidence that Duffy-negative heterozygosity reduces erythrocyte susceptibility to Pv infection. As this reduction was not associated with greater susceptibility to Pf malaria, our in vivo observations provide evidence that Pv-targeted control measures can be developed safely

Novel Association Strategy with Copy Number Variation for Identifying New Risk Loci of Human Diseases

Copy number variations (CNV) are important causal genetic variations for human disease; however, the lack of a statistical model has impeded the systematic testing of CNVs associated with disease in large-scale cohort.Here, we developed a novel integrated strategy to test CNV-association in genome-wide case-control studies. We converted the single-nucleotide polymorphism (SNP) signal to copy number states using a well-trained hidden Markov model. We mapped the susceptible CNV-loci through SNP site-specific testing to cope with the physiological complexity of CNVs. We also ensured the credibility of the associated CNVs through further window-based CNV-pattern clustering. Genome-wide data with seven diseases were used to test our strategy and, in total, we identified 36 new susceptible loci that are associated with CNVs for the seven diseases: 5 with bipolar disorder, 4 with coronary artery disease, 1 with Crohn's disease, 7 with hypertension, 9 with rheumatoid arthritis, 7 with type 1 diabetes and 3 with type 2 diabetes. Fifteen of these identified loci were validated through genotype-association and physiological function from previous studies, which provide further confidence for our results. Notably, the genes associated with bipolar disorder converged in the phosphoinositide/calcium signaling, a well-known affected pathway in bipolar disorder, which further supports that CNVs have impact on bipolar disorder.Our results demonstrated the effectiveness and robustness of our CNV-association analysis and provided an alternative avenue for discovering new associated loci of human diseases

Linkage Group Selection: Towards Identifying Genes Controlling Strain Specific Protective Immunity in Malaria

Protective immunity against blood infections of malaria is partly specific to the genotype, or strain, of the parasites. The target antigens of Strain Specific Protective Immunity are expected, therefore, to be antigenically and genetically distinct in different lines of parasite. Here we describe the use of a genetic approach, Linkage Group Selection, to locate the target(s) of Strain Specific Protective Immunity in the rodent malaria parasite Plasmodium chabaudi chabaudi. In a previous such analysis using the progeny of a genetic cross between P. c. chabaudi lines AS-pyr1 and CB, a location on P. c. chabaudi chromosome 8 containing the gene for merozoite surface protein-1, a known candidate antigen for Strain Specific Protective Immunity, was strongly selected. P. c. chabaudi apical membrane antigen-1, another candidate for Strain Specific Protective Immunity, could not have been evaluated in this cross as AS-pyr1 and CB are identical within the cell surface domain of this protein. Here we use Linkage Group Selection analysis of Strain Specific Protective Immunity in a cross between P. c. chabaudi lines CB-pyr10 and AJ, in which merozoite surface protein-1 and apical membrane antigen-1 are both genetically distinct. In this analysis strain specific immune selection acted strongly on the region of P. c. chabaudi chromosome 8 encoding merozoite surface protein-1 and, less strongly, on the P. c. chabaudi chromosome 9 region encoding apical membrane antigen-1. The evidence from these two independent studies indicates that Strain Specific Protective Immunity in P. c. chabaudi in mice is mainly determined by a narrow region of the P. c. chabaudi genome containing the gene for the P. c. chabaudi merozoite surface protein-1 protein. Other regions, including that containing the gene for P. c. chabaudi apical membrane antigen-1, may be more weakly associated with Strain Specific Protective Immunity in these parasites

Genome-Wide Population-Based Association Study of Extremely Overweight Young Adults – The GOYA Study

Background: Thirty-two common variants associated with body mass index (BMI) have been identified in genome-wide association studies, explaining ~1.45% of BMI variation in general population cohorts. We performed a genome-wide association study in a sample of young adults enriched for extremely overweight individuals. We aimed to identify new loci associated with BMI and to ascertain whether using an extreme sampling design would identify the variants known to be associated with BMI in general populations. Methodology/Principal Findings: From two large Danish cohorts we selected all extremely overweight young men and women (n = 2,633), and equal numbers of population-based controls (n = 2,740, drawn randomly from the same populations as the extremes, representing ~212,000 individuals). We followed up novel (at the time of the study) association signals (

Analysis of Mitochondrial DNA Sequences in Childhood Encephalomyopathies Reveals New Disease-Associated Variants

BACKGROUND: Mitochondrial encephalomyopathies are a heterogeneous group of clinical disorders generally caused due to mutations in either mitochondrial DNA (mtDNA) or nuclear genes encoding oxidative phosphorylation (OXPHOS). We analyzed the mtDNA sequences from a group of 23 pediatric patients with clinical and morphological features of mitochondrial encephalopathies and tried to establish a relationship of identified variants with the disease. METHODOLOGY/PRINCIPLE FINDINGS: Complete mitochondrial genomes were amplified by PCR and sequenced by automated DNA sequencing. Sequencing data was analyzed by SeqScape software and also confirmed by BLASTn program. Nucleotide sequences were compared with the revised Cambridge reference sequence (CRS) and sequences present in mitochondrial databases. The data obtained shows that a number of known and novel mtDNA variants were associated with the disease. Most of the non-synonymous variants were heteroplasmic (A4136G, A9194G and T11916A) suggesting their possibility of being pathogenic in nature. Some of the missense variants although homoplasmic were showing changes in highly conserved amino acids (T3394C, T3866C, and G9804A) and were previously identified with diseased conditions. Similarly, two other variants found in tRNA genes (G5783A and C8309T) could alter the secondary structure of Cys-tRNA and Lys-tRNA. Most of the variants occurred in single cases; however, a few occurred in more than one case (e.g. G5783A and A10149T). CONCLUSIONS AND SIGNIFICANCE: The mtDNA variants identified in this study could be the possible cause of mitochondrial encephalomyopathies with childhood onset in the patient group. Our study further strengthens the pathogenic score of known variants previously reported as provisionally pathogenic in mitochondrial diseases. The novel variants found in the present study can be potential candidates for further investigations to establish the relationship between their incidence and role in expressing the disease phenotype. This study will be useful in genetic diagnosis and counseling of mitochondrial diseases in India as well as worldwide

Similarity in Recombination Rate Estimates Highly Correlates with Genetic Differentiation in Humans

Recombination varies greatly among species, as illustrated by the poor conservation of the recombination landscape between humans and chimpanzees. Thus, shorter evolutionary time frames are needed to understand the evolution of recombination. Here, we analyze its recent evolution in humans. We calculated the recombination rates between adjacent pairs of 636,933 common single-nucleotide polymorphism loci in 28 worldwide human populations and analyzed them in relation to genetic distances between populations. We found a strong and highly significant correlation between similarity in the recombination rates corrected for effective population size and genetic differentiation between populations. This correlation is observed at the genome-wide level, but also for each chromosome and when genetic distances and recombination similarities are calculated independently from different parts of the genome. Moreover, and more relevant, this relationship is robustly maintained when considering presence/absence of recombination hotspots. Simulations show that this correlation cannot be explained by biases in the inference of recombination rates caused by haplotype sharing among similar populations. This result indicates a rapid pace of evolution of recombination, within the time span of differentiation of modern humans

De Novo Transcriptome Sequencing in Anopheles funestus Using Illumina RNA-Seq Technology

BACKGROUND: Anopheles funestus is one of the primary vectors of human malaria, which causes a million deaths each year in sub-Saharan Africa. Few scientific resources are available to facilitate studies of this mosquito species and relatively little is known about its basic biology and evolution, making development and implementation of novel disease control efforts more difficult. The An. funestus genome has not been sequenced, so in order to facilitate genome-scale experimental biology, we have sequenced the adult female transcriptome of An. funestus from a newly founded colony in Burkina Faso, West Africa, using the Illumina GAIIx next generation sequencing platform. METHODOLOGY/PRINCIPAL FINDINGS: We assembled short Illumina reads de novo using a novel approach involving iterative de novo assemblies and "target-based" contig clustering. We then selected a conservative set of 15,527 contigs through comparisons to four Dipteran transcriptomes as well as multiple functional and conserved protein domain databases. Comparison to the Anopheles gambiae immune system identified 339 contigs as putative immune genes, thus identifying a large portion of the immune system that can form the basis for subsequent studies of this important malaria vector. We identified 5,434 1:1 orthologues between An. funestus and An. gambiae and found that among these 1:1 orthologues, the protein sequence of those with putative immune function were significantly more diverged than the transcriptome as a whole. Short read alignments to the contig set revealed almost 367,000 genetic polymorphisms segregating in the An. funestus colony and demonstrated the utility of the assembled transcriptome for use in RNA-seq based measurements of gene expression. CONCLUSIONS/SIGNIFICANCE: We developed a pipeline that makes de novo transcriptome sequencing possible in virtually any organism at a very reasonable cost ($6,300 in sequencing costs in our case). We anticipate that our approach could be used to develop genomic resources in a diversity of systems for which full genome sequence is currently unavailable. Our An. funestus contig set and analytical results provide a valuable resource for future studies in this non-model, but epidemiologically critical, vector insect

Austro-Asiatic Tribes of Northeast India Provide Hitherto Missing Genetic Link between South and Southeast Asia

Northeast India, the only region which currently forms a land bridge between the Indian subcontinent and Southeast Asia, has been proposed as an important corridor for the initial peopling of East Asia. Given that the Austro-Asiatic linguistic family is considered to be the oldest and spoken by certain tribes in India, Northeast India and entire Southeast Asia, we expect that populations of this family from Northeast India should provide the signatures of genetic link between Indian and Southeast Asian populations. In order to test this hypothesis, we analyzed mtDNA and Y-Chromosome SNP and STR data of the eight groups of the Austro-Asiatic Khasi from Northeast India and the neighboring Garo and compared with that of other relevant Asian populations. The results suggest that the Austro-Asiatic Khasi tribes of Northeast India represent a genetic continuity between the populations of South and Southeast Asia, thereby advocating that northeast India could have been a major corridor for the movement of populations from India to East/Southeast Asia