Genomic and Phenomic Study of Mammary Pathogenic Escherichia coli

Escherichia coli is a major etiological agent of intra-mammary infections (IMI) in cows, leading to acute mastitis and causing great economic losses in dairy production worldwide. Particular strains cause persistent IMI, leading to recurrent mastitis. Virulence factors of mammary pathogenic E. coli (MPEC) involved pathogenesis of mastitis as well as those differentiating strains causing acute or persistent mastitis are largely unknown. This study aimed to identify virulence markers in MPEC through whole genome and phenome comparative analysis. MPEC strains causing acute (VL2874 and P4) or persistent (VL2732) mastitis were compared to an environmental strain (K71) and to the genomes of strains representing different E. coli pathotypes. Intra-mammary challenge in mice confirmed experimentally that the strains studied here have different pathogenic potential, and that the environmental strain K71 is non-pathogenic in the mammary gland. Analysis of whole genome sequences and predicted proteomes revealed high similarity among MPEC, whereas MPEC significantly differed from the non-mammary pathogenic strain K71, and from E. coli genomes from other pathotypes. Functional features identified in MPEC genomes and lacking in the non-mammary pathogenic strain were associated with synthesis of lipopolysaccharide and other membrane antigens, ferric-dicitrate iron acquisition and sugars metabolism. Features associated with cytotoxicity or intra-cellular survival were found specifically in the genomes of strains from severe and acute (VL2874) or persistent (VL2732) mastitis, respectively. MPEC genomes were relatively similar to strain K-12, which was subsequently shown here to be possibly pathogenic in the mammary gland. Phenome analysis showed that the persistent MPEC was the most versatile in terms of nutrients metabolized and acute MPEC the least. Among phenotypes unique to MPEC compared to the non-mammary pathogenic strain were uric acid and D-serine metabolism. This study reveals virulence factors and phenotypic characteristics of MPEC that may play a role in pathogenesis of E. coli mastitis.     

Pathway Thermodynamics Highlights Kinetic Obstacles in Central Metabolism

In metabolism research, thermodynamics is usually used to determine the directionality of a reaction or the feasibility of a pathway. However, the relationship between thermodynamic potentials and fluxes is not limited to questions of directionality: thermodynamics also affects the kinetics of reactions through the flux-force relationship, which states that the logarithm of the ratio between the forward and reverse fluxes is directly proportional to the change in Gibbs energy due to a reaction (Δ_rG′). Accordingly, if an enzyme catalyzes a reaction with a Δ_rG′ of -5.7 kJ/mol then the forward flux will be roughly ten times the reverse flux. As Δ_rG′ approaches equilibrium (Δ_rG′ = 0 kJ/mol), exponentially more enzyme counterproductively catalyzes the reverse reaction, reducing the net rate at which the reaction proceeds. Thus, the enzyme level required to achieve a given flux increases dramatically near equilibrium. Here, we develop a framework for quantifying the degree to which pathways suffer these thermodynamic limitations on flux. For each pathway, we calculate a single thermodynamically-derived metric (the Max-min Driving Force, MDF), which enables objective ranking of pathways by the degree to which their flux is constrained by low thermodynamic driving force. Our framework accounts for the effect of pH, ionic strength and metabolite concentration ranges and allows us to quantify how alterations to the pathway structure affect the pathway''s thermodynamics. Applying this methodology to pathways of central metabolism sheds light on some of their features, including metabolic bypasses (e.g., fermentation pathways bypassing substrate-level phosphorylation), substrate channeling (e.g., of oxaloacetate from malate dehydrogenase to citrate synthase), and use of alternative cofactors (e.g., quinone as an electron acceptor instead of NAD). The methods presented here place another arrow in metabolic engineers'' quiver, providing a simple means of evaluating the thermodynamic and kinetic quality of different pathway chemistries that produce the same molecules.     

Searching for Collective Behavior in a Large Network of Sensory Neurons

