Understanding the genetic basis of phenotype variability in individuals with neurocognitive disorders

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

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.

Neuropsychiatric disease (NPD) clustering in families with Autism Spectrum Disorder (ASD): genetic, epigenetic and environmental issues

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder characterized by deficits in social communication/interaction and by unusual repetitive and restricted behaviors and interests. ASD often co-occurs in the same families with other neuropsychiatric diseases (NPD), such as intellectual disability, schizophrenia, epilepsy, depression and attention deficit hyperactivity disorder. Genetic factors have an important role in ASD etiology. Multiple copy number variants (CNVs) and single nucleotide variants (SNVs) in candidate genes have been associated with an increased risk to develop ASD. Nevertheless, recent heritability estimates and the high genotypic and phenotypic heterogeneity characteristic of ASD indicate a role of environmental and epigenetic factors, such as long noncoding RNA (lncRNA) and microRNA (miRNA), as modulators of genetic expression and further clinical presentation. Both miRNA and lncRNA are functional RNA molecules that are transcribed from DNA but not translated into proteins, instead they act as powerful regulators of gene expression. While miRNA are small noncoding RNAs with 22-25 nucleotides in length that act at the post-transcriptional level of gene expression, the lncRNA are bigger molecules (>200 nucleotides in length) that are capped, spliced, and polyadenylated, similar to messenger RNA. Although few lncRNA were well characterized until date, there is a great evidence that they are implicated in several levels of gene expression (transcription/post-transcription/post-translation, organization of protein complexes, cell– cell signaling as well as recombination) as shown in figure 1.

Data supporting the co-expression of PDHA1 gene and of its paralogue PDHA2 in somatic cells of a family

This article presents a dataset proving the simultaneous presence of a 5′UTR-truncated PDHA1 mRNA and a full-length PDHA2 mRNA in the somatic cells of a PDC-deficient female patient and all members of her immediate family (parents and brother). We have designed a large set of primer pairs in order to perform detailed RT-PCR assays allowing the clear identification of both PDHA1 and PDHA2 mRNA species in somatic cells. In addition, two different experimental approaches were used to elucidate the copy number of PDHA1 gene in the patient and her mother. The interpretation and discussion of these data, along with further extensive experiments concerning the origin of this altered gene expression and its potential therapeutic consequences, can be found in “Complex genetic findings in a female patient with pyruvate dehydrogenase complex deficiency: null mutations in the PDHX gene associated with unusual expression of the testis-specific PDHA2 gene in her somatic cells” (A. Pinheiro, M.J. Silva, C. Florindo, et al., 2016).

The clinical significance of 15q11.2 BP1-BP2 duplications: - Where do we stand?

The 15q11.2-q13 region has been well characterized, being associated with a range of syndromatic copy number variants (CNVs), and comprises five established break points sites (BP1 to BP5). While the clinical effect for BP1-BP3, BP2-BP3 and BP4-BP5 CNVs is well established, the same cannot be said for BP1-BP2 CNVs. Recently the 15q11.2 BP1-BP2 deletion has been reviewed, emerging as a microdeletion syndrome with low penetrance and variable expressivity being the CNV frequently inherited from a healthy parent. This microdeletion is considered to be a risk factor for several neurodevelopment disorders. For the reciprocal duplication the picture has been less conclusive. Aiming for a better understanding of the clinical significance of this CNV, we collected patients with intellectual disability and/or other clinical features, referred for microarray testing, gathering clinical details for the ones with the duplication. Data was collected from two genetic laboratories. With a total of 1545 patients, we identified eleven carrying the duplication at 15q11.2 BP1-BP2. It was possible to assess inheritance in only four cases, all inherited from a healthy parent. All patients presented intellectual disability,and facial dysmorphism was the second most common feature observed. Microcephaly, autism, congenital abnormalities, dystonia and cataplexy where reported individually. The magnitude of the effect of 15q11.2 duplication remains elusive, and the outcome unclear, posing a major challenge to genetic counseling. Nevertheless, we expect the collection of more of these cases will establish this gain, as it happened with the reciprocal deletion, as a microduplication syndrome with low penetrance and variable expressivity.

A 1.77 Mb deletion in 3p26.3 encompassing CNTN6 and CNTN4 genes: case report

