Genetic structure of the genus Leptospira by multilocus enzyme electrophoresis

Thirty strains from the 11 species of the genus Leptospira were studied by multilocus enzyme electrophoresis at 12 enzyme loci, all of which were polymorphic. The mean number of alleles per locus was 6.5. Twenty-five electrophoretic types were distinguished. Grouping of the strains by cluster analysis was in general agreement with species delineation as determined by DNA-DNA hybridization, except for the strains of Leptospira meyeri and Leptospira inadai, which were scattered throughout the genus, reflecting previously recognized taxonomic uncertainties. Analysis of the clonality within Leptospira interrogans sensu stricto indicated that this population was relatively heterogeneous and a lack of gene linkage disequilibrium could not be excluded. There was a genetic discrimination between the pathogenic species and the saprophytic ones. The phenotypically intermediate species (L. inadai and Leptospira fainei) were also genetically separated and were probably closer to the saprophytes than to the pathogens.

Comparaison des coûts du diagnostic moléculaire des rétinites pigmentaires entre le séquençage selon Sanger et un séquençage à haut débit

La rétinite pigmentaire (RP) et l'amaurose congénitale de Leber (LCA) sont deux maladies héréditaires classées dans le groupe des rétinopathies pigmentaires. Plus de 100 gènes ou loci ont étés identifiés dans les RP (comptant pour 60% des patients) et 14 gènes pour les LCA (responsables de 70% des cas). A elles-deux et en ne prenant que les plus fréquentes, ces maladies représentent 948 exons à analyser lorsqu'on recherche une mutation chez un patient. La recherche de mutations au moyen du séquençage génomique classique (séquençage selon Sanger) de tous les gènes de ces deux maladies implique des délais d'analyse importants et des coûts très élevés. La méthode de séquençage à haut débit (avec le séquenceur Roche GS Junior) permet grâce à une puce à ADN le séquençage simultané des 948 exons. Le but de mon travail de Master est de comparer ces deux approches afin de déterminer celle qui est la plus économique et la plus efficiente en temps. Pour cela, j'ai d'abord établi la liste de tous les gènes impliqués dans les RP et LCA, puis identifié tous les exons ainsi que les promoteurs et séquences 3' non traduites. J'ai ensuite calculé le coût théorique d'une analyse de tous les gènes avec chacune des méthodes. J'ai également estimé les coûts à facturer à l'assurance concernant le séquençage à haut débit sur la base des coûts facturables à l'assurance de la méthode Sanger et des bénéfices du laboratoire. Le séquençage des 948 exons par le séquençage à haut débit (avec le GS Junior) représente la technique de séquençage la plus économique et la plus efficiente en temps et constitue donc la méthode de choix dans le screening diagnostic des gènes impliqués dans les RP et LCA. Cette méthode est plus rapide, les réactifs et la machine sont moins coûteux et la laborantine peut analyser un nombre plus important d'exons en un temps moindre, donc elle coûtera moins cher au laboratoire. Cette méthode est donc d'un grand intérêt pour les patients, les assurances et le laboratoire. Cette nouvelle technique de séquençage soulève de nouvelles interrogations telles que la décision de savoir quelle information doit être donnée aux médecins, aux assurances et aux patients. Interrogations auxquelles il devient de plus en plus pressant de répondre.

Functional analysis of a conserved DNA methyltransferase in caulobacter crescentus

