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.

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

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.

Centennial variations in sunspot number, open solar flux and streamer belt width: 1. correction of the sunspot number record since 1874

We analyse the widely-used international/ Zürich sunspot number record, R, with a view to quantifying a suspected calibration discontinuity around 1945 (which has been termed the “Waldmeier discontinuity” [Svalgaard, 2011]). We compare R against the composite sunspot group data from the Royal Greenwich Observatory (RGO) network and the Solar Optical Observing Network (SOON), using both the number of sunspot groups, N{sub}G{\sub}, and the total area of the sunspots, A{sub}G{\sub}. In addition, we compare R with the recently developed interdiurnal variability geomagnetic indices IDV and IDV(1d). In all four cases, linearity of the relationship with R is not assumed and care is taken to ensure that the relationship of each with R is the same before and after the putative calibration change. It is shown the probability that a correction is not needed is of order 10{sup}−8{\sup} and that R is indeed too low before 1945. The optimum correction to R for values before 1945 is found to be 11.6%, 11.7%, 10.3% and 7.9% using A{sub}G{\sub}, N{sub)G{\sub}, IDV, and IDV(1d), respectively. The optimum value obtained by combining the sunspot group data is 11.6% with an uncertainty range 8.1-14.8% at the 2σ level. The geomagnetic indices provide an independent yet less stringent test but do give values that fall within the 2σ uncertainty band with optimum values are slightly lower than from the sunspot group data. The probability of the correction needed being as large as 20%, as advocated by Svalgaard [2011], is shown to be 1.6 × 10{sup}−5{\sup}.

Centennial variations in sunspot number, open solar flux and streamer belt width: 2. Comparison with the geomagnetic data

We investigate the relationship between interdiurnal variation geomagnetic activity indices, IDV and IDV(1d), corrected sunspot number, R{sub}C{\sub}, and the group sunspot number R{sub}G{\sub}. R{sub}C{\sub} uses corrections for both the “Waldmeier discontinuity”, as derived in Paper 1 [Lockwood et al., 2014c], and the “Wolf discontinuity” revealed by Leussu et al. [2013]. We show that the simple correlation of the geomagnetic indices with R{sub}C{\sub}{sup}n{\sup} or R{sub}G{\sub}{sup}n{\sup} masks a considerable solar cycle variation. Using IDV(1d) or IDV to predict or evaluate the sunspot numbers, the errors are almost halved by allowing for the fact that the relationship varies over the solar cycle. The results indicate that differences between R{sub}C{\sub} and R{sub}G{\sub} have a variety of causes and are highly unlikely to be attributable to errors in either R{sub}G{\sub} alone, as has recently been assumed. Because it is not known if R{sub}C{\sub} or R{sub}G{\sub} is a better predictor of open flux emergence before 1874, a simple sunspot number composite is suggested which, like R{sub}G{\sub}, enables modelling of the open solar flux for 1610 onwards in Paper 3, but maintains the characteristics of R{sub}C{\sub}.

Centennial variations in sunspot number, open solar flux and streamer belt width: 3. Modelling

From the variation of near-Earth interplanetary conditions, reconstructed for the mid-19th century to the present day using historic geomagnetic activity observations, Lockwood and Owens [2014] have suggested that Earth remains within a broadened streamer belt during solar cycles when the Open Solar Flux (OSF) is low. From this they propose that the Earth was immersed in almost constant slow solar wind during the Maunder minimum (c. 1650-1710). In this paper, we extend continuity modelling of the OSF to predict the streamer belt width using both group sunspot numbers and corrected international sunspot numbers to quantify the emergence rate of new OSF. The results support the idea that the solar wind at Earth was persistently slow during the Maunder minimum because the streamer belt was broad.