Coprophagia in early life tunes expression of immune genes after weaning in rabbit ileum

Coprophagia by suckling rabbits, i.e. ingestion of feces from their mother, reduces mortality after weaning. We hypothesized that this beneficial effect of coprophagia is immune-mediated at the intestinal level. Therefore, this study investigated immune development after weaning by analyzing the ileal transcriptome at day 35 and 49 in rabbits with differential access to coprophagia in early life. Rabbit pups had access between day 1 and 15 to (i) no feces (NF) or (ii) feces from unrelated does (Foreign Feces, FF) or (iii) feces from unrelated does treated with antibiotics (FFab). 350 genes were differentially expressed between day 35 and day 49 in suckling rabbits with access to coprophagia. These genes coded for antimicrobial peptides, a mucin, cytokines and chemokines, pattern recognition receptors, proteins involved in immunoglobulin A secretion and in interferon signaling pathway. Strikingly, prevention of coprophagia or access to feces from antibiotic-treated does in early life blunted immune development between day 35 et 49 in the ileum of rabbits. Thus, coprophagia might be crucial for the maturation of intestinal immunity in rabbits and could explain why this behavior improves survival.

Sex and fetal genome influence gene expression in pig endometrium at the end of gestation

BACKGROUND

A fine balance of feto-maternal resource allocation is required to support pregnancy, which depends on interactions between maternal and fetal genetic potential, maternal nutrition and environment, endometrial and placental functions. In particular, some imprinted genes have a role in regulating maternal-fetal nutrient exchange, but few have been documented in the endometrium. The aim of this study is to describe the expression of 42 genes, with parental expression, in the endometrium comparing two extreme breeds: Large White (LW); Meishan (MS) with contrasting neonatal mortality and maturity at two days of gestation (D90-D110). We investigated their potential contribution to fetal maturation exploring genes-fetal phenotypes relationships. Last, we hypothesized that the fetal genome and sex influence their endometrial expression. For this purpose, pure and reciprocally crossbred fetuses were produced using LW and MS breeds. Thus, in the same uterus, endometrial samples were associated with its purebred or crossbred fetuses.

RESULTS

Among the 22 differentially expressed genes (DEGs), 14 DEGs were differentially regulated between the two days of gestation. More gestational changes were described in LW (11 DEGs) than in MS (2 DEGs). Nine DEGs were differentially regulated between the two extreme breeds, highlighting differences in the regulation of endometrial angiogenesis, nutrient transport and energy metabolism. We identified DEGs that showed high correlations with indicators of fetal maturation, such as ponderal index at D90 and fetal blood fructose level and placental weight at D110. We pointed out for the first time the influence of fetal sex and genome on endometrial expression at D90, highlighting AMPD3, CITED1 and H19 genes. We demonstrated that fetal sex affects the expression of five imprinted genes in LW endometrium. Fetal genome influenced the expression of four genes in LW endometrium but not in MS endometrium. Interestingly, both fetal sex and fetal genome interact to influence endometrial gene expression.

CONCLUSIONS

These data provide evidence for some sexual dimorphism in the pregnant endometrium and for the contribution of the fetal genome to feto-maternal interactions at the end of gestation. They suggest that the paternal genome may contribute significantly to piglet survival, especially in crossbreeding production systems.

Multi-tissue metabolic and transcriptomic responses to a short-term heat stress in swine

BACKGROUND

Heat stress (HS) is an increasing threat for pig production with a wide range of impacts. When submitted to high temperatures, pigs will use a variety of strategies to alleviate the effect of HS. While systemic adaptations are well known, tissue-specific changes remain poorly understood. In this study, thirty-two pigs were submitted to a 5-day HS at 32 °C.

RESULTS

Transcriptomic and metabolomic analyses were performed on several tissues. The results revealed differentially expressed genes and metabolites in different tissues. Specifically, 481, 1774, 71, 1572, 17, 164, and 169 genes were differentially expressed in muscle, adipose tissue, liver, blood, thyroid, pituitary, and adrenal glands, respectively. Regulatory glands (pituitary, thyroid, and adrenal) had a lower number of regulated genes, perhaps indicating an earlier sensitivity to HS. In addition, 7, 8, 2, and 8 metabolites were differentially produced in muscle, liver, plasma, and urine, respectively. The study also focused on the oxidative stress pathway in muscle and liver by performing a correlation analysis between genes and metabolites.

CONCLUSIONS

This study has identified various adaptation mechanisms in swine that enable them to cope with heat stress (HS). These mechanisms include a global decrease in energetic metabolism, as well as changes in metabolic precursors that are linked with protein and lipid catabolism and anabolism. Notably, the adaptation mechanisms differ significantly between regulatory (pituitary, thyroid and adrenal glands) and effector tissues (muscle, adipose tissue, liver and blood). Our findings provide new insights into the comprehension of HS adaptation mechanisms in swine.

Asterics: a simple tool for the ExploRation and Integration of omiCS data

BACKGROUND

The rapid development of omics acquisition techniques has induced the production of a large volume of heterogeneous and multi-level omics datasets, which require specific and sometimes complex analyses to obtain relevant biological information. Here, we present ASTERICS (version 2.5), a publicly available web interface for the analyses of omics datasets.

RESULTS

ASTERICS is designed to make both standard and complex exploratory and integration analysis workflows easily available to biologists and to provide high quality interactive plots. Special care has been taken to provide a comprehensive documentation of the implemented analyses and to guide users toward sound analysis choices regarding some specific omics data. Data and analyses are organized in a comprehensive graphical workflow within ASTERICS workspace to facilitate the understanding of successive data editions and analyses leading to a given result.

CONCLUSION

ASTERICS provides an easy to use platform for omics data exploration and integration. The modular organization of its open source code makes it easy to incorporate new workflows and analyses by external contributors. ASTERICS is available at https://asterics.miat.inrae.fr and can also be deployed using provided docker images.

Plasma 1H-NMR metabolic and amino acid profiles of newborn piglets from two lines divergently selected for residual feed intake

Together with environmental factors, physiological maturity at birth is a major determinant for neonatal survival and postnatal development in mammalian species. Maturity at birth is the outcome of complex mechanisms of intra-uterine development and maturation during the end of gestation. In pig production, piglet preweaning mortality averages 20% of the litter and thus, maturity is a major welfare and economic concern. Here, we used both targeted and untargeted metabolomic approaches to provide a deeper understanding of the maturity in a model of lines of pigs divergently selected on residual feed intake (RFI), previously shown to have contrasted signs of maturity at birth. Analyses were conducted on plasma metabolome of piglets at birth and integrated with other phenotypic characteristics associated to maturity. We confirmed proline and myo-inositol, previously described for their association with delayed growth, as potential markers of maturity. Urea cycle and energy metabolism were found more regulated in piglets from high and low RFI lines, respectively, suggesting a better thermoregulation ability for the low RFI (with higher feed efficiency) piglets.

