A Bayesian Dirichlet-Multinomial Test for Cross-Group Differences

Testing for differences within data sets is an important issue across various applications. Our work is primarily motivated by the analysis of microbiomial composition, which has been increasingly relevant and important with the rise of DNA sequencing. We first review classical frequentist tests that are commonly used in tackling such problems. We then propose a Bayesian Dirichlet-multinomial framework for modeling the metagenomic data and for testing underlying differences between the samples. A parametric Dirichlet-multinomial model uses an intuitive hierarchical structure that allows for flexibility in characterizing both the within-group variation and the cross-group difference and provides very interpretable parameters. A computational method for evaluating the marginal likelihoods under the null and alternative hypotheses is also given. Through simulations, we show that our Bayesian model performs competitively against frequentist counterparts. We illustrate the method through analyzing metagenomic applications using the Human Microbiome Project data.

Behavioural Genetics in the Postgenomics Era

There is growing evidence that the complexity of higher organisms does not correlate with the ‘complexity’ of the genome (the human genome contains fewer protein coding genes than corn, and many genes are preserved across species). Rather, complexity is associated with the complexity of the pathways and processes whereby the cell utilises the deoxyribonucleic acid molecule, and much else, in the process of phenotype formation. These pro- cesses include the activity of the epigenome, noncoding ribonucleic acids, alternative splicing and post-transla- tional modifications. Not accidentally, all of these pro- cesses appear to be of particular importance for the human brain, the most complex organ in nature. Because these processes can be highly environmentally reactive, they are a key to understanding behavioural plasticity and highlight the importance of the developmental process in explaining behavioural outcomes.

RECQ5 promotes recombination and mutagenesis at targeted nicks through disruption of RAD51 filaments

Thesis (Ph.D.)--University of Washington, 2016-08

Pharmacogenetics, pharmacogenomics and airway disease

The availability of a draft sequence for the human genome will revolutionise research into airway disease. This review deals with two of the most important areas impinging on the treatment of patients: pharmacogenetics and pharmacogenomics. Considerable inter-individual variation exists at the DNA level in targets for medication, and variability in response to treatment may, in part, be determined by this genetic variation. Increased knowledge about the human genome might also permit the identification of novel therapeutic targets by expression profiling at the RNA (genomics) or protein (proteomics) level. This review describes recent advances in pharmacogenetics and pharmacogenomics with regard to airway disease

Doença enxerto-versus hospedeiro: determinação da expressão de genes e alterações de vias associadas ao sistema imunológico usando metodos da biofísica

O transplante de medula óssea (TMO) é um procedimento terapêutico importante em casos relacionados à pacientes com leucemia ou linfoma. Em decorrência desse processo, uma reação conhecida como doença enxerto-versus-hospedeiro (GVHD) pode ocorrer em pacientes susceptíveis como conseqüência da presença de células imunocompetentes do doador. Entretanto, não existe um modelo para descrever completamente as ações relacionadas ao mecanismo imunológico da GVHD desde a fase que inicializa a doença até a fase efetora. O Objetivo geral deste estudo é a investigação da resposta imunológica considerando-se o sistema HLA (antígenos leucocitários humano) em pacientes que desenvolveram a GVHD em decorrência do TMO. O National Cancer Institute (NCI) – Pathway interaction Database e Reactome foram usados como bases de dados com o objetivo de se estudar a expressão de genes e vias relacionados às Classes I e II do sistema HLA (antígenos leucocitários humano). O estudo considerou a mudança de expressão de genes relacionados às 17 vias do sistema imunológico com potencialidade para se expressar em pacientes que desenvolveram a GVHD associada à TMO. Dados referentes aos transcriptomas foram obtidos utilizando-se a plataforma GPL570 Affymetrix Genoma Humano U133 Plus. A atividade relativa foi usada para determinar as alterações das vias em amostras de GVHD em relação ao controle. As análises foram realizadas utilizando-se o software Via Complex e Bioconductor. Observou-se aumento significativo da expressão de genes ralacionados às vias do sistema imune adaptativo, antígenos associados às Classe I e II do HLA, fosforilação de CD3 e CD247, sinalização dos receptores de células T em CD4+ nativas e ativação de NF-kapa β nas células B. Também observou-se alterações significativas na mudança de expressão dos genes associados às vias relacionadas à super família de moléculas B7:CD28\CTLA-4 quando comparadas ao controle. Isso pode indicar a necessidade de geração de um segundo sinal co-estimulador em GVHD, acionado pelas moléculas dessa super família. O aumento da expressão do gene CD69 nas amostras experimentais caracteriza a ativação celular e, portanto, a sinalização de estímulos em GVHD. Os achados obtidos neste estudo contribuem para melhor elucidar o mecanismo imunopatogênico associado à GVHD. P

