Evidence for time dependency of molecular rate estimates

Long-term changes in the genetic composition of a population occur by the fixation of new mutations, a process known as substitution. The rate at which mutations arise in a population and the rate at which they are fixed are expected to be equal under neutral conditions (Kimura, 1968). Between the appearance of a new mutation and its eventual fate of fixation or loss, there will be a period in which it exists as a transient polymorphism in the population (Kimura and Ohta, 1971). If the majority of mutations are deleterious (and nonlethal), the fixation probabilities of these transient polymorphisms are reduced and the mutation rate will exceed the substitution rate (Kimura, 1983). Consequently, different apparent rates may be observed on different time scales of the molecular evolutionary process (Penny, 2005; Penny and Holmes, 2001). The substitution rate of the mitochondrial protein-coding genes of birds and mammals has been traditionally recognized to be about 0.01 substitutions/site/million years (Myr) (Brown et al., 1979; Ho, 2007; Irwin et al., 1991; Shields and Wilson, 1987), with the noncoding D-loop evolving several times more quickly (e.g., Pesole et al., 1992; Quinn, 1992). Over the past decade, there has been mounting evidence that instantaneous mutation rates substantially exceed substitution rates, in a range of organisms (e.g., Denver et al., 2000; Howell et al., 2003; Lambert et al., 2002; Mao et al., 2006; Mumm et al., 1997; Parsons et al., 1997; Santos et al., 2005). The immediate reaction to the first of these findings was that the polymorphisms generated by the elevated mutation rate are short-lived, perhaps extending back only a few hundred years (Gibbons, 1998; Macaulay et al., 1997). That is, purifying selection was thought to remove these polymorphisms very rapidly.

Tuberous sclerosis complex: clinical features, diagnosis, and prevalence within Northern Ireland

Tuberous sclerosis complex (TSC) is an autosomal dominant disorder caused by mutations in the TSC1 and TSC2 genes on chromosomes 9 and 16 respectively. Diagnosis is based on clinical features but can be difficult as a result of variable phenotypic expression. With the advantage of mutation analysis in making a diagnosis of TSC, and improved identification of the associated clinical features, there have been few new data on its prevalence and on the proportion of cases due to new mutations. We have performed a retrospective epidemiological study on the prevalence of TSC, the clinical features attributed to it, and the availability of mutational analysis. We identified 73 known patients with TSC (5 deceased): 39 were female and 34 male. Ages ranged from 10 months to 69 years, with a mean age of 27 years 11 months (SD 16y 10mo). The point prevalence of TSC in our study was estimated at I out of 24 956 on the prevalence day (30 April 2004). The majority of patients (42.5%) were diagnosed at less than 15 months of age; 25% were not given a diagnosis on first developing symptoms. In all, 93.2% had epilepsy and 71.2% had a learning disability.* A mutation was identified in 95.8% of those tested (26% TSC1 and 74% TSC2). TSC2 mutations were correlated with a more severe phenotype. The new mutation rate was calculated at 64%. We conclude that the prevalence of TSC is higher than previously calculated. We recommend that all children with epilepsy be assessed for features of TSC. Larger studies will be required to assess the prevalence of mutations in each gene, and genotype-phenotype correlation.

Pathogenèse moléculaire de la neuropathie sensitive et motrice héréditaire avec agénésie du corps calleux

La neuropathie sensitive et motrice héréditaire avec agénésie du corps calleux (NSMH/ACC) se traduit par une atteinte neurodégénérative sévère associée à des anomalies développementales dans le système nerveux central et du retard mental. Bien que rare dans le monde, ce désordre autosomique récessif est particulièrement fréquent dans la population Québécoise du Canada Français du fait d’un effet fondateur. L’unique étude réalisée sur la mutation québécoise du gène qui code pour le co-transporteur de potassiumchlore 3 (KCC3) a montré qu’il y a une perte de fonction de la protéine. Cependant, la maladie est également retrouvée hors du Québec et il reste encore à élucider les pathomécanismes mis en jeu. Nous avons donc séquencé les 26 exons du gène KCC3 chez des individus recrutés dans le monde entier et suspectés d’être atteints de la maladie. Nous avons ainsi identifié trois nouvelles mutations. L’étude fonctionnelle de ces mutations nous a confirmé la perte de fonction systématique des co-transporteurs mutés. Puisque l’inactivation de KCC3 se produit majoritairement via l’élimination de segments peptidiques en C-terminus, nous avons concentré notre attention sur l’identification des interactions qui s’y produisent. À l’aide d’approches double hybride, pull-down et immunomarquage, nous avons déterminé que KCC3 interagit avec la créatine kinase CK-B et que cette interaction est perturbée par les mutations tronquantes. De plus, l’utilisation d’un inhibiteur de créatine kinase inactive KCC3, ce qui démontre qu’il existe bien un lien fonctionnel et pathologique entre KCC3 et ses partenaires C-terminaux. Nous avons aussi identifié des anomalies majeures de localisation membranaire des KCC3 mutés. Que KCC3 soit tronqué ou pleine longueur, sa distribution subcellulaire est affectée dans des cellules en culture, dans les ovocytes de Xenopes et dans des échantillons de cerveau de patients. La perte d’interaction entre KCC3 et CK-B et/ou les défauts de transit intracellulaire de KCC3 sont donc les mécanismes pathologiques majeurs de la NSMH/ACC.

