Cross-alteration of murine skin and tick microbiome concomitant with pathogen transmission after Ixodes ricinus bite.

BACKGROUND: Ticks are major vectors of diseases affecting humans such as Lyme disease or domestic animals such as anaplasmosis. Cross-alteration of the vertebrate host skin microbiome and the tick microbiome may be essential during the process of tick feeding and for the mechanism of pathogen transmission. However, it has been poorly investigated. METHODS: We used mice bitten by field-collected ticks (nymphs and adult ticks) in different experimental conditions to investigate, by 16S rRNA gene metabarcoding, the impact of blood feeding on both the mouse skin microbiome and the tick microbiome. We also investigated by PCR and 16S rRNA gene metabarcoding, the diversity of microorganisms transmitted to the host during the process of tick bite at the skin interface and the dissemination of the pathogen in host tissues (blood, heart, and spleen). RESULTS: Most of the commensal bacteria present in the skin of control mice were replaced during the blood-feeding process by bacteria originating from the ticks. The microbiome of the ticks was also impacted by the blood feeding. Several pathogens including tick-borne pathogens (Borrelia/Borreliella, Anaplasma, Neoehrlichia, Rickettsia) and opportunistic bacteria (Williamsia) were transmitted to the skin microbiome and some of them disseminated to the blood or spleen of the mice. In the different experiments of this study, skin microbiome alteration and Borrelia/Borreliella transmission were different depending on the tick stages (nymphs or adult female ticks). CONCLUSIONS: Host skin microbiome at the bite site was deeply impacted by the tick bite, to an extent which suggests a role in the tick feeding, in the pathogen transmission, and a potentially important impact on the skin physiopathology. The diversified taxonomic profiles of the tick microbiome were also modified by the blood feeding. Video Abstract.

Low-Diversity Microbiota in Apical Periodontitis and High Blood Pressure Are Signatures of the Severity of Apical Lesions in Humans.

(1) Background: In developed countries, the prevalence of apical periodontitis (AP) varies from 20% to 50% for reasons that could be associated with the apical periodontitis microbiota ecology. (2) Methods: We performed a clinical study in the Odontology department of Toulouse hospital in France, to sequence the 16S rRNA gene of AP microbiota and collect clinical parameters from 94 patients. Forty-four patients were characterized with a PAI (periapical index of AP severity) score lower or equal to 3, while the others had superior scores (n = 50). (3) Results: The low diversity of granuloma microbiota is associated with the highest severity (PAI = 5) of periapical lesions (Odds Ratio 4.592, IC 95% [1.6329; 14.0728]; p = 0.001; notably, a lower relative abundance of Burkholderiaceae and a higher relative abundance of Pseudomonas and Prevotella). We also identified that high blood pressure (HBP) is associated with the increase in PAI scores. (4) Conclusions: Our data show that a low diversity of bacterial ecology of the AP is associated with severe PAI scores, suggesting a causal mechanism. Furthermore, a second risk factor was blood pressure associated with the severity of apical periodontitis.

Prolonged SARS-CoV-2 RNA virus shedding and lymphopenia are hallmarks of COVID-19 in cancer patients with poor prognosis.

Patients with cancer are at higher risk of severe coronavirus infectious disease 2019 (COVID-19), but the mechanisms underlying virus-host interactions during cancer therapies remain elusive. When comparing nasopharyngeal swabs from cancer and noncancer patients for RT-qPCR cycle thresholds measuring acute respiratory syndrome coronavirus-2 (SARS-CoV-2) in 1063 patients (58% with cancer), we found that malignant disease favors the magnitude and duration of viral RNA shedding concomitant with prolonged serum elevations of type 1 IFN that anticorrelated with anti-RBD IgG antibodies. Cancer patients with a prolonged SARS-CoV-2 RNA detection exhibited the typical immunopathology of severe COVID-19 at the early phase of infection including circulation of immature neutrophils, depletion of nonconventional monocytes, and a general lymphopenia that, however, was accompanied by a rise in plasmablasts, activated follicular T-helper cells, and non-naive Granzyme B+FasL+, EomeshighTCF-1high, PD-1+CD8+ Tc1 cells. Virus-induced lymphopenia worsened cancer-associated lymphocyte loss, and low lymphocyte counts correlated with chronic SARS-CoV-2 RNA shedding, COVID-19 severity, and a higher risk of cancer-related death in the first and second surge of the pandemic. Lymphocyte loss correlated with significant changes in metabolites from the polyamine and biliary salt pathways as well as increased blood DNA from Enterobacteriaceae and Micrococcaceae gut family members in long-term viral carriers. We surmise that cancer therapies may exacerbate the paradoxical association between lymphopenia and COVID-19-related immunopathology, and that the prevention of COVID-19-induced lymphocyte loss may reduce cancer-associated death.