Increased frequency of the autism broader phenotype in mothers transmitting etiological CNVs to sons affected by Autism Spectrum Disorder (ASD)

Autism Spectrum Disorder (ASD) is a frequent and complex neurodevelopmental disorder, characterized by impairments in social communication and repetitive behaviors and with a high male to female ratio: ~4:1. Genetic factors, including rare Copy Number Variants (CNVs), have a substantial impact in ASD risk1, and are associated with specific phenotypic manifestations2. Recent studies reported that rare inherited CNVs are enriched in mothers of ASD children compared with mothers of controls and are preferentially transmitted from mothers to ASD children suggesting a sex bias in CNV transmission; further, the imbalanced transmission of small pathogenic CNVs from unaffected mothers to their sons with ASD has been described3,4. An increased prevalence of autism-like personality traits is found in unaffected relatives of ASD children, suggesting a genetic liability of a broader autism phenotype (BAP)5. The BAP in parents of autistic children can be assessed by the Social Responsiveness Scale (SRS)6 and Broad Autism Phenotype Questionnaire (BAPQ)7 reports. The SRS is 65-item questionnaire to identify sub-clinical social impairments and interpersonal behaviour in individuals. The BAPQ is a 36-item questionnaire measures social aloofness, rigid personality, and pragmatic language deficits in both parents and children.

Epidemiology and genetic architecture of blood pressure: a family based study of Generation Scotland

