Decay of linkage disequilibrium within genes across HGDP-CEPH human samples: most population isolates do not show increased LD

Background: It is well known that the pattern of linkage disequilibrium varies between human populations, with remarkable geographical stratification. Indirect association studies routinely exploit linkage disequilibrium around genes, particularly in isolated populations where it is assumed to be higher. Here, we explore both the amount and the decay of linkage disequilibrium with physical distance along 211 gene regions, most of them related to complex diseases, across 39 HGDP-CEPH population samples, focusing particularly on the populations defined as isolates. Within each gene region and population we use r2 between all possible single nucleotide polymorphism (SNP) pairs as a measure of linkage disequilibrium and focus on the proportion of SNP pairs with r2 greater than 0.8.Results: Although the average r2 was found to be significantly different both between and within continental regions, a much higher proportion of r2 variance could be attributed to differences between continental regions (2.8% vs. 0.5%, respectively). Similarly, while the proportion of SNP pairs with r2 > 0.8 was significantly different across continents for all distance classes, it was generally much more homogenous within continents, except in the case of Africa and the Americas. The only isolated populations with consistently higher LD in all distance classes with respect to their continent are the Kalash (Central South Asia) and the Surui (America). Moreover, isolated populations showed only slightly higher proportions of SNP pairs with r2 > 0.8 per gene region than non-isolated populations in the same continent. Thus, the number of SNPs in isolated populations that need to be genotyped may be only slightly less than in non-isolates. Conclusion: The "isolated population" label by itself does not guarantee a greater genotyping efficiency in association studies, and properties other than increased linkage disequilibrium may make these populations interesting in genetic epidemiology.

African signatures of recent positive selection in human FOXI1

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene underhuman-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences.Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans,20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39worldwide human populations from the HGDP-CEPH diversity panel.Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant inEuropeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsbstatistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recentepisode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on theputatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies stronglycorrelated with the absolute geographical latitude of the populations sampled.Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate mightbe related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediatedwater-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.

Impact of antenna correlation on the capacity of multiantenna channels

This paper applies random matrix theory to obtain analytical characterizations of the capacity of correlated multiantenna channels. The analysis is not restricted to the popular separable correlation model, but rather it embraces a more general representation that subsumesmost of the channel models that have been treated in the literature. For arbitrary signal-to-noise ratios (SNR), the characterization is conducted in the regime of large numbers of antennas. For the low- and high-SNR regions, in turn, we uncover compact capacity expansions that are valid for arbitrary numbers of antennas and that shed insight on how antenna correlation impacts the tradeoffs between power, bandwidth and rate.

High-SNR power offset in multiantenna communication

The analysis of the multiantenna capacity in the high-SNR regime has hitherto focused on the high-SNR slope (or maximum multiplexing gain), which quantifies the multiplicative increase as function of the number of antennas. This traditional characterization is unable to assess the impact of prominent channel features since, for a majority of channels, the slope equals the minimum of the number of transmit and receive antennas. Furthermore, a characterization based solely on the slope captures only the scaling but it has no notion of the power required for a certain capacity. This paper advocates a more refined characterization whereby, as function of SNRjdB, the high-SNR capacity is expanded as an affine function where the impact of channel features such as antenna correlation, unfaded components, etc, resides in the zero-order term or power offset. The power offset, for which we find insightful closed-form expressions, is shown to play a chief role for SNR levels of practical interest.

Capacity-achieving input covariance for single-user multi-antenna channels

We characterize the capacity-achieving input covariance for multi-antenna channels known instantaneously at the receiver and in distribution at the transmitter. Our characterization, valid for arbitrary numbers of antennas, encompasses both the eigenvectors and the eigenvalues. The eigenvectors are found for zero-mean channels with arbitrary fading profiles and a wide range of correlation and keyhole structures. For the eigenvalues, in turn, we present necessary and sufficient conditions as well as an iterative algorithm that exhibits remarkable properties: universal applicability, robustness and rapid convergence. In addition, we identify channel structures for which an isotropic input achieves capacity.

Capacity-approaching rate function for layered multiantenna architectures