Intestinal Dysbiosis Amplifies Acetaminophen-Induced Acute Liver Injury.

BACKGROUND & AIMS: Acute liver failure (ALF) represents an unmet medical need in Western countries. Although the link between intestinal dysbiosis and chronic liver disease is well-established, there is little evidence for a functional role of gut-liver interaction during ALF. Here we hypothesized that intestinal dysbiosis may affect ALF. METHODS: To test this hypothesis, we assessed the association of proton pump inhibitor (PPI) or long-term antibiotics (ABx) intake, which have both been linked to intestinal dysbiosis, and occurrence of ALF in the 500,000 participants of the UK BioBank population-based cohort. For functional studies, male Nlrp6-/- mice were used as a dysbiotic mouse model and injected with a sublethal dose of acetaminophen (APAP) or lipopolysaccharide (LPS) to induce ALF. RESULTS: Multivariate Cox regression analyses revealed a significantly increased risk (odds ratio, 2.3-3) for developing ALF in UK BioBank participants with PPI or ABx. Similarly, dysbiotic Nlrp6-/- mice displayed exacerbated APAP- and LPS-induced liver injury, which was linked to significantly reduced gut and liver tissue microbiota diversity and correlated with increased intestinal permeability at baseline. Fecal microbiota transfer (FMT) from Nlrp6-/- mice into wild-type (WT) mice augmented liver injury on APAP treatment in recipient WT mice, resembling the inflammatory phenotype of Nlrp6-/- mice. Specifically, FMT skewed monocyte polarization in WT mice toward a Ly6Chi inflammatory phenotype, suggesting a critical function of these cells as sensors of gut-derived signals orchestrating the inflammatory response. CONCLUSIONS: Our data show an important yet unknown function of intestinal microbiota during ALF. Intestinal dysbiosis was transferrable to healthy WT mice via FMT and aggravated liver injury. Our study highlights intestinal microbiota as a targetable risk factor for ALF.

Sources of variation of DNA methylation in rainbow trout: combined effects of temperature and genetic background.

Phenotypic plasticity is a key component of the ability of organisms to respond to changing environmental conditions. In this study, we aimed to study the establishment of DNA methylation marks in response to an environmental stress in rainbow trout and to assess whether these marks depend on the genetic background. The environmental stress chosen here was temperature, a known induction factor of epigenetic marks in fish. To disentangle the role of epigenetic mechanisms such as DNA methylation in generating phenotypic variations, nine rainbow trout isogenic lines with no genetic variability within a line were used. For each line, half of the eggs were incubated at standard temperature (11°C) and the other half at high temperature (16°C), from eyed-stage to hatching. In order to gain a first insight into the establishment of DNA methylation marks in response to an early temperature regime (control 11°C vs. heated 16°C), we have studied the expression of 8 dnmt3 (DNA methyltransferase) genes, potentially involved in de novo methylation, and analysed global DNA methylation in the different rainbow trout isogenic lines using LUMA (LUminometric Methylation Assay). Finally, finer investigation of genome-wide methylation patterns was performed using EpiRADseq, a reduced-representation library approach based on the ddRADseq (Double Digest Restriction Associated DNA) protocol, for six rainbow trout isogenic lines. We have demonstrated that thermal history during embryonic development alters patterns of DNA methylation, but to a greater or lesser extent depending on the genetic background.

Insertions and Duplications in the Polyproline Region of the Hepatitis E Virus.

