Multilevel modelling for inference of genetic regulatory networks

Time-course experiments with microarrays are often used to study dynamic biological systems and genetic regulatory networks (GRNs) that model how genes influence each other in cell-level development of organisms. The inference for GRNs provides important insights into the fundamental biological processes such as growth and is useful in disease diagnosis and genomic drug design. Due to the experimental design, multilevel data hierarchies are often present in time-course gene expression data. Most existing methods, however, ignore the dependency of the expression measurements over time and the correlation among gene expression profiles. Such independence assumptions violate regulatory interactions and can result in overlooking certain important subject effects and lead to spurious inference for regulatory networks or mechanisms. In this paper, a multilevel mixed-effects model is adopted to incorporate data hierarchies in the analysis of time-course data, where temporal and subject effects are both assumed to be random. The method starts with the clustering of genes by fitting the mixture model within the multilevel random-effects model framework using the expectation-maximization (EM) algorithm. The network of regulatory interactions is then determined by searching for regulatory control elements (activators and inhibitors) shared by the clusters of co-expressed genes, based on a time-lagged correlation coefficients measurement. The method is applied to two real time-course datasets from the budding yeast (Saccharomyces cerevisiae) genome. It is shown that the proposed method provides clusters of cell-cycle regulated genes that are supported by existing gene function annotations, and hence enables inference on regulatory interactions for the genetic network.

A novel requirement for C. elegans Alix/ALX-1 in RME-1-mediated membrane transport

BACKGROUND: Alix/Bro1p family proteins have recently been identified as important components of multivesicular endosomes (MVEs) and are involved in the sorting of endocytosed integral membrane proteins, interacting with components of the ESCRT complex, the unconventional phospholipid LBPA, and other known endocytosis regulators. During infection, Alix can be co-opted by enveloped retroviruses, including HIV, providing an important function during virus budding from the plasma membrane. In addition, Alix is associated with the actin cytoskeleton and might regulate cytoskeletal dynamics. RESULTS: Here we demonstrate a novel physical interaction between the only apparent Alix/Bro1p family protein in C. elegans, ALX-1, and a key regulator of receptor recycling from endosomes to the plasma membrane, called RME-1. The analysis of alx-1 mutants indicates that ALX-1 is required for the endocytic recycling of specific basolateral cargo in the C. elegans intestine, a pathway previously defined by the analysis of rme-1 mutants. The expression of truncated human Alix in HeLa cells disrupts the recycling of major histocompatibility complex class I, a known Ehd1/RME-1-dependent transport step, suggesting the phylogenetic conservation of this function. We show that the interaction of ALX-1 with RME-1 in C. elegans, mediated by RME-1/YPSL and ALX-1/NPF motifs, is required for this recycling process. In the C. elegans intestine, ALX-1 localizes to both recycling endosomes and MVEs, but the ALX-1/RME-1 interaction appears to be dispensable for ALX-1 function in MVEs and/or late endosomes. CONCLUSIONS: This work provides the first demonstration of a requirement for an Alix/Bro1p family member in the endocytic recycling pathway in association with the recycling regulator RME-1.

Specialist contact lens fitting

