Priming of a DNA metabarcoding approach for species identification and inventory in marine macrobenthic communities

Dissertação de mestrado em Molecular Genetics

SUCLG2 identified as both a determinator of CSF Aβ1-42 levels and an attenuator of cognitive decline in Alzheimer's disease.

Cerebrospinal fluid amyloid-beta 1-42 (Aβ1-42) and phosphorylated Tau at position 181 (pTau181) are biomarkers of Alzheimer's disease (AD). We performed an analysis and meta-analysis of genome-wide association study data on Aβ1-42 and pTau181 in AD dementia patients followed by independent replication. An association was found between Aβ1-42 level and a single-nucleotide polymorphism in SUCLG2 (rs62256378) (P = 2.5×10(-12)). An interaction between APOE genotype and rs62256378 was detected (P = 9.5 × 10(-5)), with the strongest effect being observed in APOE-ε4 noncarriers. Clinically, rs62256378 was associated with rate of cognitive decline in AD dementia patients (P = 3.1 × 10(-3)). Functional microglia experiments showed that SUCLG2 was involved in clearance of Aβ1-42.

Methylation profile of TP53 regulatory pathway and mtDNA alterations in breast cancer patients lacking TP53 mutations.

The present study investigated promoter hypermethylation of TP53 regulatory pathways providing a potential link between epigenetic changes and mitochondrial DNA (mtDNA) alterations in breast cancer patients lacking a TP53 mutation. The possibility of using the cancer-specific alterations in serum samples as a blood-based test was also explored. Triple-matched samples (cancerous tissues, matched adjacent normal tissues and serum samples) from breast cancer patients were screened for TP53 mutations, and the promoter methylation profile of P14(ARF), MDM2, TP53 and PTEN genes was analyzed as well as mtDNA alterations, including D-loop mutations and mtDNA content. In the studied cohort, no mutation was found in TP53 (DNA-binding domain). Comparison of P14(ARF) and PTEN methylation patterns showed significant hypermethylation levels in tumor tissues (P < 0.05 and <0.01, respectively) whereas the TP53 tumor suppressor gene was not hypermethylated (P < 0.511). The proportion of PTEN methylation was significantly higher in serum than in the normal tissues and it has a significant correlation to tumor tissues (P < 0.05). mtDNA analysis revealed 36.36% somatic and 90.91% germline mutations in the D-loop region and also significant mtDNA depletion in tumor tissues (P < 0.01). In addition, the mtDNA content in matched serum was significantly lower than in the normal tissues (P < 0.05). These data can provide an insight into the management of a therapeutic approach based on the reversal of epigenetic silencing of the crucial genes involved in regulatory pathways of the tumor suppressor TP53. Additionally, release of significant aberrant methylated PTEN in matched serum samples might represent a promising biomarker for breast cancer.

A genome-wide association study of early menopause and the combined impact of identified variants.

Early menopause (EM) affects up to 10% of the female population, reducing reproductive lifespan considerably. Currently, it constitutes the leading cause of infertility in the western world, affecting mainly those women who postpone their first pregnancy beyond the age of 30 years. The genetic aetiology of EM is largely unknown in the majority of cases. We have undertaken a meta-analysis of genome-wide association studies (GWASs) in 3493 EM cases and 13 598 controls from 10 independent studies. No novel genetic variants were discovered, but the 17 variants previously associated with normal age at natural menopause as a quantitative trait (QT) were also associated with EM and primary ovarian insufficiency (POI). Thus, EM has a genetic aetiology which overlaps variation in normal age at menopause and is at least partly explained by the additive effects of the same polygenic variants. The combined effect of the common variants captured by the single nucleotide polymorphism arrays was estimated to account for ∼30% of the variance in EM. The association between the combined 17 variants and the risk of EM was greater than the best validated non-genetic risk factor, smoking.

Caffeine intake and CYP1A2 variants associated with high caffeine intake protect non-smokers from hypertension.

The 15q24.1 locus, including CYP1A2, is associated with blood pressure (BP). The CYP1A2 rs762551 C allele is associated with lower CYP1A2 enzyme activity. CYP1A2 metabolizes caffeine and is induced by smoking. The association of caffeine consumption with hypertension remains controversial. We explored the effects of CYP1A2 variants and CYP1A2 enzyme activity on BP, focusing on caffeine as the potential mediator of CYP1A2 effects. Four observational (n = 16 719) and one quasi-experimental studies (n = 106) including European adults were conducted. Outcome measures were BP, caffeine intake, CYP1A2 activity and polymorphisms rs762551, rs1133323 and rs1378942. CYP1A2 variants were associated with hypertension in non-smokers, but not in smokers (CYP1A2-smoking interaction P = 0.01). Odds ratios (95% CIs) for hypertension for rs762551 CC, CA and AA genotypes were 1 (reference), 0.78 (0.59-1.02) and 0.66 (0.50-0.86), respectively, P = 0.004. Results were similar for the other variants. Higher CYP1A2 activity was linearly associated with lower BP after quitting smoking (P = 0.049 and P = 0.02 for systolic and diastolic BP, respectively), but not while smoking. In non-smokers, the CYP1A2 variants were associated with higher reported caffeine intake, which in turn was associated with lower odds of hypertension and lower BP (P = 0.01). In Mendelian randomization analyses using rs1133323 as instrument, each cup of caffeinated beverage was negatively associated with systolic BP [-9.57 (-16.22, -2.91) mmHg]. The associations of CYP1A2 variants with BP were modified by reported caffeine intake. These observational and quasi-experimental results strongly support a causal role of CYP1A2 in BP control via caffeine intake.

