Consequences of genetic erosion on fitness and phenotypic plasticity in European tree frog populations (Hyla arborea).

The detrimental effects of genetic erosion on small isolated populations are widely recognized contrary to their interactions with environmental changes. The ability of genotypes to plastically respond to variability is probably essential for the persistence of these populations. Genetic erosion impact may be exacerbated if inbreeding affects plastic responses or if their maintenance were at higher phenotypic costs. To understand the interplay 'genetic erosion-fitness-phenotypic plasticity', we experimentally compared, in different environments, the larval performances and plastic responses to predation of European tree frogs (Hyla arborea) from isolated and connected populations. Tadpoles from isolated populations were less performant, but the traits affected were environmental dependant. Heterosis observed in crosses between isolated populations allowed attributing their low fitness to inbreeding. Phenotypic plasticity can be maintained in the face of genetic erosion as inducible defences in response to predator were identical in all populations. However, the higher survival and developmental costs for isolated populations in harsh conditions may lead to an additional fitness loss for isolated populations.

A comparative phenotypic and genomic analysis of C57BL/6J and C57BL/6N mouse strains.

BACKGROUND: The mouse inbred line C57BL/6J is widely used in mouse genetics and its genome has been incorporated into many genetic reference populations. More recently large initiatives such as the International Knockout Mouse Consortium (IKMC) are using the C57BL/6N mouse strain to generate null alleles for all mouse genes. Hence both strains are now widely used in mouse genetics studies. Here we perform a comprehensive genomic and phenotypic analysis of the two strains to identify differences that may influence their underlying genetic mechanisms. RESULTS: We undertake genome sequence comparisons of C57BL/6J and C57BL/6N to identify SNPs, indels and structural variants, with a focus on identifying all coding variants. We annotate 34 SNPs and 2 indels that distinguish C57BL/6J and C57BL/6N coding sequences, as well as 15 structural variants that overlap a gene. In parallel we assess the comparative phenotypes of the two inbred lines utilizing the EMPReSSslim phenotyping pipeline, a broad based assessment encompassing diverse biological systems. We perform additional secondary phenotyping assessments to explore other phenotype domains and to elaborate phenotype differences identified in the primary assessment. We uncover significant phenotypic differences between the two lines, replicated across multiple centers, in a number of physiological, biochemical and behavioral systems. CONCLUSIONS: Comparison of C57BL/6J and C57BL/6N demonstrates a range of phenotypic differences that have the potential to impact upon penetrance and expressivity of mutational effects in these strains. Moreover, the sequence variants we identify provide a set of candidate genes for the phenotypic differences observed between the two strains.

Phenotypic plasticity in the Caribbean sponge Callyspongia vaginalis (Porifera: Haplosclerida)

Sponge morphological plasticity has been a long-standing source of taxonomic difficulty. In the Caribbean, several morphotypes of the sponge Callyspongia vaginalis have been observed. To determine the taxonomic status of three of these morphotypes and their relationship with the congeneric species C. plicifera and C. fallax, we compared the spicule composition, spongin fiber skeleton and sequenced fragments of the mitochondrial genes 16S and COI and nuclear genes 28S and 18S ribosomal RNA. Phylogenetic analyses with ribosomal markers 18S and 28S rRNA confirmed the position of our sequences within the Callyspongiidae. None of the genetic markers provided evidence for consistent differentiation among the three morphotypes of C. vaginalis and C. fallax, and only C. plicifera stood as a distinct species. The 16S mtDNA gene was the most variable molecular marker for this group, presenting a nucleotide variability (π = 0.024) higher than that reported for COI. Unlike recent studies for other sponge genera, our results indicate that species in the genus Callyspongia maintain a high degree of phenotypic plasticity, and that morphological characteristics may not reflect reproductive boundaries in C. vaginalis.

When genetic load does not correlate with phenotypic spectrum: lessons from the GnRH receptor (GNRHR).

