Sleep deprivation effects on circadian clock gene expression in the cerebral cortex parallel electroencephalographic differences among mouse strains.

Sleep deprivation (SD) results in increased electroencephalographic (EEG) delta power during subsequent non-rapid eye movement sleep (NREMS) and is associated with changes in the expression of circadian clock-related genes in the cerebral cortex. The increase of NREMS delta power as a function of previous wake duration varies among inbred mouse strains. We sought to determine whether SD-dependent changes in circadian clock gene expression parallel this strain difference described previously at the EEG level. The effects of enforced wakefulness of incremental durations of up to 6 h on the expression of circadian clock genes (bmal1, clock, cry1, cry2, csnk1epsilon, npas2, per1, and per2) were assessed in AKR/J, C57BL/6J, and DBA/2J mice, three strains that exhibit distinct EEG responses to SD. Cortical expression of clock genes subsequent to SD was proportional to the increase in delta power that occurs in inbred strains: the strain that exhibits the most robust EEG response to SD (AKR/J) exhibited dramatic increases in expression of bmal1, clock, cry2, csnkIepsilon, and npas2, whereas the strain with the least robust response to SD (DBA/2) exhibited either no change or a decrease in expression of these genes and cry1. The effect of SD on circadian clock gene expression was maintained in mice in which both of the cryptochrome genes were genetically inactivated. cry1 and cry2 appear to be redundant in sleep regulation as elimination of either of these genes did not result in a significant deficit in sleep homeostasis. These data demonstrate transcriptional regulatory correlates to previously described strain differences at the EEG level and raise the possibility that genetic differences underlying circadian clock gene expression may drive the EEG differences among these strains.

A variant in the gene FUT9 is associated with susceptibility to placental malaria infection

Malaria in pregnancy forms a substantial part of the worldwide burden of malaria, with an estimated annual death toll of up to 200,000 infants, as well as increased maternal morbidity and mortality. Studies of genetic susceptibility to malaria have so far focused on infant malaria, with only a few studies investigating the genetic basis of placental malaria, focusing only on a limited number of candidate genes. The aim of this study therefore was to identify novel host genetic factors involved in placental malaria infection. To this end we carried out a nested case-control study on 180 Mozambican pregnant women with placental malaria infection, and 180 controls within an intervention trial of malaria prevention. We genotyped 880 SNPs in a set of 64 functionally related genes involved in glycosylation and innate immunity. A SNP located in the gene FUT9, rs3811070, was significantly associated with placental malaria infection (OR = 2.31, permutation p-value = 0.028). Haplotypic analysis revealed a similarly strong association of a common haplotype of four SNPs including rs3811070. FUT9 codes for a fucosyl-transferase that is catalyzing the last step in the biosynthesis of the Lewis-x antigen, which forms part of the Lewis blood group-related antigens. These results therefore suggest an involvement of this antigen in the pathogenesis of placental malaria infection.

Linguistic and maternal genetic diversity are not correlated in native mexicans

Mesoamerica, defined as the broad linguistic and cultural area from middle southern Mexico to Costa Rica, might have played a pivotal role during the colonization of theAmerican continent. It has been suggested that the Mesoamerican isthmus could have played an important role in severely restricting prehistorically gene flow between North and SouthAmerica. Although the Native American component has been already described in admixedMexican populations, few studies have been carried out in native Mexican populations. In thisstudy we present mitochondrial DNA (mtDNA) sequence data for the first hypervariable region (HVR-I) in 477 unrelated individuals belonging to eleven different native populations from Mexico. Almost all the Native Mexican mtDNAs could be classified into the four pan-Amerindian haplogroups (A2, B2, C1 and D1); only three of them could be allocated to the rare Native American lineage D4h3. Their haplogroup phylogenies are clearly star-like, as expected from relatively young populations that have experienced diverse episodes of genetic drift (e.g. extensive isolation, genetic drift and founder effects) and posterior population expansions. In agreement with this observation is the fact that Native Mexican populations show a high degree of heterogeneity in their patterns of haplogroup frequencies. HaplogroupX2a was absent in our samples, supporting previous observations where this clade was only detected in the American northernmost areas. The search for identical sequences in the American continent shows that, although Native Mexican populations seem to show a closer relationship to North American populations, they cannot be related to a single geographical region within the continent. Finally, we did not find significant population structure on the maternal lineages when considering the four main and distinct linguistic groups represented in our Mexican samples (Oto-Manguean, Uto-Aztecan, Tarascan, and Mayan), suggesting that genetic divergence predates linguistic diversification in Mexico.