CcrM is a DNA methyltransferase that methylates the adenine in GANTC motifs in the chromo-some of the bacterial model Caulobacter crescentus. The loss of the CcrM homolog is lethal in C. crescentus and in several other species of Alphaproteobacteria. In this research, we used different experimental and bioinformatic approaches to determine why CcrM is so critical to the physiology of C. crescentus. We first showed that CcrM is a resident orphan DNA methyltransferase in non-Rickettsiales Alphaproteobacteria and that its gene is strictly conserved in this clade (with only one ex¬ception among the genomes sequenced so far). In C. crescentus, cells depleted in CcrM in rich medium quickly lose viability and present an elongated phenotype characteristic of an im¬pairment in cell division. Using minimal medium instead of rich medium as selective and main¬tenance substrate, we could generate a AccrM mutant that presents a viability comparable to the wild type strain and only mild morphological defects. On the basis of a transcriptomic ap¬proach, we determined that several genes essential for cell division were downregulated in the AccrM strain in minimal medium. We offered decisive arguments to support that the efficient transcription of two of these genes, ftsZ and mipZ, coding respectively for the Z-ring forming GTPase FtsZ and an inhibitor of FtsZ polymerization needed for the correct positioning of the Z- ring at mid-cell, requires the methylation of an adenine in a conserved GANTC motif located in their core promoter region. We propose a model, according to which the genome of C. crescentus encodes a transcriptional activator that requires a methylated adenine in a GANTC context to bind to DNA and suggest that this transcriptional regulator might be the global cell-cycle regulator GcrA. In addition, combining a classic genetic approach and in vitro evolution experiments, we showed that the mortality and cell division defects of the AccrM strain in rich medium are mainly due to limiting intracellular levels of the FtsZ protein. We also studied the dynamics of GANTC methylation in C. crescentus using the SMRT technol¬ogy developed by Pacific Biosciences. Our findings support the commonly accepted model, accord¬ing to which the methylation state of GANTC motifs varies during the cell cycle of C. crescentus: before the initiation of DNA replication, the GANTC motifs are fully-methylated (methylated on both strands); when the DNA gets replicated, the GANTC motifs become hemi-methylated (methyl¬ated on one strand only) and this occurs at different times during replication for different loci along the chromosome depending on their position relative to the origin of replication; the GANTC mo¬tifs are only remethylated after DNA replication has finished as a consequence of the massive and short-lived expression of CcrM in predivisional cells. About 30 GANTC motifs in the C. crescentus chromosome were found to be undermethylated in most of the bacterial population; these might be protected from CcrM activity by DNA binding proteins and some of them could be involved in methylation-based bistable transcriptional switches. - CcrM est une ADN méthyltransférase qui méthyle les adénines dans le contexte GANTC dans le génome de la bactérie modèle Caulobacter crescentus. La perte de l'homologue de CcrM chez C. crescentus et chez plusieurs autres espèces d'Alphaproteobactéries est létale. Dans le courant de cette recherche, nous tentons de déterminer pourquoi la protéine CcrM est cruciale pour la survie de C. crescentus. Nous démontrons d'abord que CcrM est une adénine méthyltransférase orpheline résidente, dont le gène fait partie du génome minimal partagé par les Alphaprotéobactéries non-Rickettsiales (à une exception près). Lorsqu'une souche de C. crescentus est privée de CcrM, sa viabilité décroît rapi¬dement et ses cellules présentent une morphologie allongée qui suggère que la division cellulaire est inhibée. Nous sommes parvenus à créer une souche AccrM en utilisant un milieu minimum, au lieu du milieu riche classiquement employé, comme milieu de sélection et de maintenance pour la souche. Lorsque nous avons étudié le transcriptome de cette souche de C. crescentus privée de CcrM, nous avons pu constater que plusieurs gènes essentiels pour le bon déroulement de la division cellulaire bactérienne étaient réprimés. En particulier, l'expression adéquate des gènes ftsZ et mipZ - qui codent, respectivement, pour FtsZ, la protéine qui constitue, au milieu de la cellule, un anneau protéique qui initie le processus de division et pour MipZ, un inhibiteur de la polymérisation de FtsZ qui est indispensable pour le bon positionnement de l'anneau FtsZ - est dépendante de la présence d'une adénine méthylée dans un motif GANTC conservé situé dans leur région promotrice. Nous présentons un modèle selon lequel le génome de C. crescentus code pour un facteur de transcription qui exige la présence d'une adénine méthylée dans un contexte GANTC pour s'attacher à l'ADN et nous suggérons qu'il pourrait s'agir du régulateur global du cycle cellulaire GcrA. En outre, nous montrons, en combinant la génétique classique et une approche basée sur l'évolution expérimentale, que la mortalité et l'inhibition de la division cellulaire caractéristiques de la souche àccrMeη milieu riche sont dues à des niveaux excessivement bas de protéine FtsZ. Nous avons aussi étudié la dynamique de la méthylation du chromosome de C. crescentus sur la base de la technologie SMRT développée par Pacific Biosciences. Nous confirmons le modèle communément accepté, qui affirme que l'état de méthylation des motifs GANTC change durant le cycle cellulaire de C. crescentus: les motifs GANTC sont complètement méthylés (méthylés sur les deux brins) avant de début de la réplication de l'ADN; ils deviennent hémi-méthylés (méthylés sur un brin seulement) une fois répliqués, ce qui arrive à différents moments durant la réplication pour différents sites le long du chromosome en fonction de leur position par rapport à l'origine de répli-cation; finalement, les motifs GANTC sont reméthylés après la fin de la réplication du chromosome lorsque la protéine CcrM est massivement, mais très transitoirement, produite. Par ailleurs, nous identifions dans le chromosome de C. crescentus environ 30 motifs GANTC qui restent en perma-nence non-méthylés dans une grande partie de la population bactérienne; ces motifs sont probable-ment protégés de l'action de CcrM par des protéines qui s'attachent à l'ADN et certains d'entre eux pourraient être impliqués dans des mécanismes de régulation générant une transcription bistable.