Selective disruption of Tcf7l2 in the pancreatic β cell impairs secretory function and lowers β cell mass.

Type 2 diabetes (T2D) is characterized by β cell dysfunction and loss. Single nucleotide polymorphisms in the T-cell factor 7-like 2 (TCF7L2) gene, associated with T2D by genome-wide association studies, lead to impaired β cell function. While deletion of the homologous murine Tcf7l2 gene throughout the developing pancreas leads to impaired glucose tolerance, deletion in the β cell in adult mice reportedly has more modest effects. To inactivate Tcf7l2 highly selectively in β cells from the earliest expression of the Ins1 gene (∼E11.5) we have therefore used a Cre recombinase introduced at the Ins1 locus. Tcfl2(fl/fl)::Ins1Cre mice display impaired oral and intraperitoneal glucose tolerance by 8 and 16 weeks, respectively, and defective responses to the GLP-1 analogue liraglutide at 8 weeks. Tcfl2(fl/fl)::Ins1Cre islets displayed defective glucose- and GLP-1-stimulated insulin secretion and the expression of both the Ins2 (∼20%) and Glp1r (∼40%) genes were significantly reduced. Glucose- and GLP-1-induced intracellular free Ca(2+) increases, and connectivity between individual β cells, were both lowered by Tcf7l2 deletion in islets from mice maintained on a high (60%) fat diet. Finally, analysis by optical projection tomography revealed ∼30% decrease in β cell mass in pancreata from Tcfl2(fl/fl)::Ins1Cre mice. These data demonstrate that Tcf7l2 plays a cell autonomous role in the control of β cell function and mass, serving as an important regulator of gene expression and islet cell coordination. The possible relevance of these findings for the action of TCF7L2 polymorphisms associated with Type 2 diabetes in man is discussed.

Characterization of Endotrypanum (Kinetoplastida: Trypanosomatidae), a Unique Parasite Infecting the Neotropical Tree Sloths (Edentata)

This article reviews current concepts of the biology of Endotrypanum spp. Data summarized here on parasite classification and taxonomic divergence found among these haemoflagellates come from our studies of molecular characterization of Endotrypanum stocks (representing an heterogenous population of reference strains and isolates from the Brazilian Amazon region) and from scientific literature. Using numerical zymotaxonomy we have demonstrated genetic diversity among these parasites. The molecular trees obtained revealed that there are, at least, three groups (distinct species?) of Endotrypanum, which are distributed in Central and South America. In concordance with this classification of the parasites there are further newer molecular data obtained using distinct markers. Moreover, comparative studies (based on the molecular genetics of the organisms) have shown the phylogenetic relationships between some Endotrypanum and related kinetoplastid lineages.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.

Multivariate QST-FST comparisons: a neutrality test for the evolution of the g matrix in structured populations.

Neutrality tests in quantitative genetics provide a statistical framework for the detection of selection on polygenic traits in wild populations. However, the existing method based on comparisons of divergence at neutral markers and quantitative traits (Q(st)-F(st)) suffers from several limitations that hinder a clear interpretation of the results with typical empirical designs. In this article, we propose a multivariate extension of this neutrality test based on empirical estimates of the among-populations (D) and within-populations (G) covariance matrices by MANOVA. A simple pattern is expected under neutrality: D = 2F(st)/(1 - F(st))G, so that neutrality implies both proportionality of the two matrices and a specific value of the proportionality coefficient. This pattern is tested using Flury's framework for matrix comparison [common principal-component (CPC) analysis], a well-known tool in G matrix evolution studies. We show the importance of using a Bartlett adjustment of the test for the small sample sizes typically found in empirical studies. We propose a dual test: (i) that the proportionality coefficient is not different from its neutral expectation [2F(st)/(1 - F(st))] and (ii) that the MANOVA estimates of mean square matrices between and among populations are proportional. These two tests combined provide a more stringent test for neutrality than the classic Q(st)-F(st) comparison and avoid several statistical problems. Extensive simulations of realistic empirical designs suggest that these tests correctly detect the expected pattern under neutrality and have enough power to efficiently detect mild to strong selection (homogeneous, heterogeneous, or mixed) when it is occurring on a set of traits. This method also provides a rigorous and quantitative framework for disentangling the effects of different selection regimes and of drift on the evolution of the G matrix. We discuss practical requirements for the proper application of our test in empirical studies and potential extensions.

Genome-wide meta-analysis of common variant differences between men and women.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10(-8)) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ~115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

Multi-ethnic fine-mapping of 14 central adiposity loci.