Hypertension is a major risk factor for cardiovascular disease and mortality, and a growing global public health concern, with up to one-third of the world’s population affected. Despite the vast amount of evidence for the benefits of blood pressure (BP) lowering accumulated to date, elevated BP is still the leading risk factor for disease and disability worldwide. It is well established that hypertension and BP are common complex traits, where multiple genetic and environmental factors contribute to BP variation. Furthermore, family and twin studies confirmed the genetic component of BP, with a heritability estimate in the range of 30-50%. Contemporary genomic tools enabling the genotyping of millions of genetic variants across the human genome in an efficient, reliable, and cost-effective manner, has transformed hypertension genetics research. This is accompanied by the presence of international consortia that have offered unprecedentedly large sample sizes for genome-wide association studies (GWASs). While GWAS for hypertension and BP have identified more than 60 loci, variants in these loci are associated with modest effects on BP and in aggregate can explain less than 3% of the variance in BP. The aims of this thesis are to study the genetic and environmental factors that influence BP and hypertension traits in the Scottish population, by performing several genetic epidemiological analyses. In the first part of this thesis, it aims to study the burden of hypertension in the Scottish population, along with assessing the familial aggregation and heritialbity of BP and hypertension traits. In the second part, it aims to validate the association of common SNPs reported in the large GWAS and to estimate the variance explained by these variants. In this thesis, comprehensive genetic epidemiology analyses were performed on Generation Scotland: Scottish Family Health Study (GS:SFHS), one of the largest population-based family design studies. The availability of clinical, biological samples, self-reported information, and medical records for study participants has allowed several assessments to be performed to evaluate factors that influence BP variation in the Scottish population. Of the 20,753 subjects genotyped in the study, a total of 18,470 individuals (grouped into 7,025 extended families) passed the stringent quality control (QC) criteria and were available for all subsequent analysis. Based on the BP-lowering treatment exposure sources, subjects were further classified into two groups. First, subjects with both a self-reported medications (SRMs) history and electronic-prescription records (EPRs; n =12,347); second, all the subjects with at least one medication history source (n =18,470). In the first group, the analysis showed a good concordance between SRMs and EPRs (kappa =71%), indicating that SRMs can be used as a surrogate to assess the exposure to BP-lowering medication in GS:SFHS participants. Although both sources suffer from some limitations, SRMs can be considered the best available source to estimate the drug exposure history in those without EPRs. The prevalence of hypertension was 40.8% with higher prevalence in men (46.3%) compared to women (35.8%). The prevalence of awareness, treatment and controlled hypertension as defined by the study definition were 25.3%, 31.2%, and 54.3%, respectively. These findings are lower than similar reported studies in other populations, with the exception of controlled hypertension prevalence, which can be considered better than other populations. Odds of hypertension were higher in men, obese or overweight individuals, people with a parental history of hypertension, and those living in the most deprived area of Scotland. On the other hand, deprivation was associated with higher odds of treatment, awareness and controlled hypertension, suggesting that people living in the most deprived area may have been receiving better quality of care, or have higher comorbidity levels requiring greater engagement with doctors. These findings highlight the need for further work to improve hypertension management in Scotland. The family design of GS:SFHS has allowed family-based analysis to be performed to assess the familial aggregation and heritability of BP and hypertension traits. The familial correlation of BP traits ranged from 0.07 to 0.20, and from 0.18 to 0.34 for parent-offspring pairs and sibling pairs, respectively. A higher correlation of BP traits was observed among first-degree relatives than other types of relative pairs. A variance-component model that was adjusted for sex, body mass index (BMI), age, and age-squared was used to estimate heritability of BP traits, which ranged from 24% to 32% with pulse pressure (PP) having the lowest estimates. The genetic correlation between BP traits showed a high correlation between systolic (SBP), diastolic (DBP) and mean arterial pressure (MAP) (G: 81% to 94%), but lower correlations with PP (G: 22% to 78%). The sibling recurrence risk ratio (λS) for hypertension and treatment were calculated as 1.60 and 2.04 respectively. These findings confirm the genetic components of BP traits in GS:SFHS, and justify further work to investigate genetic determinants of BP. Genetic variants reported in the recent large GWAS of BP traits were selected for genotyping in GS:SFHS using a custom designed TaqMan® OpenArray®. The genotyping plate included 44 single nucleotide polymorphisms (SNPs) that have been previously reported to be associated with BP or hypertension at genome-wide significance level. A linear mixed model that is adjusted for age, age-squared, sex, and BMI was used to test for the association between the genetic variants and BP traits. Of the 43 variants that passed the QC, 11 variants showed statistically significant association with at least one BP trait. The phenotypic variance explained by these variant for the four BP traits were 1.4%, 1.5%, 1.6%, and 0.8% for SBP, DBP, MAP, and PP, respectively. The association of genetic risk score (GRS) that were constructed from selected variants has showed a positive association with BP level and hypertension prevalence, with an average effect of one mmHg increase with each 0.80 unit increases in the GRS across the different BP traits. The impact of BP-lowering medication on the genetic association study for BP traits has been established, with typical practice of adding a fixed value (i.e. 15/10 mmHg) to the measured BP values to adjust for BP treatment. Using the subset of participants with the two treatment exposure sources (i.e. SRMs and EPRs), the influence of using either source to justify the addition of fixed values in SNP association signal was analysed. BP phenotypes derived from EPRs were considered the true phenotypes, and those derived from SRMs were considered less accurate, with some phenotypic noise. Comparing SNPs association signals between the four BP traits in the two model derived from the different adjustments showed that MAP was the least impacted by the phenotypic noise. This was suggested by identifying the same overlapped significant SNPs for the two models in the case of MAP, while other BP traits had some discrepancy between the two sources

Tissue Proteomes: Quantitative Mass Spectrometry of Murine Liver and Ovarian Endometrioma