Recombinant strains of hepatitis E virus (HEV) with insertions of human genomic fragments or HEV sequence duplications in the sequence encoding the polyproline region (PPR) were previously described in chronically infected patients. Such genomic rearrangements confer a replicative advantage in vitro but little is known about their frequency, location, or origin. As the sequences of only a few virus genomes are available, we analyzed the complete genomes of 114 HEV genotype 3 strains from immunocompromised (n = 85) and immunocompetent (n = 29) patients using the single molecular real-time sequencing method to determine the frequency, location, and origin of inserted genomic fragments, plus the proportions of variants with genomic rearrangements in each virus quasispecies. We also examined the amino acid compositions and post-translational modifications conferred by these rearrangements by comparing them to sequences without human gene insertions or HEV gene duplications. We found genomic rearrangements in 7/114 (6.1%) complete genome sequences (4 HEV-3f, 1 HEV-3e, 1 HEV-3 h, and 1 HEV-3chi-new), all from immunocompromised patients, and 3/7 were found at the acute phase of infection. Six of the seven patients harbored virus-host recombinant variants, including one patient with two different recombinant variants. We also detected three recombinant variants with genome duplications of the PPR or PPR + X domains in a single patient. All the genomic rearrangements (seven human fragment insertions of varying origins and three HEV genome duplications) occurred in the PPR. The sequences with genomic rearrangements had specific characteristics: increased net load (p < 0.001) and more ubiquitination (p < 0.001), phosphorylation (p < 0.001), and acetylation (p < 0.001) sites. The human fragment insertions and HEV genome duplications had slightly different characteristics. We believe this is the first description of HEV strains with genomic rearrangements in patients at the acute phase of infection; perhaps these strains are directly transmitted. Clearly, genomic rearrangements produce a greater net load with duplications and insertions having different features. Further studies are needed to clarify the mechanisms by which such modifications influence HEV replication.

One complete and three draft genome sequences of four Brochothrix thermosphacta strains, CD 337, TAP 175, BSAS1 3 and EBP 3070.

Brochothrix thermosphacta is one of the dominant bacterial species associated with spoilage of chilled meat and seafood products through the production of various metabolites responsible for off-odors. However, metabolic pathways leading to meat and seafood spoilage are not all well known. The production of spoiling molecules seems to depend both on strains and on food matrix. Several B. thermosphacta genome sequences have been reported, all issued from meat isolates. Here, we report four genome sequences, one complete and three as drafts. The four B. thermosphacta strains CD 337, TAP 175, BSAS1 3, and EBP 3070 were isolated from different ecological niches (seafood or meat products either spoiled or not and bovine slaughterhouse). These strains known as phenotypically and genetically different were selected to represent intraspecies diversity. CD 337 genome is 2,594,337 bp long, complete and circular, containing 2593 protein coding sequences and 28 RNA genes. TAP 175, BSAS1 3, and EBP 3070 genomes are arranged in 57, 83, and 71 contigs, containing 2515, 2668, and 2611 protein-coding sequences, respectively. These genomes were compared with two other B. thermosphacta complete genome sequences. The main genome content differences between strains are phages, plasmids, restriction/modification systems, and cell surface functions, suggesting a similar metabolic potential but a different niche adaptation capacity.

The complete methylome of an entomopathogenic bacterium reveals the existence of loci with unmethylated Adenines.

DNA methylation can serve to control diverse phenomena in eukaryotes and prokaryotes, including gene regulation leading to cell differentiation. In bacteria, DNA methylomes (i.e., methylation state of each base of the whole genome) have been described for several species, but methylome profile variation during the lifecycle has rarely been studied, and only in a few model organisms. Moreover, major phenotypic changes have been reported in several bacterial strains with a deregulated methyltransferase, but the corresponding methylome has rarely been described. Here we report the first methylome description of an entomopathogenic bacterium, Photorhabdus luminescens. Eight motifs displaying a high rate of methylation (>94%) were identified. The methylome was strikingly stable over course of growth, but also in a subpopulation responsible for a critical step in the bacterium's lifecycle: successful survival and proliferation in insects. The rare unmethylated GATC motifs were preferentially located in putative promoter regions, and most of them were methylated after Dam methyltransferase overexpression, suggesting that DNA methylation is involved in gene regulation. Our findings bring key insight into bacterial methylomes and encourage further research to decipher the role of loci protected from DNA methylation in gene regulation.