Whole-genome sequencing of cryopreserved resources from French Large White pigs at two distinct sampling times reveals strong signatures of convergent and divergent selection between the dam and sire lines

BACKGROUND

Numerous genomic scans for positive selection have been performed in livestock species within the last decade, but often a detailed characterization of the detected regions (gene or trait under selection, timing of selection events) is lacking. Cryopreserved resources stored in reproductive or DNA gene banks offer a great opportunity to improve this characterization by providing direct access to recent allele frequency dynamics, thereby differentiating between signatures from recent breeding objectives and those related to more ancient selection constraints. Improved characterization can also be achieved by using next-generation sequencing data, which helps narrowing the size of the detected regions while reducing the number of associated candidate genes.

METHODS

We estimated genetic diversity and detected signatures of recent selection in French Large White pigs by sequencing the genomes of 36 animals from three distinct cryopreserved samples: two recent samples from dam (LWD) and sire (LWS) lines, which had diverged from 1995 and were selected under partly different objectives, and an older sample from 1977 prior to the divergence.

RESULTS

French LWD and LWS lines have lost approximately 5% of the SNPs that segregated in the 1977 ancestral population. Thirty-eight genomic regions under recent selection were detected in these lines and the corresponding selection events were further classified as convergent between lines (18 regions), divergent between lines (10 regions), specific to the dam line (6 regions) or specific to the sire line (4 regions). Several biological functions were found to be significantly enriched among the genes included in these regions: body size, body weight and growth regardless of the category, early life survival and calcium metabolism more specifically in the signatures in the dam line and lipid and glycogen metabolism more specifically in the signatures in the sire line. Recent selection on IGF2 was confirmed and several other regions were linked to a single candidate gene (ARHGAP10, BMPR1B, GNA14, KATNA1, LPIN1, PKP1, PTH, SEMA3E or ZC3HAV1, among others).

CONCLUSIONS

These results illustrate that sequencing the genome of animals at several recent time points generates considerable insight into the traits, genes and variants under recent selection in a population. This approach could be applied to other livestock populations, e.g. by exploiting the rich biological resources stored in cryobanks.

Joint Automatic Metabolite Identification and Quantification of a Set of 1H NMR Spectra

Metabolomics is a promising approach to characterize phenotypes or to identify biomarkers. It is also easily accessible through NMR, which can provide a comprehensive understanding of the metabolome of any living organisms. However, the analysis of ¹H NMR spectrum remains difficult, mainly due to the different problems encountered to perform automatic identification and quantification of metabolites in a reproducible way. In addition, methods that perform automatic identification and quantification of metabolites are often designed to process one given complex mixture spectrum at a time. Hence, when a set of complex mixture spectra coming from the same experiment has to be processed, the approach is simply repeated independently for every spectrum, despite their resemblance. Here, we present new methods that are the first to either align spectra or to identify and quantify metabolites by integrating information coming from several complex spectra of the same experiment. The performances of these new methods are then evaluated on both simulated and real datasets. The results show an improvement in the metabolite identification and in the accuracy of metabolite quantifications, especially when the concentration is low. This joint procedure is available in version 2.0 of ASICS package.

Thermoregulation at birth differs between piglets from two genetic lines divergent for residual feed intake

Thermoregulation is essential to piglets' neonatal survival. This study used infrared thermography (IRT) to assess thermoregulation abilities of piglets from two lines divergent for residual feed intake (RFI). At birth, morphology (weight, length, width and circumference), vigour (respiration, mobility and vocalisation), and rectal temperature were recorded from piglets of the 11th generation of the low RFI (LRFI, more efficient; n = 34) and the high RFI (HRFI, less efficient; n = 28) lines. Infrared thermography images were taken at 8, 15, 30 and 60 min post partum. Temperatures of the ear base and tip, and of the back (i.e. shoulders to rumps) were extracted (Thermacam Researcher Pro 2.0) and analysed with linear mixed models (SAS 9.4). Piglets had different average hourly weight gain (HRFI = 7.1 ± 1.3 g/h, LRFI = 3.6 ± 1.3 g/h; P < 0,001) but did not differ in morphology or vigour. All temperatures increased overtime. At birth, piglets' rectal temperature was correlated with the initial temperature of the ear base and the maximum back temperature (0.37 and 0.33, respectively; P < 0.05). High residual feed intake piglets had lower ear tip temperatures than LRFI piglets at 15 (24.7 ± 0.37 °C vs. 26.3 ± 0.36 °C, respectively; F1, 63.5 = 9.11, P < 0.005) and 30 min post partum (26.2 ± 0.47 °C vs. 27.6 ± 0.44 °C, respectively; F1, 66.9 = 4.52, P < 0.05). Moreover, thermal pattern of the ear tip differed between the two genetic lines. In conclusion, IRT allowed non-invasive assessment of piglets' thermoregulation abilities and indicated an influence of genetic selection for RFI on neonatal thermoregulation abilities.

ASICS: an R package for a whole analysis workflow of 1D 1H NMR spectra

MOTIVATION

In metabolomics, the detection of new biomarkers from Nuclear Magnetic Resonance (NMR) spectra is a promising approach. However, this analysis remains difficult due to the lack of a whole workflow that handles spectra pre-processing, automatic identification and quantification of metabolites and statistical analyses, in a reproducible way.

RESULTS

We present ASICS, an R package that contains a complete workflow to analyse spectra from NMR experiments. It contains an automatic approach to identify and quantify metabolites in a complex mixture spectrum and uses the results of the quantification in untargeted and targeted statistical analyses. ASICS was shown to improve the precision of quantification in comparison to existing methods on two independent datasets. In addition, ASICS successfully recovered most metabolites that were found important to explain a two level condition describing the samples by a manual and expert analysis based on bucketing. It also found new relevant metabolites involved in metabolic pathways related to risk factors associated with the condition.

AVAILABILITY AND IMPLEMENTATION

ASICS is distributed as an R package, available on Bioconductor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Proteomic analysis of adipose tissue during the last weeks of gestation in pure and crossbred Large White or Meishan fetuses gestated by sows of either breed

Background

The degree of adipose tissue development at birth may influence neonatal survival and subsequent health outcomes. Despite their lower birth weights, piglets from Meishan sows (a fat breed with excellent maternal ability) have a higher survival rate than piglets from Large White sows (a lean breed). To identify the main pathways involved in subcutaneous adipose tissue maturation during the last month of gestation, we compared the proteome and the expression levels of some genes at d 90 and d 110 of gestation in purebred and crossbred Large White or Meishan fetuses gestated by sows of either breed.

Results

