Neanderthal admixture in current human populations

In the present body of work two primary subjects have been addressed, both individually and in their correspondence, namely (1) the potential for Neanderthals to have contributed to the Modern Human population, and (2) the genetic diversity of one of the most prehistorically impactful human popuations, the Armenians. The first subject is addressed by assessing 1000 mutations in 384 current humans, particularly for those mutations which appear to derive from the Neanderthal lineage. Additionally, the validity of the Neanderthal sequences themselves is evaluated through alignment analysis of fragementary DNA derived from the Vindija Cave sample. Armenian genetic diversity is analyzed through the autosomal short tandem repeats, y-chromsome single nucleotide polymorphisms, and y-chromosome short tandem repeats. The diversity found indicates that Armenians are a diverse group which has been genetically influenced by the various migrations and invasions which have entered their historic lands. Further, we find evidence that Armenians may be closely associated with the peopling of Europe.

Overweight/Obesity and HIV Disease Progression in HIV+ Adults in Botswana

Studies indicate that overweight and obesity protect against HIV-disease progression in antiretroviral therapy (ART)-naïve patients. We examined retrospectively the relationship of overweight/obesity with HIV-disease progression in ART-naïve HIV+ adults in Botswana in a case-control study with 18-month follow-up, which included 217 participants, 139 with BMI 18.0-24.9 kg/m2 and 78 with BMI ≥25 kg/m2. Archived plasma samples were used to determine inflammatory markers: leptin and bacterial endotoxin lipopolysaccharide (LPS), and genotype single nucleotide polymorphisms (SNPs) of the Fat Mass and Obesity Associated Gene (FTO). At baseline, BMI was inversely associated with risk for AIDS-defining conditions (HR=0.218; 95%CI=0.068, 0.701, P=0.011), and higher fat mass was associated with reduced risk of the combined outcome of CD4+cell count ≤250/µL and AIDS-defining conditions, whichever occurred earlier (HR=0.918; 95%CI=0.847, 0.994, P=0.036) over 18 months, adjusting for age, gender, marriage, children, and baseline CD4+cell count and HIV-viral load. FTO-SNP rs17817449 was associated with BMI (OR=1.082; 95%CI=1.001, 1.169; P=0.047). Fat mass was associated with the risk alleles of rs1121980 (OR=1.065; 95%CI=1.009, 1.125, P=0.021), rs8050136 (OR=1.078; 95%CI=1.021, 1.140; P=0.007), and rs17817449 (OR=1.086; 95%CI=1.031, 1.145; P=0.002), controlling for age, gender, tribe, total energy intake, and activity. There were no associations of SNPs with markers of disease progression. Leptin levels were positively associated with BMI (β=1.764; 95%CI=0.788, 2.739; P=0.022) and fat mass (β=0.112; 95%CI=0.090, 0.135; P<0.001), but inversely with viral load (β=-0.305; 95%CI=-0.579, -.031; P=0.030). LPS levels were inversely associated with BMI (OR=0.790, 95%CI=0.630, 0.990; P=0.041), and fat mass (OR=0.852, 95%CI=0.757, 0.958; P=0.007) and directly with viral load (OR=2.608, 95%CI=1.111, 6.124; P=0.028), adjusting for age, gender, smoking and %fat mass. In this cohort, overweight/obesity predicted slower HIV-disease progression. Obesity may confer an advantage in maintaining fat stores to support the overactive immune system. FTO-SNPs may contribute to the variation in fat mass; however, they were not associated with HIV-disease progression. Our findings suggest that the obesity paradox may be explained by the association of increased LPS with lower BMI and higher viral load; while viral load decreased with increasing leptin levels. Studies in African populations are needed to clarify whether genetic variation and inflammation mediate the obesity paradox in HIV-disease progression.

Avaliação de polimorfismos funcionais nos genes de reparo XRCC1, APEX1, XPD e XPF em carcinomas de células escamosas orais