The simultaneous use of multiple transmit and receive antennas can unleash very large capacity increases in rich multipath environments. Although such capacities can be approached by layered multi-antenna architectures with per-antenna rate control, the need for short-term feedback arises as a potential impediment, in particular as the number of antennas—and thus the number of rates to be controlled—increases. What we show, however, is that the need for short-term feedback in fact vanishes as the number of antennas and/or the diversity order increases. Specifically, the rate supported by each transmit antenna becomes deterministic and a sole function of the signal-to-noise, the ratio of transmit and receive antennas, and the decoding order, all of which are either fixed or slowly varying. More generally, we illustrate -through this specific derivation— the relevance of some established random CDMA results to the single-user multi-antenna problem.

Kinesin spindle protein SiRNA slows tumor progression.

The kinesin spindle protein (KSP), a member of the kinesin superfamily of microtubule-based motors, plays a critical role in mitosis as it mediates centrosome separation and bipolar spindle assembly and maintenance. Inhibition of KSP function leads to cell cycle arrest at mitosis with the formation of monoastral microtubule arrays, and ultimately, to cell death. Several KSP inhibitors are currently being studied in clinical trials and provide new opportunities for the development of novel anticancer therapeutics. RNA interference (RNAi) may represent a powerful strategy to interfere with key molecular pathways involved in cancer. In this study, we have established an efficient method for intratumoral delivery of siRNA. We evaluated short interfering RNA (siRNA) duplexes targeting luciferase as surrogate marker or KSP sequence. To examine the potential feasibility of RNAi therapy, the siRNA was transfected into pre-established lesions by means of intratumor electro-transfer of RNA therapeutics (IERT). This technology allowed cell permeation of the nucleic acids and to efficiently knock down gene expression, albeit transiently. The KSP-specific siRNA drastically reduced outgrowth of subcutaneous melanoma and ovarian cancer lesions. Our results show that intratumoral electro-transfer of siRNA is feasible and KSP-specific siRNA may provide a novel strategy for therapeutic intervention. J. Cell. Physiol. 228: 58-64, 2013. © 2012 Wiley Periodicals, Inc.

Impact of genetic SLC28 transporter and ITPA variants on ribavirin serum level, hemoglobin drop and therapeutic response in patients with HCV infection.

BACKGROUND & AIMS: In the last decade, pegylated interferon-α (PegIFN-α) plus ribavirin (RBV) was the standard treatment of chronic hepatitis C for genotype 1, and it remains the standard for genotypes 2 and 3. Recent studies reported associations between RBV-induced anemia and genetic polymorphisms of concentrative nucleoside transporters such as CNT3 (encoded by SLC28A3) and inosine triphosphatase (encoded by ITPA). We aimed at studying genetic determinants of RBV kinetics, efficacy and treatment-associated anemia. METHODS: We included 216 patients from two Swiss study cohorts (61% HCV genotype 1, 39% genotypes 2 or 3). Patients were analyzed for SLC28A2 single nucleotide polymorphism (SNP) rs11854484, SLC28A3 rs56350726, and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101, and followed for treatment-associated hemoglobin changes and sustained virological response (SVR). In 67 patients, RBV serum levels were additionally measured during treatment. RESULTS: Patients with SLC28A2 rs11854484 genotype TT had higher dosage- and body weight-adjusted RBV levels than those with genotypes TC or CC (p=0.02 and p=0.06 at weeks 4 and 8, respectively). ITPA SNP rs1127354 was associated with hemoglobin drop ≥3 g/dl during treatment, in genotype (relative risk (RR)=2.1, 95% CI 1.3-3.5) as well as allelic analyses (RR=2.0, 95%CI 1.2-3.4). SLC28A3 rs56350726 was associated with SVR in genotype (RR=2.2; 95% CI 1.1-4.3) as well as allelic analyses (RR=2.0, 95% CI 1.1-3.4). CONCLUSIONS: The newly identified association between RBV serum levels and SLC28A2 rs11854484 genotype, as well as the replicated association of ITPA and SLC28A3 genetic polymorphisms with RBV-induced anemia and treatment response, may support individualized treatment of chronic hepatitis C and warrant further investigation in larger studies.

A polymorphism in the leptin gene promoter is associated with anemia in patients with HIV disease.

