Automated SNP Detection in Expressed Sequence Tags: Statistical Considerations and Application to Maritime Pine Sequences

We developed an automated pipeline for the detection of single nucleotide polymorphisms (SNPs) in expressed sequence tag (EST) data sets, by combining three DNA sequence analysis programs: Phred, Phrap and PolyBayes. This application requires access to the individual electrophoregram traces. First, a reference set of 65 SNPs was obtained from the sequencing of 30 gametes in 13 maritime pine (Pinus pinaster Ait.) gene fragments (6671 bp), resulting in a frequency of 1 SNP every 102.6 bp. Second, parameters of the three programs were optimized in order to retrieve as many true SNPs, while keeping the rate of false positive as low as possible. Overall, the efficiency of detection of true SNPs was 83.1%. However, this rate varied largely as a function of the rare SNP allele frequency: down to 41% for rare SNP alleles (frequency < 10%), up to 98% for allele frequencies above 10%. Third, the detection method was applied to the 18498 assembled maritime pine (Pinus pinaster Ait.) ESTs, allowing to identify a total of 1400 candidate SNPs, in contigs containing between 4 and 20 sequence reads. These genetic resources, described for the first time in a forest tree species, were made available at http://www.pierroton.inra/genetics/Pinesnps. We also derived an analytical expression for the SNP detection probability as a function of the SNP allele frequency, the number of haploid genomes used to generate the EST sequence database, and the sample size of the contigs considered for SNP detection. The frequency of the SNP allele was shown to be the main factor influencing the probability of SNP detection.     

A profibrotic function of IL-12p40 in experimental pulmonary fibrosis

The p40 subunit of IL-12 (IL-12p40), but not the heterodimeric form IL-12p70, is secreted during the development of silica-induced lung fibrosis in C57BL/6 mice. To delineate the contribution of IL-12p40 to the lung inflammatory and fibrotic processes, we compared the pulmonary responses with silica particles of IL-12p35-deficient mice (IL-12p35(-/-), able to produce IL-12p40) and IL-12p40-deficient mice (IL-12p40(-/-)). IL-12p35(-/-) and IL-12p40(-/-) animals developed strikingly contrasting responses to silica in comparison with wild-type C57BL/6 mice. Although the IL-12p40(-/-) mice exhibited limited inflammatory and fibrotic reactions, the IL-12p35(-/-) mice presented a robust and well-developed pulmonary inflammation and fibrosis. Furthermore, the silica-induced increase in lung IL-12p40 content was significantly higher in IL-12p35(-/-) mice than in wild-type controls, and was associated with extensive lung fibrosis and pulmonary macrophage infiltration. The contrasting responses observed between these two IL-12 subunit-deficient murine strains were not accompanied by a strict type 1 or type 2 polarization as estimated by the measurements of lung IFN-gamma/IgG2a and IL-4/IgG1 content. In vitro proliferation, type I collagen expression, as well as myofibroblast differentiation of purified pulmonary fibroblasts were not affected by treatment with exogenous rIL-12p40. In vivo, supplementation with rIL-12p40 restored the impaired pulmonary fibrotic response and macrophage accumulation in silica-treated IL-12p40(-/-) mice, and also promoted fibrosis and macrophage influx in wild-type mice. Together, our data suggest that IL-12p40 plays an important role in silica-induced pulmonary inflammation and fibrosis, possibly by exacerbating macrophage recruitment.     

Bidirectional integrative regulation of Cav1.2 calcium channel by microRNA miR-103: role in pain

Chronic pain states are characterized by long-term sensitization of spinal cord neurons that relay nociceptive information to the brain. Among the mechanisms involved, up-regulation of Cav1.2-comprising L-type calcium channel (Cav1.2-LTC) in spinal dorsal horn have a crucial role in chronic neuropathic pain. Here, we address a mechanism of translational regulation of this calcium channel. Translational regulation by microRNAs is a key factor in the expression and function of eukaryotic genomes. Because perfect matching to target sequence is not required for inhibition, theoretically, microRNAs could regulate simultaneously multiple mRNAs. We show here that a single microRNA, miR-103, simultaneously regulates the expression of the three subunits forming Cav1.2-LTC in a novel integrative regulation. This regulation is bidirectional since knocking-down or over-expressing miR-103, respectively, up- or down-regulate the level of Cav1.2-LTC translation. Functionally, we show that miR-103 knockdown in naive rats results in hypersensitivity to pain. Moreover, we demonstrate that miR-103 is down-regulated in neuropathic animals and that miR-103 intrathecal applications successfully relieve pain, identifying miR-103 as a novel possible therapeutic target in neuropathic chronic pain.     