A total of 52 proteins in fetal subcutaneous adipose tissue were identified as differentially expressed over the course of gestation. Many proteins involved in energy metabolism were more abundant, whereas some proteins participating in cytoskeleton organization were reduced in abundance on d 110 compared with d 90. Irrespective of age, 24 proteins differed in abundance between fetal genotypes, and an interaction effect between fetal age and genotype was observed for 13 proteins. The abundance levels of proteins known to be responsive to nutrient levels such as aldolase and fatty acid binding proteins, as well as the expression levels of FASN, a key lipogenic enzyme, and MLXIPL, a pivotal transcriptional mediator of glucose-related stimulation of lipogenic genes, were elevated in the adipose tissue of pure and crossbred fetuses from Meishan sows. These data suggested that the adipose tissue of these fetuses had superior metabolic functionality, whatever their paternal genes. Conversely, proteins participating in redox homeostasis and apoptotic cell clearance had a lower abundance in Meishan than in Large White fetuses. Time-course differences in adipose tissue protein abundance were revealed between fetal genotypes for a few secreted proteins participating in responses to organic substances, such as alpha-2-HS-glycoprotein, transferrin and albumin.

Conclusions

These results underline the importance of not only fetal age but also maternal intrauterine environment in the regulation of several proteins in subcutaneous adipose tissue. These proteins may be used to estimate the maturity grade of piglet neonates.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0244-2) contains supplementary material, which is available to authorized users.

Transgenerational epigenetic inheritance in birds

While it has been shown that epigenetics accounts for a portion of the variability of complex traits linked to interactions with the environment, the real contribution of epigenetics to phenotypic variation remains to be assessed. In recent years, a growing number of studies have revealed that epigenetic modifications can be transmitted across generations in several animal species. Numerous studies have demonstrated inter- or multi-generational effects of changing environment in birds, but very few studies have been published showing epigenetic transgenerational inheritance in these species. In this review, we mention work conducted in parent-to-offspring transmission analyses in bird species, with a focus on the impact of early stressors on behaviour. We then present recent advances in transgenerational epigenetics in birds, which involve germline linked non-Mendelian inheritance, underline the advantages and drawbacks of working on birds in this field and comment on future directions of transgenerational studies in bird species.

Integrated Analysis of Proteomic and Transcriptomic Data Highlights Late Fetal Muscle Maturation Process

In pigs, the perinatal period is the most critical time for survival. Piglet maturation, which occurs at the end of gestation, is an important determinant of early survival. Skeletal muscle plays a key role in adaptation to extra-uterine life, e.g. motor function and thermoregulation. Progeny from two breeds with extreme neonatal mortality rates were analyzed at 90 and 110 days of gestation (dg). The Large White breed is a highly selected breed for lean growth and exhibits a high rate of neonatal mortality, whereas the Meishan breed is fatter and more robust and has a low neonatal mortality. Our aim was to identify molecular signatures underlying late fetal longissimus muscle development. First, integrated analysis was used to explore relationships between co-expression network models built from a proteomic data set (bi-dimensional electrophoresis) and biological phenotypes. Second, correlations with a transcriptomic data set (microarrays) were investigated to combine different layers of expression with a focus on transcriptional regulation. Muscle glycogen content and myosin heavy chain polymorphism were good descriptors of muscle maturity and were used for further data integration analysis. Using 89 identified unique proteins, network inference, correlation with biological phenotypes and functional enrichment revealed that mitochondrial oxidative metabolism was a key determinant of neonatal muscle maturity. Some proteins, including ATP5A1 and CKMT2, were important nodes in the network related to muscle metabolism. Transcriptomic data suggest that overexpression of mitochondrial PCK2 was involved in the greater glycogen content of Meishan fetuses at 110 dg. GPD1, an enzyme involved in the mitochondrial oxidation of cytosolic NADH, was overexpressed in Meishan. Thirty-one proteins exhibited a positive correlation between mRNA and protein levels in both extreme fetal genotypes, suggesting transcriptional regulation. Gene ontology enrichment and Ingenuity analyses identified PPARGC1A and ESR1 as possible transcriptional factors positively involved in late fetal muscle maturation.

Time course study of the response to LPS targeting the pig immune gene networks

Background

Stress is a generic term used to describe non-specific responses of the body to all kinds of challenges. A very large variability in the response can be observed across individuals, depending on numerous conditioning factors like genetics, early influences and life history. As a result, there is a wide range of individual vulnerability and resilience to stress, also called robustness. The importance of robustness-related traits in breeding strategies is increasing progressively towards the production of animals with a high level of production under a wide range of climatic conditions and management systems, together with a lower environmental impact and a high level of animal welfare. The present study aims at describing blood transcriptomic, hormonal, and metabolic responses of pigs to a systemic challenge using lipopolysaccharide (LPS). The objective is to analyze the individual variation of the biological responses in relation to the activity of the HPA axis measured by the levels of plasma cortisol after LPS and ACTH in 120 juvenile Large White (LW) pigs. The kinetics of the response was measured with biological variables and whole blood gene expression at 4 time points. A multilevel statistical analysis was used to take into account the longitudinal aspect of the data.

Results

Cortisol level reaches its peak 4 h after LPS injection. The characteristic changes of white blood cell count to LPS were observed, with a decrease of total count, maximal at t=+4 h, and the mirror changes in the respective proportions of lymphocytes and granulocytes. The lymphocytes / granulocytes ratio was maximal at t=+1 h. An integrative statistical approach was used and provided a set of candidate genes for kinetic studies and ongoing complementary studies focused on the LPS-stimulated inflammatory response.

Conclusions

The present study demonstrates the specific biomarkers indicative of an inflammation in swine. Furthermore, these stress responses persist for prolonged periods of time and at significant expression levels, making them good candidate markers for evaluating the efficacy of anti-inflammatory drugs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4363-5) contains supplementary material, which is available to authorized users.

1HNMR-Based metabolomic profiling method to develop plasma biomarkers for sensitivity to chronic heat stress in growing pigs

The negative impact of heat stress (HS) on the production performances in pig faming is of particular concern. Novel diagnostic methods are needed to predict the robustness of pigs to HS. Our study aimed to assess the reliability of blood metabolome to predict the sensitivity to chronic HS of 10 F1 (Large White × Creole) sire families (SF) reared in temperate (TEMP) and in tropical (TROP) regions (n = 56±5 offsprings/region/SF). Live body weight (BW) and rectal temperature (RT) were recorded at 23 weeks of age. Average daily feed intake (AFDI) and average daily gain were calculated from weeks 11 to 23 of age, together with feed conversion ratio. Plasma blood metabolome profiles were obtained by Nuclear Magnetic Resonance spectroscopy (^1HNMR) from blood samples collected at week 23 in TEMP. The sensitivity to hot climatic conditions of each SF was estimated by computing a composite index of sensitivity (I_sens) derived from a linear combination of t statistics applied to familial BW, ADFI and RT in TEMP and TROP climates. A model of prediction of sensitivity was established with sparse Partial Least Square Discriminant Analysis (sPLS-DA) between the two most robust SF (n = 102) and the two most sensitive ones (n = 121) using individual metabolomic profiles measured in TEMP. The sPLS-DA selected 29 buckets that enabled 78% of prediction accuracy by cross-validation. On the basis of this training, we predicted the proportion of sensitive pigs within the 6 remaining families (n = 337). This proportion was defined as the predicted membership of families to the sensitive category. The positive correlation between this proportion and I_sens (r = 0.97, P < 0.01) suggests that plasma metabolome can be used to predict the sensitivity of pigs to hot climate.