What is meant by the term ‘specialist contact lens fitting’? Or put another way, what would be considered non-specialist contact lens fitting? Is there such a thing as routine contact lens fitting? Soft or silicone hydrogel fitting for daily wear would probably be considered as routine contact lens fitting, but would extended or flexible wear remain in the same category or would they be considered a specialist fit? Different eras will classify different products as being ‘specialist’. Certainly twenty years ago soft toric contact lenses were considered as being speciality lenses but today would be thought of as routine lenses. Conversely, gas permeable lenses were thought of as mainstream twenty years ago but now are considered as speciality lenses. Although this would not be the same globally, as in some countries (such as Netherlands, France and Japan) gas permeable lens fitting remains popular and is not on the decline as in other countries (Canada, Australia and Sweden) [1]. Bandage soft lenses applied after surface laser refractive procedures would be considered as therapeutic lenses but in reality they are just plano thin hydrogel lenses worn constantly for 3–4 days to allow the underlying epithelium to convalesce and are then removed [2]. Some patients find that wearing hydrogel lenses during periods when they suffer from seasonal allergies actually improves their ocular comfort as the contact lens acts as a barrier to the allergen [3] and [4]. Scleral lenses have long been considered speciality lenses, apart from a time when they were the only lenses available but at that time all contact lens work would have been considered speciality practice! Nowadays we see the advent of mini-scleral designs and we see large diameter gas permeable lenses too. It is possible that these lenses increase the popularity of gas permeable lenses again and they become more main stream. So it would seem that the lines between routine and speciality contact lens fitting are not clear. Whether a lens is classed a specialist fit or not would depend on the lens type, why it was fitted, where in the world the fitting was being done and even the era in which it was fitted. This begs the question as to what would be considered entry level knowledge in contact lens fitting. This may not be an issue for most BCLA members or CLAE readers but certainly would be for bodies such as the College of Optometrists (UK) or the Association of British Dispensing Opticians when they are planning the final registration examinations for budding practitioners or when planning the level of higher level qualifications such as College Certificates or Diplomas. Similarly for training institutions when they are planning their course content. This becomes even trickier when trying to devise a qualification that spans across many countries, like the European Diploma in Optometry and Optics. How do we know if the training and examination level is correct? One way would be to analyse things when they go wrong and if patterns of malpractice are seen then maybe that could be used as an indicator to more training being needed. There were 162 Fitness to Practice Hearing at the General Optical Council between 2001 and 2010. Forty-seven of these were clinically related case, 39 fraud related, and 76 others. Of the clinical ones only 3 were contact lens related. So it would appear that as whole, in the profession, contact lens clinical skills are not being questioned too often (although it seems a few of us can’t keep our hands out the cookie jar!).

A laser induced breakdown spectroscopy quantitative analysis method based on the robust least squares support vector machine regression model

Data fluctuation in multiple measurements of Laser Induced Breakdown Spectroscopy (LIBS) greatly affects the accuracy of quantitative analysis. A new LIBS quantitative analysis method based on the Robust Least Squares Support Vector Machine (RLS-SVM) regression model is proposed. The usual way to enhance the analysis accuracy is to improve the quality and consistency of the emission signal, such as by averaging the spectral signals or spectrum standardization over a number of laser shots. The proposed method focuses more on how to enhance the robustness of the quantitative analysis regression model. The proposed RLS-SVM regression model originates from the Weighted Least Squares Support Vector Machine (WLS-SVM) but has an improved segmented weighting function and residual error calculation according to the statistical distribution of measured spectral data. Through the improved segmented weighting function, the information on the spectral data in the normal distribution will be retained in the regression model while the information on the outliers will be restrained or removed. Copper elemental concentration analysis experiments of 16 certified standard brass samples were carried out. The average value of relative standard deviation obtained from the RLS-SVM model was 3.06% and the root mean square error was 1.537%. The experimental results showed that the proposed method achieved better prediction accuracy and better modeling robustness compared with the quantitative analysis methods based on Partial Least Squares (PLS) regression, standard Support Vector Machine (SVM) and WLS-SVM. It was also demonstrated that the improved weighting function had better comprehensive performance in model robustness and convergence speed, compared with the four known weighting functions.

Characterization of microvesicles released from human red blood cells