CONTEXT: A broad spectrum of GnRH-deficient phenotypes has been identified in individuals with both mono- and biallelic GNRHR mutations. OBJECTIVE: The objective of the study was to determine the correlation between the severity of the reproductive phenotype(s) and the number and functional severity of rare sequence variants in GNRHR. SUBJECTS: Eight hundred sixty-three probands with different forms of GnRH deficiency, 46 family members and 422 controls were screened for GNRHR mutations. The 70 subjects (32 patients and 38 family members) harboring mutations were divided into four groups (G1-G4) based on the functional severity of the mutations (complete or partial loss of function) and the number of affected alleles (monoallelic or biallelic) with mutations, and these classes were mapped on their clinical phenotypes. RESULTS: The prevalence of heterozygous rare sequence variants in GNRHR was significantly higher in probands vs. controls (P < 0.01). Among the G1-G3 groups (homozygous subjects with successively decreasing severity and number of mutations), the hypogonadotropic phenotype related to their genetic load. In contrast, subjects in G4, with only monoallelic mutations, demonstrated a greater diversity of clinical phenotypes. CONCLUSIONS: In patients with GnRH deficiency and biallelic mutations in GNRHR, genetic burden defined by severity and dose is associated with clinical phenotype. In contrast, for patients with monoallelic GNRHR mutations this correlation does not hold. Taken together, these data indicate that as-yet-unidentified genetic and/or environmental factors may combine with singly mutated GNRHR alleles to produce reproductive phenotypes.

Therapeutic drug monitoring and metabolites profil analysis of antiviral drugs