Effect of oral calcium carbonate on aortic calcification in apolipoprotein E-deficient (apoE-/-) mice with chronic renal failure.

BACKGROUND: In chronic kidney disease (CKD) patients, the intake of calcium-based phosphate binders is associated with a marked progression of coronary artery and aortic calcification, in contrast to patients receiving calcium-free phosphate binders. The aim of this study was to reexamine the role of calcium carbonate in vascular calcification and to analyse its effect on aortic calcification-related gene expression in chronic renal failure (CRF). METHODS: Mice deficient in apolipoprotein E underwent either sham operation or subtotal nephrectomy to create CRF. They were then randomly assigned to one of the three following groups: a control non-CRF group and a CRF group fed on standard diet, and a CRF group fed on calcium carbonate enriched diet, for a period of 8 weeks. Aortic atherosclerotic plaque and calcification were evaluated using quantitative morphologic image processing. Aortic gene and protein expression was examined using immunohistochemistry and Q-PCR methods. RESULTS: Calcium carbonate supplementation was effective in decreasing serum phosphorus but was associated with a higher serum calcium concentration. Compared with standard diet, calcium carbonate enriched diet unexpectedly induced a significant decrease of both plaque (p<0.05) and non-plaque-associated calcification surface (p<0.05) in CRF mice. It also increased osteopontin (OPN) protein expression in atherosclerotic lesion areas of aortic root. There was also a numerical increase in OPN and osteoprotegerin gene expression in total thoracic aorta but the difference did not reach the level of significance. Finally, calcium carbonate did not change the severity of atherosclerotic lesions. CONCLUSION: In this experimental model of CRF, calcium carbonate supplementation did not accelerate but instead decreased vascular calcification. If our observation can be extrapolated to humans, it appears to question the contention that calcium carbonate supplementation, at least when given in moderate amounts, necessarily enhances vascular calcification. It is also compatible with the hypothesis of a preponderant role of phosphorus over that of calcium in promoting vascular calcification in CRF.

Genetically programmed autoinducer destruction reduces virulence gene expression and swarming motility in Pseudomonas aeruginosa PAO1.

Virulence in the opportunistic human pathogen Pseudomonas aeruginosa is controlled by cell density via diffusible signalling molecules ('autoinducers') of the N-acylhomoserine lactone (AHL) type. Two Bacillus sp. isolates (A23 and A24) with AHL-degrading activity were identified among a large collection of rhizosphere bacteria. From isolate A24 a gene was cloned which was similar to the aiiA gene, encoding an AHL lactonase in another Bacillus strain. Expression of the aiiA homologue from isolate A24 in P. aeruginosa PAO1 reduced the amount of the quorum sensing signal N-oxododecanoyl-L-homoserine lactone and completely prevented the accumulation of the second AHL signal, N-butyryl-L-homoserine lactone. This strongly reduced AHL content correlated with a markedly decreased expression and production of several virulence factors and cytotoxic compounds such as elastase, rhamnolipids, hydrogen cyanide and pyocyanin, and strongly reduced swarming. However, no effect was observed on flagellar swimming or on twitching motility, and aiiA expression did not affect bacterial adhesion to a polyvinylchloride surface. In conclusion, introduction of an AHL degradation gene into P. aeruginosa could block cell-cell communication and exoproduct formation, but failed to interfere with surface colonization.

Effects of epigenetic regulatory elements on telomeric gene expression