Base excision repair (BER) and nucleotide excision repair (NER) pathways play critical role in maintaining genome integrity. Polymorphisms in BER and NER genes which modulate the DNA repair capacity may affect the susceptibility and prognosis of oral cancer. This study was conducted with genomic DNA from 92 patients with oral squamous cell carcinomas (OSCC) and 130 controls. The cases were followed up to explore the associations between BER and NER genes polymorphisms and the risk and prognosis of OSCC. Four single-nucleotide polymorphisms (SNPs) in XRCC1 (rs25487), APEX1 (rs1130409), XPD (rs13181) and XPF (rs1799797) genes were tested by polymerase chain reaction – quantitative real time method. The GraphPad Prism version 6.0.1 statistical software was applied for statistical analysis of association. Odds ratio (OR), hazard ratio (HR), and their 95 % confidence intervals (CIs) were calculated by logistic regression. Kaplan-Meier curve and Cox proportional hazard model were used for prognostic analysis. The presence of polymorphic variants in XRCC1, APEX1, XPD and XPF genes were not associated with an increased risk of OSCC. Gene-environment interactions with smoking were not significant for any polymorphism. The presence of polymorphic variants of the XPD gene in association with alcohol consumption conferred an increased risk of 1.86 (95% CI: 0.86 – 4.01, p=0.03) for OSCC. Only APEX1 was associated with decreased specific survival (HR 3.94, 95% CI: 1.31 – 11.88, p=0.01). These results suggest an interaction between polymorphic variants of the XPF gene and alcohol consumption. Additionally APEX1 may represent a prognostic marker for OSCC.

Valor prognóstico de polimorfismos nos genes de reparo do DNA XRCC3 E RAD51 em pacientes com carcinoma epidermóide oral e de orofaringe

Faults in the genes responsible for repairs to the DNA can influence the onset of cancer or affect the response to treatment. This research evaluated the frequency of three single nucleotide polymorphisms (SNPs) in two repair genes DNA RAD51 172g> T (rs1801321), RAD51 135G> C (rs1801320) and XRCC3 T241M (rs861539) in individuals without cancer (n = 130) and patients with oral squamous cell carcinoma (OSC) and carcinoma oropharyngeal squamous (ORSC) (n = 126) and investigated possible relationships of these findings with clinical and pathological data and clinical outcomes: tumor response to radiotherapy and chemotherapy, disease-free survival, and overall survival. It was found that the allele and genotype frequencies were in equilibrium Hard-Weinberg equilibrium. The presence of at least one polymorphic allele in XRCC3 (rs861539) gene is associated with histological grade (WHO) higher (p = 0.007). We observed a higher recurrence rate trend (p = 0.08) and more advanced stage (p = 0.08) in the group that had at least one polymorphic allele of RAD51 gene (rs1801321). The presence of the analyzed SNPs not proved to be a risk factor for the development of CEO or CEOR; however, when combined with smoking or drinking, increased the risk of developing cancer from three to one hundred and fifty times. The tumor response to radiotherapy and chemotherapy was similar in patients with and without SNPs. No polymorphism showed statistical significance in relation to recurrence-free survival or overall survival. We conclude that the presence of at least one polymorphic allele of the SNPs rs861539 in XRCC3 gene, rs1801320 and rs1801321 in the RAD51 gene increase the risk of development of OSC and ORSC, when associated with the habit of drinking or smoking. Polymorphisms studied in XRCC3 and RAD51 genes are not associated with response to radiation therapy, relapse-free survival or overall survival.

Associations between smoking behavior-related alleles and the risk of melanoma.

Several studies have reported that cigarette smoking is inversely associated with the risk of melanoma. This study further tested whether incorporating genetic factors will provide another level of evaluation of mechanisms underlying the association between smoking and risk of melanoma. We investigated the association between SNPs selected from genome-wide association studies (GWAS) on smoking behaviors and risk of melanoma using 2,298 melanoma cases and 6,654 controls. Among 16 SNPs, three (rs16969968 [A], rs1051730 [A] and rs2036534 [C] in the 15q25.1 region) reached significance for association with melanoma risk in men (0.01 < = P values < = 0.02; 0.85 < = Odds Ratios (ORs) <= 1.20). There was association between the genetic scores based on the number of smoking behavior-risk alleles and melanoma risk with P-trend = 0.005 among HPFS. Further association with smoking behaviors indicating those three SNPs (rs16969968 [A], rs1051730 [A] and rs2036534 [C]) significantly associated with number of cigarettes smoked per day, CPD, with P = 0.009, 0.011 and 0.001 respectively. The SNPs rs215605 in the PDE1C gene and rs6265 in the BDNF gene significantly interacted with smoking status on melanoma risk (interaction P = 0.005 and P = 0.003 respectively). Our study suggests that smoking behavior-related SNPs are likely to play a role in melanoma development and the potential public health importance of polymorphisms in the CHRNA5-A3-B4 gene cluster. Further larger studies are warranted to validate the findings.