Chromosome microarray analysis is a powerful diagnostic tool and is being used as a first-line approach to detect chromosome imbalances associated with intellectual disability, dysmorphic features and congenital abnormalities. This test enables the identification of new copy number variants (CNVs) and their association with new microdeletion/microduplication syndromes in patients previously without diagnosis. We report the case of a 7 year-old female with moderate intellectual disability, severe speech delay and auto and hetero aggressivity with a previous 45,XX,der(13;14)mat karyotype performed at a younger age. Affymetrix CytoScan 750K chromosome microarray analysis was performed detecting a 1.77 Mb deletion at 3p26.3, encompassing 2 OMIM genes, CNTN6 and CNTN4. These genes play an important role in the formation, maintenance, and plasticity of functional neuronal networks. Deletions or mutations in CNTN4 gene have been implicated in intellectual disability and learning disabilities. Disruptions or deletions in the CNTN6 gene have been associated with development delay and other neurodevelopmental disorders. The haploinsufficiency of these genes has been suggested to participate to the typical clinical features of 3p deletion syndrome. Nevertheless inheritance from a healthy parent has been reported, suggesting incomplete penetrance and variable phenotype for this CNV. We compare our patient with other similar reported cases, adding additional value to the phenotype-genotype correlation of deletions in this region.

From genes to radioresistance in head and neck squamous cell carcinoma

Head and Neck Cancers (HNC) are a group of tumours located in the upper aero-digestive tract. Head and Neck Squamous Cell Carcinoma (HNSCC) represent about 90% of all HNC cases. It has been considered the sixth most malignant tumour worldwide and, despite clinical and technological advances, the five-year survival rate has not improved much in the last years. Nowadays, HNSCC is well established as a heterogeneous disease and that its development is due to accumulation of genetic events. Apart from the majority of the patients being diagnosed in an advanced stage, HNSCC is also a disease with poor therapeutic outcome. One of the therapeutic approaches is radiotherapy. However, this approach has different drawbacks like the radioresistance acquired by some tumour cells, leading to a worse prognosis. A major knowledge in radiation biology is imperative to improve this type of treatment and avoid late toxicities, maintaining patient quality of life in the subsequent years after treatment. Then, identification of genetic markers associated to radiotherapy response in patients and possible alterations in cells after radiotherapy are essential steps towards an improved diagnosis, higher survival rate and a better life quality. Not much is known about the radiation effects on cells, so, the principal aim of this study was to contribute to a more extensive knowledge about radiation treatment in HNSCC. For this, two commercial cell lines, HSC-3 and BICR-10, were used and characterized resorting to karyotyping, aCGH and MS-MLPA. These cell lines were submitted to different doses of irradiation and the resulting genetic and methylation alterations were evaluated. Our results showed a great difference in radiation response between the two cell lines, allowing the conclusion that HSC-3 was much more radiosensitive than BICR-10. Bearing this in mind, analysis of cell death, cell cycle and DNA damages was performed to try to elucidate the motifs behind this difference. The characterization of both cell lines allowed the confirmation that HSC-3 was derived from a metastatic tumour and the hypothesis that BICR-10 was derived from a dysplasia. Furthermore, this pilot study enabled the suggestion of some genetic and epigenetic alterations that cells suffer after radiation treatment. Additionally, it also allowed the association of some genetic characteristics that could be related to the differences in radiation response observable in this two cell lines. Taken together all of our results contribute to a better understanding of radiation effects on HNSCC allowing one further step towards the prediction of patients’ outcome, better choice of treatment approaches and ultimately a better quality of life.

Exploration génétique de l'hypothyroïdie congénitale par dysgénésie thyroïdienne