Background/Aims: Extracellular vesicles (EVs) are spherical fragments of cell membrane released from various cell types under physiological as well as pathological conditions. Based on their size and origin, EVs are classified as exosome, microvesicles (MVs) and apoptotic bodies. Recently, the release of MVs from human red blood cells (RBCs) under different conditions has been reported. MVs are released by outward budding and fission of the plasma membrane. However, the outward budding process itself, the release of MVs and the physical properties of these MVs have not been well investigated. The aim of this study is to investigate the formation process, isolation and characterization of MVs released from RBCs under conditions of stimulating Ca2+ uptake and activation of protein kinase C. Methods: Experiments were performed based on single cell fluorescence imaging, fluorescence activated cell sorter/flow cytometer (FACS), scanning electron microscopy (SEM), atomic force microscopy (AFM) and dynamic light scattering (DLS). The released MVs were collected by differential centrifugation and characterized in both their size and zeta potential. Results: Treatment of RBCs with 4-bromo-A23187 (positive control), lysophosphatidic acid (LPA), or phorbol-12 myristate-13 acetate (PMA) in the presence of 2 mM extracellular Ca2+ led to an alteration of cell volume and cell morphology. In stimulated RBCs, exposure of phosphatidylserine (PS) and formation of MVs were observed by using annexin V-FITC. The shedding of MVs was also observed in the case of PMA treatment in the absence of Ca2+, especially under the transmitted bright field illumination. By using SEM, AFM and DLS the morphology and size of stimulated RBCs, MVs were characterized. The sizes of the two populations of MVs were 205.8 ± 51.4 nm and 125.6 ± 31.4 nm, respectively. Adhesion of stimulated RBCs and MVs was observed. The zeta potential of MVs was determined in the range from - 40 mV to - 10 mV depended on the solutions and buffers used. Conclusion: An increase of intracellular Ca2+ or an activation of protein kinase C leads to the formation and release of MVs in human RBCs.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Biometria de frutos e sementes, crescimento e produção de mudas de Campomanesia adamantium Blume e Campomanesia pubescens O. Berg

Campomanesia adamantium and Campomanesia pubescens are morphologically similar, occur in the same regions of the Cerrado, are difficult to differentiate and exhibit naturally slow growth and development. The objectives of this study were to analyze fruit and seed biometric data, emergence capacity and seedling growth characteristics. Fruit from both species was collected and used to measure biometric data of the fruit and seeds and to set up two emergence and two seedling growth experiments. C. adamantium fruit is round and has wider, pale yellow seeds while C. pubescens fruit is ellipsoidal to pyriform with yellow gold seeds. C. adamantium’s greater fresh fruit mass and biometric variability favors selection of promising material for commercialization. Seed drying reduced the speed and rate of seedling emergence of C. adamantium, but had no effect on C. pubescens. Leaf number, height, shoot dry weight and the ratio of shoot dry weight to root dry weight were greater in C. adamantium than in the slower growing C. pubescens. Increases in substrate volume favor seedling development. Slow-release fertilizer application at 1g per 115 cm3 substrate increased leaf formation, plant height and shoot dry weight of seedlings of both species and lateral budding in C. adamantium. Until 120 days after transplant, lateral budding was not observed in C. pubescens seedlings. For all traits evaluated in this experiment, C. adamantium seedling growth was greater than that of C. pubescens.

Mechanisms of Cdc42 Polarization in Yeast

Polarization is important for the function and morphology of many different cell types. The keys regulators of polarity in eukaryotes are the Rho-family GTPases. In the budding yeast Saccharomyces cerevisiae, which must polarize in order to bud and to mate, the master regulator is the highly conserved Rho GTPase, Cdc42. During polarity establishment, active Cdc42 accumulates at a site on the plasma membrane characterizing the “front” of the cell where the bud will emerge. The orientation of polarization is guided by upstream cues that dictate the site of Cdc42 clustering. However, in the absence of upstream cues, yeast can still polarize in a random direction during symmetry breaking. Symmetry breaking suggests cells possess an autocatalytic polarization mechanism that can amplify stochastic fluctuations of polarity proteins through a positive feedback mechanism.

Two different positive feedback mechanisms have been proposed to polarize Cdc42 in budding yeast. One model posits that Cdc42 activation must be localized to a site at the plasma membrane. Another model posits that Cdc42 delivery must be localized to a particular site at the plasma membrane. Although both mechanisms could work in parallel to polarize Cdc42, it is unclear which mechanism is critical to polarity establishment. We directly tested the predictions of the two positive feedback models using genetics and live microscopy. We found that localized Cdc42 activation is necessary for polarity establishment.