Risk Prediction for Epithelial Ovarian Cancer in 11 United States-Based Case-Control Studies: Incorporation of Epidemiologic Risk Factors and 17 Confirmed Genetic Loci.

Previously developed models for predicting absolute risk of invasive epithelial ovarian cancer have included a limited number of risk factors and have had low discriminatory power (area under the receiver operating characteristic curve (AUC) < 0.60). Because of this, we developed and internally validated a relative risk prediction model that incorporates 17 established epidemiologic risk factors and 17 genome-wide significant single nucleotide polymorphisms (SNPs) using data from 11 case-control studies in the United States (5,793 cases; 9,512 controls) from the Ovarian Cancer Association Consortium (data accrued from 1992 to 2010). We developed a hierarchical logistic regression model for predicting case-control status that included imputation of missing data. We randomly divided the data into an 80% training sample and used the remaining 20% for model evaluation. The AUC for the full model was 0.664. A reduced model without SNPs performed similarly (AUC = 0.649). Both models performed better than a baseline model that included age and study site only (AUC = 0.563). The best predictive power was obtained in the full model among women younger than 50 years of age (AUC = 0.714); however, the addition of SNPs increased the AUC the most for women older than 50 years of age (AUC = 0.638 vs. 0.616). Adapting this improved model to estimate absolute risk and evaluating it in prospective data sets is warranted.

Association between the oxytocin receptor (OXTR) gene and mesolimbic responses to rewards.

BACKGROUND: There has been significant progress in identifying genes that confer risk for autism spectrum disorders (ASDs). However, the heterogeneity of symptom presentation in ASDs impedes the detection of ASD risk genes. One approach to understanding genetic influences on ASD symptom expression is to evaluate relations between variants of ASD candidate genes and neural endophenotypes in unaffected samples. Allelic variations in the oxytocin receptor (OXTR) gene confer small but significant risk for ASDs for which the underlying mechanisms may involve associations between variability in oxytocin signaling pathways and neural response to rewards. The purpose of this preliminary study was to investigate the influence of allelic variability in the OXTR gene on neural responses to monetary rewards in healthy adults using functional magnetic resonance imaging (fMRI). METHODS: The moderating effects of three single nucleotide polymorphisms (SNPs) (rs1042778, rs2268493 and rs237887) of the OXTR gene on mesolimbic responses to rewards were evaluated using a monetary incentive delay fMRI task. RESULTS: T homozygotes of the rs2268493 SNP demonstrated relatively decreased activation in mesolimbic reward circuitry (including the nucleus accumbens, amygdala, insula, thalamus and prefrontal cortical regions) during the anticipation of rewards but not during the outcome phase of the task. Allelic variation of the rs1042778 and rs237887 SNPs did not moderate mesolimbic activation during either reward anticipation or outcomes. CONCLUSIONS: This preliminary study suggests that the OXTR SNP rs2268493, which has been previously identified as an ASD risk gene, moderates mesolimbic responses during reward anticipation. Given previous findings of decreased mesolimbic activation during reward anticipation in ASD, the present results suggest that OXTR may confer ASD risk via influences on the neural systems that support reward anticipation.

Peroxisome Proliferator-Activated Receptor-γ Coactivator 1-α (PPARGC1A) Genetic Associations with Type 2 Diabetes in Three Ethnicities