A human genome contains more than 20 000 protein-encoding genes. A human proteome, instead, has been estimated to be much more complex and dynamic. The most powerful tool to study proteins today is mass spectrometry (MS). MS based proteomics is based on the measurement of the masses of charged peptide ions in a gas-phase. The peptide amino acid sequence can be deduced, and matching proteins can be found, using software to correlate MS-data with sequence database information. Quantitative proteomics allow the estimation of the absolute or relative abundance of a certain protein in a sample. The label-free quantification methods use the intrinsic MS-peptide signals in the calculation of the quantitative values enabling the comparison of peptide signals from numerous patient samples. In this work, a quantitative MS methodology was established to study aromatase overexpressing (AROM+) male mouse liver and ovarian endometriosis tissue samples. The workflow of label-free quantitative proteomics was optimized in terms of sensitivity and robustness, allowing the quantification of 1500 proteins with a low coefficient of variance in both sample types. Additionally, five statistical methods were evaluated for the use with label-free quantitative proteomics data. The proteome data was integrated with other omics datasets, such as mRNA microarray and metabolite data sets. As a result, an altered lipid metabolism in liver was discovered in male AROM+ mice. The results suggest a reduced beta oxidation of long chain phospholipids in the liver and increased levels of pro-inflammatory fatty acids in the circulation in these mice. Conversely, in the endometriosis tissues, a set of proteins highly specific for ovarian endometrioma were discovered, many of which were under the regulation of the growth factor TGF-β1. This finding supports subsequent biomarker verification in a larger number of endometriosis patient samples.

Transcriptoma diferencial entre células-tronco mesenquimais humanas jovens e senescentes

Human mesenchymal stem cells (MSC) are powerful sources for cell therapy in regenerative medicine. The long time cultivation can result in replicative senescence or can be related to the emergence of chromosomal alterations responsible for the acquisition of tumorigenesis features in vitro. In this study, for the first time, the expression profile of MSC with a paracentric chromosomal inversion (MSC/inv) was compared to normal karyotype (MSC/n) in early and late passages. Furthermore, we compared the transcriptome of each MSC in early passages with late passages. MSC used in this study were obtained from the umbilical vein of three donors, two MSC/n and one MSC/inv. After their cryopreservation, they have been expanded in vitro until reached senescence. Total RNA was extracted using the RNeasy mini kit (Qiagen) and marked with the GeneChip ® 3 IVT Express Kit (Affymetrix Inc.). Subsequently, the fragmented aRNA was hybridized on the microarranjo Affymetrix Human Genome U133 Plus 2.0 arrays (Affymetrix Inc.). The statistical analysis of differential gene expression was performed between groups MSC by the Partek Genomic Suite software, version 6.4 (Partek Inc.). Was considered statistically significant differences in expression to p-value Bonferroni correction ˂.01. Only signals with fold change ˃ 3.0 were included in the list of differentially expressed. Differences in gene expression data obtained from microarrays were confirmed by Real Time RT-PCR. For the interpretation of biological expression data were used: IPA (Ingenuity Systems) for analysis enrichment functions, the STRING 9.0 for construction of network interactions; Cytoscape 2.8 to the network visualization and analysis bottlenecks with the aid of the GraphPad Prism 5.0 software. BiNGO Cytoscape pluggin was used to access overrepresentation of Gene Ontology categories in Biological Networks. The comparison between senescent and young at each group of MSC has shown that there is a difference in the expression parttern, being higher in the senescent MSC/inv group. The results also showed difference in expression profiles between the MSC/inv versus MSC/n, being greater when they are senescent. New networks were identified for genes related to the response of two of MSC over cultivation time. Were also identified genes that can coordinate functional categories over represented at networks, such as CXCL12, SFRP1, xvi EGF, SPP1, MMP1 e THBS1. The biological interpretation of these data suggests that the population of MSC/inv has different constitutional characteristics, related to their potential for differentiation, proliferation and response to stimuli, responsible for a distinct process of replicative senescence in MSC/inv compared to MSC/n. The genes identified in this study are candidates for biomarkers of cellular senescence in MSC, but their functional relevance in this process should be evaluated in additional in vitro and/or in vivo assays

Novel Mutation in the ATP-Binding Cassette Transporter A3 (ABCA3) Encoding Gene Causes Respiratory Distress Syndrome in A Term Newborn in Southwest Iran