Identification des bases génétiques des myopathies à multi-minicores avec ou sans cardiomyopathie

Thèse réalisée en cotutelle avec l'Université Pierre et Marie Curie, Paris 6(UPMC, Paris, France).

Études génétiques de familles récessives d’ataxies et de paraplégies spastiques

Au cours des dernières années, la génétique a subi une progression phénoménale suite au développement de nouvelles technologies de séquençage. En effet, le séquençage de l’exome entier chez des familles a permis l’identification de nouveaux gènes impliqués pour plusieurs maladies. La neurologie a d’ailleurs bénéficié de ces avancées et plusieurs gènes ont été mis en évidence comme causatifs pour différents désordres neurologiques. Dans ce travail il sera question de deux désordres du mouvement pour lequel nous avons utilisés des technologies de séquençage traditionnelles, en l’occurrence le séquençage par Sanger, ainsi que de nouvelles technologies pour le séquençage de l’exome entier afin d’identifier de nouveaux gènes causatifs. Le premier désordre du mouvement qui sera décrit est l’ataxie, où ne seront abordées que les ataxies de cause génétiques, à transmission récessive. Le premier chapitre relatera les nouvelles mutations qui ont été trouvées chez des canadiens-français souffrant de l’ataxie de Beauce. Il sera aussi question de nouvelles mutations retrouvées dans deux autres populations, confirmant l’implication du gène SYNE1 dans les cas d’ataxie cérébelleuse à travers le monde. Le second chapitre fera la démonstration qu’il est souhaitable d’utiliser le séquençage de l’exome entier dans le but de poser un diagnostic clinique. En effet, il a été possible de trouver la cause génétique d’une famille comportant deux membres atteints d’atrophie congénitale du cervelet, où le symptôme prédominant est l’ataxie. Le séquençage de l’exome a permis la mise en évidence de mutations dans le gène PMM2, déjà connues pour cause le syndrome des glycoprotéines déficientes en hydrates de carbone. Dans un second temps, il sera question d’un autre désordre du mouvement la paraplégie spastique familiale (PSF). Le chapitre 3 relatera les mutations trouvées dans le gène CYP7B1 dans notre cohorte de patients PSF.

Polimorfismos del gen Butirilcolinesterasa responsables de reacciones adversas en pacientes consumidores de “cocaína”