Diversity of hepatitis E virus genotype 3.

Hepatitis E virus genotype 3 (HEV-3) can lead to chronic infection in immunocompromised patients, and ribavirin is the treatment of choice. Recently, mutations in the polymerase gene have been associated with ribavirin failure but their frequency before treatment according to HEV-3 subtypes has not been studied on a large data set. We used single-molecule real-time sequencing technology to sequence 115 new complete genomes of HEV-3 infecting French patients. We analyzed phylogenetic relationships, the length of the polyproline region, and mutations in the HEV polymerase gene. Eighty-five (74%) were in the clade HEV-3efg, 28 (24%) in HEV-3chi clade, and 2 (2%) in HEV-3ra clade. Using automated partitioning of maximum likelihood phylogenetic trees, complete genomes were classified into subtypes. Polyproline region length differs within HEV-3 clades (from 189 to 315 nt). Investigating mutations in the polymerase gene, distinct polymorphisms between HEV-3 subtypes were found (G1634R in 95% of HEV-3e, G1634K in 56% of HEV-3ra, and V1479I in all HEV-3efg, clade HEV-3ra, and HEV-3k strains). Subtype-specific polymorphisms in the HEV-3 polymerase have been identified. Our study provides new complete genome sequences of HEV-3 that could be useful for comparing strains circulating in humans and the animal reservoir.

Genomic repeats, misassembly and reannotation: a case study with long-read resequencing of Porphyromonas gingivalis reference strains.

BACKGROUND: Without knowledge of their genomic sequences, it is impossible to make functional models of the bacteria that make up human and animal microbiota. Unfortunately, the vast majority of publicly available genomes are only working drafts, an incompleteness that causes numerous problems and constitutes a major obstacle to genotypic and phenotypic interpretation. In this work, we began with an example from the class Bacteroidia in the phylum Bacteroidetes, which is preponderant among human orodigestive microbiota. We successfully identify the genetic loci responsible for assembly breaks and misassemblies and demonstrate the importance and usefulness of long-read sequencing and curated reannotation. RESULTS: We showed that the fragmentation in Bacteroidia draft genomes assembled from massively parallel sequencing linearly correlates with genomic repeats of the same or greater size than the reads. We also demonstrated that some of these repeats, especially the long ones, correspond to misassembled loci in three reference Porphyromonas gingivalis genomes marked as circularized (thus complete or finished). We prove that even at modest coverage (30X), long-read resequencing together with PCR contiguity verification (rrn operons and an integrative and conjugative element or ICE) can be used to identify and correct the wrongly combined or assembled regions. Finally, although time-consuming and labor-intensive, consistent manual biocuration of three P. gingivalis strains allowed us to compare and correct the existing genomic annotations, resulting in a more accurate interpretation of the genomic differences among these strains. CONCLUSIONS: In this study, we demonstrate the usefulness and importance of long-read sequencing in verifying published genomes (even when complete) and generating assemblies for new bacterial strains/species with high genomic plasticity. We also show that when combined with biological validation processes and diligent biocurated annotation, this strategy helps reduce the propagation of errors in shared databases, thus limiting false conclusions based on incomplete or misleading information.

The Complete Genome Sequence of the Fish Pathogen Tenacibaculum maritimum Provides Insights into Virulence Mechanisms.

Tenacibaculum maritimum is a devastating bacterial pathogen of wild and farmed marine fish with a broad host range and a worldwide distribution. We report here the complete genome sequence of the T. maritimum type strain NCIMB 2154T. The genome consists of a 3,435,971-base pair circular chromosome with 2,866 predicted protein-coding genes. Genes encoding the biosynthesis of exopolysaccharides, the type IX secretion system, iron uptake systems, adhesins, hemolysins, proteases, and glycoside hydrolases were identified. They are likely involved in the virulence process including immune escape, invasion, colonization, destruction of host tissues, and nutrient scavenging. Among the predicted virulence factors, type IX secretion-mediated and cell-surface exposed proteins were identified including an atypical sialidase, a sphingomyelinase and a chondroitin AC lyase which activities were demonstrated in vitro.