Introduction: ABCA3 glycoprotein belongs to the ATP-binding cassette (ABC) superfamily of transporters, which utilize the energy derived from hydrolysis of ATP for the translocation of a wide variety of substrates across the plasma membrane. Mutations in the ABCA3 gene are knowingly causative for fatal surfactant deficiency, particularly respiratory distress syndrome (RDS) in term babies. Case Presentation: In this study, Sanger sequencing of the whole ABCA3 gene (NCBI NM_001089) was performed in a neonatal boy with severe RDS. A homozygous mutation has been identified in the patient. Parents were heterozygous for the same missense mutation GGA > AGA at position 202 in exon 6 of the ABCA3 gene (c.604G > A; p.G202R). Furthermore, 70 normal individuals have been analyzed for the mentioned change with negative results. Conclusions: Regarding Human Genome Mutation Database (HGMD) and other literature recherche, the detected change is a novel mutation and has not been reported before. Bioinformatics mutation predicting tools prefer it as pathogenic.

Comparative survey of the relative impact of mRNA features on local ribosome profiling read density

Ribosome profiling (Ribo-seq), a promising technology for exploring ribosome decoding rates, is characterized by the presence of infrequent high peaks in ribosome footprint density and by long alignment gaps. Here, to reduce the impact of data heterogeneity we introduce a simple normalization method, Ribo-seq Unit Step Transformation (RUST). RUST is robust and outperforms other normalization techniques in the presence of heterogeneous noise. We illustrate how RUST can be used for identifying mRNA sequence features that affect ribosome footprint densities globally. We show that a few parameters extracted with RUST are sufficient for predicting experimental densities with high accuracy. Importantly the application of RUST to 30 publicly available Ribo-seq data sets revealed a substantial variation in sequence determinants of ribosome footprint frequencies, questioning the reliability of Ribo-seq as an accurate representation of local ribosome densities without prior quality control. This emphasizes our incomplete understanding of how protocol parameters affect ribosome footprint densities.

Microbioma oral humano

O microbioma oral humano é constituído por um vasto conjunto de microrganismos presentes na cavidade oral. Analisando a cavidade oral podemos verificar que nela existem mais de 700 espécies de bactérias responsáveis pelo domínio de parte do microbioma humano, tornando-a um importante local de estudo. É um dos habitats com maior diversidade no corpo humano onde esses microrganismos se apresentam de forma organizada e estruturada. Estes habitats estão intimamente relacionados com o desenvolvimento do sistema imunitário e com a proteção contra agentes patogénicos. O microbioma oral é único e específico em cada indivíduo, sofrendo variações em indivíduos diferentes. Na origem da diversidade do microbioma oral estão associados fatores como genética, dieta e localização geográfica, tendo também grande importância a localização anatómica e a idade do indivíduo. O Projeto Microbioma Humano surgiu com a finalidade de identificar diversos microrganismos presentes no ser humano, bem como compreender os principais fatores responsáveis pelas suas alterações. O estudo do microbioma oral tem sido possível graças a novas técnicas moleculares, que ajudaram a ultrapassar certas limitações de cultivo de determinas espécies bacterianas. O estudo do microbioma, das interações entre as comunidades microbianas e a sua relação com o hospedeiro são a chave para a prevenção de certas doenças orais infeciosas como a cárie dentária e a doença periodontal.

Protein markers of B. xylophilus Steiner & Buhrer, 1934 (Nickle, 1970) populations using quantitative proteomics and character compatibility.

The Pine Wood Nematode (PWN) Bursaphelenchus xylophilus is a severe forest pathogen in countries where it has been introduced and is considered a worldwide quarantine organism. In this study, protein markers for differentiating populations of this nematode were identified by studying differences among four selected Iberian and one American population. These populations were compared by quantitative proteomics (iTRAQ). From a total of 2860 proteins identified using the public database from the B. xylophilus genome project, 216 were unambiguous and significantly differentially regulated in the studied populations. Comparisons of their pairwise ratio were statistically treated and supported in order to convert them into discrete character states, suggesting that 141 proteins were not informative as population specific markers. Application of the Character Compatibility methodology on the remaining 75 proteins (belonging to families with different biological functions) excludes 27 which are incompatible among them. Considering only the compatible proteins, the method selects a subset of 30 specific unique protein markers which allowed the compared classification of the Iberian isolates. This approach makes it easier search for diagnostic tools and phylogenetic inference within species and populations of a pathogen exhibiting a high level of genetic diversity.