While this explains how active Cdc42 localizes to a particular site at the plasma membrane, it does not address how Cdc42 concentrates at that site. Several different mechanisms have been proposed to concentrate Cdc42. The GDI can extract Cdc42 from membranes and selective mobilize GDP-Cdc42 in the cytoplasm. It was proposed that selectively mobilizing GDP-Cdc42 in combination with local activation could locally concentrate total Cdc42 at the polarity site. Although the GDI is important for rapid Cdc42 accumulation at the polarity site, it is not essential to Cdc42 concentration. It was proposed that delivery of Cdc42 by actin-mediated vesicle can act as a backup pathway to concentrate Cdc42. However, we found no evidence for an actin-dependent concentrating pathway. Live microscopy experiments reveal that prenylated proteins are not restricted to membranes, and can enter the cytoplasm. We found that the GDI-independent concentrating pathway still requires Cdc42 to exchange between the plasma membrane and the cytoplasm, which is supported by computational modeling. In the absence of the GDI, we found that Cdc42 GAP became essential for polarization. We propose that the GAP limits GTP-Cdc42 leak into the cytoplasm, which would be prohibitive to Cdc42 polarization.

From Polarity to Morphogenesis PAK Behaviors and Mechanism for Bud Sensing in Morphogenesis Checkpoint

Bud formation by Saccharomyces cerevisiae is a fundamental process for yeast proliferation. Bud emergence is initiated by the polarization of the cytoskeleton, leading to local secretory vesicle delivery and gulcan synthase activity. The master regulator of polarity establishment is a small Rho-family GTPase – Cdc42. Cdc42 forms a clustered patch at the incipient budding site in late G1 and mediates downstream events which lead to bud emergence. Cdc42 promotes morphogenesis via its various effectors. PAKs (p21-activated kinases) are important Cdc42 effectors which mediate actin cytoskeleton polarization and septin filament assembly. The PAKs Cla4 and Ste20 share common binding domains for GTP-Cdc42 and they are partially redundant in function. However, we found that Cla4 and Ste20 behaved differently during the polarization and this depended on their different membrane interaction domains. Also, Cla4 and Ste20 compete for a limited number of binding sites at the polarity patch during bud emergence. These results suggest that PAKs may be differentially regulated during polarity establishment.

Morphogenesis of yeast must be coordinated with the nuclear cycle to enable successful proliferation. Many environmental stresses temporarily disrupt bud formation, and in such circumstances, the morphogenesis checkpoint halts nuclear division until bud formation can resume. Bud emergence is essential for degradation of the mitotic inhibitor, Swe1. Swe1 is localized to the septin cytoskeleton at the bud neck by the Swe1-binding protein Hsl7. Neck localization of Swe1 is required for Swe1 degradation. Although septins form a ring at the presumptive bud site prior to bud emergence, Hsl7 is not recruited to the septins until after bud emergence, suggesting that septins and/or Hsl7 respond to a “bud sensor”. Here we show that recruitment of Hsl7 to the septin ring depends on a combination of two septin-binding kinases: Hsl1 and Elm1. We elucidate which domains of these kinases are needed, and show that artificial targeting of those domains suffices to recruit Hsl7 to septin rings even in unbudded cells. Moreover, recruitment of Elm1 is responsive to bud emergence. Our findings suggest that Elm1 plays a key role in sensing bud emergence.

Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi

© Medina et al.Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.

Genetic Analysis of Mitotic Recombination in Saccharomyces cerevisiae

Mitotic genome instability can occur during the repair of double-strand breaks (DSBs) in DNA, which arise from endogenous and exogenous sources. Studying the mechanisms of DNA repair in the budding yeast, Saccharomyces cerevisiae has shown that Homologous Recombination (HR) is a vital repair mechanism for DSBs. HR can result in a crossover event, in which the broken molecule reciprocally exchanges information with a homologous repair template. The current model of double-strand break repair (DSBR) also allows for a tract of information to non-reciprocally transfer from the template molecule to the broken molecule. These “gene conversion” events can vary in size and can occur in conjunction with a crossover event or in isolation. The frequency and size of gene conversions in isolation and gene conversions associated with crossing over has been a source of debate due to the variation in systems used to detect gene conversions and the context in which the gene conversions are measured.