Evidence for polygenic adaptation to pathogens in the human genome.

Most approaches aiming at finding genes involved in adaptive events have focused on the detection of outlier loci, which resulted in the discovery of individually "significant" genes with strong effects. However, a collection of small effect mutations could have a large effect on a given biological pathway that includes many genes, and such a polygenic mode of adaptation has not been systematically investigated in humans. We propose here to evidence polygenic selection by detecting signals of adaptation at the pathway or gene set level instead of analyzing single independent genes. Using a gene-set enrichment test to identify genome-wide signals of adaptation among human populations, we find that most pathways globally enriched for signals of positive selection are either directly or indirectly involved in immune response. We also find evidence for long-distance genotypic linkage disequilibrium, suggesting functional epistatic interactions between members of the same pathway. Our results show that past interactions with pathogens have elicited widespread and coordinated genomic responses, and suggest that adaptation to pathogens can be considered as a primary example of polygenic selection.

Common variants associated with plasma triglycerides and risk for coronary artery disease.

Triglycerides are transported in plasma by specific triglyceride-rich lipoproteins; in epidemiological studies, increased triglyceride levels correlate with higher risk for coronary artery disease (CAD). However, it is unclear whether this association reflects causal processes. We used 185 common variants recently mapped for plasma lipids (P < 5 × 10(-8) for each) to examine the role of triglycerides in risk for CAD. First, we highlight loci associated with both low-density lipoprotein cholesterol (LDL-C) and triglyceride levels, and we show that the direction and magnitude of the associations with both traits are factors in determining CAD risk. Second, we consider loci with only a strong association with triglycerides and show that these loci are also associated with CAD. Finally, in a model accounting for effects on LDL-C and/or high-density lipoprotein cholesterol (HDL-C) levels, the strength of a polymorphism's effect on triglyceride levels is correlated with the magnitude of its effect on CAD risk. These results suggest that triglyceride-rich lipoproteins causally influence risk for CAD.

A step-by-step tutorial to use HierFstat to analyse populations hierarchically structured at multiple levels.

The populations of parasites and infectious agents are most of the time structured in complex hierarchy that lies beyond the classical nested design described by Wright's F-statistics (F(IS), F(ST) and F(IT)). In this note we propose a user-friendly step-by-step notice for using recent software (HierFstat) that computes and test fixation indices for any hierarchical structure. We add some tricks and tips for some special data kind (haploid, single locus), some other procedure (bootstrap over loci) and how to handle crossed factors.

Effects of long-term averaging of quantitative blood pressure traits on the detection of genetic associations.

Blood pressure (BP) is a heritable, quantitative trait with intraindividual variability and susceptibility to measurement error. Genetic studies of BP generally use single-visit measurements and thus cannot remove variability occurring over months or years. We leveraged the idea that averaging BP measured across time would improve phenotypic accuracy and thereby increase statistical power to detect genetic associations. We studied systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP) averaged over multiple years in 46,629 individuals of European ancestry. We identified 39 trait-variant associations across 19 independent loci (p < 5 × 10(-8)); five associations (in four loci) uniquely identified by our LTA analyses included those of SBP and MAP at 2p23 (rs1275988, near KCNK3), DBP at 2q11.2 (rs7599598, in FER1L5), and PP at 6p21 (rs10948071, near CRIP3) and 7p13 (rs2949837, near IGFBP3). Replication analyses conducted in cohorts with single-visit BP data showed positive replication of associations and a nominal association (p < 0.05). We estimated a 20% gain in statistical power with long-term average (LTA) as compared to single-visit BP association studies. Using LTA analysis, we identified genetic loci influencing BP. LTA might be one way of increasing the power of genetic associations for continuous traits in extant samples for other phenotypes that are measured serially over time.