Transgene expression in eukaryotic cells strongly depends on the locus of integration in the host genome and often results in limited transcription level because of unfavorable chromatin structure at the integration site. Epigenetic regulators are DNA sequences which are believed to act on the chromatin structure and may protect transgenes from this so-called position effect. Despite being extensively used to increase transgene expression, the mechanism of action of many of these elements remains largely unknown. Here we evaluated different epigenetic regulatory DNA elements for their ability to protect transgene transcription at telomeres, a defined chromatin environment associated to low or inconsistent expression caused by the Telomere Position Effect (TPE). For the assessment of the effects of epigenetic regulators at telomeres, a novel dual reporter system had to be designed. Telomeric integration of the newly-developed dual reporter system carrying different epigenetic regulators showed that MARs (Matrix Attachment Regions), a UCOE (Ubiquitous Chromatin-Opening Element) or the chicken cHS4 insulator have strong barrier activity which prevented TPE from spreading toward the centromere, resulting in stable and in some cases increased expression of a telomeric-distal reporter gene. In addition, MARs and STAR element 40 resulted in an increase of cells expressing the telomeric-proximal reporter gene, suggesting also an anti-silencing effect. Chromatin immunoprecipitation assays revealed that at telomeres MARs actively promote the deposition of euchromatic histone marks, especially acetylation of both histone H3 and H4, which might be involved in MARs' barrier and transcriptional activator activities. Differently, the chromatin in proximity of the UCOE element was depleted of several repressive chromatin marks, such as trimethylated lysine 9 and lysine 27 on histone H3 and trimethylated lysine 20 of histone H4, possibly favoring the preservation of an open chromatin structure at the integration site. We conclude that epigenetic regulatory elements that may be used to enhance and sustain transgene expression have all a specific epigenetic signature which might be at the basis of their mechanism of action, and that a combination of different classes of epigenetic regulators might be advantageous when high levels of protein expression are required. - L'expression des transgènes dans les cellules eucaryotes est fortement influencée par leur site d'intégration dans le génome. En effet, une structure chromatinienne défavorable au niveau du site d'intégration peut fortement limiter le degré d'expression d'un transgène. Il existe toutefois des séquences d'ADN qui, en agissant sur la structure de la chromatine, permettent de limiter cet effet de position et, par conséquent, de promouvoir l'expression soutenue d'un transgène. Ces éléments génomiques, connus comme régulateurs épigénétiques, sont largement utilisés dans plusieurs domaines où une expression élevée et soutenue est requise, malgré un mode de fonctionnement parfois méconnu. Dans cette étude, j'ai évalué la capacité de différents régulateurs épigénétiques à protéger la transcription de transgènes au niveau des télomères, régions chromatiniennes bien définies qui ont été associées à un fort effet de silençage, connu comme «effet de position télomérique». Pour cela, un nouveau système à deux gènes rapporteurs a été développé. Lorsque des MARs (Matrix Attachment Regions, séquences d'ADN pouvant s'associer à la matrice nucléaire), un UCOE (Ubiquitous Chromatin-Opening Element, élément pouvant ouvrir la chromatine) ou l'isolateur génétique cHS4 (dérivé du locus de la β-globine de poulet) sont placés entre les deux gènes rapporteurs, une forte activité barrière bloquant la propagation de la chromatine répressive télomérique est observée, résultant en un plus grand nombre de cellules exprimant le gène télomérique-distal. D'autre part, une augmentation du nombre de cellules exprimant le gène télomérique-proximal, observée en présence des éléments MAR et STAR 40 (Stabilizing Anti-Repressor element 40, un élément pouvant prévenir la répression génique), suggère aussi un faible effet anti-silençage pour ces éléments. Des expériences d'immunoprécipitation de la chromatine démontrent qu'au télomère, les MARs favorisent l'assemblage de marqueurs de la chromatine active, surtout l'acétylation des histones H3 et H4, qui pourraient être à la base de l'activité barrière et de celle d'activateur transcriptionel. Différemment, la chromatine à proximité de l'élément UCOE est particulièrement pauvre en marqueurs de la chromatine silencieuse, comme la trimethylation des lysines 9 et 27 de l'histone H3, ainsi que la trimethylation de la lysine 20 de l'histone H4. Cela suggère que UCOE pourrait préserver une structure chromatinienne ouverte au site d'intégration, favorisant l'expression des gènes à sa proximité. En conclusion, les régulateurs épigénétiques analysés lors de cette étude ont tous montré une signature épigénétique spécifique qui pourrait être à la base de leurs mécanismes de fonctionnement, suggérant aussi qu'une utilisation d'éléments épigénétiques de classe différente dans un même vecteur d'expression pourrait être avantageuse lorsque de hauts et soutenus niveaux d'expression sont nécessaires.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin.

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

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.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics.

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Biotic Factors Affecting Expression of the 2,4-Diacetylphloroglucinol Biosynthesis Gene phlA in Pseudomonas fluorescens Biocontrol Strain CHA0 in the Rhizosphere.