Comparing the intestinal transcriptome of Meishan and Large White piglets during late fetal development reveals genes involved in glucose and lipid metabolism and immunity as valuable clues of intestinal maturity

Background

Maturity of intestinal functions is critical for neonatal health and survival, but comprehensive description of mechanisms underlying intestinal maturation that occur during late gestation still remain poorly characterized. The aim of this study was to investigate biological processes specifically involved in intestinal maturation by comparing fetal jejunal transcriptomes of two representative porcine breeds (Large White, LW; Meishan, MS) with contrasting neonatal vitality and maturity, at two key time points during late gestation (gestational days 90 and 110). MS and LW sows inseminated with mixed semen (from breed LW and MS) gave birth to both purebred and crossbred fetuses. We hypothesized that part of the differences in neonatal maturity between the two breeds results from distinct developmental profiles of the fetal intestine during late gestation. Reciprocal crossed fetuses were used to analyze the effect of parental genome. Transcriptomic data and 23 phenotypic variables known to be associated with maturity trait were integrated using multivariate analysis with expectation of identifying relevant genes-phenotypic variable relationships involved in intestinal maturation.

Results

A moderate maternal genotype effect, but no paternal genotype effect, was observed on offspring intestinal maturation. Four hundred and four differentially expressed probes, corresponding to 274 differentially expressed genes (DEGs), more specifically involved in the maturation process were further studied. In day 110-MS fetuses, Ingenuity® functional enrichment analysis revealed that 46% of DEGs were involved in glucose and lipid metabolism, cell proliferation, vasculogenesis and hormone synthesis compared to day 90-MS fetuses. Expression of genes involved in immune pathways including phagocytosis, inflammation and defense processes was changed in day 110-LW compared to day 90-LW fetuses (corresponding to 13% of DEGs). The transcriptional regulator PPARGC1A was predicted to be an important regulator of differentially expressed genes in MS. Fetal blood fructose level, intestinal lactase activity and villous height were the best predicted phenotypic variables with probes mostly involved in lipid metabolism, carbohydrate metabolism and cellular movement biological pathways.

Conclusions

Collectively, our findings indicate that the neonatal maturity of pig intestine may rely on functional development of glucose and lipid metabolisms, immune phagocyte differentiation and inflammatory pathways. This process may partially be governed by PPARGC1A.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4001-2) contains supplementary material, which is available to authorized users.

Handling missing rows in multi-omics data integration: multiple imputation in multiple factor analysis framework

Background

In omics data integration studies, it is common, for a variety of reasons, for some individuals to not be present in all data tables. Missing row values are challenging to deal with because most statistical methods cannot be directly applied to incomplete datasets. To overcome this issue, we propose a multiple imputation (MI) approach in a multivariate framework. In this study, we focus on multiple factor analysis (MFA) as a tool to compare and integrate multiple layers of information. MI involves filling the missing rows with plausible values, resulting in M completed datasets. MFA is then applied to each completed dataset to produce M different configurations (the matrices of coordinates of individuals). Finally, the M configurations are combined to yield a single consensus solution.

Results

We assessed the performance of our method, named MI-MFA, on two real omics datasets. Incomplete artificial datasets with different patterns of missingness were created from these data. The MI-MFA results were compared with two other approaches i.e., regularized iterative MFA (RI-MFA) and mean variable imputation (MVI-MFA). For each configuration resulting from these three strategies, the suitability of the solution was determined against the true MFA configuration obtained from the original data and a comprehensive graphical comparison showing how the MI-, RI- or MVI-MFA configurations diverge from the true configuration was produced. Two approaches i.e., confidence ellipses and convex hulls, to visualize and assess the uncertainty due to missing values were also described. We showed how the areas of ellipses and convex hulls increased with the number of missing individuals. A free and easy-to-use code was proposed to implement the MI-MFA method in the R statistical environment.

Conclusions

We believe that MI-MFA provides a useful and attractive method for estimating the coordinates of individuals on the first MFA components despite missing rows. MI-MFA configurations were close to the true configuration even when many individuals were missing in several data tables. This method takes into account the uncertainty of MI-MFA configurations induced by the missing rows, thereby allowing the reliability of the results to be evaluated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1273-5) contains supplementary material, which is available to authorized users.

Exploring transcriptomic diversity in muscle revealed that cellular signaling pathways mainly differentiate five Western porcine breeds

Background

Among transcriptomic studies, those comparing species or populations can increase our understanding of the impact of the evolutionary forces on the differentiation of populations. A particular situation is the one of short evolution time with breeds of a domesticated species that underwent strong selective pressures. In this study, the gene expression diversity across five pig breeds has been explored in muscle. Samples came from: 24 Duroc, 33 Landrace, 41 Large White dam line, 10 Large White sire line and 39 Piétrain. From these animals, 147 muscle samples obtained at slaughter were analyzed using the porcine Agilent 44 K v1 microarray.

Results

A total of 12,358 genes were identified as expressed in muscle after normalization and 1,703 genes were declared differential for at least one breed (FDR < 0.001). The functional analysis highlighted that gene expression diversity is mainly linked to cellular signaling pathways such as the PI3K (phosphoinositide 3-kinase) pathway. The PI3K pathway is known to be involved in the control of development of the skeletal muscle mass by affecting extracellular matrix - receptor interactions, regulation of actin cytoskeleton pathways and some metabolic functions. This study also highlighted 228 spots (171 unique genes) that differentiate the breeds from each other. A common subgroup of 15 genes selected by three statistical methods was able to differentiate Duroc, Large White and Piétrain breeds.

Conclusions

This study on transcriptomic differentiation across Western pig breeds highlighted a global picture: mainly signaling pathways were affected. This result is consistent with the selection objective of increasing muscle mass. These transcriptional changes may indicate selection pressure or simply breed differences which may be driven by human selection. Further work aiming at comparing genetic and transcriptomic diversities would further increase our understanding of the consequences of human impact on livestock species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2259-9) contains supplementary material, which is available to authorized users.

Time course of the response to ACTH in pig: biological and transcriptomic study

Background

HPA axis plays a major role in physiological homeostasis. It is also involved in stress and adaptive response to the environment. In farm animals in general and specifically in pigs, breeding strategies have highly favored production traits such as lean growth rate, feed efficiency and prolificacy at the cost of robustness. On the hypothesis that the HPA axis could contribute to the trade-off between robustness and production traits, we have designed this experiment to explore individual variation in the biological response to the main stress hormone, cortisol, in pigs. We used ACTH injections to trigger production of cortisol in 120 juvenile Large White (LW) pigs from 28 litters and the kinetics of the response was measured with biological variables and whole blood gene expression at 4 time points. A multilevel statistical analysis was used to take into account the longitudinal aspect of the data.

Results