In Chapter 2, I use an unbiased system that measures the frequency and size of gene conversion events, as well as the association of gene conversion events with crossing over between homologs in diploid yeast. We show mitotic gene conversions occur at a rate of 1.3x10-6 per cell division, are either large (median 54.0kb) or small (median 6.4kb), and are associated with crossing over 43% of the time.

DSBs can arise from endogenous cellular processes such as replication and transcription. Two important RNA/DNA hybrids are involved in replication and transcription: R-loops, which form when an RNA transcript base pairs with the DNA template and displaces the non-template DNA strand, and ribonucleotides embedded into DNA (rNMPs), which arise when replicative polymerase errors insert ribonucleotide instead of deoxyribonucleotide triphosphates. RNaseH1 (encoded by RNH1) and RNaseH2 (whose catalytic subunit is encoded by RNH201) both recognize and degrade the RNA in within R-loops while RNaseH2 alone recognizes, nicks, and initiates removal of rNMPs embedded into DNA. Due to their redundant abilities to act on RNA:DNA hybrids, aberrant removal of rNMPs from DNA has been thought to lead to genome instability in an rnh201Δ background.

In Chapter 3, I characterize (1) non-selective genome-wide homologous recombination events and (2) crossing over on chromosome IV in mutants defective in RNaseH1, RNaseH2, or RNaseH1 and RNaseH2. Using a mutant DNA polymerase that incorporates 4-fold fewer rNMPs than wild type, I demonstrate that the primary recombinogenic lesion in the RNaseH2-defective genome is not rNMPs, but rather R-loops. This work suggests different in-vivo roles for RNaseH1 and RNaseH2 in resolving R-loops in yeast and is consistent with R-loops, not rNMPs, being the the likely source of pathology in Aicardi-Goutières Syndrome patients defective in RNaseH2.

Using Coding to Improve Localization and Backscatter Communication Performance in Low-Power Sensor Networks

Backscatter communication is an emerging wireless technology that recently has gained an increase in attention from both academic and industry circles. The key innovation of the technology is the ability of ultra-low power devices to utilize nearby existing radio signals to communicate. As there is no need to generate their own energetic radio signal, the devices can benefit from a simple design, are very inexpensive and are extremely energy efficient compared with traditional wireless communication. These benefits have made backscatter communication a desirable candidate for distributed wireless sensor network applications with energy constraints.

The backscatter channel presents a unique set of challenges. Unlike a conventional one-way communication (in which the information source is also the energy source), the backscatter channel experiences strong self-interference and spread Doppler clutter that mask the information-bearing (modulated) signal scattered from the device. Both of these sources of interference arise from the scattering of the transmitted signal off of objects, both stationary and moving, in the environment. Additionally, the measurement of the location of the backscatter device is negatively affected by both the clutter and the modulation of the signal return.

This work proposes a channel coding framework for the backscatter channel consisting of a bi-static transmitter/receiver pair and a quasi-cooperative transponder. It proposes to use run-length limited coding to mitigate the background self-interference and spread-Doppler clutter with only a small decrease in communication rate. The proposed method applies to both binary phase-shift keying (BPSK) and quadrature-amplitude modulation (QAM) scheme and provides an increase in rate by up to a factor of two compared with previous methods.

Additionally, this work analyzes the use of frequency modulation and bi-phase waveform coding for the transmitted (interrogating) waveform for high precision range estimation of the transponder location. Compared to previous methods, optimal lower range sidelobes are achieved. Moreover, since both the transmitted (interrogating) waveform coding and transponder communication coding result in instantaneous phase modulation of the signal, cross-interference between localization and communication tasks exists. Phase discriminating algorithm is proposed to make it possible to separate the waveform coding from the communication coding, upon reception, and achieve localization with increased signal energy by up to 3 dB compared with previous reported results.

The joint communication-localization framework also enables a low-complexity receiver design because the same radio is used both for localization and communication.

Simulations comparing the performance of different codes corroborate the theoretical results and offer possible trade-off between information rate and clutter mitigation as well as a trade-off between choice of waveform-channel coding pairs. Experimental results from a brass-board microwave system in an indoor environment are also presented and discussed.

