Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

Quantitative Ultrasound of the Heel and Fracture Risk Assessment: an Updated Meta-Analysis

Objectives: Quantitative ultrasound (QUS) is an attractive method for assessing fracture risk because it is portable, inexpensive, without ionizing radiation, and available in areas of the world where DXA is not readily accessible or affordable. However, the diversity of QUS scanners and variability of fracture outcomes measured in different studies is an important obstacle to widespread utilisation of QUS for fracture risk assessment. We aimed in this review to assess the predictive power of heel QUS for fractures considering different characteristics of the association (QUS parameters and fracture outcomes measured, QUS devices, study populations, and independence from DXA-measured bone density).Materials/Methods : We conducted an inverse-variance randomeffects meta-analysis of prospective studies with heel QUS measures at baseline and fracture outcomes in their follow-up. Relative risks (RR) per standard deviation (SD) of different QUS parameters (broadband ultrasound attenuation [BUA], speed of sound &SOS;, stiffness index &SI;, and quantitative ultrasound index [QUI]) for various fracture outcomes (hip, vertebral, any clinical, any osteoporotic, and major osteoporotic fractures) were reported based on study questions.Results : 21 studies including 55,164 women and 13,742 men were included with a total follow-up of 279,124 person-years. All four QUS parameters were associated with risk of different fractures. For instance, RR of hip fracture for 1 SD decrease of BUA was 1.69 (95% CI 1.43-2.00), SOS was 1.96 (95% CI 1.64-2.34), SI was 2.26 (95%CI 1.71-2.99), and QUI was 1.99 (95% CI 1.49-2.67). Validated devices from different manufacturers predicted fracture risks with a similar performance (meta-regression p-values>0.05 for difference of devices). There was no sign of publication bias among the studies. QUS measures predicted fracture with a similar performance in men and women. Meta-analysis of studies with QUS measures adjusted for hip DXA showed a significant and independent association with fracture risk (RR/SD for BUA =1.34 [95%CI 1.22-1.49]).Conclusions : This study confirms that QUS of the heel using validated devices predicts risk of different fracture outcomes in elderly men and women. Further research and international collaborations are needed for standardisation of QUS parameters across various manufacturers and inclusion of QUS in fracture risk assessment tools. Disclosure of Interest : None declared.

Genetic variability of hull-less barley accessions based on molecular and quantitative data

The objective of this work was to characterize and quantify the genetic, molecular, and agronomic variability of hull-less barley genotypes, for the selection of parents and identification of genotypes adapted to the irrigated production system in the Brazilian Cerrado. Eighteen hull-less barley accessions were evaluated, and three covered barley accessions served as reference. The characterization was based on 157 RAPD molecular markers and ten agronomic traits. Genetic distance matrices were obtained based on molecular markers and quantitative traits. Graphic grouping and dispersion analyses were performed. Genetic, molecular, and agronomic variability was high among genotypes. Ethiopian accessions were genetically more similar, and the Brazilian ones were genetically more distant. For agronomic traits, two more consistent groupings were obtained, one with the most two-rowed materials, and the other with six-rowed materials. The more diverging materials were the two-rowed CI 13453, CN Cerrado 5, CN Cerrado 1, and CN Cerrado 2. The PI 356466, CN Cerrado 1, PI 370799, and CI 13453 genotypes show agronomic traits of interest and, as genetically different genotypes, they are indicated for crossing, in breeding programs.

Quantitative morphometry analysis of the fetal brain using clinical MR imaging