Cortisol level reached its peak 1 h after ACTH injection. White blood cell composition was modified with a decrease of lymphocytes and monocytes and an increase of granulocytes (FDR<0.05). Basal level of cortisol was correlated with birth and weaning weights. Microarray analysis identified 65 unique genes of which expression responded to the injection of ACTH (adjusted P<0.05). These genes were classified into 4 clusters with distinctive kinetics in response to ACTH injection. The first cluster identified genes strongly correlated to cortisol and previously reported as being regulated by glucocorticoids. In particular, DDIT4, DUSP1, FKBP5, IL7R, NFKBIA, PER1, RGS2 and RHOB were shown to be connected to each other by the glucocorticoid receptor NR3C1. Most of the differentially expressed genes that encode transcription factors have not been described yet as being important in transcription networks involved in stress response. Their co-expression may mean co-regulation and they could thus provide new patterns of biomarkers of the individual sensitivity to cortisol.

Conclusions

We identified 65 genes as biological markers of HPA axis activation at the gene expression level. These genes might be candidates for a better understanding of the molecular mechanisms of the stress response.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2118-8) contains supplementary material, which is available to authorized users.

Sperm DNA methylation analysis in swine reveals conserved and species-specific methylation patterns and highlights an altered methylation at the GNAS locus in infertile boars

Male infertility is an increasing health issue in today's society for both human and livestock populations. In livestock, male infertility slows the improvement of animal selection programs and agricultural productivity. There is increasing evidence that epigenetic marks play an important role in the production of good-quality sperm. We therefore screened for specific or common epigenetic signatures of livestock infertility. To do so, we compared DNA methylation level in sperm DNA from fertile and infertile boars. We evaluated first the global level of sperm DNA methylation and found no difference between the two groups of boars. We then selected 42 loci of interest, most of them known to be imprinted in human or mice, and assessed the imprinting status of five of them not previously described in swine tissues: WT1, CNTN3, IMPACT, QPCT, and GRB10. DNA methylation level was then quantified in fertile and infertile boars at these 42 loci. Results from fertile boars indicated that the methylation level of the selected loci is highly conserved between pig, human, and mice, with a few exceptions, including the POU5F1 (OCT4) promoter and RTL1. Comparison between fertile and infertile boars revealed that one imprinted region, the GNAS locus, shows an increase in sperm DNA methylation in three out of eight infertile boars with low semen quality. This increase in DNA methylation is associated with an altered expression of the genes belonging to the GNAS locus, suggesting a new role for GNAS in the proper formation of functional gametes.

Muscle transcriptomic investigation of late fetal development identifies candidate genes for piglet maturity

Background

In pigs, the perinatal period is the most critical time for survival. Piglet maturation, which occurs at the end of gestation, leads to a state of full development after birth. Therefore, maturity is an important determinant of early survival. Skeletal muscle plays a key role in adaptation to extra-uterine life, e.g. glycogen storage and thermoregulation. In this study, we performed microarray analysis to identify the genes and biological processes involved in piglet muscle maturity. Progeny from two breeds with extreme muscle maturity phenotypes were analyzed at two time points during gestation (gestational days 90 and 110). The Large White (LW) breed is a selected breed with an increased rate of mortality at birth, whereas the Meishan (MS) breed produces piglets with extremely low mortality at birth. The impact of the parental genome was analyzed with reciprocal crossed fetuses.

Results

Microarray analysis identified 12,326 differentially expressed probes for gestational age and genotype. Such a high number reflects an important transcriptomic change that occurs between 90 and 110 days of gestation. 2,000 probes, corresponding to 1,120 unique annotated genes, involved more particularly in the maturation process were further studied. Functional enrichment and graph inference studies underlined genes involved in muscular development around 90 days of gestation, and genes involved in metabolic functions, such as gluconeogenesis, around 110 days of gestation. Moreover, a difference in the expression of key genes, e.g. PCK2, LDHA or PGK1, was detected between MS and LW just before birth. Reciprocal crossing analysis resulted in the identification of 472 genes with an expression preferentially regulated by one parental genome. Most of these genes (366) were regulated by the paternal genome. Among these paternally regulated genes, some known imprinted genes, such as MAGEL2 or IGF2, were identified and could have a key role in the maturation process.

Conclusion

These results reveal the biological mechanisms that regulate muscle maturity in piglets. Maturity is also under the conflicting regulation of the parental genomes. Crucial genes, which could explain the biological differences in maturity observed between LW and MS breeds, were identified. These genes could be excellent candidates for a key role in the maturity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-797) contains supplementary material, which is available to authorized users.

The structure of a gene co-expression network reveals biological functions underlying eQTLs

What are the commonalities between genes, whose expression level is partially controlled by eQTL, especially with regard to biological functions? Moreover, how are these genes related to a phenotype of interest? These issues are particularly difficult to address when the genome annotation is incomplete, as is the case for mammalian species. Moreover, the direct link between gene expression and a phenotype of interest may be weak, and thus difficult to handle. In this framework, the use of a co-expression network has proven useful: it is a robust approach for modeling a complex system of genetic regulations, and to infer knowledge for yet unknown genes. In this article, a case study was conducted with a mammalian species. It showed that the use of a co-expression network based on partial correlation, combined with a relevant clustering of nodes, leads to an enrichment of biological functions of around 83%. Moreover, the use of a spatial statistics approach allowed us to superimpose additional information related to a phenotype; this lead to highlighting specific genes or gene clusters that are related to the network structure and the phenotype. Three main results are worth noting: first, key genes were highlighted as a potential focus for forthcoming biological experiments; second, a set of biological functions, which support a list of genes under partial eQTL control, was set up by an overview of the global structure of the gene expression network; third, pH was found correlated with gene clusters, and then with related biological functions, as a result of a spatial analysis of the network topology.

Phenotypic prediction based on metabolomic data for growing pigs from three main European breeds

Predicting phenotypes is a statistical and biotechnical challenge, both in medicine (predicting an illness) and animal breeding (predicting the carcass economical value on a young living animal). High-throughput fine phenotyping is possible using metabolomics, which describes the global metabolic status of an individual, and is the closest to the terminal phenotype. The purpose of this work was to quantify the prediction power of metabolomic profiles for commonly used production phenotypes from a single blood sample from growing pigs. Several statistical approaches were investigated and compared on the basis of cross validation: raw data vs. signal preprocessing (wavelet transformation), with a single-feature selection method. The best results in terms of prediction accuracy were obtained when data were preprocessed using wavelet transformations on the Daubechies basis. The phenotypes related to meat quality were not well predicted because the blood sample was taken some time before slaughter, and slaughter is known to have a strong influence on these traits. By contrast, phenotypes of potential economic interest (e.g., lean meat percentage and ADFI) were well predicted (R(2) = 0.7; P < 0.0001) using metabolomic data.

Genetic variability of transcript abundance in pig peri-mortem skeletal muscle: eQTL localized genes involved in stress response, cell death, muscle disorders and metabolism

Background