The Genetic Investigation of Anthropometric Traits (GIANT) consortium identified 14 loci in European Ancestry (EA) individuals associated with waist-to-hip ratio (WHR) adjusted for body mass index. These loci are wide and narrowing the signals remains necessary. Twelve of 14 loci identified in GIANT EA samples retained strong associations with WHR in our joint EA/individuals of African Ancestry (AA) analysis (log-Bayes factor >6.1). Trans-ethnic analyses at five loci (TBX15-WARS2, LYPLAL1, ADAMTS9, LY86 and ITPR2-SSPN) substantially narrowed the signals to smaller sets of variants, some of which are in regions that have evidence of regulatory activity. By leveraging varying linkage disequilibrium structures across different populations, single-nucleotide polymorphisms (SNPs) with strong signals and narrower credible sets from trans-ethnic meta-analysis of central obesity provide more precise localizations of potential functional variants and suggest a possible regulatory role. Meta-analysis results for WHR were obtained from 77 167 EA participants from GIANT and 23 564 AA participants from the African Ancestry Anthropometry Genetics Consortium. For fine mapping we interrogated SNPs within ± 250 kb flanking regions of 14 previously reported index SNPs from loci discovered in EA populations by performing trans-ethnic meta-analysis of results from the EA and AA meta-analyses. We applied a Bayesian approach that leverages allelic heterogeneity across populations to combine meta-analysis results and aids in fine-mapping shared variants at these locations. We annotated variants using information from the ENCODE Consortium and Roadmap Epigenomics Project to prioritize variants for possible functionality.

Persistent pulmonary arterial hypertension in the newborn (PPHN): a frequent manifestation of TMEM70 defective patients.

INTRODUCTION: Mutations in the TMEM70 are the most common cause of nuclear ATP synthase deficiency resulting in a distinctive phenotype characterized by severe neonatal hypotonia, hypertrophic cardiomyopathy (HCMP), facial dysmorphism, severe lactic acidosis, hyperammonemia and 3-methylglutaconic aciduria (3-MGA). METHODS AND RESULTS: We collected 9 patients with genetically confirmed TMEM70 defect from 8 different families. Six were homozygous for the c.317-2A>G mutation, 2 were compound heterozygous for mutations c.317-2A>G and c.628A>C and 1 was homozygous for the novel c.701A>C mutation. Generalized hypotonia, lactic acidosis, hyperammonemia and 3-MGA were present in all since birth. Five patients presented acute respiratory distress at birth requiring intubation and ventilatory support. HCMP was detected in 5 newborns and appeared a few months later in 3 additional children. Five patients showed a severe and persistent neonatal pulmonary hypertension (PPHN) requiring Nitric Oxide (NO) and/or sildenafil administration combined in 2 cases with high-frequency oscillatory (HFO) ventilation. In 3 of these patients, echocardiography detected signs of HCMP at birth. CONCLUSIONS: PPHN is a life-threatening poorly understood condition with bad prognosis if untreated. Pulmonary hypertension has rarely been reported in mitochondrial disorders and, so far, it has been described in association with TMEM70 deficiency only in one patient. This report further expands the clinical and genetic spectrum of the syndrome indicating PPHN as a frequent and life-threatening complication regardless of the type of mutation. Moreover, in these children PPHN appears even in the absence of an overt cardiomyopathy, thus representing an early sign and a clue for diagnosis.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy.

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

The transmembrane serine protease (TMPRSS3) mutated in deafness DFNB8/10 activates the epithelial sodium channel (ENaC) in vitro.

TMPRSS3 encodes a transmembrane serine protease that contains both LDLRA and SRCR domains and is mutated in non-syndromic autosomal recessive deafness (DFNB8/10). To study its function, we cloned the mouse ortholog which maps to Mmu17, which is structurally similar to the human gene and encodes a polypeptide with 88% identity to the human protein. RT-PCR and RNA in situ hybridization on rat and mouse cochlea revealed that Tmprss3 is expressed in the spiral ganglion, the cells supporting the organ of Corti and the stria vascularis. RT-PCR on mouse tissues showed expression in the thymus, stomach, testis and E19 embryos. Transient expression of wild-type or tagged TMPRSS3 protein showed a primary localization in the endoplasmic reticulum. The epithelial amiloride-sensitive sodium channel (ENaC), which is expressed in many sodium-reabsorbing tissues including the inner ear and is regulated by membrane-bound channel activating serine proteases (CAPs), is a potential substrate of TMPRSS3. In the Xenopus oocyte expression system, proteolytic processing of TMPRSS3 was associated with increased ENaC mediated currents. In contrast, 6 TMPRSS3 mutants (D103G, R109W, C194F, W251C, P404L, C407R) causing deafness and a mutant in the catalytic triad of TMPRSS3 (S401A), failed to undergo proteolytic cleavage and activate ENaC. These data indicate that important signaling pathways in the inner ear are controlled by proteolytic cleavage and suggest: (i) the existence of an auto-catalytic processing by which TMPRSS3 would become active, and (ii) that ENaC could be a substrate of TMPRSS3 in the inner ear.