In vivo fetal magnetic resonance imaging provides aunique approach for the study of early human braindevelopment [1]. In utero cerebral morphometry couldpotentially be used as a marker of the cerebralmaturation and help to distinguish between normal andabnormal development in ambiguous situations. However,this quantitative approach is a major challenge becauseof the movement of the fetus inside the amniotic cavity,the poor spatial resolution provided by very fast MRIsequences and the partial volume effect. Extensiveefforts are made to deal with the reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution [2,3,4]. Frameworkswere developed for the segmentation of specific regionsof the fetal brain such as posterior fossa, brainstem orgerminal matrix [5,6], or for the entire brain tissue[7,8], applying the Expectation-Maximization MarkovRandom Field (EM-MRF) framework. However, many of theseprevious works focused on the young fetus (i.e. before 24weeks) and use anatomical atlas priors to segment thedifferent tissue or regions. As most of the gyraldevelopment takes place after the 24th week, acomprehensive and clinically meaningful study of thefetal brain should not dismiss the third trimester ofgestation. To cope with the rapidly changing appearanceof the developing brain, some authors proposed a dynamicatlas [8]. To our opinion, this approach however faces arisk of circularity: each brain will be analyzed /deformed using the template of its biological age,potentially biasing the effective developmental delay.Here, we expand our previous work [9] to proposepost-processing pipeline without prior that allow acomprehensive set of morphometric measurement devoted toclinical application. Data set & Methods: Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences (TR 7000 ms, TE 180 ms, FOV 40 x 40 cm,slice thickness 5.4mm, in plane spatial resolution1.09mm). For each fetus, 6 axial volumes shifted by 1 mmwere acquired under motherâeuro?s sedation (about 1min pervolume). First, each volume is segmentedsemi-automatically using region-growing algorithms toextract fetal brain from surrounding maternal tissues.Inhomogeneity intensity correction [10] and linearintensity normalization are then performed. Brain tissues(CSF, GM and WM) are then segmented based on thelow-resolution volumes as presented in [9]. Ahigh-resolution image with isotropic voxel size of 1.09mm is created as proposed in [2] and using B-splines forthe scattered data interpolation [11]. Basal gangliasegmentation is performed using a levet setimplementation on the high-resolution volume [12]. Theresulting white matter image is then binarized and givenas an input in FreeSurfer software(http://surfer.nmr.mgh.harvard.edu) to providetopologically accurate three-dimensional reconstructionsof the fetal brain according to the local intensitygradient. References: [1] Guibaud, Prenatal Diagnosis29(4) (2009). [2] Rousseau, Acad. Rad. 13(9), 2006. [3]Jiang, IEEE TMI 2007. [4] Warfield IADB, MICCAI 2009. [5]Claude, IEEE Trans. Bio. Eng. 51(4) 2004. [6] Habas,MICCAI 2008. [7] Bertelsen, ISMRM 2009. [8] Habas,Neuroimage 53(2) 2010. [9] Bach Cuadra, IADB, MICCAI2009. [10] Styner, IEEE TMI 19(39 (2000). [11] Lee, IEEETrans. Visual. And Comp. Graph. 3(3), 1997. [12] BachCuadra, ISMRM 2010.

Quantitative multiparameter assays to measure the effect of adjuvants on human antigen-specific CD8 T-cell responses.

Large numbers and functionally competent T cells are required to protect from diseases for which antibody-based vaccines have consistently failed (1), which is the case for many chronic viral infections and solid tumors. Therefore, therapeutic vaccines aim at the induction of strong antigen-specific T-cell responses. Novel adjuvants have considerably improved the capacity of synthetic vaccines to activate T cells, but more research is necessary to identify optimal compositions of potent vaccine formulations. Consequently, there is a great need to develop accurate methods for the efficient identification of antigen-specific T cells and the assessment of their functional characteristics directly ex vivo. In this regard, hundreds of clinical vaccination trials have been implemented during the last 15 years, and monitoring techniques become more and more standardized.

Prediction of hip fracture risk by quantitative ultrasound in more than 7000 Swiss women > or =70 years of age: comparison of three technologically different bone ultrasound devices in the SEMOF study.

To compare the prediction of hip fracture risk of several bone ultrasounds (QUS), 7062 Swiss women > or =70 years of age were measured with three QUSs (two of the heel, one of the phalanges). Heel QUSs were both predictive of hip fracture risk, whereas the phalanges QUS was not. INTRODUCTION: As the number of hip fracture is expected to increase during these next decades, it is important to develop strategies to detect subjects at risk. Quantitative bone ultrasound (QUS), an ionizing radiation-free method, which is transportable, could be interesting for this purpose. MATERIALS AND METHODS: The Swiss Evaluation of the Methods of Measurement of Osteoporotic Fracture Risk (SEMOF) study is a multicenter cohort study, which compared three QUSs for the assessment of hip fracture risk in a sample of 7609 elderly ambulatory women > or =70 years of age. Two QUSs measured the heel (Achilles+; GE-Lunar and Sahara; Hologic), and one measured the heel (DBM Sonic 1200; IGEA). The Cox proportional hazards regression was used to estimate the hazard of the first hip fracture, adjusted for age, BMI, and center, and the area under the ROC curves were calculated to compare the devices and their parameters. RESULTS: From the 7609 women who were included in the study, 7062 women 75.2 +/- 3.1 (SD) years of age were prospectively followed for 2.9 +/- 0.8 years. Eighty women reported a hip fracture. A decrease by 1 SD of the QUS variables corresponded to an increase of the hip fracture risk from 2.3 (95% CI, 1.7, 3.1) to 2.6 (95% CI, 1.9, 3.4) for the three variables of Achilles+ and from 2.2 (95% CI, 1.7, 3.0) to 2.4 (95% CI, 1.8, 3.2) for the three variables of Sahara. Risk gradients did not differ significantly among the variables of the two heel QUS devices. On the other hand, the phalanges QUS (DBM Sonic 1200) was not predictive of hip fracture risk, with an adjusted hazard risk of 1.2 (95% CI, 0.9, 1.5), even after reanalysis of the digitalized data and using different cut-off levels (1700 or 1570 m/s). CONCLUSIONS: In this elderly women population, heel QUS devices were both predictive of hip fracture risk, whereas the phalanges QUS device was not.

Equilibration and disequilibration between monazite and garnet: indication from phase-composition and quantitative texture analysis

The integration of information which can be gained from accessory [i.e. age (t)] and rock-forming minerals [i.e. temperature (T) and pressure (P)] requires a more profound understanding of the equilibration kinetics during metamorphic processes. This paper presents an approach comparing conventional P-T estimate from equilibrated assemblages of rock-forming minerals with temperature data derived from yttrium-garnet-monazite (YGM) and yttrium-garnet-xenotime (YGX) geothermometry. Such a comparison provides an initial indication on differences between equilibration of major and trace elements. Regarding this purpose, two migmatites, two polycyclic and one monocyclic gneiss from the Central Alps (Switzerland, northern Italy) were investigated. While the polycyclic samples exhibit trace-element equilibration between monazite and garnet grains assigned to the same metamorphic event, there are relics of monazite and garnet obviously surviving independent of their textural position. These observations suggest that surface processes dominate transport processes during equilibration of those samples. The monocyclic gneiss, on the contrary, displays rare isolated monazite with equilibration of all elements, despite comparably large transport distances. With a nearly linear crystal-size distribution of the garnet grain population, growth kinetics, related to the major elements, were likely surface-controlled in this sample. In contrast to these completely equilibrated examples, the migmatites indicate disequilibrium between garnet and monazite with a change in REE patterns on garnet transects. The cause for this disequilibrium may be related to a potential disequilibrium initiated by a changing bulk chemistry during melt segregation. While migmatite environments are expected to support high transport rates (i.e. high temperatures and melt presence), the evolution of equilibration in migmatites is additionaly related to change in chemistry. As a key finding, surface-controlled equilibration kinetics seem to dominate transport-controlled processes in the investigated samples. This may be decisive information towards the understanding of age data derived from monazite.

Analytical methods for the quantitative determination of selective serotonin reuptake inhibitors for therapeutic drug monitoring purposes in patients.

Five selective serotonin reuptake inhibitors (SSRIs) have been introduced recently: citalopram, fluoxetine, fluvoxamine, paroxetine and sertraline. Although no therapeutic window has been defined for SSRIs, in contrast to tricyclic antidepressants, analytical methods for therapeutic drug monitoring of SSRIs are useful in several instances. SSRIs differ widely in their chemical structure and in their metabolism. The fact that some of them have N-demethylated metabolites, which are also SSRIs, requires that methods be available which allow therapeutic drug monitoring of the parent compounds and of these active metabolites. most procedures are based on prepurification of the SSRIs by liquid-liquid extraction before they are submitted to separation by chromatographic procedures (high-performance liquid chromatography, gas chromatography, thin layer chromatography) and detection by various detectors (UV, fluorescence, electrochemical detector, nitrogen-phosphorus detector, mass spectrometry). This literature review shows that most methods allow quantitative determination of SSRIs in plasma, in the lower ng/ml range, and that they are, therefore, suitable for therapeutic drug monitoring purposes of this category of drugs.

Genomic determinants of the efficiency of internal ribosomal entry sites of viral and cellular origin

Variation in cellular gene expression levels has been shown to be inherited. Expression is controlled at transcriptional and post-transcriptional levels. Internal ribosome entry sites (IRES) are used by viruses to bypass inhibition of cap-dependent translation, and by eukaryotic cells to control translation under conditions when protein synthesis is inhibited. We aimed at identifying genomic determinants of variability in IRES-mediated translation of viral [Encephalomyocarditis virus (EMCV)] and cellular IRES [X-linked inhibitor-of-apoptosis (XIAP) and c-myc]. Bicistronic lentiviral constructs expressing two fluorescent reporters were used to transduce laboratory and B lymphoblastoid cell lines [15 CEPH pedigrees (n = 205) and 50 unrelated individuals]. IRES efficiency varied according to cell type and among individuals. Control of IRES activity has a significant genetic component (h(2) of 0.47 and 0.36 for EMCV and XIAP, respectively). Quantitative linkage analysis identified a suggestive locus (LOD 2.35) on chromosome 18q21.2, and genome-wide association analysis revealed of a cluster of SNPs on chromosome 3, intronic to the FHIT gene, marginally associated (P = 5.9E-7) with XIAP IRES function. This study illustrates the in vitro generation of intermediate phenotypes by using cell lines for the evaluation of genetic determinants of control of elements such as IRES.

Connexin 43 expression in human hypertrophied heart due to pressure and volume overload.

Left ventricular hypertrophy (LVH) is due to pressure overload or mechanical stretch and is thought to be associated with remodeling of gap-junctions. We investigated whether the expression of connexin 43 (Cx43) is altered in humans in response to different degrees of LVH. The expression of Cx43 was analyzed by quantitative polymerase chain reaction, Western blot analysis and immunohistochemistry on left ventricular biopsies from patients undergoing aortic or mitral valve replacement. Three groups were analyzed: patients with aortic stenosis with severe LVH (n=9) versus only mild LVH (n=7), and patients with LVH caused by mitral regurgitation (n=5). Cx43 mRNA expression and protein expression were similar in the three groups studied. Furthermore, immunohistochemistry revealed no change in Cx43 distribution. We can conclude that when compared with mild LVH or with LVH due to volume overload, severe LVH due to chronic pressure overload is not accompanied by detectable changes of Cx43 expression or spatial distribution.

The protective role of transferrin in Müller glial cells after iron-induced toxicity.

PURPOSE: Transferrin (Tf) expression is enhanced by aging and inflammation in humans. We investigated the role of transferrin in glial protection. METHODS: We generated transgenic mice (Tg) carrying the complete human transferrin gene on a C57Bl/6J genetic background. We studied human (hTf) and mouse (mTf) transferrin localization in Tg and wild-type (WT) C57Bl/6J mice using immunochemistry with specific antibodies. Müller glial (MG) cells were cultured from explants and characterized using cellular retinaldehyde binding protein (CRALBP) and vimentin antibodies. They were further subcultured for study. We incubated cells with FeCl(3)-nitrilotriacetate to test for the iron-induced stress response; viability was determined by direct counting and measurement of lactate dehydrogenase (LDH) activity. Tf expression was determined by reverse transcriptase-quantitative PCR with human- or mouse-specific probes. hTf and mTf in the medium were assayed by ELISA or radioimmunoassay (RIA), respectively. RESULTS: mTf was mainly localized in retinal pigment epithelium and ganglion cell layers in retina sections of both mouse lines. hTf was abundant in MG cells. The distribution of mTf and hTf mRNA was consistent with these findings. mTf and hTf were secreted into the medium of MG cell primary cultures. Cells from Tg mice secreted hTf at a particularly high level. However, both WT and Tg cell cultures lose their ability to secrete Tf after a few passages. Tg MG cells secreting hTf were more resistant to iron-induced stress toxicity than those no longer secreted hTf. Similarly, exogenous human apo-Tf, but not human holo-Tf, conferred resistance to iron-induced stress on MG cells from WT mice. CONCLUSIONS: hTf localization in MG cells from Tg mice was reminiscent of that reported for aged human retina and age-related macular degeneration, both conditions associated with iron deposition. The role of hTf in protection against toxicity in Tg MG cells probably involves an adaptive mechanism developed in neural retina to control iron-induced stress.

"Soap-Bubble" visualization and quantitative analysis of 3D coronary magnetic resonance angiograms.

In order to compare coronary magnetic resonance angiography (MRA) data obtained with different scanning methodologies, adequate visualization and presentation of the coronary MRA data need to be ensured. Furthermore, an objective quantitative comparison between images acquired with different scanning methods is desirable. To address this need, a software tool ("Soap-Bubble") that facilitates visualization and quantitative comparison of 3D volume targeted coronary MRA data was developed. In the present implementation, the user interactively specifies a curved subvolume (enclosed in the 3D coronary MRA data set) that closely encompasses the coronary arterial segments. With a 3D Delaunay triangulation and a parallel projection, this enables the simultaneous display of multiple coronary segments in one 2D representation. For objective quantitative analysis, frequently explored quantitative parameters such as signal-to-noise ratio (SNR); contrast-to-noise ratio (CNR); and vessel length, sharpness, and diameter can be assessed. The present tool supports visualization and objective, quantitative comparisons of coronary MRA data obtained with different scanning methods. The first results obtained in healthy adults and in patients with coronary artery disease are presented.

The Munc13 proteins differentially regulate readily releasable pool dynamics and calcium-dependent recovery at a central synapse.

The Munc13 gene family encodes molecules located at the synaptic active zone that regulate the reliability of synapses to encode information over a wide range of frequencies in response to action potentials. In the CNS, proteins of the Munc13 family are critical in regulating neurotransmitter release and synaptic plasticity. Although Munc13-1 is essential for synaptic transmission, it is paradoxical that Munc13-2 and Munc13-3 are functionally dispensable at some synapses, although their loss in other synapses leads to increases in frequency-dependent facilitation. We addressed this issue at the calyx of Held synapse, a giant glutamatergic synapse that we found to express all these Munc13 isoforms. We studied their roles in the regulation of synaptic transmission and their impact on the reliability of information transfer. Through detailed electrophysiological analyses of Munc13-2, Munc13-3, and Munc13-2-3 knock-out and wild-type mice, we report that the combined loss of Munc13-2 and Munc13-3 led to an increase in the rate of calcium-dependent recovery and a change in kinetics of release of the readily releasable pool. Furthermore, viral-mediated overexpression of a dominant-negative form of Munc13-1 at the calyx demonstrated that these effects are Munc13-1 dependent. Quantitative immunohistochemistry using Munc13-fluorescent protein knock-in mice revealed that Munc13-1 is the most highly expressed Munc13 isoform at the calyx and the only one highly colocalized with Bassoon at the active zone. Based on these data, we conclude that Munc13-2 and Munc13-3 isoforms limit the ability of Munc13-1 to regulate calcium-dependent replenishment of readily releasable pool and slow pool to fast pool conversion in central synapses.

Population genetics of Capercaillie (Tetrao urogallus) in the Jura and the Pyrenees: a non-invasive approach to avian conservation genetics

RÉSUMÉ Le Grand tétras est un galliforme de montagne apparenté au faisan et au tétras lyre. Il est distribué de manière continue à travers la toundra et les montagnes de moyenne altitude en Europe de l'ouest. Toutefois, les populations d'Europe de l'ouest ont subi un déclin constant au cours des derniers siècles. Les causes de ce déclin sont probablement liées à l'activité humaine, telle .que l'élevage ou le tourisme, qui ont engendré une modification et une fragmentation de l'habitat de l'espèce. Malheureusement, les populations soumises à de forts déclins démographiques peuvent subir des effets génétiques (augmentation de la consanguinité et perte de diversité génétique) pouvant diminuer leur potentiel de reproduction et conduire irrémédiablement à l'extinction. Cette thèse présente les analyses conduites dans le but d'estimer l'impact du déclin démographique des populations de Grand tétras sur l'étendue et la distribution de leur variabilité génétique dans le Jura et dans les Pyrénées. Du fait de la législation locale protégeant les tétraonidés en général, mais également en raison de la biologie très cryptique du Grand tétras, l'ensemble des analyses de cette étude a été réalisé à partir de matériel génétique extrait des fientes (ou échantillonnage génétique non invasif). Dans la première partie de l'étude, je détaille les protocoles d'extraction. d'ADN et d'amplification par PCR modifiés à partir des protocoles classiques utilisant des échantillons conventionnels, riches en ADN. L'utilisation d'ADN fécal impose des contraintes dues à la mauvaise qualité et à la faible quantité du matériel génétique à disposition dans les fientes. Ces contraintes ont pu être partiellement contournées en réalisant des répétitions multiples du génotypage afin d'obtenir un degré de fiabilité suffisante. J'ai également analysé les causes de la dégradation de l'ADN dans les excréments. Parmi les causes les plus communes, telles que l'activité bactérienne, l'hydrolyse spontanée et la dégradation enzymatique par les DNases libres, c'est ce dernier facteur qui apparaît comme étant la cause majeure et la plus rapide responsable de la dégradation de la qualité des échantillons. La rapidité de l'action enzymatique suggère que les plans d'échantillonnages de excréments sur le terrain pourraient être optimisés en les réalisant dans des conditions climatiques froides et sèches, favorisant ainsi l'inhibition des DNases. La seconde partie de la thèse est une étude par simulation visant à déterminer la capacité du logiciel Structure à identifier les structures génétiques complexes et hiérarchiques fréquemment rencontrées dans les populations naturelles, et ce en utilisant différents types de marqueurs génétiques. Les troisième et quatrième parties de cette thèse décrivent le statut génétique des populations résiduelles du Jura et des Pyrénées à partir de l'analyse de 11 loci microsatellites. Nous n'avons pas pu mettre en évidence dans les deux populations des effets liés à la consanguinité ou à la réduction de la diversité génétique. De plus, la différenciation génétique entre les patches d'habitats favorables reste modérée et corrélée à la distance géographique, ce qui suggère que la dispersion d'individus entre les patches a été importante au moins pendant ces dernières générations. La comparaison des paramètres de la diversité génétique avec ceux d'autres populations de Grand tétras, ou d'autres espèces proches, indique que la population du Jura a retenu une proportion importante de sa diversité originelle. Ces résultats suggèrent que le déclin récent des populations a jusqu'ici eu un impact modéré sur les facteurs génétiques et que ces populations semblent avoir conservé le potentiel génétique nécessaire à leur survie à long terme. Finalement, en cinquième partie, l'analyse de l'apparentement entre les mâles qui participent à la parade sur les places de chant (leks) indique que ces derniers sont distribués en agrégats de manière non aléatoire, préférentiellement entre individus apparentés. De plus, la corrélation entre les distances génétique et géographique entre les leks est en accord avec les motifs d'isolement par la distance mis en évidence à d'autres niveaux hiérarchiques (entre patches d'habitat et populations), ainsi qu'avec les études menées sur d'autres espèces ayant choisi ce même système de reproduction. En conclusion, cette première étude basée uniquement sur de l'ADN nucléaire aviaire extrait à partir de fèces a fourni des informations nouvelles qui n'auraient pas pu être obtenues par une méthode d'observation sur le terrain ou d'échantillonnage génétique classique. Aucun oiseau n'a été dérangé ou capturé, et les résultats sont comparables à d'autres études concernant des espèces proches. Néanmoins, la taille de ces populations approche des niveaux au-dessous desquels la survie à long terme est fortement incertaine. La persistance de la diversité génétique pour les prochaines générations reste en conséquence liée à la survie des adultes et à une reprise du succès de la reproduction. ABSTRACT Capercaillie (Tetrao urogallus) is a large grouse that is continuously distributed across the tundra and the mid-high mountains of Western Europe. However, the populations in Western Europe have been showing a constant decline during the last decades. The causes for this decline are possibly related to human activities, such as cattle breeding and tourism that have both led to habitat modification and fragmentation. Unfortunately, populations that have undergone drastic demographic bottlenecks often go through genetic processes of inbreeding and loss of diversity that decrease their fitness and eventually lead to extinction. This thesis presents the investigations conducted to estimate the impact of the demographic decline of capercaillie populations on the extent and distribution of their genetic variability in the Jura and in the Pyrenees mountains. Because grouse are protected by wildlife legislation, and also because of the cryptic behaviour of capercaillie, all DNA material used in this study was extracted from faeces (non-invasive genetic sampling). In the first part of my thesis, I detail the protocols of DNA extraction and PCR amplification adapted from classical methods using conventional DNA-rich samples. The use of faecal DNA imposes specific constraints due to the low quantity and the highly degraded genetic material available. These constraints are partially overcome by performing multiple genotyping repetitions to obtain sufficient reliability. I also investigate the causes of DNA degradation in faeces. Among the main degraders, namely bacterial activity, spontaneous hydrolysis, and free-¬DNase activities, the latter was pointed out as the most important according to our experiments. These enzymes degrade DNA very rapidly, and, as a consequence, faeces sampling schemes must be planned preferably in cold and dry weather conditions, allowing for enzyme activity inhibition. The second part of the thesis is a simulation study aiming to assess the capacity of the software Structure to detect population structure in hierarchical models relevant to situations encountered in wild populations, using several genetic markers. The methods implemented in Structure appear efficient in detecting the highest hierarchical structure. The third and fourth parts of the thesis describe the population genetics status of the remaining Jura and Pyrenees populations using 11 microsatellite loci. In either of these populations, no inbreeding nor reduced genetic diversity was detected. Furthermore, the genetic differentiation between patches defined by habitat suitability remains moderate and correlated with geographical distance, suggesting that significant dispersion between patches was at work at least until the last generations. The comparison of diversity indicators with other species or other populations of capercaillie indicate that population in the Jura has retained a large part of its original genetic diversity. These results suggest that the recent decline has had so forth a moderate impact on• genetic factors and that these populations might have retained the potential for long term survival, if the decline is stopped. Finally, in the fifth part, the analysis of relatedness between males participating in the reproduction parade, or lek, indicate that capercaillie males, like has been shown for some other grouse species, gather on leks• among individuals that are more related than the average of the population. This pattern appears to be due to both population structure and kin-association. As a conclusion, this first study relying exclusively on nuclear DNA extracted from faeces has provided novel information that was not available through field observation or classical genetic sampling. No bird has been captured or disturbed, and the results are consistent with other studies of closely related species. However, the size of these populations is approaching thresholds below which long-term survival is unlikely. The persistence of genetic diversity for the forthcoming generations remains therefore bond to adult survival and to the increase of reproduction success.