La butirilcolinesterasa humana (BChE; EC 3.1.1.8) es una enzima polimórfica sintetizada en el hígado y en el tejido adiposo, ampliamente distribuida en el organismo y encargada de hidrolizar algunos ésteres de colina como la procaína, ésteres alifáticos como el ácido acetilsalicílico, fármacos como la metilprednisolona, el mivacurium y la succinilcolina y drogas de uso y/o abuso como la heroína y la cocaína. Es codificada por el gen BCHE (OMIM 177400), habiéndose identificado más de 100 variantes, algunas no estudiadas plenamente, además de la forma más frecuente, llamada usual o silvestre. Diferentes polimorfismos del gen BCHE se han relacionado con la síntesis de enzimas con niveles variados de actividad catalítica. Las bases moleculares de algunas de esas variantes genéticas han sido reportadas, entre las que se encuentra las variantes Atípica (A), fluoruro-resistente del tipo 1 y 2 (F-1 y F-2), silente (S), Kalow (K), James (J) y Hammersmith (H). En este estudio, en un grupo de pacientes se aplicó el instrumento validado Lifetime Severity Index for Cocaine Use Disorder (LSI-C) para evaluar la gravedad del consumo de “cocaína” a lo largo de la vida. Además, se determinaron Polimorfismos de Nucleótido Simple (SNPs) en el gen BCHE conocidos como responsables de reacciones adversas en pacientes consumidores de “cocaína” mediante secuenciación del gen y se predijo el efecto delos SNPs sobre la función y la estructura de la proteína, mediante el uso de herramientas bio-informáticas. El instrumento LSI-C ofreció resultados en cuatro dimensiones: consumo a lo largo de la vida, consumo reciente, dependencia psicológica e intento de abandono del consumo. Los estudios de análisis molecular permitieron observar dos SNPs codificantes (cSNPs) no sinónimos en el 27.3% de la muestra, c.293A>G (p.Asp98Gly) y c.1699G>A (p.Ala567Thr), localizados en los exones 2 y 4, que corresponden, desde el punto de vista funcional, a la variante Atípica (A) [dbSNP: rs1799807] y a la variante Kalow (K) [dbSNP: rs1803274] de la enzima BChE, respectivamente. Los estudios de predicción In silico establecieron para el SNP p.Asp98Gly un carácter patogénico, mientras que para el SNP p.Ala567Thr, mostraron un comportamiento neutro. El análisis de los resultados permite proponer la existencia de una relación entre polimorfismos o variantes genéticas responsables de una baja actividad catalítica y/o baja concentración plasmática de la enzima BChE y algunas de las reacciones adversas ocurridas en pacientes consumidores de cocaína.

Mutational meltdown in primary endosymbionts: selection limits Muller's ratchet

Background Primary bacterial endosymbionts of insects (p-endosymbionts) are thought to be undergoing the process of Muller's ratchet where they accrue slightly deleterious mutations due to genetic drift in small populations with negligible recombination rates. If this process were to go unchecked over time, theory predicts mutational meltdown and eventual extinction. Although genome degradation is common among p-endosymbionts, we do not observe widespread p-endosymbiont extinction, suggesting that Muller's ratchet may be slowed or even stopped over time. For example, selection may act to slow the effects of Muller's ratchet by removing slightly deleterious mutations before they go to fixation thereby causing a decrease in nucleotide substitutions rates in older p-endosymbiont lineages. Methodology/Principal Findings To determine whether selection is slowing the effects of Muller's ratchet, we determined the age of the Candidatus Riesia/sucking louse assemblage and analyzed the nucleotide substitution rates of several p-endosymbiont lineages that differ in the length of time that they have been associated with their insect hosts. We find that Riesia is the youngest p-endosymbiont known to date, and has been associated with its louse hosts for only 13–25 My. Further, it is the fastest evolving p-endosymbiont with substitution rates of 19–34% per 50 My. When comparing Riesia to other insect p-endosymbionts, we find that nucleotide substitution rates decrease dramatically as the age of endosymbiosis increases. Conclusions/Significance A decrease in nucleotide substitution rates over time suggests that selection may be limiting the effects of Muller's ratchet by removing individuals with the highest mutational loads and decreasing the rate at which new mutations become fixed. This countering effect of selection could slow the overall rate of endosymbiont extinction.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Detecção de cinco novas mutações em pacientes brasileiros com gangliosidose GM1

A Gangliosidose GM1 é um Erro Inato do Metabolismo (EIM) causado pela deficiência da enzima B-galactosidase ácida. Essa doença é caracterizada pelo acúmulo de metabólitos não degradados, principalmente gangliosídeo GM1, nos lisossomos de vários tipos celulares. Baseado na idade de início e na atividade residual da enzima, a Gangliosidose GM1 é classificada em três diferentes tipos: infantil, juvenil e adulto. O gene da B-galactosidase ácida (GLB1, GeneBank M27507) está situado no cromossomo 3 e possui mais de 60 kb, contendo 16 exons. Cerca de 50 mutações associadas à doença estão descritas na literatura. No sul do Brasil, há uma alta freqüência dessa doença (1:17.000 nascidos vivos). Neste trabalho, vinte pacientes diagnosticados no Hospital de Clínicas de Porto Alegre (Brasil) tiveram o gene GLB1 investigado por SSCP (Single Strand Conformational Polymorphism) usando DNA extraído de sangue periférico. Através desta triagem foram encontradas 52 alterações de mobilidade do DNA, indicando a presença de mutações. As amostras relativas aos exons 2 e 15 foram submetidas a sequenciamento direto com seqüenciador ABI31O(Applied Biosystens) utilizado kit BigDye 3.1. Cinco novas mutações no gene GLB1 (F63Y, R38G, Y36S, Y64F e R59C) e duas mutações já descritas (R59H e 1622-1627insG) foram encontradas. Este trabalho possibilitou a genotipagem completa de 6 pacientes e parcial de 5, e direcionou a investigação de mutações, contribuindo diretamente no diagnóstico da enfermidade e permitindo a realização de estudos de correlação genótipo/fenótipo destes pacientes.