The genetics of transcript-level variation is an exciting field that has recently given rise to many studies. Genetical genomics studies have mainly focused on cell lines, blood cells or adipose tissues, from human clinical samples or mice inbred lines. Few eQTL studies have focused on animal tissues sampled from outbred populations to reflect natural genetic variation of gene expression levels in animals. In this work, we analyzed gene expression in a whole tissue, pig skeletal muscle sampled from individuals from a half sib F2 family shortly after slaughtering.

Results

QTL detection on transcriptome measurements was performed on a family structured population. The analysis identified 335 eQTLs affecting the expression of 272 transcripts. The ontologic annotation of these eQTLs revealed an over-representation of genes encoding proteins involved in processes that are expected to be induced during muscle development and metabolism, cell morphology, assembly and organization and also in stress response and apoptosis. A gene functional network approach was used to evidence existing biological relationships between all the genes whose expression levels are influenced by eQTLs. eQTLs localization revealed a significant clustered organization of about half the genes located on segments of chromosome 1, 2, 10, 13, 16, and 18. Finally, the combined expression and genetic approaches pointed to putative cis-drivers of gene expression programs in skeletal muscle as COQ4 (SSC1), LOC100513192 (SSC18) where both the gene transcription unit and the eQTL affecting its expression level were shown to be localized in the same genomic region. This suggests cis-causing genetic polymorphims affecting gene expression levels, with (e.g. COQ4) or without (e.g. LOC100513192) potential pleiotropic effects that affect the expression of other genes (cluster of trans-eQTLs).

Conclusion

Genetic analysis of transcription levels revealed dependence among molecular phenotypes as being affected by variation at the same loci. We observed the genetic variation of molecular phenotypes in a specific situation of cellular stress thus contributing to a better description of muscle physiologic response. In turn, this suggests that large amounts of genetic variation, mediated through transcriptional networks, can drive transient cell response phenotypes and contribute to organismal adaptative potential.

Transcriptome profiling of sheep granulosa cells and oocytes during early follicular development obtained by laser capture microdissection

Background

Successful achievement of early folliculogenesis is crucial for female reproductive function. The process is finely regulated by cell-cell interactions and by the coordinated expression of genes in both the oocyte and in granulosa cells. Despite many studies, little is known about the cell-specific gene expression driving early folliculogenesis. The very small size of these follicles and the mixture of types of follicles within the developing ovary make the experimental study of isolated follicular components very difficult.

The recently developed laser capture microdissection (LCM) technique coupled with microarray experiments is a promising way to address the molecular profile of pure cell populations. However, one main challenge was to preserve the RNA quality during the isolation of single cells or groups of cells and also to obtain sufficient amounts of RNA.

Using a new LCM method, we describe here the separate expression profiles of oocytes and follicular cells during the first stages of sheep folliculogenesis.

Results

We developed a new tissue fixation protocol ensuring efficient single cell capture and RNA integrity during the microdissection procedure. Enrichment in specific cell types was controlled by qRT-PCR analysis of known genes: six oocyte-specific genes (SOHLH2, MAEL, MATER, VASA, GDF9, BMP15) and three granulosa cell-specific genes (KL, GATA4, AMH).

A global gene expression profile for each follicular compartment during early developmental stages was identified here for the first time, using a bovine Affymetrix chip. Most notably, the granulosa cell dataset is unique to date. The comparison of oocyte vs. follicular cell transcriptomes revealed 1050 transcripts specific to the granulosa cell and 759 specific to the oocyte.

Functional analyses allowed the characterization of the three main cellular events involved in early folliculogenesis and confirmed the relevance and potential of LCM-derived RNA.

Conclusions

The ovary is a complex mixture of different cell types. Distinct cell populations need therefore to be analyzed for a better understanding of their potential interactions. LCM and microarray analysis allowed us to identify novel gene expression patterns in follicular cells at different stages and in oocyte populations.

Expression levels of 25 genes in liver and testis located in a QTL region for androstenone on SSC7q1.2

A quantitative trait locus (QTL) for boar fat androstenone levels has been identified near the SSC7 centromere in a Large White × Meishan cross. Backcrosses were produced to isolate the Chinese haplotype in a European genetic background. The expression of 25 genes from the QTL region was studied in the testes and livers of 5-month-old backcross boars, with the aim of identifying the causal gene. Using Fluidigm, a new high-throughput technology, the expression of 25 genes was measured in a single real-time PCR experiment. This study found six significantly down-regulated genes (C6ORF106, C6ORF81, CLPS, SLC26A8, SRPK1 and MAPK14) in the testes of MS-LW backcross boars. However, according to current knowledge, none of the genes appear to be related to androstenone metabolism. In the livers, none of the genes were significantly up- or down-regulated, including TEAD3, which was previously designated as a possible candidate to explain this QTL.

Pathway results from the chicken data set using GOTM, Pathway Studio and Ingenuity softwares

Background

As presented in the introduction paper, three sets of differentially regulated genes were found after the analysis of the chicken infection data set from EADGENE. Different methods were used to interpret these results.

Results

GOTM, Pathway Studio and Ingenuity softwares were used to investigate the three lists of genes. The three softwares allowed the analysis of the data and highlighted different networks. However, only one set of genes, showing a differential expression between primary and secondary response gave significant biological interpretation.

Conclusion

Combining these databases that were developed independently on different annotation sources supplies a useful tool for a global biological interpretation of microarray data, even if they may contain some imperfections (e.g. gene not or not well annotated).

Methods for interpreting lists of affected genes obtained in a DNA microarray experiment

Background

The aim of this paper was to describe and compare the methods used and the results obtained by the participants in a joint EADGENE (European Animal Disease Genomic Network of Excellence) and SABRE (Cutting Edge Genomics for Sustainable Animal Breeding) workshop focusing on post analysis of microarray data. The participating groups were provided with identical lists of microarray probes, including test statistics for three different contrasts, and the normalised log-ratios for each array, to be used as the starting point for interpreting the affected probes. The data originated from a microarray experiment conducted to study the host reactions in broilers occurring shortly after a secondary challenge with either a homologous or heterologous species of Eimeria.

Results

Several conceptually different analytical approaches, using both commercial and public available software, were applied by the participating groups. The following tools were used: Ingenuity Pathway Analysis, MAPPFinder, LIMMA, GOstats, GOEAST, GOTM, Globaltest, TopGO, ArrayUnlock, Pathway Studio, GIST and AnnotationDbi. The main focus of the approaches was to utilise the relation between probes/genes and their gene ontology and pathways to interpret the affected probes/genes. The lack of a well-annotated chicken genome did though limit the possibilities to fully explore the tools. The main results from these analyses showed that the biological interpretation is highly dependent on the statistical method used but that some common biological conclusions could be reached.

Conclusion

It is highly recommended to test different analytical methods on the same data set and compare the results to obtain a reliable biological interpretation of the affected genes in a DNA microarray experiment.

Using microarrays to identify positional candidate genes for QTL: the case study of ACTH response in pigs

BACKGROUND