Genetic heterogeneity, lifestyle factors, gene-gene or gene-environment interactions are the determinants of T2D which puts Hispanics and populations with African ancestry at higher risk of developing T2D. In this dissertation, the genetic associations of PPARGC1A polymorphisms with T2D and its related phenotypes (metabolic markers) in Haitian Americans (cases=110, controls=116), African Americans (cases=120, controls=124) and Cuban Americans (cases=160, controls=181) of South Florida were explored. Five single nucleotide polymorphisms of gene PPARGC1A were evaluated in each ethnicity for their disease association. In Haitian Americans, rs7656250 (OR= 0.22, pp=0.03) had significant protective association with T2D but had risk association in African Americans for rs7656250 (OR=1.02, p=0.96) and rs4235308 (OR=2.53, p=0.03). We found that in Haitian American females, both rs7656250 (OR=0.23, pp=0.03) had protective association with T2D. In African American females, rs7656250 (OR=1.14, p=0.78) had risk association whereas in males, it had significant protective effect (OR=0.37, p=0.04). However, the risk association exhibited by rs4235308 was stronger in African American females (OR=2.69, p=0.03) than males (OR=1.16, p=0.72). In Cuban Americans, only rs7656250 showed significant risk association with T2D (OR=6.87, p=0.02) which was stronger in females alone (OR=7.67, p=0.01). We also observed significant differences among correlations of PPARGC1A SNPs and T2D phenotypes. Positive correlation was observed for log Hs-CRP with rs3774907 (pp=0.03) in Cuban Americans respectively. Correlation of log A1C with rs7656250 (p=0.02) was positive in Cuban Americans while it was negative for rs3774907 in Haitian Americans (ppPPARGC1A correlations with T2D and its phenotypes among the three ethnicities studied (ii) the associations of PPARGC1A SNPs showed significant effect modification by sex. The findings suggest that variations in effects of PPARGC1A gene polymorphisms among three ethnicities and between sexes may have biomedical implications for the development of T2D as well as the phenotypes related to T2D.

The differential impact of oxytocin receptor (OXTR) genotypes on the risk of autism spectrum disorder (ASD) and resulting social communication deficits

Background: Autism spectrum disorder (ASD) is multifactorial and is likely the result of complex interactions between multiple environmental and genetic factors. Recently, it has been suggested that each symptom cluster of the disorder, such as poor social communication, may be mediated by different genetic influences. Genes in the oxytocin pathway, which mediates social behaviours in humans, have been studied with single nucleotide polymorphisms (SNPs) in the oxytocin receptor gene (OXTR) being implicated in ASD. This thesis examines the presence of different oxytocin receptor genotypes, and their associations with ASD and resulting social communication deficits. Methods: The relationship between four OXTR variants and ASD was evaluated in 607 ASD simplex (SPX) families. Cases were compared to their unaffected siblings using a conditional logistic approach. Odds ratios and associated 95 percent confidence intervals were obtained. A second sample of 235 individuals with a diagnosis of ASD was examined to evaluate whether these four OXTR variants were associated with social communication scores on the Autism Diagnostic Interview – Revised (ADI-R). Parameter estimates and associated 95 percent confidence intervals were generated using a linear regression approach. Multiple testing issues were addressed using false discovery adjustments. Results: The rs53576 AG genotype was significantly associated with a lower risk of ASD (OR = 0.707, 95% CI: 0.512-0.975). A single genotype (AG) provided by the rs2254298 marker was found to be significantly associated with higher social communication scores (Parameter estimate = 1.833, SE = 0.762, p = 0.0171). This association was also seen in a Caucasian only and mothers as the respondent samples. No association was significant following false discovery rate adjustments. Conclusion: The findings from these studies provide limited support for the role of OXTR SNPs in ASD, especially in social communication skills. The clinical significance of these associations remains unknown, however, it is likely that these associations do not play a role in the severity of symptoms associated with ASD. Rather, they may be important in the appearance of social deficits due to the rs2254298 markers association with enlarged amygdalas.

Evidence for three genetic loci involved in both anorexia nervosa risk and variation of body mass index

The maintenance of normal body weight is disrupted in patients with anorexia nervosa (AN) for prolonged periods of time. Prior to the onset of AN, premorbid body mass index (BMI) spans the entire range from underweight to obese. After recovery, patients have reduced rates of overweight and obesity. As such, loci involved in body weight regulation may also be relevant for AN and vice versa. Our primary analysis comprised a cross-trait analysis of the 1000 single-nucleotide polymorphisms (SNPs) with the lowest P-values in a genome-wide association meta-analysis (GWAMA) of AN (GCAN) for evidence of association in the largest published GWAMA for BMI (GIANT). Subsequently we performed sex-stratified analyses for these 1000 SNPs. Functional ex vivo studies on four genes ensued. Lastly, a look-up of GWAMA-derived BMI-related loci was performed in the AN GWAMA. We detected significant associations (P-values <5 × 10(-5), Bonferroni-corrected P<0.05) for nine SNP alleles at three independent loci. Interestingly, all AN susceptibility alleles were consistently associated with increased BMI. None of the genes (chr. 10: CTBP2, chr. 19: CCNE1, chr. 2: CARF and NBEAL1; the latter is a region with high linkage disequilibrium) nearest to these SNPs has previously been associated with AN or obesity. Sex-stratified analyses revealed that the strongest BMI signal originated predominantly from females (chr. 10 rs1561589; Poverall: 2.47 × 10(-06)/Pfemales: 3.45 × 10(-07)/Pmales: 0.043). Functional ex vivo studies in mice revealed reduced hypothalamic expression of Ctbp2 and Nbeal1 after fasting. Hypothalamic expression of Ctbp2 was increased in diet-induced obese (DIO) mice as compared with age-matched lean controls. We observed no evidence for associations for the look-up of BMI-related loci in the AN GWAMA. A cross-trait analysis of AN and BMI loci revealed variants at three chromosomal loci with potential joint impact. The chromosome 10 locus is particularly promising given that the association with obesity was primarily driven by females. In addition, the detected altered hypothalamic expression patterns of Ctbp2 and Nbeal1 as a result of fasting and DIO implicate these genes in weight regulation.

Increasing the power of genome wide association studies in natural populations using repeated measures : evaluation and implementation

1. Genomewide association studies (GWAS) enable detailed dissections of the genetic basis for organisms' ability to adapt to a changing environment. In long-term studies of natural populations, individuals are often marked at one point in their life and then repeatedly recaptured. It is therefore essential that a method for GWAS includes the process of repeated sampling. In a GWAS, the effects of thousands of single-nucleotide polymorphisms (SNPs) need to be fitted and any model development is constrained by the computational requirements. A method is therefore required that can fit a highly hierarchical model and at the same time is computationally fast enough to be useful. 2. Our method fits fixed SNP effects in a linear mixed model that can include both random polygenic effects and permanent environmental effects. In this way, the model can correct for population structure and model repeated measures. The covariance structure of the linear mixed model is first estimated and subsequently used in a generalized least squares setting to fit the SNP effects. The method was evaluated in a simulation study based on observed genotypes from a long-term study of collared flycatchers in Sweden. 3. The method we present here was successful in estimating permanent environmental effects from simulated repeated measures data. Additionally, we found that especially for variable phenotypes having large variation between years, the repeated measurements model has a substantial increase in power compared to a model using average phenotypes as a response. 4. The method is available in the R package RepeatABEL. It increases the power in GWAS having repeated measures, especially for long-term studies of natural populations, and the R implementation is expected to facilitate modelling of longitudinal data for studies of both animal and human populations.

Multicentric Genome-Wide Association Study for Primary Spontaneous Pneumothorax

Despite elevated incidence and recurrence rates for Primary Spontaneous Pneumothorax (PSP), little is known about its etiology, and the genetics of idiopathic PSP remains unexplored. To identify genetic variants contributing to sporadic PSP risk, we conducted the first PSP genome-wide association study. Two replicate pools of 92 Portuguese PSP cases and of 129 age- and sex-matched controls were allelotyped in triplicate on the Affymetrix Human SNP Array 6.0 arrays. Markers passing quality control were ranked by relative allele score difference between cases and controls (|RASdiff|), by a novel cluster method and by a combined Z-test. 101 single nucleotide polymorphisms (SNPs) were selected using these three approaches for technical validation by individual genotyping in the discovery dataset. 87 out of 94 successfully tested SNPs were nominally associated in the discovery dataset. Replication of the 87 technically validated SNPs was then carried out in an independent replication dataset of 100 Portuguese cases and 425 controls. The intergenic rs4733649 SNP in chromosome 8 (between LINC00824 and LINC00977) was associated with PSP in the discovery (P = 4.07E-03, ORC[95% CI] = 1.88[1.22-2.89]), replication (P = 1.50E-02, ORC[95% CI] = 1.50[1.08-2.09]) and combined datasets (P = 8.61E-05, ORC[95% CI] = 1.65[1.29-2.13]). This study identified for the first time one genetic risk factor for sporadic PSP, but future studies are warranted to further confirm this finding in other populations and uncover its functional role in PSP pathogenesis.

Association between the il28b rs12979860 polymorphism and therapy response in patients infected with genotype 1 of hepatitis c virus in Central Brazil

v. 45, n.2, p.152-160, abr/.jun. 2016.

Genome Wide Association Studies of Phagocytosis and the Cellular Immune Response in Drosophila melanogaster.

Phagocytosis of bacteria by specialized blood cells, known as hemocytes, is a vital component of Drosophila cellular immunity. To identify novel genes that mediate the cellular response to bacteria, we conducted three separate genetic screens using the Drosophila Genetic Reference Panel (DGRP). Adult DGRP lines were tested for the ability of their hemocytes to phagocytose the Gram-positive bacteria Staphylococcus aureus or the Gram-negative bacteria Escherichia coli. The DGRP lines were also screened for the ability of their hemocytes to clear S. aureus infection through the process of phagosome maturation. Genome-wide association analyses were performed to identify potentially relevant single nucleotide polymorphisms (SNPs) associated with the cellular immune phenotypes. The S. aureus phagosome maturation screen identified SNPs near or in 528 candidate genes, many of which have no known role in immunity. Three genes, dpr10, fred, and CG42673, were identified whose loss-of-function in blood cells significantly impaired the innate immune response to S. aureus. The DGRP S. aureus screens identified variants in the gene, Ataxin 2 Binding Protein-1 (A2bp1) as important for the cellular immune response to S. aureus. A2bp1 belongs to the highly conserved Fox-1 family of RNA-binding proteins. Genetic studies revealed that A2bp1 transcript levels must be tightly controlled for hemocytes to successfully phagocytose S. aureus. The transcriptome of infected and uninfected hemocytes from wild type and A2bp1 mutant flies was analyzed and it was found that A2bp1 negatively regulates the expression of the Immunoglobulin-superfamily member Down syndrome adhesion molecule 4 (Dscam4). Silencing of A2bp1 and Dscam4 in hemocytes rescues the fly’s immune response to S. aureus indicating that Dscam4 negatively regulates S. aureus phagocytosis. Overall, we present an examination of the cellular immune response to bacteria with the aim of identifying and characterizing roles for novel mediators of innate immunity in Drosophila. By screening panel of lines in which all genetic variants are known, we successfully identified a large set of candidate genes that could provide a basis for future studies of Drosophila cellular immunity. Finally, we describe a novel, immune-specific role for the highly conserved Fox-1 family member, A2bp1.

Bioanalytical applications of microfluidic devices

The first part of the thesis describes a new patterning technique--microfluidic contact printing--that combines several of the desirable aspects of microcontact printing and microfluidic patterning and addresses some of their important limitations through the integration of a track-etched polycarbonate (PCTE) membrane. Using this technique, biomolecules (e.g., peptides, polysaccharides, and proteins) were printed in high fidelity on a receptor modified polyacrylamide hydrogel substrate. The patterns obtained can be controlled through modifications of channel design and secondary programming via selective membrane wetting. The protocols support the printing of multiple reagents without registration steps and fast recycle times. The second part describes a non-enzymatic, isothermal method to discriminate single nucleotide polymorphisms (SNPs). SNP discrimination using alkaline dehybridization has long been neglected because the pH range in which thermodynamic discrimination can be done is quite narrow. We found, however, that SNPs can be discriminated by the kinetic differences exhibited in the dehybridization of PM and MM DNA duplexes in an alkaline solution using fluorescence microscopy. We combined this method with multifunctional encoded hydrogel particle array (fabricated by stop-flow lithography) to achieve fast kinetics and high versatility. This approach may serve as an effective alternative to temperature-based method for analyzing unamplified genomic DNA in point-of-care diagnostic.