Molecular taxonomy of the Formica rufa group (red wood ants) (Hymenoptera: Formicidae): a new cryptic species in the Swiss Alps?

Because of their beneficial impact on forest ecosystems, European red wood ants (Formica rufa group) are protected by law in many European countries and are considered to be among the most reliable bioindicators of forest stability. However, their taxonomy has been much debated and, unfortunately, it is too often neglected. This happens mainly because the morphology-based method for species delimitation requires lots of time and experience. We therefore employed 9 microsatellite loci and mitochondrial DNA (COI gene) to verify the power of genetic markers for red wood ant species delimitation and to investigate the cryptic diversity of these ants within the Eastern Swiss Alps. We analyzed 83 nests belonging to all red wood ant species that occur in the Swiss National Park area. Genetic data indicated that these species represent different genetic pools. Moreover, results showed that Formica aquilonia YARROW, 1955 and F. paralugubris SEIFERT, 1996 often hybridize within the Park, confirming that these two species are genetically very close and could have diverged only recently. Nevertheless, microsatellites also revealed that one entire population, located in the Minger Valley and morphologically identified as F. lugubris ZETTERSTEDT, 1838, is genetically different to all other analyzed F. lugubris populations found within the same area and to other red wood ant species. These findings, confirmed by mitochondrial DNA analyses, suggest the existence of a new cryptic species within the Eastern Swiss Alps. This putative cryptic species has been provisionally named F. lugubris-A2. These results have a great importance for future conservation plans, monitoring and evolutionary studies on these protected ants.

Systematic molecular and cytogenetic screening of 100 patients with marfanoid syndromes and intellectual disability.

The association of marfanoid habitus (MH) and intellectual disability (ID) has been reported in the literature, with overlapping presentations and genetic heterogeneity. A hundred patients (71 males and 29 females) with a MH and ID were recruited. Custom-designed 244K array-CGH (Agilent®; Agilent Technologies Inc., Santa Clara, CA) and MED12, ZDHHC9, UPF3B, FBN1, TGFBR1 and TGFBR2 sequencing analyses were performed. Eighty patients could be classified as isolated MH and ID: 12 chromosomal imbalances, 1 FBN1 mutation and 1 possibly pathogenic MED12 mutation were found (17%). Twenty patients could be classified as ID with other extra-skeletal features of the Marfan syndrome (MFS) spectrum: 4 pathogenic FBN1 mutations and 4 chromosomal imbalances were found (2 patients with both FBN1 mutation and chromosomal rearrangement) (29%). These results suggest either that there are more loci with genes yet to be discovered or that MH can also be a relatively non-specific feature of patients with ID. The search for aortic complications is mandatory even if MH is associated with ID since FBN1 mutations or rearrangements were found in some patients. The excess of males is in favour of the involvement of other X-linked genes. Although it was impossible to make a diagnosis in 80% of patients, these results will improve genetic counselling in families.

GENCODE: the reference human genome annotation for The ENCODE Project.

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

Genetic causes of transitions from sexual reproduction to asexuality in plants and animals.

The persistence of sexual reproduction in the face of competition from asexual invaders is more likely if asexual lineages are produced infrequently or have low fitness. The generation rate and success of new asexual lineages will be influenced by the proximate mechanisms underlying transitions to asexuality. As such, characterization of these mechanisms can help explain the distribution of reproductive modes among natural populations. Here, we synthesize the literature addressing proximate causes of transitions from sexual to asexual reproduction in plants and animals. In cyclical and facultatively asexual taxa, individual mutations can cause obligate asexuality. The evolution of asexuality in obligately sexual groups is more complex, requiring the simultaneous acquisition of two traits generally controlled by different genetic factors: unreduced gamete formation and spontaneous development of unfertilized gametes. At least three 'pre-adaptations' could favour transitions to obligate asexuality in obligate sexuals. First, linkage among loci affecting separate key components of asexuality facilitates its spread, with evidence for these linkage blocks in plants. Second, asexuality should evolve more readily in haplodiploids; support for this hypothesis comes from two examples where a single locus causes transitions to asexuality. Third, standing genetic variation for the production of unreduced gametes could facilitate transitions to asexuality, but whether the ability to produce unreduced gametes contributes to the evolution of obligate asexuality remains unclear. We close by reviewing the associations between asexuality, hybridization and polyploidy, and argue that current data suggest that hybridization is more likely to play a causal role in transitions to asexuality than polyploidy.