To study factors associated with anemia and its effect on survival in HIV-infected persons treated with modern combined antiretroviral therapy (cART), we characterized the prevalence of anemia in the Veterans Aging Cohort Study (VACS) and used a candidate gene approach to identify proinflammatory gene single nucleotide polymorphisms (SNPs) associated with anemia in HIV disease. The study comprised 1597 HIV(+) and 865 HIV(-) VACS subjects with DNA, blood, and annotated clinical data available for analysis. Anemia was defined according to World Health Organization criteria (hemoglobin < 13 g/dL and < 12 g/dL in men and women, respectively). The prevalence of anemia in HIV(+) and HIV(-) subjects was 23.1% and 12.9%, respectively. Independent of HIV status, anemia was present in 23.4% and 8% in blacks and whites, respectively. Analysis of our candidate genes revealed that the leptin -2548 G/A SNP was associated with anemia in HIV(+), but not HIV(-), patients, with the AA and AG genotypes significantly predicting anemia (P < .003 and P < .039, respectively, logistic regression). This association was replicated in an independent cohort of HIV(+) women. Our study provides novel insight into the association between genetic variability in the leptin gene and anemia in HIV(+) individuals.

CNOT3 is a modifier of PRPF31 mutations in retinitis pigmentosa with incomplete penetrance.

Heterozygous mutations in the PRPF31 gene cause autosomal dominant retinitis pigmentosa (adRP), a hereditary disorder leading to progressive blindness. In some cases, such mutations display incomplete penetrance, implying that certain carriers develop retinal degeneration while others have no symptoms at all. Asymptomatic carriers are protected from the disease by a higher than average expression of the PRPF31 allele that is not mutated, mainly through the action of an unknown modifier gene mapping to chromosome 19q13.4. We investigated a large family with adRP segregating an 11-bp deletion in PRPF31. The analysis of cell lines derived from asymptomatic and affected individuals revealed that the expression of only one gene among a number of candidates within the 19q13.4 interval significantly correlated with that of PRPF31, both at the mRNA and protein levels, and according to an inverse relationship. This gene was CNOT3, encoding a subunit of the Ccr4-not transcription complex. In cultured cells, siRNA-mediated silencing of CNOT3 provoked an increase in PRPF31 expression, confirming a repressive nature of CNOT3 on PRPF31. Furthermore, chromatin immunoprecipitation revealed that CNOT3 directly binds to a specific PRPF31 promoter sequence, while next-generation sequencing of the CNOT3 genomic region indicated that its variable expression is associated with a common intronic SNP. In conclusion, we identify CNOT3 as the main modifier gene determining penetrance of PRPF31 mutations, via a mechanism of transcriptional repression. In asymptomatic carriers CNOT3 is expressed at low levels, allowing higher amounts of wild-type PRPF31 transcripts to be produced and preventing manifestation of retinal degeneration.

Stratification by smoking status reveals an association of CHRNA5-A3-B4 genotype with body mass index in never smokers.

We previously used a single nucleotide polymorphism (SNP) in the CHRNA5-A3-B4 gene cluster associated with heaviness of smoking within smokers to confirm the causal effect of smoking in reducing body mass index (BMI) in a Mendelian randomisation analysis. While seeking to extend these findings in a larger sample we found that this SNP is associated with 0.74% lower body mass index (BMI) per minor allele in current smokers (95% CI -0.97 to -0.51, P = 2.00 × 10(-10)), but also unexpectedly found that it was associated with 0.35% higher BMI in never smokers (95% CI +0.18 to +0.52, P = 6.38 × 10(-5)). An interaction test confirmed that these estimates differed from each other (P = 4.95 × 10(-13)). This difference in effects suggests the variant influences BMI both via pathways unrelated to smoking, and via the weight-reducing effects of smoking. It would therefore be essentially undetectable in an unstratified genome-wide association study of BMI, given the opposite association with BMI in never and current smokers. This demonstrates that novel associations may be obscured by hidden population sub-structure. Stratification on well-characterized environmental factors known to impact on health outcomes may therefore reveal novel genetic associations.

Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Sequence determinants of a microtubule tip localization signal (MtLS).