Maximum entropy models are the least structured probability distributions that exactly reproduce a chosen set of statistics measured in an interacting network. Here we use this principle to construct probabilistic models which describe the correlated spiking activity of populations of up to 120 neurons in the salamander retina as it responds to natural movies. Already in groups as small as 10 neurons, interactions between spikes can no longer be regarded as small perturbations in an otherwise independent system; for 40 or more neurons pairwise interactions need to be supplemented by a global interaction that controls the distribution of synchrony in the population. Here we show that such “K-pairwise” models—being systematic extensions of the previously used pairwise Ising models—provide an excellent account of the data. We explore the properties of the neural vocabulary by: 1) estimating its entropy, which constrains the population''s capacity to represent visual information; 2) classifying activity patterns into a small set of metastable collective modes; 3) showing that the neural codeword ensembles are extremely inhomogenous; 4) demonstrating that the state of individual neurons is highly predictable from the rest of the population, allowing the capacity for error correction.     

Identification of proteins that form specific complexes with the highly conserved protein Translin in Schizosaccharomyces pombe

Translin is a single-stranded DNA and RNA binding protein that has a high affinity for G-rich sequences. TRAX is a Translin paralog that associates with Translin. Both Translin and TRAX were highly conserved in eukaryotes. The nucleic acid binding form of Translin is a barrel-shaped homo-octamer. A Translin–TRAX hetero-octamer having a similar structure also binds nucleic acids. Previous reports suggested that Translin may be involved in chromosomal translocations, telomere metabolism and the control of mRNA transport and translation. More recent studies have indicated that Translin–TRAX hetero-octamers are involved in RNA silencing. To gain a further insight into the functions of Translin, we have undertaken to systematically search for proteins with which it forms specific complexes in living cells. Here we report the results of such a search conducted in the fission yeast Schizosaccharomyces pombe, a suitable model system. This search was carried out by affinity purification and immuno-precipitation techniques, combined with differential labeling of the intracellular proteins with the stable isotopes ¹⁵N and ¹⁴N. We identified for the first time two proteins containing an RNA Recognition Motif (RRM), which are specifically associated with the yeast Translin: (1) the pre-mRNA-splicing factor srp1 that belongs to the highly conserved SR family of proteins and (2) vip1, a protein conserved in fungi. Our data also support the presence of RNA in these intracellular complexes. Our experimental approach should be generally applicable to studies of weak intracellular protein–protein interactions and provides a clear distinction between false positive vs. truly interacting proteins.     

A case study of virilizing adrenal tumor in an adolescent female elite tennis player--insight into the use of anabolic steroids in young athletes

A 14-year-old Caucasian girl was referred to the endocrine clinic for evaluation of voice deepening, facial hirsutism, and acne starting 2 years previously. She had been a competitive tennis player since age 7 years, practicing for 4-6 hours daily. On physical examination she was noticed to have a masculine appearance with mild facial acne and moderate hirsutism. Tanner stage was 1 for breast tissue and 5 for pubic hair. Her androgen levels (testosterone, androstenedione, dehydroepiandrosterone sulfate) were extremely elevated. Adrenal ultrasonography revealed a round left 4.6 × 5.3-cm adrenal mass. Laparoscopic left adrenalectomy was performed. The histologic findings were compatible with a benign adrenocortical tumor. Postoperatively, androgen levels dropped to within the normal range. Breast development proceeded normally, menarche occurred 2 months after tumor resection, and menses has been regular since then. Muscle strength of the dominant and nondominant upper and lower extremities was measured 1 month before surgery and 1 year later, using an isokinetic dynamometer (Biodex Systems II, Biodex, Shirley, NY, USA). There was no significant decrease in overall muscle strength after removal of the virilizing tumor and the marked drop in circulating androgens. In addition, the patient maintained her age category, number 1, national tennis ranking. The results suggest that even extremely high levels of tumor-related circulating androgens had no evident effect on muscle strength and competitive performance in a female adolescent tennis player. The lack of beneficial effect on performance in adolescents, combined with the potentially hazardous side effects of anabolic steroids, suggests that teenage athletes should avoid their use.     

Colon stem cell and crypt dynamics exposed by cell lineage reconstruction