The widespread use of combination antiretroviral therapy (ARVs) has considerably improved the prognosis of patients infected with HIV. Conversely, considerable advances have been recently realized for the therapy of hepatitis C infection with the recent advent of potent new anti-HCV drugs that allow an increasing rate HCV infection cure. Despite their overall efficacy, a significant number of patients do not achieve or maintain adequate clinical response, defined as an undetectable viral load for HIV, and a sustained virological response (or cure) in HCV infection. Treatment failure therefore still remains an important issue besides drugs toxicities and viral resistance which is not uncommon in a significant percentage of patients who do not reach adequate virological suppression. The reasons of variability in drug response are multifactorial and apart from viral genetics, other factors such as environmental factors, drug- drug interactions, and imperfect compliance may have profound impact on antiviral drugs' clinical response. The possibility of measuring plasma concentration of antiviral drugs enables to guide antiviral drug therapy and ensure optimal drug exposure. The overall objective of this research was to widen up the current knowledge on pharmacokinetic and pharmacogenetic factors that influence the clinical response and toxicity of current and newly approved antiretroviral and anti-HCV drugs. To that endeavour, analytical methods using liquid chromatography coupled with tandem mass spectrometry have been developed and validated for the precise and accurate measurement of new antiretroviral and anti-HCV drugs . These assays have been applied for the TDM of ARVs and anti-HCV in patients infected with either HIV or HCV respectively, and co-infected with HIV- HCV. A pharmacokinetic population model was developed to characterize inter and intra-patient variability of rilpivirine, the latest marketed Non Nucleoside Reverse transcriptase (NNRTI) Inhibitor of HIVand to identify genetic and non genetic covariates influencing rilpivirine exposure. None of the factors investigated so far showed however any influence of RPV clearance. Importantly, we have found that the standard daily dosage regimen (25 mg QD) proposed for rilpivirine results in concentrations below the proposed therapeutic target in about 40% of patients. In these conditions, virologie escape is a potential risk that remains to be further investigated, notably via the TDM approach that can be a useful tool to identify patients who are at risk for being exposed to less than optimal levels of rilpivirine in plasma. Besides the last generation NNRTI rilpivirine, we have studied efavirenz, the major NNRTI clinically used so far. Namely for efavirenz, we aimed at identifying a potential new marker of toxicity that may be incriminated for the neuropsychological sides effects and hence discontinuation of efavirenz therapy. To that endeavour, a comprehensive analysis of phase I and phase II metabolites profiles has been performed in plasma, CSF and in urine from patients under efavirenz therapy. We have found that phase II metabolites of EFV constitute the major species circulating in blood, sometimes exceeding the levels of the parent drug efavirenz. Moreover we have identified a new metabolite of efavirenz in humans, namely the 8-OH-EFV- sulfate which is present at high concentrations in all body compartments from patients under efavirenz therapy. These investigations may open the way to possible alternate phenotypic markers of efavirenz toxicity. Finally, the specific influence of P-glycoprotein on the cellular disposition of a series ARVs (NNRTIs and Pis] has been studies in in vitro cell systems using the siRNA silencing approach. -- Depuis l'introduction de la thérapie antirétrovirale (ARVs) la morbidité et la mortalité liées au VIH ont considérablement diminué. En parallèle le traitement contre le virus de l'hépatite C (VHC) a connu récemment d'énormes progrès avec l'arrivée de nouveaux médicaments puissants, ce qui a permis une augmentation considérable de la guérison de l'infection par le VHC. En dépit de l'efficacité de ces traitements antiviraux, les échecs thérapeutiques ainsi que les effets secondaires des traitements restent un problème important. Une réponse imparfaite ou la toxicité du traitement est certainement multifactorielle. Le suivi thérapeutique des médicaments [Therapeutic Drug Monitoring TDM) à travers la mesure des concentrations plasmatiques constitue une approche importante pour guider le traitement médicamenteux et de s'assurer que les patients sont exposés à des concentrations optimales des médicaments dans le sang, et puissent tirer tout le bénéfice potentiel du traitement. L'objectif global de cette thèse était d'étudier les facteurs pharmacocinétiques et pharmacogénétiques qui influencent l'exposition des médicaments antiviraux (ARVs et anti- VHC) récemment approuvés. A cet effet, des méthodes de quantification des concentrations plasmatiques des médicaments antirétroviraux, anti-VHC ainsi que pour certains métabolites ont été développées et validées en utilisant la Chromatographie liquide couplée à la spectrométrie de masse tandem. Ces méthodes ont été utilisées pour le TDM des ARVs et pour les agents anti-VHC chez les patients infectés par le VIH, et le VHC, respectivement, mais aussi chez les patients co-infectés par le VIH-VHC. Un modèle de pharmacocinétique de population a été développé pour caractériser la variabilité inter-et intra-patient du médicament rilpivirine, un inhibiteur non nucléosidique de la transcriptase de VIH et d'identifier les variables génétiques et non génétiques influençant l'exposition au médicament. Aucun des facteurs étudiés n'a montré d'influence notable sur la clairance de la rilpivirine. Toutefois, la concentration résiduelle extrapolée selon le modèle de pharmacocinétique de population qui a été développé, a montré qu'une grande proportion des patients présente des concentrations minimales inférieures à la cible thérapeutique proposée. Dans ce contexte, la relation entre les concentrations minimales et l'échappement virologique nécessite une surveillance étroite des taux sanguins des patients recevant de la rilpivirine. A cet effet, le suivi thérapeutique est un outil important pour l'identification des patients à risque soient sous-exposés à lai rilpivirine. Pour identifier de nouveaux marqueurs de la toxicité qui pourraient induire l'arrêt du traitement, le profil des métabolites de phase I et de phase II a été étudié dans différentes matrices [plasma, LCR et urine) provenant de patients recevant de l'efavirenz. Les métabolites de phase II, qui n'avaient à ce jour jamais été investigués, constituent les principales espèces présentes dans les matrices étudiées. Au cours de ces investigations, un nouveau métabolite 8- OH-EFV-sulfate a été identifié chez l'homme, et ce dernier est. présent à des concentrations importantes. L'influence de certains facteurs pharmacogénétique des patients sur le profil des métabolites a été étudiée et ouvre la voie à de possibles nouveaux marqueurs phénotypiques alternatifs qui pourraient possiblement mieux prédire la toxicité associée au traitement par l'efavirenz. Finalement, nous nous sommes intéressés à étudier dans un modèle in vitro certains facteurs, comme la P-glycoprotéine, qui influencent la disposition cellulaire de certains médicaments antirétroviraux, en utilisant l'approche par la technologie du siRNA permettant de bloquer sélectivement l'expression du gène de cette protéine d'efflux des médicaments. -- Depuis l'introduction de la thérapie antiretrovirale (ARVs] la morbidité et la mortalité liées au VIH ont considérablement diminué. En parallèle le traitement contre le virus de l'hépatite C (VHC) a connu récemment d'énormes progrès avec l'arrivée de nouveaux médicaments puissants, ce qui a permis une augmentation considérable de la guérison de l'infection par le VHC. En dépit de l'efficacité de ces traitements antiviraux, les échecs thérapeutiques ainsi que les effets secondaires des traitements restent un problème important. Il a pu être démontré que la concentration de médicament présente dans l'organisme est corrélée avec l'efficacité clinique pour la plupart des médicaments agissant contre le VIH et contre le VHC. Les médicaments antiviraux sont généralement donnés à une posologie fixe et standardisée, à tous les patients, il existe cependant une importante variabilité entre les concentrations sanguines mesurées chez les individus. Cette variabilité peut être expliquée par plusieurs facteurs démographiques, environnementaux ou génétiques. Dans ce contexte, le suivi des concentrations sanguines (ou Therapeutic Drug Monitoring, TDM) permet de contrôler que les patients soient exposés à des concentrations suffisantes (pour bloquer la réplication du virus dans l'organisme) et éviter des concentrations excessives, ce qui peut entraîner l'apparition d'intolérence au traitement. Le but de ce travail de thèse est d'améliorer la compréhension des facteurs pharmacologiques et génétiques qui peuvent influencer l'efficacité et/ou la toxicité des médicaments antiviraux, dans le but d'améliorer le suivi des patients. A cet effet, des méthodes de dosage très sensibles et ont été mises au point pour permettre de quantifier les médicaments antiviraux dans le sang et dans d'autres liquides biologiques. Ces méthodes de dosage sont maintenant utilisées d'une part dans le cadre de la prise en charge des patients en routine et d'autre part pour diverses études cliniques chez les patients infectés soit par le HIV, le HCV ou bien coinfectés par les deux virus. Une partie de ce travail a été consacrée à l'investigation des différents facteurs démographiques, génétiques et environnementaux qui pourraient l'influencer la réponse clinique à la rilpivirine, un nouveau médicament contre le VIH. Toutefois, parmi tous les facteurs étudiés à ce jour, aucun n'a permis d'expliquer la variabilité de l'exposition à la rilpivirine chez les patients. On a pu cependant observer qu'à la posologie standard recommandée, un pourcentage relativement élevé de patients pourrait présenter des concentrations inférieures à la concentration sanguine minimale actuellement proposée. Il est donc utile de surveiller étroitement les concentrations de rilpivirine chez les patients pour identifier sans délai ceux qui risquent d'être sous-exposés. Dans l'organisme, le médicament subit diverses transformations (métabolisme) par des enzymes, notamment dans le foie, il est transporté dans les cellules et tissus par des protéines qui modulent sa concentration au site de son action pharmacologique. A cet effet, différents composés (métabolites) produits dans l'organisme après l'administration d'efavirenz, un autre médicament anti-VIH, ont été étudiés. En conclusion, nous nous sommes intéressés à la fois aux facteurs pharmacologiques et génétiques des traitements antiviraux, une approche qui s'inscrit dans l'optique d'une stratégie globale de prise en charge du patient. Dans ce contexte, le suivi des concentrations sanguines de médicaments constitue une des facettes du domaine émergent de la Médecine Personnalisée qui vise à maximiser le bénéfice thérapeutique et le profil de tolérance des médicaments antiviraux