Microarray studies can supplement QTL studies by suggesting potential candidate genes in the QTL regions, which by themselves are too large to provide a limited selection of candidate genes. Here we provide a case study where we explore ways to integrate QTL data and microarray data for the pig, which has only a partial genome sequence. We outline various procedures to localize differentially expressed genes on the pig genome and link this with information on published QTL. The starting point is a set of 237 differentially expressed cDNA clones in adrenal tissue from two pig breeds, before and after treatment with adrenocorticotropic hormone (ACTH).

RESULTS

Different approaches to localize the differentially expressed (DE) genes to the pig genome showed different levels of success and a clear lack of concordance for some genes between the various approaches. For a focused analysis on 12 genes, overlapping QTL from the public domain were presented. Also, differentially expressed genes underlying QTL for ACTH response were described. Using the latest version of the draft sequence, the differentially expressed genes were mapped to the pig genome. This enabled co-location of DE genes and previously studied QTL regions, but the draft genome sequence is still incomplete and will contain many errors. A further step to explore links between DE genes and QTL at the pathway level was largely unsuccessful due to the lack of annotation of the pig genome. This could be improved by further comparative mapping analyses but this would be time consuming.

CONCLUSION

This paper provides a case study for the integration of QTL data and microarray data for a species with limited genome sequence information and annotation. The results illustrate the challenges that must be addressed but also provide a roadmap for future work that is applicable to other non-model species.

Characterization of porcine ASB6 gene and transcripts-comparison of mammalian transcripts

A member of the porcine Ankyrin repeat and suppressor of cytokine signaling (SOCS) Box protein family (ASB), designed as ASB6, was sequenced and the genomic organization of the six exons was determined. We present here a detailed analysis of ASB6 transcripts in pigs. We demonstrate the existence of an alternative transcript resulting from intron retention. This secondary transcript, if functional, encodes a protein without SOCS box. A comparison of mammalian ASB6 transcripts is performed to demonstrate the importance of transcripts encoding for a truncated ASB6 protein.

A muscle transcriptome analysis identifies positional candidate genes for a complex trait in pig

Muscle tenderness is an important complex trait for meat quality and thus for genetic improvement through animal breeding. However, the physiological or genetic control of tenderness development in muscle is still poorly understood. In this work, using transcriptome analysis, we found a relationship between gene expression variability and tenderness. Muscle (longissimus dorsi) samples from 30 F(2) pigs were characterized by Warner-Bratzler Shear Force (WBSF) on cooked meat as a measurement of muscle tenderness. Gene expression levels were measured using microarrays for 17 muscle samples selected to represent a range of WBSF values. Using a linear regression model, we determined that samples with WBSF values above 30 N could be effectively analysed for genes exhibiting a significant association of their expression level on shear force (false discovery rate <0.05). These genes were shown to be involved in three functional networks: cell cycle, energy metabolism and muscle development. Twenty-two genes were mapped on the pig genome and 12 were found to be located in regions previously reported to contain quantitative trait loci (QTL) affecting pig meat tenderness (chromosomes 2, 6 and 13). Some genes appear therefore as positional candidate genes for QTL.

Gene array and real time PCR analysis of the adrenal sensitivity to adrenocorticotropic hormone in pig

Background

Variability in hypothalamic-pituitary-adrenal (HPA) axis activity has been shown to be influenced by genetic factors and related to great metabolic differences such as obesity. The aim of this study was to investigate molecular bases of genetic variability of the adrenal sensitivity to ACTH, a major source of variability, in Meishan (MS) and Large White (LW) pigs, MS being reported to exhibit higher basal cortisol levels, response to ACTH and fatness than LW. A pig cDNA microarray was used to identify changes in gene expression in basal conditions and in response to ACTH stimulation.

Results

Genotype and/or ACTH affected the expression of 211 genes related to transcription, cell growth/maintenance, signal transduction, cell structure/adhesion/extra cellular matrix and protein kinase/phosphatase activity. No change in the expression of known key regulator proteins of the ACTH signaling pathway or of steroidogenic enzymes was found. However, Mdh2, Sdha, Suclg2, genes involved in the tricarboxylic acid (TCA) pathway, were over-expressed in MS pigs. Higher TCA cycle activity in MS than in LW may thus result in higher steroidogenic activity and thus explain the typically higher cortisol levels in MS compared to LW. Moreover, up-regulation of Star and Ldlr genes in MS and/or in response to ACTH suggest that differences in the adrenal function between MS and LW may also involve mechanisms requisite for cholesterol supply to steroidogenesis.

Conclusion

The present study provides new potential candidate genes to explain genetic variations in the adrenal sensitivity to ACTH and better understand relationship between HPA axis activity and obesity.

Longitudinal Analysis of Gene Expression in Porcine Skeletal Muscle After Post-Injection Local Injury

PURPOSE

The purpose of this study is to describe the time course of gene expression in a skeletal muscle local injury induced by an intramuscular (IM) injection, and to compare the dynamics of gene expression with pathological events.

MATERIALS AND METHODS

Ten piglets received 4 IM injections of propylene glycol in the longissimus dorsi muscles 6 h, 2, 7, and 21 days before euthanasia, where control and injected muscle sites were sampled for RNA isolation and microscopic examination. The hybridization of nylon cDNA microarrays was carried out with radioactive probes obtained from the muscle RNA.

RESULTS

153 genes were found under- or over-expressed at least once among the investigated time-conditions. The eight most discriminant genes were also identified: Two genes (GTP-binding protein RAD and Ankyrin repeat domain protein) were over-expressed at 6 h and six genes between 2 and 21 days (Osteonectin, Fibronectin, Matrix metalloproteinase-2, Collagen alpha 1(I) chain, Collagen alpha 2(I) chain, and Thymosin beta-4). Necrosis, inflammation and regeneration were observed through both the dynamics of gene expression profiles and through the microscopic examinations.

CONCLUSION

Our data demonstrate that several pathways are involved in post-injection muscle injury, and that necrosis, inflammation and regeneration are not sequential but occur in parallel.

Determination of PRKAG1 coding sequence and mapping of PRKAG1 and PRKAG2 relatively to porcine back fat thickness QTL

PRKAG1, PRKAG2 and PRKAG3 encode three isoforms of AMP-activated protein kinase gamma chain. A major effect on meat quality and a medium effect on back fat thickness of the RN- mutation in the PRKAG3 gene has previously been reported. We have now mapped PRKAG1 and PRKAG2 at expected locations on SSC5 and SSC18 by analysis of radiation hybrids (IMpRH panel). PRKAG2 has been mapped in a region where no quantitative trait loci (QTL) has been reported. PRKAG1 has been mapped close to (but probably outside) a region containing a QTL influencing fatness traits. We have determined the full coding sequence of PRKAG1. No missense mutation was identified when comparing the coding sequence of one Meishan and one Large White boars. Further work is, however, required to determine if a polymorphism in PRKAG1 could be responsible for a part of the variability observed on fatness traits.