ABSTRACT Production of the polyketide antimicrobial metabolite 2,4-diacetyl-phloroglucinol (DAPG) is a key factor in the biocontrol activity of Pseudomonas fluorescens CHA0. Strain CHA0 carrying a translational phlA'-'lacZ fusion was used to monitor expression of the phl biosynthetic genes in vitro and in the rhizosphere. Expression of the reporter gene accurately reflected actual production of DAPG in vitro and in planta as determined by direct extraction of the antimicrobial compound. In a gnotobiotic system containing a clay and sand-based artificial soil, reporter gene expression was significantly greater in the rhizospheres of two monocots (maize and wheat) compared with gene expression in the rhizospheres of two dicots (bean and cucumber). We observed this host genotype effect on bacterial gene expression also at the level of cultivars. Significant differences were found among six additional maize cultivars tested under gnotobiotic conditions. There was no difference between transgenic maize expressing the Bacillus thuringiensis insecticidal gene cry1Ab and the near-isogenic parent line. Plant age had a significant impact on gene expression. Using maize as a model, expression of the phlA'-'lacZ reporter gene peaked at 24 h after planting of pregerminated seedlings, and dropped to a fourth of that value within 48 h, remaining at that level throughout 22 days of plant growth. Root infection by Pythium ultimum stimulated bacterial gene expression on both cucumber and maize, and this was independent of differences in rhizosphere colonization on these host plants. To our knowledge, this is the first comprehensive evaluation of how biotic factors that commonly confront bacterial inoculants in agricultural systems (host genotype, host age, and pathogen infection) modulate the expression of key biocontrol genes for disease suppression.

Asymmetric and differential gene introgression at a contact zone between two highly divergent lineages of field voles (Microtus agrestis).

Secondary contact zones have the potential to shed light on the mode and rate at which reproductive isolation accumulates during allopatric speciation. We investigated the population genetics of a contact zone between two highly divergent lineages of field voles (Microtus agrestis) in the Swiss Jura mountains. To shed light on the processes underlying introgression, we used maternally, paternally, and bi-parentally inherited markers. Though the two lineages maintained a strong genetic structure, we found some hybrids and evidence of gene flow. The extent of introgression varied with the mode of inheritance, being highest for mtDNA and absent for the Y chromosome. In addition, introgression was asymmetric, occurring only from the Northern to the Southern lineage. Both patterns seem parsimoniously explained by neutral processes linked to differences in effective sizes and sex-biased dispersal rates. The lineage with lower effective population size was also the more introgressed, and the mode-of-inheritance effect correlated with the male-biased dispersal rate of microtine rodents. We cannot exclude, however, that Haldane's effect contributed to the latter, as we found a marginally significant deficit in males (the heterogametic sex) among hybrids. We propose a possible demographic scenario to account for the patterns documented, and empirical extensions to further investigate this contact zone.

Phylogeographical structure, postglacial recolonization and barriers to gene flow in the distinctive Valais chromosome race of the common shrew (Sorex araneus).

Using one male-inherited and eight biparentally inherited microsatellite markers, we investigate the population genetic structure of the Valais chromosome race of the common shrew (Sorex araneus) in the Central Alps of Europe. Unexpectedly, the Y-chromosome microsatellite suggests nearly complete absence of male gene flow among populations from the St-Bernard and Simplon regions (Switzerland). Autosomal markers also show significant genetic structuring among these two geographical areas. Isolation by distance is significant and possible barriers to gene flow exist in the study area. Two different approaches are used to better understand the geographical patterns and the causes of this structuring. Using a principal component analysis for which testing procedure exists, and partial Mantel tests, we show that the St-Bernard pass does not represent a significant barrier to gene flow although it culminates at 2469 m, close to the highest altitudinal record for this species. Similar results are found for the Simplon pass, indicating that both passes represented potential postglacial recolonization routes into Switzerland from Italian refugia after the last Pleistocene glaciations. In contrast with the weak effect of these mountain passes, the Rhône valley lowlands significantly reduce gene flow in this species. Natural obstacles (the large Rhône river) and unsuitable habitats (dry slopes) are both present in the valley. Moreover, anthropogenic changes to landscape structures are likely to have strongly reduced available habitats for this shrew in the lowlands, thereby promoting genetic differentiation of populations found on opposite sides of the Rhône valley.

Association between variants of the leptin receptor gene (LEPR) and overweight: a systematic review and an analysis of the CoLaus study.