Microtubule plus-end-tracking proteins (+TIPs) specifically localize to the growing plus-ends of microtubules to regulate microtubule dynamics and functions. A large group of +TIPs contain a short linear motif, SXIP, which is essential for them to bind to end-binding proteins (EBs) and target microtubule ends. The SXIP sequence site thus acts as a widespread microtubule tip localization signal (MtLS). Here we have analyzed the sequence-function relationship of a canonical MtLS. Using synthetic peptide arrays on membrane supports, we identified the residue preferences at each amino acid position of the SXIP motif and its surrounding sequence with respect to EB binding. We further developed an assay based on fluorescence polarization to assess the mechanism of the EB-SXIP interaction and to correlate EB binding and microtubule tip tracking of MtLS sequences from different +TIPs. Finally, we investigated the role of phosphorylation in regulating the EB-SXIP interaction. Together, our results define the sequence determinants of a canonical MtLS and provide the experimental data for bioinformatics approaches to carry out genome-wide predictions of novel +TIPs in multiple organisms.

Polymorphisms of the " macrophage migration inhibitory factor " gene in infectious diseases

Summary The proinflammatory cytokine macrophage migration inhibitory factor (MIF) has emerged as a central mediator of inflammation and innate immune defense against infections. MIF has been shown to play an important role in the pathogenesis of infectious diseases like sepsis, tuberculosis and autoimmune inflammatory diseases, such as arthritis, inflammatory bowel disease and asthma. Two functional polymorphisms of the MIF gene promoter, a five to eight CATT repeat microsatellite at position -794 and a G/C SNP at position -173, have been associated with increased susceptibility to or severity of autoimmune inflammatory diseases like arthritis, colitis and atopy. The aim of this thesis was to define whether, and if so by which mechanisms, MIF gene polymorphisms influence the susceptibility to or the outcome of one of the most severe and one of the most prevalent infectious diseases: meningococcal sepsis and tuberculosis, respectively. The results of the comparison between 1106 patients suffering from severe meningococcal infections and 434 healthy volunteers showed that carriers of the CATT5-5 genotype were protected from meningococcemia. A transmission disequilibrium test involving 106 families confirmed this association. At baseline and after stimulation with Neisseria meningitidis, the CATT5 MIF promoter drove lower transcriptional activity than the CATT6 or CATT7 alleles in human monocytic cells and whole blood of CATT5-5 healthy individuals tended to produce less MIF than whole blood of CATT6-6 individuals. Beyond, we describe several new MIF gene polymorphisms in Africans. Genotyping the CATT microsatellite and the -173*G/C SNP revealed great genetic diversity in six African ethnic groups. Comparing 471 African tuberculosis cases and 932 matched healthy controls, we observed ethnicity dependent associations of the -173*G/C and the CATT5-8 with susceptibility to or severity of tuberculosis, but confirmation in larger cohorts ìs needed. In conclusion, we report that homozygous carriage of a low expression allele of the MIF gene protects from meningococcal disease. These results support the concept that analyses of MIF genotypes in patients with sepsis may help to classify patients into risk categories and to identify those patients who may benefit from anti-MIF therapeutic strategies.

Association study of 182 candidate genes in anorexia nervosa.

We performed association studies with 5,151 SNPs that were judged as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN) based on location under reported linkage peaks, previous results in the literature (182 candidate genes), brain expression, biological plausibility, and estrogen responsivity. We employed a case-control design that tested each SNP individually as well as haplotypes derived from these SNPs in 1,085 case individuals with AN diagnoses and 677 control individuals. We also performed separate association analyses using three increasingly restrictive case definitions for AN: all individuals with any subtype of AN (All AN: n = 1,085); individuals with AN with no binge eating behavior (AN with No Binge Eating: n = 687); and individuals with the restricting subtype of AN (Restricting AN: n = 421). After accounting for multiple comparisons, there were no statistically significant associations for any individual SNP or haplotype block with any definition of illness. These results underscore the importance of large samples to yield appropriate power to detect genotypic differences in individuals with AN and also motivate complementary approaches involving Genome-Wide Association (GWA) studies, Copy Number Variation (CNV) analyses, sequencing-based rare variant discovery assays, and pathway-based analysis in order to make up for deficiencies in traditional candidate gene approaches to AN.