L'hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT, ectopie dans plus de 80 %) a une prévalence de 1 cas sur 4000 naissances vivantes. L’HCDT est la conséquence d'une défaillance de la thyroïde embryonnaire à se différencier, à se maintenir ou à migrer vers sa localisation anatomique (partie antérieure du cou), qui aboutit à une absence totale de la thyroïde (athyréose) ou à une ectopie thyroïdienne (linguale ou sublinguale). Les HCDT sont principalement non-syndromiques (soit 98% des cas sont non-familiale), ont un taux de discordance de 92% chez les jumeaux monozygotes, et ont une prédominance féminine et ethnique (i.e., Caucasienne). La majorité des cas d’HCDT n’a pas de cause connue, mais est associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdie). Des mutations germinales dans les facteurs de transcription liés à la thyroïde (NKX2.1, FOXE1, PAX8, NKX2.5) ont été identifiées dans seulement 3% des patients atteints d’HCDT sporadiques et l’analyse de liaisons exclue ces gènes dans les rares familles multiplex avec HCDT. Nous supposons que le manque de transmission familiale claire d’HCDT peut résulter de la nécessité d’au moins deux « hits » génétiques différents dans des gènes importants pour le développement thyroïdien. Pour répondre au mieux nos questions de recherche, nous avons utilisé deux approches différentes: 1) une approche gène candidat, FOXE1, seul gène impliqué dans l’ectopie dans le modèle murin et 2) une approche en utilisant les techniques de séquençage de nouvelle génération (NGS) afin de trouver des variants génétiques pouvant expliquer cette pathologie au sein d’une cohorte de patients avec HCDT. Pour la première approche, une étude cas-contrôles a été réalisée sur le promoteur de FOXE1. Il a récemment été découvert qu’une région du promoteur de FOXE1 est différentiellement méthylée au niveau de deux dinucléotides CpG consécutifs, définissant une zone cruciale de contrôle de l’expression de FOXE1. L’analyse d’association basée sur les haplotypes a révélé qu’un haplotype (Hap1: ACCCCCCdel1C) est associé avec le HCDT chez les Caucasiens (p = 5x10-03). Une réduction significative de l’activité luciférase est observée pour Hap1 (réduction de 68%, p<0.001) comparé au promoteur WT de FOXE1. Une réduction de 50% de l’expression de FOXE1 dans une lignée de cellules thyroïdienne humaine est suffisante pour réduire significativement la migration cellulaire (réduction de 55%, p<0.05). Un autre haplotype (Hap2: ACCCCCCC) est observé moins fréquemment chez les Afro-Américain comparés aux Caucasiens (p = 1.7x10-03) et Hap2 diminue l’activité luciférase (réduction de 26%, p<0.001). Deux haplotypes distincts sont trouvés fréquemment dans les contrôles Africains (Black-African descents). Le premier haplotype (Hap3: GTCCCAAC) est fréquent (30.2%) chez les contrôles Afro-Américains comparés aux contrôles Caucasiens (6.3%; p = 2.59 x 10-9) tandis que le second haplotype (Hap4: GTCCGCAC) est trouvé exclusivement chez les contrôles Afro-Américains (9.4%) et est absent chez les contrôles Caucasiens (P = 2.59 x 10-6). Pour la deuxième approche, le séquençage de l’exome de l’ADN leucocytaire entre les jumeaux MZ discordants n’a révélé aucune différence. D'où l'intérêt du projet de séquençage de l’ADN et l’ARN de thyroïdes ectopiques et orthotopiques dans lesquelles de l'expression monoallélique aléatoire dans a été observée, ce qui pourrait expliquer comment une mutation monoallélique peut avoir des conséquences pathogéniques. Finalement, le séquençage de l’exome d’une cohorte de 36 cas atteints d’HCDT a permis d’identifier de nouveaux variants probablement pathogéniques dans les gènes récurrents RYR3, SSPO, IKBKE et TNXB. Ces quatre gènes sont impliqués dans l’adhésion focale (jouant un rôle dans la migration cellulaire), suggérant un rôle direct dans les défauts de migration de la thyroïde. Les essais de migration montrent une forte diminution (au moins 60% à 5h) de la migration des cellules thyroïdiennes infectées par shRNA comparés au shCtrl dans 2 de ces gènes. Des zebrafish KO (-/- et +/-) pour ces nouveaux gènes seront réalisés afin d’évaluer leur impact sur l’embryologie de la thyroïde.

Drinking water disinfection by-products, genetic polymorphisms, and birth outcomes in a european mother-child cohort study

International audience

Reassessing the Evolutionary History of the 17q21 Inversion Polymorphism

A polymorphic inversion that lies on chromosome 17q21 comprises two major haplotype families (H1 and H2) that not only differ in orientation but also in copy-number. Although the processes driving the spread of the inversion-associated lineage (H2) in humans remain unclear, a selective advantage has been proposed for one of its subtypes. Here, we genotyped a large panel of individuals from previously overlooked populations using a custom array with a unique panel of H2-specific single nucleotide polymorphisms and found a patchy distribution of H2 haplotypes in Africa, with North Africans displaying a higher frequency of inverted subtypes, when compared with Sub-Saharan groups. Interestingly, North African H2s were found to be closer to "non-African" chromosomes further supporting that these populations may have diverged more recently from groups outside Africa. Our results uncovered higher diversity within the H2 family than previously described, weakening the hypothesis of a strong selective sweep on all inverted chromosomes and suggesting a rather complex evolutionary history at this locus.