(Table 1) Magnetic susceptibility, lithologic composition and stable isotope record for ODP Hole 132-810C

Site 810 was drilled atop Shatsky Rise during Ocean Drilling Program (ODP) Leg 132. The principal objective at Site 810 was to drill interbedded cherts and chalks of Mesozoic age using the diamond coring system (DCS). The objective was not achieved because of difficulties in setting up the reentry cone on the seafloor; however, a shortened section of Cretaceous-Cenozoic nannofossil ooze was recovered with the advanced piston corer (APC). Although the section is interrupted by hiatuses, the upper 50 m carry detailed information relating to biogenic productivity, water chemistry, and eolian input during the Pliocene and Pleistocene. Four holes were drilled at Site 810. Hole 810A consists of a single mud-line core for an ongoing ODP geriatric study. The second hole (Hole 810B) was washed to 60 mbsf (without core recovery) to provide information required for setting the 16-in. casing attached to the reentry cone. Hole 810C penetrated 136.1 mbsf, mostly with the APC, with a total recovery of 143.81 m of nannofossil ooze. A reentry cone was placed over Hole 810D but no casing was successfully suspended in the hole and no sediment was cored. This data report presents the results of shore-based high-resolution analyses of carbonate and oxygen isotopic variations in the upper 50 m of the section at Site 810 and compares these variations with the shipboard determinations of magnetic susceptibility and GRAPE bulk density from the multisensor track.

Effects of pH on asexual reproduction and statolith formation of the scyphozoan, Aurelia labiata, 2010

Although anthropogenic infuences such as global warming, overfishing, and eutrophication may contribute to jellyfish blooms, little is known about the effects of ocean acidification on jellyfish. Most medusae form statoliths of calcium sulfate hemihydrate that are components of their balance organs (statocysts). This study was designed to test the effects of pH (7.9, within the average current range, 7.5, expected by 2100, and 7.2, expected by 2300) combined with two temperatures (9 and 15°C) on asexual reproduction and statolith formation of the moon jellyfish, Aurelia labiata. Polyp survival was 100% after 122 d in seawater in all six temperature and pH combinations. Because few polyps at 9°C strobilated, and temperature effects on budding were consistent with published results, we did not analyze data from those three treatments further. At 15°C, there were no significant effects of pH on the numbers of ephyrae or buds produced per polyp or on the numbers of statoliths per statocyst; however, statolith size was signi?cantly smaller in ephyrae released from polyps reared at low pH. Our results indicate that A. labiata polyps are quite tolerant of low pH, surviving and reproducing asexually even at the lowest tested pH; however, the effects of small statoliths on ephyra fitness are unknown. Future research on the behavior of ephyrae with small statoliths would further our understanding of how ocean acidi?cation may affect jellyfish survival in nature.

Polymer Process Partitioning: Extruding Plastic Pipe

When plastic pipe is solidified, it proceeds through a long cooling chamber. Inside this chamber, inside the hollow extrudate, the plastic is molten, and this inner surface solidifies last. Sag, the flow due to the self-weight of the molten plastic, then happens in this cooling chamber, and sometimes, thickened regions (called knuckles) arise in the lower quadrants, especially of large diameter thickwalled pipes. To compensate for sag, engineers normally shift the die centerpiece downward. This thesis focuses on the consequences of this decentering. Specifically, when the molten polymer is viscoelastic, as is normally the case, a downward lateral force is exerted on the mandrel. Die eccentricity also affects the downstream axial force on the mandrel. These forces govern how rigidly the mandrel must be attached (normally, on a spider die). We attack this flow problem in eccentric cylindrical coordinates, using the Oldroyd 8-constant constitutive model framework. Specifically, we revise the method of Jones (1964), called polymer process partitioning. We estimate both axial and lateral forces. We develop a corresponding map to help plastics engineers predict the extrudate shape, including extrudate knuckles. From the mass balance over the postdie region, we then predict the shape of the extrudate entering the cooling chamber. We further include expressions for the stresses in the extruded polymer melt. We include detailed dimensional worked examples to show process engineers how to use our results to design pipe dies, and especially to suppress extrudate knuckling.