Accurate Identification and Quantification of DNA Species by Next-Generation Sequencing in Adeno-Associated Viral Vectors Produced in Insect Cells.

Recombinant adeno-associated viral (rAAV) vectors have proven excellent tools for the treatment of many genetic diseases and other complex diseases. However, the illegitimate encapsidation of DNA contaminants within viral particles constitutes a major safety concern for rAAV-based therapies. Moreover, the development of rAAV vectors for early-phase clinical trials has revealed the limited accuracy of the analytical tools used to characterize these new and complex drugs. Although most published data concerning residual DNA in rAAV preparations have been generated by quantitative PCR, we have developed a novel single-strand virus sequencing (SSV-Seq) method for quantification of DNA contaminants in AAV vectors produced in mammalian cells by next-generation sequencing (NGS). Here, we describe the adaptation of SSV-Seq for the accurate identification and quantification of DNA species in rAAV stocks produced in insect cells. We found that baculoviral DNA was the most abundant contaminant, representing less than 2.1% of NGS reads regardless of serotype (2, 8, or rh10). Sf9 producer cell DNA was detected at low frequency (≤0.03%) in rAAV lots. Advanced computational analyses revealed that (1) baculoviral sequences close to the inverted terminal repeats preferentially underwent illegitimate encapsidation, and (2) single-nucleotide variants were absent from the rAAV genome. The high-throughput sequencing protocol described here enables effective DNA quality control of rAAV vectors produced in insect cells, and is adapted to conform with regulatory agency safety requirements.

Complete mitochondrial genome sequence of the freshwater diatom Asterionella formosa.

We report the complete mitochondrial genome sequence of the freshwater diatom Asterionella formosa. The large 61.9 kb circular sequence encodes 34 proteins and 25 tRNAs that are universally conserved in other sequenced diatoms. We fully resolved a unique 24 kb region containing highly conserved repeated sequence units, possibly collocating with an origin of replication.

A Panel of Embryonic Stem Cell Lines Reveals the Variety and Dynamic of Pluripotent States in Rabbits.

Conventional rabbit embryonic stem cell (ESC) lines are derived from the inner cell mass (ICM) of pre-implantation embryos using methods and culture conditions that are established for primate ESCs. In this study, we explored the capacity of the rabbit ICM to give rise to ESC lines using conditions similar to those utilized to generate naive ESCs in mice. On single-cell dissociation and culture in fibroblast growth factor 2 (FGF2)-free, serum-supplemented medium, rabbit ICMs gave rise to ESC lines lacking the DNA-damage checkpoint in the G1 phase like mouse ESCs, and with a pluripotency gene expression profile closer to the rabbit ICM/epiblast profiles. These cell lines can be converted to FGF2-dependent ESCs after culture in conventional conditions. They can also colonize the rabbit pre-implantation embryo. These results indicate that rabbit epiblast cells can be coaxed toward different types of pluripotent stem cells and reveal the dynamics of pluripotent states in rabbit ESCs.

Genome-Wide Study of Structural Variants in Bovine Holstein, Montbéliarde and Normande Dairy Breeds.

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.

Interaction of macrocyclic lactones with P-glycoprotein: structure-affinity relationship.

P-glycoprotein (P-gp) is involved in the ATP-dependant cellular efflux of a large number of drugs including ivermectin, a macrocyclic lactone (ML) endectocide, widely used in livestock and human antiparasitic therapy. The interactions of P-gp with ivermectin and other MLs were studied. In a first approach, the ability of ivermectin (IVM), eprinomectin (EPR), abamectin (ABA), doramectin (DOR), selamectin (SEL), or moxidectin (MOX) to inhibit the rhodamine123 efflux was measured in recombinant cells overexpressing P-gp. Then, the influence of these compounds on the P-gp ATPase activity was tested on membrane vesicles prepared from fibroblasts overexpressing P-gp. All the MLs tested increased the intracellular rhodamine123. However, the potency of MOX to inhibit P-gp function was 10 times lower than the other MLs. They all inhibited the basal and decreased the verapamil-stimulated P-gp ATPase activity. But SEL and MOX were less potent than the other MLs when competing with verapamil. According to the structural specificity of SEL and MOX, we conclude that the integrity of the sugar moiety is determinant to achieve the optimal interaction of macrocyclic lactones with P-gp. The structure-affinity relationship for interaction with P-gp is important information for improving ML bioavailability and reversal of multidrug resistance (MDR).