Phenotypic and molecular assessment of seven patients with 6p25 deletion syndrome: relevance to ocular dysgenesis and hearing impairment.

BACKGROUND: Thirty-nine patients have been described with deletions involving chromosome 6p25. However, relatively few of these deletions have had molecular characterization. Common phenotypes of 6p25 deletion syndrome patients include hydrocephalus, hearing loss, and ocular, craniofacial, skeletal, cardiac, and renal malformations. Molecular characterization of deletions can identify genes that are responsible for these phenotypes. METHODS: We report the clinical phenotype of seven patients with terminal deletions of chromosome 6p25 and compare them to previously reported patients. Molecular characterization of the deletions was performed using polymorphic marker analysis to determine the extents of the deletions in these seven 6p25 deletion syndrome patients. RESULTS: Our results, and previous data, show that ocular dysgenesis and hearing impairment are the two most highly penetrant phenotypes of the 6p25 deletion syndrome. While deletion of the forkhead box C1 gene (FOXC1) probably underlies the ocular dysgenesis, no gene in this region is known to be involved in hearing impairment. CONCLUSIONS: Ocular dysgenesis and hearing impairment are the two most common phenotypes of 6p25 deletion syndrome. We conclude that a locus for dominant hearing loss is present at 6p25 and that this locus is restricted to a region distal to D6S1617. Molecular characterization of more 6p25 deletion patients will aid in refinement of this locus and the identification of a gene involved in dominant hearing loss.