Advanced neuroimaging methodologies to improve connectivity detection in normal and abnormal language brain networks

The language connectome was in-vivo investigated using multimodal non-invasive quantitative MRI. In PPA patients (n=18) recruited by the IRCCS ISNB, Bologna, cortical thickness measures showed a predominant reduction on the left hemisphere (p<0.005) with respect to matched healthy controls (HC) (n=18), and an accuracy of 86.1% in discrimination from Alzheimer’s disease patients (n=18). The left temporal and para-hippocampal gyri significantly correlated (p<0.01) with language fluency. In PPA patients (n=31) recruited by the Northwestern University Chicago, DTI measures were longitudinally evaluated (2-years follow-up) under the supervision of Prof. M. Catani, King’s College London. Significant differences with matched HC (n=27) were found, tract-localized at baseline and widespread in the follow-up. Language assessment scores correlated with arcuate (AF) and uncinate (UF) fasciculi DTI measures. In left-ischemic stroke patients (n=16) recruited by the NatBrainLab, King’s College London, language recovery was longitudinally evaluated (6-months follow-up). Using arterial spin labelling imaging a significant correlation (p<0.01) between language recovery and cerebral blood flow asymmetry, was found in the middle cerebral artery perfusion, towards the right. In HC (n=29) recruited by the DIBINEM Functional MR Unit, University of Bologna, an along-tract algorithm was developed suitable for different tractography methods, using the Laplacian operator. A higher left superior temporal gyrus and precentral operculum AF connectivity was found (Talozzi L et al., 2018), and lateralized UF projections towards the left dorsal orbital cortex. In HC (n=50) recruited in the Human Connectome Project, a new tractography-driven approach was developed for left association fibres, using a principal component analysis. The first component discriminated cortical areas typically connected by the AF, suggesting a good discrimination of cortical areas sharing a similar connectivity pattern. The evaluation of morphological, microstructural and metabolic measures could be used as in-vivo biomarkers to monitor language impairment related to neurodegeneration or as surrogate of cognitive rehabilitation/interventional treatment efficacy.

Harnessing the power of Drosophila model to obtain new insights on MYC and ABCC1 cellular functions: from neuroblastoma to autism spectrum disorder

MYCN amplification is a genetic hallmark of the childhood tumour neuroblastoma. MYCN-MAX dimers activate the expression of genes promoting cell proliferation. Moreover, MYCN seems to transcriptionally repress cell differentiation even in absence of MAX. We adopted the Drosophila eye as model to investigate the effect of high MYC to MAX expression ratio on cells. We found that dMyc overexpression in eye cell precursors inhibits cell differentiation and induces the ectopic expression of Antennapedia (the wing Hox gene). The further increase of MYC/MAX ratio results in an eye-to-wing homeotic transformation. Notably, dMyc overexpression phenotype is suppressed by low levels of transcriptional co-repressors and MYCN associates to the promoter of Deformed (the eye Hox gene) in proximity to repressive sites. Hence, we envisage that, in presence of high MYC/MAX ratio, the “free MYC” might inhibit Deformed expression, leading in turn to the ectopic expression of Antennapedia. This suggests that MYCN might reinforce its oncogenic role by affecting the physiological homeotic program. Furthermore, poor neuroblastoma outcome associates with a high level of the MRP1 protein, encoded by the ABCC1 gene and known to promote drug efflux in cancer cells. Intriguingly, this correlation persists regardless of chemotherapy and ABCC1 overexpression enhances neuroblastoma cell motility. We found that Drosophila dMRP contributes to the adhesion between the dorsal and ventral epithelia of the wing by inhibiting the function of integrin receptors, well known regulators of cell adhesion and migration. Besides, integrins play a crucial role during synaptogenesis and ABCC1 locus is included in a copy number variable region of the human genome (16p13.11) involved in neuropsychiatric diseases. Interestingly, we found that the altered dMRP/MRP1 level affects nervous system development in Drosophila embryos. These preliminary findings point out novel ABCC1 functions possibly defining ABCC1 contribution to neuroblastoma and to the pathogenicity of 16p13.11 deletion/duplication