FROM GENES TO GENOMES: LOCAL ADAPTATION AND ADAPTIVE POTENTIAL IN TWO ARCTIC SEABIRDS

By investigating the mechanisms underlying the evolution and the maintenance of local adaptations we can help predict how species will adapt to future environmental change. In this thesis I investigate local adaptation and adaptive potential in thick-billed and common murres (Uria lomvia and U. aalge), two arctic seabirds of international conservation concern. Thanks to the recent development of new genomic methods, I address three major themes that are relevant for both the development of evolutionary theory and conservation: 1) the role of gene flow in the origin and maintenance of adaptation; 2) levels and distribution of standing genetic variation, and their contribution to adaptive potential; and 3) the genomic mechanisms maintaining an adaptive dimorphism within a single interbreeding population. First, I review the literature on genomics of local adaptation with gene flow and find that adaptation can be maintained despite gene flow, that gene flow itself can promote adaptation, and that genetic architecture is important in the origin and maintenance of local adaptations. Second, I genotype genome-wide markers and toll-like receptor genes (TLRs) to investigate local adaptation and adaptive potential in thick-billed murres. Thick-billed murres do not show signatures of local adaptation to their breeding grounds, but outlier loci group birds according to their non-breeding distributions, suggesting that selection and/or demographic connectivity in the winter may explain patterns of differentiation in this species. Genetic variation at TLRs does not decrease with increasing latitude as predicted, but tests of selection and measures of genetic diversity suggest differences in local selective regimes at most genes. Thick-billed murres show high levels of standing genetic variation and their adaptive potential will mostly depend on rate and magnitude of environmental change. Finally, I improve and annotate the assembly of the highly heterozygous genome of the thick-billed murre. Using this assembly as a reference, I perform whole genome analyses to investigate the genomic basis of an adaptive dimorphism in Atlantic common murres. I show for the first time that a 60 kb complex copy number variant in a non-coding region maintains differences in plumage and cold adaptation despite high gene flow.

Evaluación del desempeño de pruebas diagnósticas de tuberculosis en pacientes VIH positivo en dos hospitales públicos de Bogotá

La tuberculosis TB es una de las principales causas de muerte en el mundo en individuos con infección por VIH. En Colombia esta coinfección soporta una carga importante en la población general convirtiéndose en un problema de salud pública. En estos pacientes las pruebas diagnósticas tienen sensibilidad inferior y la enfermedad evoluciona con mayor frecuencia hacia formas diseminadas y rápidamente progresivas y su diagnóstico oportuno representa un reto en Salud. El objetivo de este proyecto es evaluar el desempeño de las pruebas diagnósticas convencionales y moleculares, para la detección de TB latente y activa pacientes con VIH, en dos hospitales públicos de Bogotá. Para TB latente se evaluó la concordancia entre las pruebas QuantiFERON-TB (QTF) y Tuberculina (PPD), sugiriendo superioridad del QTF sobre la PPD. Se evaluaron tres pruebas diagnósticas por su sensibilidad y especificidad, baciloscopia (BK), GenoType®MTBDR plus (Genotype) y PCR IS6110 teniendo como estándar de oro el cultivo. Los resultados de sensibilidad (S) y especificidad (E) de cada prueba con una prevalencia del 19,4 % de TB pulmonar y extrapulmonar en los pacientes que participaron del estudio fue: BK S: 64% E: 99,1%; Genotype S: 77,8% E: 94,5%; PCRIS6110 S: 73% E: 95,5%, de la misma forma se determinaron los valores predictivos positivos y negativos (VPP y VPN) BK: 88,9% y 94,8%, Genotype S: 77,8% E: 94,5%; PCRIS6110 S: 90% y 95,7%. Se concluyó bajo análisis de curva ROC que las pruebas muestran un rendimiento diagnóstico similar por separado en el diagnóstico de TB en pacientes con VIH, aumentando su rendimiento diagnostico cuando se combinan

Registro de pacientes con distrofinopatías en Colombia