Transcriptional modulations by RXR agonists are only partially subordinated to PPARalpha signaling and attest additional, organ-specific, molecular cross-talks.

Nuclear hormone receptors (NR) are important transcriptional regulators of numerous genes involved in diverse pathophysiological and therapeutic functions. Following ligand activation, class II NR share the ability to heterodimerize with the retinoid X receptor (RXR). It is established that RXR activators, rexinoids, transactivate several peroxisome proliferator-activated receptor alpha (PPARalpha) target genes in a PPARalpha-dependent manner. We hypothesized that, once activated, RXR might signal through quiescent NR other than PPARalpha, in an organ-specific manner. To study this putative phenomenon in vivo, we developed an array of 120 genes relevant to the class II NR field. The genes were selected using both published data and high-density screenings performed on RXR or PPARalpha agonist-treated mice. Wild-type C57BL/6J and PPARalpha-deficient mice were treated with fenofibrate (PPARalpha activator) or LGD1069 (RXR activator). Using our customized array, we studied the hepatic, cardiac, and renal expression of this panel of 120 genes and compared them in both murine genotypes. The results obtained from this study confirmed the ability of an RXR agonist to modulate PPARalpha-restricted target genes in the liver and the kidney. Furthermore, we show that various organ-specific regulations occurring in both genotypes (PPARalpha +/+ or -/-) are highly indicative of the ability of RXR to recruit other class II NR pathways. Further development of this molecular tool may lead to a better understanding of the permissiveness of class II nuclear receptor dimers in vivo.

MDR1-deficient genotype in Collie dogs hypersensitive to the P-glycoprotein substrate ivermectin.

Multidrug resistance (MDR) phenotypes in cancer cells are associated with overexpression of the drug carrier P-glycoprotein. The antiparasitic drug ivermectin, one of its substrates, abnormally accumulates in the brain of transgenic mice lacking the P-glycoprotein, resulting in neurotoxicity. Similarly, an enhanced sensitivity to ivermectin has been reported in certain dogs of the Collie breed. To explore the basis of this phenotype, we analyzed the canine P-glycoprotein-encoding MDR1 gene, and we report the first characterization of the cDNA for wild-type (Beagle) P-glycoprotein. The corresponding transcripts from ivermectin-sensitive Collies revealed a homozygous 4-bp exonic deletion. We established, by genetic testings, that the MDR1 frame shift is predictable. Accordingly, no P-glycoprotein was detected in the homozygote-deficient dogs. In conclusion, we characterized a unique case of naturally occurring gene invalidation. This provides a putative novel model that remains to be exploited in the field of human therapeutics and that might significantly affect tissue distribution and drug bioavailability studies.

Comparative effect of fenofibrate on hepatic desaturases in wild-type and peroxisome proliferator-activated receptor alpha-deficient mice.

In this study is presented the effect of fenofibrate, a prototypical peroxisome proliferator of the fibrate class, on wild-type and peroxisome proliferator-activated receptor alpha (PPARalpha)-/- mouse liver FA profile, desaturase mRNA levels, and activities. We established that, following peroxisome proliferator exposure, the hepatic FA profile was greatly modified. These modifications in hepatic FA content required the expression of PPARalpha, as they are suppressed in transgenic mice deficient in this nuclear receptor. Following peroxisome proliferator exposure, delta6- and delta5-desaturase mRNA levels and activities were increased in wild-type but not in PPARalpha-deficient mouse liver. These results suggest the involvement of PPARalpha in the control of hepatic delta6- and delta5-desaturases in mice. Their roles in minimizing long-chain PUFA depletion in the liver during peroxisome proliferator exposure are discussed.