MODY 2: Mutation identification and molecular ancestry in a Brazilian family

Maturity Onset Diabetes of the Young (MODY) is a heterogeneous group of genetic diseases characterized by a primary defect in insulin secretion and hyperglycemia, non-ketotic disease, monogenic autosomal dominant mode of inheritance, age at onset less than 25. years, and lack of auto-antibodies. It accounts for 2-5% of all cases of non-type 1 diabetes. MODY subtype 2 is caused by mutations in the glucokinase (GCK) gene. In this study, we sequenced the GCK gene of two volunteers with clinical diagnosis for MODY2 and we were able to identify four mutations including one for a premature stop codon (c.76C>T). Based on these results, we have developed a specific PCR-RFLP assay to detect this mutation and tested 122 related volunteers from the same family. This mutation in the GCK gene was detected in 21 additional subjects who also had the clinical features of this genetic disease. In conclusion, we identified new GCK gene mutations in a Brazilian family of Italian descendance, with one due to a premature stop codon located in the second exon of the gene. We also developed a specific assay that is fast, cheap and reliable to detect this mutation. Finally, we built a molecular ancestry model based on our results for the migration of individuals carrying this genetic mutation from Northern Italy to Brazil. © 2012 Elsevier B.V.

Mutation Spectrum in the Large GTPase Dynamin 2, and Genotype-Phenotype Correlation in Autosomal Dominant Centronuclear Myopathy

Centronuclear myopathy (CNM) is a genetically heterogeneous disorder associated with general skeletal muscle weakness, type I fiber predominance and atrophy, and abnormally centralized nuclei. Autosomal dominant CNM is due to mutations in the large GTPase dynamin 2 (DNM2), a mechanochemical enzyme regulating cytoskeleton and membrane trafficking in cells. To date, 40 families with CNM-related DNM2 mutations have been described, and here we report 60 additional families encompassing a broad genotypic and phenotypic spectrum. In total, 18 different mutations are reported in 100 families and our cohort harbors nine known and four new mutations, including the first splice-site mutation. Genotype-phenotype correlation hypotheses are drawn from the published and new data, and allow an efficient screening strategy for molecular diagnosis. In addition to CNM, dissimilar DNM2 mutations are associated with Charcot-Marie-Tooth (CMT) peripheral neuropathy (CMTD1B and CMT2M), suggesting a tissue-specific impact of the mutations. In this study, we discuss the possible clinical overlap of CNM and CMT, and the biological significance of the respective mutations based on the known functions of dynamin 2 and its protein structure. Defects in membrane trafficking due to DNM2 mutations potentially represent a common pathological mechanism in CNM and CMT. Hum Mutat 33: 949-959, 2012. (C) 2012 Wiley Periodicals, Inc.

Long-term follow-up of patients with Bartter syndrome type I and II

BACKGROUND: Little information is available on a long-term follow-up in Bartter syndrome type I and II. METHODS: Clinical presentation, treatment and long-term follow-up (5.0-21, median 11 years) were evaluated in 15 Italian patients with homozygous (n = 7) or compound heterozygous (n = 8) mutations in the SLC12A1 (n = 10) or KCNJ1 (n = 5) genes. RESULTS: Thirteen new mutations were identified. The 15 children were born pre-term with a normal for gestational age body weight. Medical treatment at the last follow-up control included supplementation with potassium in 13, non-steroidal anti-inflammatory agents in 12 and gastroprotective drugs in five patients. At last follow-up, body weight and height were within normal ranges in the patients. Glomerular filtration rate was <90 mL/min/1.73 m(2) in four patients (one of them with a pathologically increased urinary protein excretion). In three patients, abdominal ultrasound detected gallstones. The group of patients with antenatal Bartter syndrome had a lower renin ratio (P < 0.05) and a higher standard deviation score (SDS) for height (P < 0.05) than a previously studied group of patients with classical Bartter syndrome. CONCLUSIONS: Patients with Bartter syndrome type I and II tend to present a satisfactory prognosis after a median follow-up of more than 10 years. Gallstones might represent a new complication of antenatal Bartter syndrome.