Stem cell dynamics in vivo are often being studied by lineage tracing methods. Our laboratory has previously developed a retrospective method for reconstructing cell lineage trees from somatic mutations accumulated in microsatellites. This method was applied here to explore different aspects of stem cell dynamics in the mouse colon without the use of stem cell markers. We first demonstrated the reliability of our method for the study of stem cells by confirming previously established facts, and then we addressed open questions. Our findings confirmed that colon crypts are monoclonal and that, throughout adulthood, the process of monoclonal conversion plays a major role in the maintenance of crypts. The absence of immortal strand mechanism in crypts stem cells was validated by the age-dependent accumulation of microsatellite mutations. In addition, we confirmed the positive correlation between physical and lineage proximity of crypts, by showing that the colon is separated into small domains that share a common ancestor. We gained new data demonstrating that colon epithelium is clustered separately from hematopoietic and other cell types, indicating that the colon is constituted of few progenitors and ruling out significant renewal of colonic epithelium from hematopoietic cells during adulthood. Overall, our study demonstrates the reliability of cell lineage reconstruction for the study of stem cell dynamics, and it further addresses open questions in colon stem cells. In addition, this method can be applied to study stem cell dynamics in other systems.     

Collective dynamics of gene expression in cell populations

The phenotypic state of the cell is commonly thought to be determined by the set of expressed genes. However, given the apparent complexity of genetic networks, it remains open what processes stabilize a particular phenotypic state. Moreover, it is not clear how unique is the mapping between the vector of expressed genes and the cell's phenotypic state. To gain insight on these issues, we study here the expression dynamics of metabolically essential genes in twin cell populations. We show that two yeast cell populations derived from a single steady-state mother population and exhibiting a similar growth phenotype in response to an environmental challenge, displayed diverse expression patterns of essential genes. The observed diversity in the mean expression between populations could not result from stochastic cell-to-cell variability, which would be averaged out in our large cell populations. Remarkably, within a population, sets of expressed genes exhibited coherent dynamics over many generations. Thus, the emerging gene expression patterns resulted from collective population dynamics. It suggests that in a wide range of biological contexts, gene expression reflects a self-organization process coupled to population-environment dynamics.     

An analytical approach differentiates between individual and collective cancer invasion

Tumour cells employ a variety of mechanisms to invade their environment and to form metastases. An important property is the ability of tumour cells to transition between individual cell invasive mode and collective mode. The switch from collective to individual cell invasion in the breast was shown recently to determine site of subsequent metastasis. Previous studies have suggested a range of invasion modes from single cells to large clusters. Here, we use a novel image analysis method to quantify and categorise invasion. We have developed a process using automated imaging for data collection, unsupervised morphological examination of breast cancer invasion using cognition network technology (CNT) to determine how many patterns of invasion can be reliably discriminated. We used Bayesian network analysis to probabilistically connect morphological variables and therefore determine that two categories of invasion are clearly distinct from one another. The Bayesian network separated individual and collective invading cell groups based on the morphological measurements, with the level of cell-cell contact the most discriminating morphological feature. Smaller invading groups were typified by smoother cellular surfaces than those invading collectively in larger groups. Interestingly, elongation was evident in all invading cell groups and was not a specific feature of single cell invasion as a surrogate of epithelial-mesenchymal transition. In conclusion, the combination of cognition network technology and Bayesian network analysis provides an insight into morphological variables associated with transition of cancer cells between invasion modes. We show that only two morphologically distinct modes of invasion exist.     

The Expression of the Beta Cell-Derived Autoimmune Ligand for the Killer Receptor Nkp46 Is Attenuated in Type 2 Diabetes

NK cells rapidly kill tumor cells, virus infected cells and even self cells. This is mediated via killer receptors, among which NKp46 (NCR1 in mice) is prominent. We have recently demonstrated that in type 1 diabetes (T1D) NK cells accumulate in the diseased pancreas and that they manifest a hyporesponsive phenotype. In addition, we found that NKp46 recognizes an unknown ligand expressed by beta cells derived from humans and mice and that blocking of NKp46 activity prevented diabetes development. Here we investigated the properties of the unknown NKp46 ligand. We show that the NKp46 ligand is mainly located in insulin granules and that it is constitutively secreted. Following glucose stimulation the NKp46 ligand translocates to the cell membrane and its secretion decreases. We further demonstrate by using several modalities that the unknown NKp46 ligand is not insulin. Finally, we studied the expression of the NKp46 ligand in type 2 diabetes (T2D) using 3 different in vivo models and 2 species; mice and gerbils. We demonstrate that the expression of the NKp46 ligand is decreased in all models of T2D studied, suggesting that NKp46 is not involved in T2D.