Production of Phosphatidylinositol 5-Phosphate by the Phosphoinositide 3-Phosphatase Myotubularin in Mammalian Cells

MTM1, the gene encoding myotubularin (MTM1), is mutated in the X-linked myotubular myopathy (XLMTM), a severe genetic muscular disorder. MTM1 is a phosphoinositide phosphatase hydrolyzing phosphatidylinositol 3-phosphate (PtdIns(3)P) in yeast and in vitro. Because this lipid is implicated in the regulation of vesicular trafficking, we used established cell lines from XLMTM patients to evaluate whether the lack of endogenous MTM1 expression could affect PtdIns(3)P labeling patterns. Our results showed that the vesicular trafficking related to early endosomes was not significantly affected in the XLMTM cell lines compared with control cells. However, in addition to PtdIns(3)P, we found that MTM1 can hydrolyze phosphatidylinositol 3,5-bisphosphate both in vitro and in mammalian cells. Using a mass assay, we demonstrated that the product generated is phosphatidylinositol 5-phosphate (PtdIns(5)P), a recently discovered phosphoinositide, the function of which is still unknown. In L6 myotubes overexpressing MTM1, hyperosmotic shock induced an increase in the mass level of PtdIns(5)P that was reduced by 50% upon overexpression of the MTM1 inactive mutant D278A. These data demonstrate for the first time a role for MTM1 in the production of PtdIns(5)P in mammalian cells, suggesting that the lack of transformation of phosphatidylinositol 3,5-bisphosphate into PtdIns(5)P might be an important component in the etiology of myotubular myopathy.

Functional redundancy in the myotubularin family

Myotubularin-related genes define a novel highly conserved family of eukaryotic proteins of at least 11 human members. The hMTM1 gene that codes for myotubularin is mutated in X-linked myotubular myopathy, a severe congenital disease. Recently, we and others have characterized myotubularin as a potent and specific phosphatidylinositol 3-phosphate 3-phosphatase. In the present study we investigated the lipid phosphatase activity and the subcellular localization of two other members of the family, hMTMR2 protein that is mutated in the demyelinating neuropathy Charcot-Marie-Tooth type 4B and the FYVE-finger containing hMTMR3 protein. Our results show that both proteins are potent phosphatidylinositol 3-phosphate 3-phosphatases either in vitro or in yeast where they interfered with vesicular trafficking. Their localization is mainly cytoplasmic, with however strong labeling of Rac-inducible plasma membrane ruffles. The fact that the ubiquitously expressed hMTM1 and hMTMR2 genes are involved in different pathologies indicates that despite their shared enzymatic activity, they are not functionally redundant, at least in certain cell types. This might be explained by subtle differences in expression and/or in recruitment and regulation at their specific site of action.

Identification by differential display of a chicken tolloid-related metalloprotease specifically expressed in the caudal notochord

While the ventralizing factor Sonic hedgehog is expressed in the entire notochord (Development 121 (1995) 2537) the latter displays distinct ventralizing activities along its rostrocaudal axis. Hence, in HH stage-10 chicken embryo, the caudal notochord exhibits floor plate inducing capacities lost by rostral regions (Development 117 (1993) 205). Therefore, we hypothesize that the caudal notochord produces some cofactors which may contribute to its ventralizing properties. In order to identify such molecules we applied the differential display strategy and isolated a secreted Tolloid-related metalloprotease displaying a regionalized expression in the notochord.

Nonselective coupling of the human mu-opioid receptor to multiple inhibitory G-protein isoforms

The human mu-opioid receptor was expressed in Saccharomyces cerevisiae. Binding of [3H]diprenorphine to yeast spheroplasts was specific and saturable (Kd = 1 nm, Bmax = 0.2-1 pmol x mg-1 of membrane proteins). Inhibition of [3H]diprenorphine binding by antagonists and agonists with varying opioid selectivities (mu, delta and kappa) occurred with the same order of potency as in mammalian tissues. Affinities of antagonists were the same with yeast spheroplasts as in reference tissues whereas those of agonists, except etorphine and buprenorphine, were 10-fold to 100-fold lower. Addition of heterotrimeric Gi,o-proteins purified from bovine brain shifted the mu-opioid receptor into a high-affinity state for agonists. Using individually purified Galpha-subunits re-associated with betagamma-dimers, we showed that alphao1, alphao2, alphai1, alphai2 and alphai3 reconstituted high-affinity agonist binding with equal efficiency. This suggests that the structural determinants of the mu-opioid receptor responsible for G-protein coupling are not able to confer a high degree of specificity towards any member of the Gi,o family. The selective effects of opioid observed in specialized tissues upon opioid stimulation may be a result of regulation of G-protein activity by cell-specific factors which should conveniently be analysed using the reconstitution assay described here.

CYP76C2, an Arabidopsis thaliana cytochrome P450 gene expressed during hypersensitive and developmental cell death

The characterisation of an Arabidopsis thaliana cytochrome P450-encoding cDNA clone, B72, preferentially expressed during the hypersensitive response (HR) provoked by the bacterial pathogen Pseudomonas syringae pathovar maculicola, is reported. The B72 cDNA clone corresponded to the CYP76C2 gene, which belongs to a small multigene family comprising four genes. HR-triggering bacteria harbouring different avirulence genes induced the accumulation of transcripts of this P450 gene. CYP76C2 gene expression was moreover associated with various processes leading to cell death such as leaf senescence, ageing of cell cultures, wounding as well as with treatment with the necrotising heavy metal salt, lead nitrate.

Genetic characterization of RRS1, a recessive locus in Arabidopsis thaliana that confers resistance to the bacterial soilborne pathogen Ralstonia solanacearum

The soilborne, vascular pathogen Ralstonia solanacearum, the causative agent of bacterial wilt, was shown to infect a range of Arabidopsis thaliana accessions. The pathogen was capable of infecting the Col-5 accession in an hrp-dependent manner, following root inoculation. Elevated bacterial population levels were found in leaves of Col-5, 4 to 5 days after root inoculation by the GMI1000 strain. Bacteria were found predominantly in the xylem vessels and spread systematically throughout the plant. The Nd-1 accession of A. thaliana was resistant to the GMI1000 strain of R. solanacearum. Bacterial concentrations detected in leaves of Nd-1, inoculated with an hrp+ strain of R. solanacearum, were only slightly higher than those detected in the susceptible accession, Col-5, following inoculation with a strain whose hrp gene cluster was deleted. Leaf inoculation of the GMI1000 strain on the resistant accession Nd-1 induced the formation of lesions in the older leaves of the rosette whereas the same strain of R. solanacearum provoked complete wilting of Col-5. Resistance to strain GMI1000 of R. solanacearum segregated as a simply inherited recessive trait in a genetic cross between Col-5 and Nd-1. F9 recombinant inbred lines generated between these two accessions were used to map a locus, RRS1, that was the major determinant of resistance between restriction fragment length polymorphism markers mi83 and mi61 on chromosome V. This region of the A. thaliana genome is known to contain many other pathogen recognition capabilities.