INTRODUCCIÓN. La distrofia muscular de Duchenne es una enfermedad neuromuscular con una herencia recesiva ligada al X que afecta a 1 de cada 3500 niños nacidos vivos. Se produce por mutaciones en el gen DMD que codifica para la distrofina. Se caracteriza por manifestaciones clínicas variables típicas de una distrofia muscular proximal progresiva. OBJETIVO. Realizar el primer registro en Colombia de los pacientes identificados con distrofinopatías, teniendo en cuenta características clínicas y paraclínicas, así como las mutaciones causales de esta patología. METODOLOGÍA Es un estudio descriptivo, transversal, de la revisión de historias clínicas de los pacientes con diagnóstico de DMD atendidos en la consulta de Genética de la Universidad del Rosario durante los años 2006 a 2015. RESULTADOS Se identificaron 99 pacientes, de los cuales 56 (56,56%) corresponden al fenotipo Duchenne y 12 (12,12%) al Becker. No fue posible clasificar a 31 pacientes (31,3%) por falta de datos clínicos. La edad de inicio de los síntomas fue en promedio de 4,41 años. Las mutaciones más frecuentes fueron las deleciones (69%), seguidas por las mutaciones puntuales(14%), las duplicaciones (11%) y por otras mutaciones (4%). CONCLUSIONES Este registro de distrofinopatías es el primero reportado en Colombia y el punto de partida para conocer la incidencia de la enfermedad, caracterización clínica y molecular de los pacientes, garantizando así el acceso oportuno a los nuevos tratamientos de medicina de precisión que permitan mejorar la calidad de vida de los pacientes y sus familias.

Development of molecular approaches to assess the diversity and abundance of marine nematodes assemblages for assessment of Good Environmental Status of the estuarine and marine ecosystems

In Europe, the concerns with the status of marine ecosystems have increased, and the Marine Directive has as main goal the achievement of Good Environmental Status (GES) of EU marine waters by 2020. Molecular tools are seen as promising and emerging approaches to improve ecosystem monitoring, and have led ecology into a new era, representing perhaps the most source of innovation in marine monitoring techniques. Benthic nematodes are considered ideal organisms to be used as biological indicator of natural and anthropogenic disturbances in aquatic ecosystems underpinning monitoring programmes on the ecological quality of marine ecosystems, very useful to assess the GES of the marine environment. dT-RFLP (directed Terminal-Restriction Fragment Length Polymorphism) allows to assess the diversity of nematode communities, but also allows studying the functioning of the ecosystem, and combined with relative real-time PCR (qPCR), provides a high-throughput semi-quantitative characterization of nematode communities. These characteristics make the two molecular tools good descriptors for the good environmental status assessment. The main aim of this study is to develop and optimize the dT-RFLP and qPCR in Mira estuary (SW coast, Portugal). A molecular phylogenetic analysis of marine and estuarine nematodes is being performed combining morphological and molecular analysis to evaluate the diversity of free-living marine nematodes in Mira estuary. After morphological identification, barcoding of 18S rDNA and COI genes are being determined for each nematode species morphologically identified. So far we generated 40 new sequences belonging to 32 different genus and 17 families, and the study has shown a good degree of concordance between traditional morphology-based identification and DNA sequences. These results will improve the assessment of marine nematode diversity and contribute to a more robust nematode taxonomy. The DNA sequences are being used to develop the dT-RFLP with the ability to easily process large sample numbers (hundreds and thousands), rather than typical of classical taxonomic or low throughput molecular analyses. A preliminary study showed that the digest enzymes used in dT-RFLP for terrestrial assemblages separated poorly the marine nematodes at taxonomic level for functional group analysis. A new digest combination was designed using the software tool DRAT (Directed Terminal Restriction Analysis Tool) to distinguished marine nematode taxa. Several solutions were provided by DRAT and tested empirically to select the solution that cuts most efficiently. A combination of three enzymes and a single digest showed to be the best solution to separate the different clusters. Parallel to this, another tool is being developed to estimate the population size (qPCR). An improvement in qPCR estimation of gene copy number using an artificial reference is being performed for marine nematodes communities to quantify the abundance. Once developed, it is proposed to validate both methodologies by determining the spatial and temporal variability of benthic nematodes assemblages across different environments. The application of these high-throughput molecular approaches for benthic nematodes will improve sample throughput and their implementation more efficient and faster as indicator of ecological status of marine ecosystems.

Genome-wide detection of CNVs and their association with meat tenderness in Nelore cattle.

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina1 High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing 6.5% of the genome. Comparing our results with previous studies revealed an overlap in 1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.