Population-based registration of phenotypic anatomic and familial data for melanoma: results from a Swiss multicentric study

CONTEXT AND OBJECTIVES: A multicentric study was set up to assess the feasibility for Swiss cancer registries of actively retrieving 3 additional variables of epidemiological and a etiological relevance for melanoma, and of potential use for the evaluation of prevention campaigns. MATERIAL AND METHODS: The skin type, family history of melanoma and precise anatomical site were retrieved for melanoma cases registered in 5 Swiss cantons (Neuchâtel, St-Gall and Appenzell, Vaud and Wallis) over 3 to 6 consecutive years (1995-2002). Data were obtained via a short questionnaire administered by the physicians - mostly dermatologists - who originally excised the lesions. As the detailed body site was routinely collected in Ticino, data from this Cancer Registry were included in the body site analysis. Relative melanoma density (RMD) was computed by the ratio of observed to expected numbers of melanomas allowing for body site surface areas, and further adjusted for site-specific melanocyte density. RESULTS: Of the 1,645 questionnaires sent, 1,420 (86.3%) were returned. The detailed cutaneous site and skin type were reliably obtained for 84.7% and 78.7% of questionnaires, and family history was known in 76% of instances. Prevalence of sun-sensitive subjects and patients with melanoma affected first-degree relatives, two target groups for early detection and surveillance campaigns were 54.1% and 3.4%, respectively. After translation into the 4th digit of the International Classification of Diseases for Oncology, the anatomical site codes from printed (original information) and pictorial support (body chart from the questionnaire) concurred for 94.6% of lesions. Discrepancies occurred mostly for lesions on the upper, outer part of the shoulder for which the clinician's textual description was "shoulder blade". This differential misclassification suggests under-estimation by about 10% of melanomas of the upper limbs and an over-estimation of 5% for truncal melanomas. Sites of highest melanoma risk were the face, the shoulder and the upper arm for sexes, the back for men and the leg for women. Three major features of this series were: (1) an unexpectedly high RMD for the face in women (6.2 vs 4.2 in men), (2) the absence of a male predominance for melanomas on the ears, and (3) for the upper limbs, a steady gradient of increasing melanoma density with increasing proximity to the trunk, regardless of sex. DISCUSSION AND CONCLUSION: The feasibility of retrieving the skin type, the precise anatomical location and family history of melanoma in a reliable manner was demonstrated thanks to the collaboration of Swiss dermatologists. Use of a schematic body drawing improves the quality of the anatomical site data and facilitate the reporting task of doctors. Age and sex patterns of RMD paralleled general indicators of sun exposure and behaviour, except for the hand (RMD=0.2). These Swiss results support some site or sun exposure specificity in the aetiology of melanoma.

The mean-variance model from the inverse of the variance-covariance matrix

En este artículo, a partir de la inversa de la matriz de varianzas y covarianzas se obtiene el modelo Esperanza-Varianza de Markowitz siguiendo un camino más corto y matemáticamente riguroso. También se obtiene la ecuación de equilibrio del CAPM de Sharpe.

Impact of incorrect assumptions on the covariance structure of random effects and/or residuals in nonlinear mixed models for repeated measures data

In this paper we analyse, using Monte Carlo simulation, the possible consequences of incorrect assumptions on the true structure of the random effects covariance matrix and the true correlation pattern of residuals, over the performance of an estimation method for nonlinear mixed models. The procedure under study is the well known linearization method due to Lindstrom and Bates (1990), implemented in the nlme library of S-Plus and R. Its performance is studied in terms of bias, mean square error (MSE), and true coverage of the associated asymptotic confidence intervals. Ignoring other criteria like the convenience of avoiding over parameterised models, it seems worst to erroneously assume some structure than do not assume any structure when this would be adequate.

Spatial weights : constructing weight-compatible exchange matrices from proximity matrices

Exchange matrices represent spatial weights as symmetric probability distributions on pairs of regions, whose margins yield regional weights, generally well-specified and known in most contexts. This contribution proposes a mechanism for constructing exchange matrices, derived from quite general symmetric proximity matrices, in such a way that the margin of the exchange matrix coincides with the regional weights. Exchange matrices generate in turn diffusive squared Euclidean dissimilarities, measuring spatial remoteness between pairs of regions. Unweighted and weighted spatial frameworks are reviewed and compared, regarding in particular their impact on permutation and normal tests of spatial autocorrelation. Applications include tests of spatial autocorrelation with diagonal weights, factorial visualization of the network of regions, multivariate generalizations of Moran's I, as well as "landscape clustering", aimed at creating regional aggregates both spatially contiguous and endowed with similar features.