BACKGROUND: Three non-synonymous single nucleotide polymorphisms (Q223R, K109R and K656N) of the leptin receptor gene (LEPR) have been tested for association with obesity-related outcomes in multiple studies, showing inconclusive results. We performed a systematic review and meta-analysis on the association of the three LEPR variants with BMI. In addition, we analysed 15 SNPs within the LEPR gene in the CoLaus study, assessing the interaction of the variants with sex. METHODOLOGY/PRINCIPAL FINDINGS: We searched electronic databases, including population-based studies that investigated the association between LEPR variants Q223R, K109R and K656N and obesity- related phenotypes in healthy, unrelated subjects. We furthermore performed meta-analyses of the genotype and allele frequencies in case-control studies. Results were stratified by SNP and by potential effect modifiers. CoLaus data were analysed by logistic and linear regressions and tested for interaction with sex. The meta-analysis of published data did not show an overall association between any of the tested LEPR variants and overweight. However, the choice of a BMI cut-off value to distinguish cases from controls was crucial to explain heterogeneity in Q223R. Differences in allele frequencies across ethnic groups are compatible with natural selection of derived alleles in Q223R and K109R and of the ancient allele in K656N in Asians. In CoLaus, the rs10128072, rs3790438 and rs3790437 variants showed interaction with sex for their association with overweight, waist circumference and fat mass in linear regressions. CONCLUSIONS: Our systematic review and analysis of primary data from the CoLaus study did not show an overall association between LEPR SNPs and overweight. Most studies were underpowered to detect small effect sizes. A potential effect modification by sex, population stratification, as well as the role of natural selection should be addressed in future genetic association studies.

Intraspecific competition reveals conditional fitness effects of single gene polymorphism at the Arabidopsis root growth regulator BRX.

Intraspecific genetic variation for morphological traits is observed in many organisms. In Arabidopsis thaliana, alleles responsible for intraspecific morphological variation are increasingly being identified. However, the fitness consequences remain unclear in most cases. Here, the fitness effects of alleles of the BRX gene are investigated. A brx loss-of-function allele, which was found in a natural accession, results in a highly branched but poorly elongated root system. Comparison between the control accession Sav-0 and an introgression of brx into this background (brxS) indicated that, surprisingly, brx loss of function did not negatively affect fitness in pure stands. However, in mixed, well-watered stands brxS performance and reproductive output decreased significantly, as the proportion of Sav-0 neighbors increased. Additional comparisons between brxS and a brxS line that was complemented by a BRX transgene confirmed a direct effect of the loss-of-function allele on plant performance, as indicated by restored competitive ability of the transgenic genotype. Further, because plant height was very similar across genotypes and because the experimental setup largely excluded shading effects, the impaired competitiveness of the brx loss-of-function genotype likely reflects below-ground competition. In summary, these data reveal conditional fitness effects of a single gene polymorphism in response to intraspecific competition in Arabidopsis.

Pregnancy outcome following maternal exposure to statins: a multicentre prospective study.

OBJECTIVE: This contribution addresses the risk associated with exposure to statins during pregnancy. DESIGN: Multicentre observational prospective controlled study. SETTING: European Network of Teratology Information Services. POPULATION: Pregnant women who contacted one of 11 participating centres, seeking advice about exposure to statins during pregnancy, or to agents known to be nonteratogenic. METHODS: Pregnancies exposed during first trimester to statins were followed up prospectively, and their outcomes were compared with a matched control group. MAIN OUTCOME MEASURES: Rates of major birth defects, live births, miscarriages, elective terminations, preterm deliveries and gestational age and birthweight at delivery. RESULTS: We collected observations from 249 exposed pregnancies and 249 controls. The difference in the rate of major birth defects between the statin-exposed and the control groups was small and statistically nonsignificant (4.1% versus 2.7% odds ratio [OR] 1.5; 95% confidence interval [95% CI] 0.5-4.5, P = 0.43). In an adjusted Cox model, the difference between miscarriage rates was also small and not significant (hazard ratio 1.36, 95% CI 0.63-2.93, P = 0.43). Premature birth was more frequent in exposed pregnancies (16.1% versus 8.5%; OR 2.1, 95% CI 1.1-3.8, P = 0.019). Nonetheless, median gestational age at birth (39 weeks, interquartile range [IQR] 37-40 versus 39 weeks, IQR 38-40, P = 0.27) and birth weight (3280 g, IQR 2835-3590 versus 3250 g, IQR 2880-3630, P = 0.95) did not differ between exposed and non-exposed pregnancies. CONCLUSIONS: This study did not detect a teratogenic effect of statins. Its statistical power remains insufficient to challenge current recommendations of treatment discontinuation during pregnancy.