CUBN is a gene locus for albuminuria.

Identification of genetic risk factors for albuminuria may alter strategies for early prevention of CKD progression, particularly among patients with diabetes. Little is known about the influence of common genetic variants on albuminuria in both general and diabetic populations. We performed a meta-analysis of data from 63,153 individuals of European ancestry with genotype information from genome-wide association studies (CKDGen Consortium) and from a large candidate gene study (CARe Consortium) to identify susceptibility loci for the quantitative trait urinary albumin-to-creatinine ratio (UACR) and the clinical diagnosis microalbuminuria. We identified an association between a missense variant (I2984V) in the CUBN gene, which encodes cubilin, and both UACR (P = 1.1 × 10(-11)) and microalbuminuria (P = 0.001). We observed similar associations among 6981 African Americans in the CARe Consortium. The associations between this variant and both UACR and microalbuminuria were significant in individuals of European ancestry regardless of diabetes status. Finally, this variant associated with a 41% increased risk for the development of persistent microalbuminuria during 20 years of follow-up among 1304 participants with type 1 diabetes in the prospective DCCT/EDIC Study. In summary, we identified a missense CUBN variant that associates with levels of albuminuria in both the general population and in individuals with diabetes.

Heterozygosity-fitness correlations among wild populations of European tree frogs (Hyla arborea) detect fixation load.

Quantifying the impacts of inbreeding and genetic drift on fitness traits in fragmented populations is becoming a major goal in conservation biology. Such impacts occur at different levels and involve different sets of loci. Genetic drift randomly fixes slightly deleterious alleles leading to different fixation load among populations. By contrast, inbreeding depression arises from highly deleterious alleles in segregation within a population and creates variation among individuals. A popular approach is to measure correlations between molecular variation and phenotypic performances. This approach has been mainly used at the individual level to detect inbreeding depression within populations and sometimes at the population level but without consideration about the genetic processes measured. For the first time, we used in this study a molecular approach considering both the interpopulation and intrapopulation level to discriminate the relative importance of inbreeding depression vs. fixation load in isolated and non-fragmented populations of European tree frog (Hyla arborea), complemented with interpopulational crosses. We demonstrated that the positive correlations observed between genetic heterozygosity and larval performances on merged data were mainly caused by co-variations in genetic diversity and fixation load among populations rather than by inbreeding depression and segregating deleterious alleles within populations. Such a method is highly relevant in a conservation perspective because, depending on how populations lose fitness (inbreeding vs. fixation load), specific management actions may be designed to improve the persistence of populations.

Restricted gene flow and subpopulation differentiation in Silene dioica

The size of breeding units and the hierarchical population structure of the dioecious perennial herb Silene dioica were investigated on four closely situated island populations in the Skeppsvik Archipelago in northern Sweden. F-statistics analyses of nine polymorphic allozyme loci revealed that plants on single islands are divided into many small breeding units, between 0.2 m2 and 6 m2. Hierarchical analyses showed that levels of differentiation among subpopulations within islands (FPL=0.080) were about twice as high as among islands (FLT=0.048). These results are discussed in the light of what is known about pollen and seed movement in the archipelago.

Forensic identification of urine samples: a comparison between nuclear and mitochondrial DNA markers.

Urine samples from 20 male volunteers of European Caucasian origin were stored at 4 degrees C over a 4-month period in order to compare the identification potential of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) markers. The amount of nDNA recovered from urines dramatically declined over time. Consequently, nDNA likelihood ratios (LRs) greater than 1,000 were obtained for 100, 70 and 55% of the urines analysed after 6, 60 and 120 days, respectively. For the mtDNA, HVI and HVII sequences were obtained for all samples tested, whatever the period considered. Nevertheless, the highest mtDNA LR of 435 was relatively low compared to its nDNA equivalent. Indeed, LRs obtained with only three nDNA loci could easily exceed this value and are quite easier to obtain. Overall, the joint use of nDNA and mtDNA markers enabled the 20 urine samples to be identified, even after the 4-month period.