Congenital disorder of glycosylation type Id (CDG Id) : phenotypic, biochemical and molecular characterization of a new patient

Rapport de synthèseLes troubles de la glycosylation (Congenital Disorders of Glycosylation, CDG) regroupent une famille de maladies multi-systémiques héréditaires causées par des défauts dans la synthèse de glycoconjugés. La glycosylation est une réaction enzymatique consistant à lier de façon covalente un glucide à une chaîne peptidique ou une protéine. Il existe deux types de glycosylation. La N-gjycosylation est l'addition de glucides aux chaînes peptidiques en croissance dès leur entrée dans la lumière du réticulum endoplasmique. Elle s'effectue sur les futures glycoprotéines membranaires et conduit à des chaînes de sucres courtes et ramifiées. La O-glycosylation est l'addition de glucides au niveau des résidus hydroxylés des acides aminés sérine et thréonine des chaînes peptidiques déjà présentes dans la lumière de l'appareil de Golgi. Elle est, dans la plupart des cas, effectuée sur îes protéoglycanes et conduit à des chaînes de sucres longues et non ramifiées. La classification des CDG repose sur le niveau de l'étape limitante de la glycosylation. Les CDG de type 1, plus fréquents, regroupent les déficits enzymatiques se situant en amont du transfert de Poligosaccharide sur la chaîne peptidique. Les CDG de type 2 regroupent ceux ayant lieu en aval de ce transfert. Parmi les nombreux différents sous-types de CDG, le CDG de type ld est causé par une anomalie de la mannosyltransferase, enzyme codée par le gène ALG3 (chromosome 3q27). Jusqu'à ce jour, six patients atteints de CDG ld ont été reportés dans la littérature. Notre travail a permis de décrire un septième patient et d'affiner les caractéristiques cliniques, biologiques, neuroradiologiques et moléculaires du CDG ld. Notre patient est notamment porteur d'une nouvelle mutation de type missense sur le gène ALG3. Tous les patients atteints de CDG ld présentent une encéphalopathie progressive avec microcéphalie, retard psychomoteur sévère et épilepsie. Une ostéopénie marquée est présente chez certains patients. Elle est parfois sous diagnostiquée et révélée uniquement lors de fracture pathologique. Les patients atteints de CDG ld présentent également des traits dysmorphiques typiques, mais aucune atteinte multi-systémique ou anomalie biologique spécifique n'est retrouvée telle que dans les autres types de CDG. Le dépistage biochimique des troubles de la glycosylation se fait par une analyse simple et peu coûteuse qui est l'analyse de la transferrine sérique par isoelectrofocusing ou par électrophorèse capillaire. Un tel dépistage devrait être effectué chez tout patient présentant une encéphalopathie d'origine indéterminée, et cela même en l'absence d'atteinte multi- systémique. Notre travail a été publié sous forme d'article de type « short report », peer-reviewed, dans le Journal of Inherited Metabolic Diseases. Le Journal est une révue spécialisée du domaine des erreirs innées du métabolisme. S'agissant d'un seul patient rapporté, l'article ne montre que très synthétiquement le travail effectué, Pour cette raison un complément à l'article avec matériel, méthodes et résultats figure ci-après et concerne la partie de recherche moléculaire de notre travail. La doctorante a non seulement encadré personnellement le patient au niveau clinique et biochimique, mais a plus particulièrement mis au point l'analyse moléculaire du gène ALG3 dans le laboratoire de Pédiatrie Moléculaire pour la première fois ; cela a impliqué l'étude du gène, le choix des oligonucleotides et l'optimisation des réactions d'amplification et séquençage.

Density matrices and momentum distributions in 3He-4He mixtures

We report variational calculations, in the hypernetted-chain (HNC)-Fermi-HNC scheme, of one-body density matrices and one-particle momentum distributions for 3He-4He mixtures described by a Jastrow correlated wave function. The 4He condensate fractions and the 3He strength poles are examined and compared with the Monte Carlo available results. The agreement has been found to be very satisfactory. Their density dependence is also studied.