Origins of Shared Genetic Variation in African Cichlids

Cichlid fishes have evolved tremendous morphological and behavioral diversity in the waters of East Africa. Within each of the Great Lakes Tanganyika, Malawi, and Victoria, the phenomena of hybridization and retention of ancestral polymorphism explain allele sharing across species. Here, we explore the sharing of single nucleotide polymorphisms (SNPs) between the major East African cichlid assemblages. A set of approximately 200 genic and nongenic SNPs was ascertained in five Lake Malawi species and genotyped in a diverse collection of 160 species from across Africa. We observed segregating polymorphism outside of the Malawi lineage for more than 50% of these loci; this holds similarly for genic versus nongenic SNPs, as well as for SNPs at putative CpG versus non-CpG sites. Bayesian and principal component analyses of genetic structure in the data demonstrate that the Lake Malawi endemic flock is not monophyletic and that river species have likely contributed significantly to Malawi genomes. Coalescent simulations support the hypothesis that river cichlids have transported polymorphism between lake assemblages. We observed strong genetic differentiation between Malawi lineages for approximately 8% of loci, with contributions from both genic and nongenic SNPs. Notably, more than half of these outlier loci between Malawi groups are polymorphic outside of the lake. Cichlid fishes have evolved diversity in Lake Malawi as new mutations combined with standing genetic variation shared across East Africa.

Evolution of the spermadhesin gene family

Spermadhesins belong to a novel family of secretory proteins of the male genital tract. They are major proteins of the seminal plasma and have been found peripherally associated to the sperm surface. So far, they have only been detected in ungulates, specifically in pig, cattle, and horse, respectively. Spermadhesins form a subgroup of the superfamily of proteins with a CUB-domain that has been found in a variety of developmentally regulated proteins. The structure and function of the spermadhesins have been investigated in the pig. They are multifunctional proteins showing a range of ligand-binding abilities, e.g. to carbohydrates, phospholipids, and protease inhibitors, suggesting that they may be involved in different steps of fertilization. We report here the genomic organization of the porcine spermadhesin gene cluster as well as a detailed comparative analysis with respect to other mammalian species. The porcine spermadhesin genes are located on SSC 14q28-q29 in a region syntenic to HSA 10q26. The pig contains five closely linked spermadhesin genes, whereas only two spermadhesin genes are present in the cattle genome. Inactive copies of spermadhesin genes are still detectable in the human, chimp, and dog genome while the corresponding region was lost from the rodent genomes of mouse and rat. Within the pig, the five spermadhesin genes contain both highly diverged and highly conserved regions. Interestingly, the pattern of divergence does not correlate with the position of the exons. Evolutionary analyses suggest that the pattern of diversity is shaped by ancestral variation, recombination, and new mutations.

Accelerated evolution and Muller's rachet in endosymbiotic bacteria.

Many bacteria live only within animal cells and infect hosts through cytoplasmic inheritance. These endosymbiotic lineages show distinctive population structure, with small population size and effectively no recombination. As a result, endosymbionts are expected to accumulate mildly deleterious mutations. If these constitute a substantial proportion of new mutations, endosymbionts will show (i) faster sequence evolution and (ii) a possible shift in base composition reflecting mutational bias. Analyses of 16S rDNA of five independently derived endosymbiont clades show, in every case, faster evolution in endosymbionts than in free-living relatives. For aphid endosymbionts (genus Buchnera), coding genes exhibit accelerated evolution and unusually low ratios of synonymous to nonsynonymous substitutions compared to ratios for the same genes for enterics. This concentration of the rate increase in nonsynonymous substitutions is expected under the hypothesis of increased fixation of deleterious mutations. Polypeptides for all Buchnera genes analyzed have accumulated amino acids with codon families rich in A+T, supporting the hypothesis that substitutions are deleterious in terms of polypeptide function. These observations are best explained as the result of Muller's ratchet within small asexual populations, combined with mutational bias. In light of this explanation, two observations reported earlier for Buchnera, the apparent loss of a repair gene and the overproduction of a chaperonin, may reflect compensatory evolution. An alternative hypothesis, involving selection on genomic base composition, is contradicted by the observation that the speedup is concentrated at nonsynonymous sites.