A 13C NMR spectrometric method for the determination of intramolecular δ13C values in fructose from plant sucrose samples

Recent developments in (13) C NMR spectrometry have allowed the determination of intramolecular (13) C/(12) C ratios with high precision. However, the analysis of carbohydrates requires their derivatization to constrain the anomeric carbon. Fructose has proved to be particularly problematic because of a byproduct occurring during derivatization and the complexity of the NMR spectrum of the derivative. Here, we describe a method to determine the intramolecular (13) C/(12) C ratios in fructose by (13) C NMR analysis of the acetyl-isopropylidene derivative. We have applied this method to measure the intramolecular (13) C/(12) C distribution in the fructosyl moiety of sucrose and have compared this with that in the glucosyl moiety. Three prominent features stand out. First, in sucrose from both C(3) and C(4) plants, the C-1 and C-2 positions of the glucosyl and fructosyl moieties are markedly different. Second, these positions in C(3) and C(4) plants show a similar profile. Third, the glucosyl and fructosyl moieties of sucrose from Crassulacean acid metabolism (CAM) metabolism have a different profile. These contrasting values can be interpreted as a result of the isotopic selectivity of enzymes that break or make covalent bonds in glucose metabolism, whereas the distinctive (13) C pattern in CAM sucrose probably indicates a substantial contribution of gluconeogenesis to glucose synthesis.     

Colonization and extinction dynamics of an annual plant metapopulation in an urban environment

The metapopulation framework considers that the spatiotemporal distribution of organisms results from a balance between the colonization and extinction of populations in a suitable and discrete habitat network. Recent spatially realistic metapopulation models have allowed patch dynamics to be investigated in natural populations but such models have rarely been applied to plants. Using a simple urban fragmented population system in which favourable habitat can be easily mapped, we studied patch dynamics in the annual plant Crepis sancta (Asteraceae). Using stochastic patch occupancy models (SPOMs) and multi‐year occupancy data we dissected extinction and colonization patterns in our system. Overall, our data were consistent with two distinct metapopulation scenarios. A metapopulation (sensu stricto) dynamic in which colonization occurs over a short distance and extinction is lowered by nearby occupied patches (rescue effect) was found in a set of patches close to the city centre, while a propagule rain model in which colonization occurs from a large external population was most consistent with data from other networks. Overall, the study highlights the importance of external seed sources in urban patch dynamics. Our analysis emphasizes the fact that plant distributions are governed not only by habitat properties but also by the intrinsic properties of colonization and dispersal of species. The metapopulation approach provides a valuable tool for understanding how colonization and extinction shape occupancy patterns in highly fragmented plant populations. Finally, this study points to the potential utility of more complex plant metapopulation models than traditionally used for analysing ecological and evolutionary processes in natural metapopulations.     

New consistent QTL in pea associated with partial resistance to Aphanomyces euteiches in multiple French and American environments

Partial resistances, often controlled by quantitative trait loci (QTL), are considered to be more durable than monogenic resistances. Therefore, a precursor to developing efficient breeding programs for polygenic resistance to pathogens should be a greater understanding of genetic diversity and stability of resistance QTL in plants. In this study, we deciphered the diversity and stability of resistance QTL to Aphanomyces euteiches in pea towards pathogen variability, environments and scoring criteria, from two new sources of partial resistance (PI?180693 and 552), effective in French and USA infested fields. Two mapping populations of 178 recombinant inbred lines each, derived from crosses between 552 or PI 180693 (partially resistant) and Baccara (susceptible), were used to identify QTL for Aphanomyces root rot resistance in controlled and in multiple French and USA field conditions using several resistance criteria. We identified a total of 135 additive-effect QTL corresponding to 23 genomic regions and 13 significant epistatic interactions associated with partial resistance to A.?euteiches in pea. Among the 23 additive-effect genomic regions identified, five were consistently detected, and showed highly stable effects towards A.?euteiches strains, environments, resistance criteria, condition tests and RIL populations studied. These results confirm the complexity of inheritance of partial resistance to A.?euteiches in pea and provide good bases for the choice of consistent QTL to use in marker-assisted selection schemes to increase current levels of resistance to A.?euteiches in pea breeding programs.     

Integration of detailed modules in a core model of body fluid homeostasis and blood pressure regulation

This paper presents a contribution to the definition of the interfaces required to perform heterogeneous model integration in the context of integrative physiology. A formalization of the model integration problem is proposed and a coupling method is presented. The extension of the classic Guyton model, a multi-organ, integrated systems model of blood pressure regulation, is used as an example of the application of the proposed method. To this end, the Guyton model has been restructured, extensive sensitivity analyses have been performed, and appropriate transformations have been applied to replace a subset of its constituting modules by integrating a pulsatile heart and an updated representation of the renin-angiotensin system. Simulation results of the extended integrated model are presented and the impacts of their integration within the original model are evaluated.     

Structural and functional analysis of the Rous Sarcoma virus negative regulator of splicing and demonstration of its activation by the 9G8 SR protein

Retroviruses require both spliced and unspliced RNAs for replication. Accumulation of Rous Sarcoma virus (RSV) unspliced RNA depends upon the negative regulator of splicing (NRS). Its 5'-part is considered as an ESE binding SR proteins. Its 3'-part contains a decoy 5'-splice site (ss), which inhibits splicing at the bona fide 5'-ss. Only the 3D structure of a small NRS fragment had been experimentally studied. Here, by chemical and enzymatic probing, we determine the 2D structure of the entire RSV NRS. Structural analysis of other avian NRSs and comparison with all sequenced avian NRSs is in favour of a phylogenetic conservation of the NRS 2D structure. By combination of approaches: (i) in vitro and in cellulo splicing assays, (ii) footprinting assays and (iii) purification and analysis of reconstituted RNP complex, we define a small NRS element retaining splicing inhibitory property. We also demonstrate the capability of the SR protein 9G8 to increase NRS activity in vitro and in cellulo. Altogether these data bring new insights on how NRS fine tune splicing activity.     

Diversity and evolution of the Hordeum murinum polyploid complex in Algeria

Population diversity and evolutionary relationships in the Hordeum murinum L. polyploid complex were explored in contrasted bioclimatic conditions from Algeria. A multidisciplinary approach based on morphological, cytogenetic, and molecular data was conducted on a large population sampling. Distribution of diploids (subsp. glaucum) and tetraploids (subsp. leporinum) revealed a strong correlation with a North-South aridity gradient. Most cytotypes exhibit regular meiosis with variable irregularities in some tetraploid populations. Morphological analyses indicate no differentiation among taxa but high variability correlated with bioclimatic parameters. Two and three different nuclear sequences (gene coding for an unspliced genomic protein kinase domain) were isolated in tetraploid and hexaploid cytotypes, respectively, among which one was identical with that found in the diploid subsp. glaucum. The tetraploids (subsp. leporinum and subsp. murinum) do not exhibit additivity for 5S and 45S rDNA loci comparative with the number observed in the related diploid (subsp. glaucum). The subgenomes in the tetraploid taxa could not be differentiated using genomic in situ hybridization (GISH). Results support an allotetraploid origin for subsp. leporinum and subsp. murinum that derives from the diploid subsp. glaucum and another unidentified diploid parent. The hexaploid (subsp. leporinum) has an allohexaploid origin involving the two genomes present in the allotetraploids and another unidentified third diploid progenitor.     

Structure of the Mycobacterium tuberculosis OmpATb protein: a model of an oligomeric channel in the mycobacterial cell wall

The pore-forming outer membrane protein OmpATb from Mycobacterium tuberculosis is a virulence factor required for acid resistance in host phagosomes. In this study, we determined the 3D structure of OmpATb by NMR in solution. We found that OmpATb is composed of two independent domains separated by a proline-rich hinge region. As expected, the high-resolution structure of the C-terminal domain (OmpATb(198-326)) revealed a module structurally related to other OmpA-like proteins from Gram-negative bacteria. The N-terminal domain of OmpATb (73-204), which is sufficient to form channels in planar lipid bilayers, exhibits a fold, which belongs to the α+β sandwich class fold. Its peculiarity is to be composed of two overlapping subdomains linked via a BON (Bacterial OsmY and Nodulation) domain initially identified in bacterial proteins predicted to interact with phospholipids. Although OmpATb(73-204) is highly water soluble, current-voltage measurements demonstrate that it is able to form conducting pores in model membranes. A HADDOCK modeling of the NMR data gathered on the major monomeric form and on the minor oligomeric populations of OmpATb(73-204) suggest that OmpATb(73-204) can form oligomeric rings able to insert into phospholipid membrane, similar to related proteins from the Type III secretion systems, which form multisubunits membrane-associated rings at the basal body of the secretion machinery.