Evolution of the leukotoxin promoter in genus Mannheimia

BACKGROUND: The Mannheimia species encompass a wide variety of bacterial lifestyles, including opportunistic pathogens and commensals of the ruminant respiratory tract, commensals of the ovine rumen, and pathogens of the ruminant integument. Here we present a scenario for the evolution of the leukotoxin promoter among representatives of the five species within genus Mannheimia. We also consider how the evolution of the leukotoxin operon fits with the evolution and maintenance of virulence. RESULTS: The alignment of the intergenic regions upstream of the leukotoxin genes showed significant sequence and positional conservation over a 225-bp stretch immediately proximal to the transcriptional start site of the lktC gene among all Mannheimia strains. However, in the course of the Mannheimia genome evolution, the acquisition of individual noncoding regions upstream of the conserved promoter region has occurred. The rate of evolution estimated branch by branch suggests that the conserved promoter may be affected to different extents by the types of natural selection that potentially operate in regulatory regions. Tandem repeats upstream of the core promoter were confined to M. haemolytica with a strong association between the sequence of the repeat units, the number of repeat units per promoter, and the phylogenetic history of this species. CONCLUSION: The mode of evolution of the intergenic regions upstream of the leukotoxin genes appears to be highly dependent on the lifestyle of the bacterium. Transition from avirulence to virulence has occurred at least once in M. haemolytica with some evolutionary success of bovine serotype A1/A6 strains. Our analysis suggests that changes in cis-regulatory systems have contributed to the derived virulence phenotype by allowing phase-variable expression of the leukotoxin protein. We propose models for how phase shifting and the associated virulence could facilitate transmission to the nasopharynx of new hosts.

Molecular genetic analysis of Dichelobacter nodosus proteases AprV2/B2, AprV5/B5 and BprV/B in clinical material from European sheep flocks

Dichelobacter nodosus, the etiological agent of ovine footrot, exists both as virulent and as benign strains, which differ in virulence mainly due to subtle differences in the three subtilisin-like proteases AprV2, AprV5 and BprV found in virulent, and AprB2, AprB5 and BprB in benign strains of D. nodosus. Our objective was a molecular genetic epidemiological analysis of the genes of these proteases by direct sequence analysis from clinical material of sheep from herds with and without history of footrot from 4 different European countries. The data reveal the two proteases known as virulent AprV2 and benign AprB2 to correlate fully to the clinical status of the individuals or the footrot history of the herd. In samples taken from affected herds, the aprV2 gene was found as a single allele whereas in samples from unaffected herds several alleles with minor modifications of the aprB2 gene were detected. The different alleles of aprB2 were related to the herds. The aprV5 and aprB5 genes were found in the form of several alleles scattered without distinction between affected and non-affected herds. However, all different alleles of aprV5 and aprB5 encode the same amino acid sequences, indicating the existence of a single protease isoenzyme 5 in both benign and virulent strains. The genes of the basic proteases BprV and BprB also exist as various alleles. However, differences found in samples from affected versus non-affected herds do not reflect the currently known epitopes that are attributed to differences in biochemical activity. The data of the study confirm the prominent role of AprV2 in the virulence of D. nodosus and shed a new light on the presence of the other protease genes and their allelic variants in clinical samples.

Frequent sensitization to Candida albicans and profilins in adult eosinophilic esophagitis

BACKGROUND: Eosinophilic esophagitis (EoE) is often associated with atopic airway and skin diseases. More than 80% of EoE patients are sensitized to aero- and/or food allergens. Immunoglobulin (Ig)E-mediated immune responses to microbes have been reported to be deleterious in connection with atopic diseases. AIM: The aim of this study was to obtain a comprehensive overview about the sensitization spectrum of adult EoE patients. METHODS: IgE in sera of 35 patients with active EoE were analyzed for reactivity to Candida albicans, as well as to a panel of recombinant and purified natural allergen components, using a microarray. RESULTS: IgE sensitization to Candida albicans was found in 43% of EoE patients. More than 80% of EoE patients were sensitized to aeroallergens and 22% to food-specific allergen components, whereas 69% of the patients exhibited specific IgE to cross-reactive allergens. Among the latter, profilins were identified as most frequent IgE cross-reactive allergen components. Interestingly, dysphagia, the main symptom of adult EoE patients following rice and/or bread ingestion, was associated with sensitization to cross-reactive allergens such as profilins, pathogenesis-related (PR) 10 and lipid transfer proteins (LTP). Intolerance toward meat rarely correlated with sensitization to animal food allergens. CONCLUSION: Candida albicans and cross-reactive plant allergen components, in particular profilins, were identified as frequent sensitizers in adult EoE patients. Specific elimination therapies are suggested to reveal their actual roles in the pathogenesis of EoE.

Analysis of Variation in Clubfoot Candidate Genes

Isolated clubfoot, a common birth defect occurring in more than 135,000 livebirths worldwide each year, is associated with significant health care and financial burdens. Clubfoot is defined by forefoot adduction, hindfoot varus, midfoot cavus and hindfoot equinus. Isolated clubfoot, which is the focus of these studies, is distinct from syndromic clubfoot because there are no other associated malformations. Population, family, twin and segregation analysis studies provide evidence that genetic and environmental factors play an etiologic role in isolated clubfoot. The studies described in this thesis were performed to define the role of genetic variation in isolated clubfoot. Interrogation of a deletion region associated with syndromic clubfoot, suggested that CASP8 and CASP10, two apoptotic genes, play a role in isolated clubfoot. To explore the role of apoptotic genes in clubfoot, SNPs spanning genes involved in the apoptotic pathway in the six chromosomal deletion regions, and limb patterning genes, HOXD and HOXA, were interrogated. SNPs in mitochondrial mediated apoptotic genes and several SNPs in HOXA and HOXD genes were modestly associated with clubfoot with the most significant SNP, rs3801776, located in the basal promoter of HOXA9. Several significant associations were found with SNPs in NFAT2 and TNIP2. Significant gene interactions were detected between SNPs in HOX and apoptotic genes. These findings suggest a model for clubfoot in which variation in one gene is not sufficient to cause the malformation but requires variation several genes to perturb protein expression sufficiently to alter muscle and foot development. These results significantly impact our knowledge base by delineating underlying mechanisms causing clubfoot.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Large scale variation in Enterococcus faecalis illustrated by the genome analysis of strain OG1RF.

BACKGROUND: Enterococcus faecalis has emerged as a major hospital pathogen. To explore its diversity, we sequenced E. faecalis strain OG1RF, which is commonly used for molecular manipulation and virulence studies. RESULTS: The 2,739,625 base pair chromosome of OG1RF was found to contain approximately 232 kilobases unique to this strain compared to V583, the only publicly available sequenced strain. Almost no mobile genetic elements were found in OG1RF. The 64 areas of divergence were classified into three categories. First, OG1RF carries 39 unique regions, including 2 CRISPR loci and a new WxL locus. Second, we found nine replacements where a sequence specific to V583 was substituted by a sequence specific to OG1RF. For example, the iol operon of OG1RF replaces a possible prophage and the vanB transposon in V583. Finally, we found 16 regions that were present in V583 but missing from OG1RF, including the proposed pathogenicity island, several probable prophages, and the cpsCDEFGHIJK capsular polysaccharide operon. OG1RF was more rapidly but less frequently lethal than V583 in the mouse peritonitis model and considerably outcompeted V583 in a murine model of urinary tract infections. CONCLUSION: E. faecalis OG1RF carries a number of unique loci compared to V583, but the almost complete lack of mobile genetic elements demonstrates that this is not a defining feature of the species. Additionally, OG1RF's effects in experimental models suggest that mediators of virulence may be diverse between different E. faecalis strains and that virulence is not dependent on the presence of mobile genetic elements.

A molecular analysis of the conditional defect in MuSVts 110 viral RNA splicing

Cells infected with MuSVts110 express a viral RNA which contains an inherent conditional defect in RNA splicing. It has been shown previously that splicing of the MuSVts110 primary transcript is essential to morphological transformation of 6m2 cells in vitro. A growth temperature of 33$\sp\circ$C is permissive for viral RNA splicing,and, consequently, 6m2 cells appear morphologically transformed at this temperature. However, 6m2 cells appear phenotypically normal when incubated at 39$\sp\circ$C, the non-permissive temperature for viral RNA splicing.^ After a shift from 39$\sp\circ$C to 33$\sp\circ$C, the coordinate splicing of previously synthesized and newly transcribed MuSVts110 RNA was achieved. By S1 nuclease analysis of total RNA isolated at various times, 5$\sp\prime$ splice site cleavage of the MuSVts110 transcript appeared to occur 60 minutes after the shift to 33$\sp\circ$C, and 30 minutes prior to detectable exon ligation. In addition, consistent with the permissive temperatures and the kinetic timeframe of viral RNA splicing after a shift to 33$\sp\circ$C, four temperature sensitive blockades to primer extension were identified 26-75 bases upstream of the 3$\sp\prime$ splice site. These blockades likely reflect four branchpoint sequences utilized in the formation of MuSVts110 lariat splicing-intermediates.^ The 54-5A4 cell line is a spontaneous revertant of 6m2 cells and appears transformed at all growth temperatures. Primer extension sequence analysis has shown that a five base deletion occurred at the 3$\sp\prime$ splice site in MuSVts110 RNA allowing the expression of a viral transforming protein in 54-5A4 in the absence of RNA splicing, whereas in the parental 6m2 cell line, a splicing event is necessary to generate a similar transforming protein. As a consequence of this deletion, splicing cannot occur and the formation of the four MuSVts110 branched-intermediates were not observed at any temperature in 54-5A4 cells. However, 5$\sp\prime$ splice site cleavage was still detected at 33$\sp\circ$C.^ Finally, we have investigated the role of the 1488 bp deletion which occurred in the generation of MuSVts110 in the activation of temperature sensitive viral RNA splicing. This deletion appears solely responsible for splice site activation. Whether intron size is the crucial factor in MuSVts110 RNA splicing or whether inhibitory sequences were removed by the deletion is currently unknown. (Abstract shortened with permission of author.) ^

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

Transcriptional analysis of liver-specific expression and differential interleukin-6 response of rat kininogen genes

Analyses of rat T1 kininogen gene/chloramphenicol acetyltransferase (T1K/CAT) constructs revealed two regions important for tissue-specific and induced regulation of T1 kininogen.^ Although the T1 kininogen gene is inducible by inflammatory cytokines, a highly homologous K kininogen gene is minimally responsive. Moreover, the basal expression of a KK/CAT construct was 5- to 7-fold higher than that of the analogous T1K/CAT construct. To examine the molecular basis of this differential regulation, a series of promoter swapping experiments was carried out. Our transfection results showed that at least two regions in the K kininogen gene are important for its high basal expression: a distal 19-bp region (C box) constituted a binding site for CCAAT/enhancer binding protein (C/EBP) family proteins and a proximal 66-bp region contained two adjacent binding sites for hepatocyte nuclear factor-3 (HNF-3). The distal HNF-3 binding site from the K kininogen promoter demonstrated a stronger affinity than that from the T1 kininogen promoter. Since C/EBP and HNF-3 are highly enriched in the liver and known to enhance transcription of liver-specific genes, differential binding affinities of these factors accounted for the higher basal expression of the K kininogen gene.^ In contrast to the K kininogen C box, the T1 kininogen C box does not bind C/EBP presumably due to their two-nucleotide divergence. This sequence divergence, however, converts it to a consensus binding sequence for two IL-6-inducible transcription factors--IL-6 response element binding protein and acute-phase response factor. To functionally determine whether C box sequences are important for their differential acute-phase response, T1 and K kininogen C boxes were swapped and analyzed after transfection into Hep3B cells. Our results showed that the T1 kininogen C box is indeed one of the IL-6 response elements in T1 kininogen promoter. Furthermore, its function can be modulated by a 5$\sp\prime$-adjacent C/EBP-binding site (B box) whose mutation significantly reduced the overall induced activity. Moreover, this B box is the target site for binding and transactivation of another IL-6 inducible transcription factor C/EBP$\delta.$ Evolutionary divergence of a few critical nucleotides can either lead to subtle changes in the binding affinities of a given transcription factor or convert a binding sequence for a constitutive factor to a site recognized by an inducible factor. (Abstract shortened by UMI.) ^

Analysis of galectin expression and identification of endogenous ligands in a human colon carcinoma cell line

The 14.5 kDa (galectin-1) and 31 kDa (galectin-3) lectins are the most well characterized members of a family of vertebrate carbohydrate-binding proteins known as the galectins. Evidence has been obtained implicating these galectins in events as diverse as cell-cell and cell-extracellular matrix interactions, growth regulation, transformation, differentiation, and programmed cell death. In the present study, sodium butyrate was found to be a potent inducer of galectin-1 in the KM12 human colon carcinoma cell line. Prior to treatment with butyrate this cell line expresses only galectin-3. These cells were utilized as an in vitro model system to study galectin expression as well as that of their endogenous ligands. The initial phase of this project involved the examination of the induction of galectin-1 by butyrate at the protein level. These studies indicated that galectin-1 induction by butyrate was relatively rapid reaching nearly maximal levels after only 24 hours. Additionally, the induction was found to be reversible upon the removal of butyrate and to precede the increase in expression of the well characterized differentiation marker, carcinoembryonic antigen (CEA). The second phase of this project involved the characterization of potential glycoprotein ligands for galectin-1 and galectin-3. This work demonstrated that the polylactosaminoglycan-containing glycoproteins laminin, CEA, and the lysosome-associated glycoproteins-1 and -2 (LAMPs-1 and -2) are capable of serving as ligands for both galectin-1 and -3. The third phase of this project involved the analysis of the induction of the galectin-1 promoter by butyrate. Through the analysis of deletion constructs transiently transfected into KM12 cells, the region of the galectin-1 promoter mediating a high level of induction by butyrate was localized primarily within a proximal portion of the promoter containing a CCAAT element and an Sp1 binding site. The CCAAT-binding activity in the KM12 nuclear extracts was subsequently dentified as NF-Y by gel shift analysis. These studies suggest that: (1) the galectins may be involved in modulating adhesive interactions in human colon carcinoma cells through the binding of several polylactosaminoglycans shown to play a role in adhesion and (2) high level induction of the galectin-1 promoter by butyrate can proceed through a discreet, proximal element containing an NF-Y-binding CCAAT box and an Sp1 site. ^

Analysis of differential gene expression in nontumorigenic glioblastoma multiforme

Loss of chromosome 10 represents the most common cytogenetic abnormality in high grade gliomas (glioblastoma multiforme). To identify genes involved in the malignant progression of human gliomas, a subtractive hybridization was performed between a tumorigenic glioblastoma cell line (LG11) and a nontumorgenic hybrid cell (LG11.3) containing an introduced chromosome 10. LG11 mRNA was subtracted from LG11.3 cDNA to produce cDNA probes enriched for sequences whose expression differs quantitatively from the parental tumorigenic cells. Both known and novel sequences were identified as a result of the subtraction. Northern blot analysis was then used to confirm differential expression of several subtracted clones. One novel clone, clone 17, identified a 2.6 kb message that showed a consistent two to four fold increase in expression in the LG11.3 nontumorigenic cells. Clone 17 (340 bp) was used successfully to screen for a near full-length version, RIG (regulated in glioma), which was 2,569 bp in size. The RIG cDNA sequence showed homology to clone 17 and to an anonymous EST (IB666), but to no previously identified genes. This screening effort also identified several independent clones representing novel sequences, most of which failed to show increased expression in the nontumorigenic GBM cells. Tissue distribution studies of RIG indicated highest levels of expression in human brain with appreciably lower levels in heart and lung. In vitro transcription and translation experiments demonstrated the ability of RIG to direct the synthesis of a 13 kD protein product. However, open reading frame analysis revealed no identify with previously described motifs or any known proteins. Using a combination of somatic cell hybrid panels and in situ hybridization, the RIG gene was mapped to chromosome 11p14-11p15. Further study of RIG and related gene products may provide insight into the negative regulation of glial oncogenesis. ^

Characterization of the human tissue transglutaminase gene promoter

Transglutaminases are a family of calcium-dependent enzymes, that catalyze the covalent cross-linking of proteins by forming $\varepsilon(\gamma$-glutamyl)lysine isopeptide bonds. In order to investigate the molecular mechanisms regulating the expression of the tissue transglutaminase gene and to determine its biological functions, the goal of this research has been to clone and characterize the human tissue transglutaminase promoter. Thirteen clones of the tissue transglutaminase gene were obtained from the screening of a human placental genomic DNA library. A 1.74 Kb fragment derived from DNA located immediately upstream of the translation start site was subcloned and sequenced. Sequence analysis of this DNA fragment revealed that it contains a TATA box (TATAA), a CAAT box (GGACAAT), and a series of potential transcription factor binding sites and hormone response elements. Four regions of significant homology, a GC-rich region, a TG-rich region, an AG-rich region, and HR1, were identified by aligning 1.8 Kb of DNA flanking the human, mouse, and guinea pig tissue transglutaminase genes.^ To measure promoter activity, we subcloned the 1.74 Kb fragment of the tissue transglutaminase gene into a luciferase reporter vector to generate transglutaminase promoter/luciferase reporter constructs. Transfection experiments showed that this DNA segment includes a functional promoter with high constitutive activity. Deletion analysis revealed that the SP1 sites or corresponding sequences contribute to this activity. We investigated the role of DNA methylation in regulating the activity of the promoter and found that in vitro methylation of tissue transglutaminase promoter/luciferase reporter constructs suppressed their basal activity. Methylation of the promoter is inversely correlated with the expression of the tissue transglutaminase gene in vivo. These results suggest that DNA methylation may be one of the mechanisms regulating the expression of the gene. The tumor suppressor gene product p53 was also shown to inhibit the activity of the promoter, suggesting that induction of the tissue transglutaminase gene is not involved in the p53-dependent programmed cell death pathway. Although retinoids regulate the expression of the tissue transglutaminase gene in vivo, retinoid-inducible activity can not be identified in 3.7 Kb of DNA 5$\sp\prime$ to the tissue transglutaminase gene.^ The structure of the 5$\sp\prime$ end of the tissue transglutaminase gene was mapped. Alignment analysis of the human tissue transglutaminase gene with other human transglutaminases showed that tissue transglutaminase is the simplest member of transglutaminase superfamily. Transglutaminase genes show a conserved core of exons and introns but diverse N-terminuses and promoters. These observations suggest that key regulatory sequences and promoter elements have been appended upstream of the core transglutaminase gene to generate the diversity of regulated expression and regulated activity characteristic of the transglutaminase gene family. ^

Functional analysis of histones H3 and H4 in Tup1 repression and in GCN5 activation

Tup1 forms a complex with Ssn6 in yeast. Ssn6-Tup1 complex is recruited via direct interactions with specific DNA binding proteins to a specific promoter region and mediates repression of several sets of genes including a-cell specific genes (asg) in $\alpha$ cells. It has been shown that repression of asgs also requires histone H4 and that Tup1 can directly interact with H3 and H4 in vitro. To address whether histone H3 is required for the repression of asgs, I have examined the effect of H3 and H4 mutations on the expression of a $\alpha$2-controlled LacZ reporter. Assay of $\beta$-glactosidase shows that mutations in either H3 or H4 cause a weak derepression of the reporter gene. Some double mutations result in a stronger derepression, while others do not. The H3 N-terminal deletion also leads to a slightly decreased expression of the reporter gene in $\alpha$ cells. Our data suggest that the N-termini of both H3 and H4 are cooperatively involved in the repression of a-cell specific genes in $\alpha$ cells, possibly through their interaction with Tup1.^ GCN5 was originally identified as a transcriptional regulator required to activate a subset of genes in yeast. Recently, it has been shown that GCN5 encodes the catalytic subunit of a nuclear histone acetyltransferase, providing the first direct link between histone acetylation and gene transcription. Recombinant Gcn5p (rGcn5p) exhibits a limited substrate specificity in vitro. However, neither the specificity of this enzyme in vivo nor the importance of particular acetylated residues to transcription or cell growth are well defined. In order to define the sites of histone acetylation mediated by Gcn5p in vivo and assess the significance of histone acetylation, more than 30 yeast strains have been constructed to bear specific H3 and/or H4 mutations in the presence or absence of GCN5 function. Our genetic data suggest that Gcn5p may have additional targets in vivo that are not identified as the targets of rGcn5p by previous studies. Western analysis using antibodies specifically recognizing particular acetylated isoforms of H3 and H4 led us to conclude that Gcn5p is necessary for full acetylation of multiple sites in both H3 and H4 in vivo. Consistent with these observations, rGcn5p still acetylates histones H3 and H4 bearing mutations either in H3 K14 or H4 K8,16, sites previously identified as the targets of acetylation by rGcn5p in H3 and H4. Our data also demonstrated that Gcn5p-mediated acetylation events are important for normal progression of the cell cycle and for transcriptional activation. Furthermore, a critical overall level of acetylation is essential for cell viability. ^

Play in Rats: Association across Contexts and Types , and Analysis of Structure.

Play has been proposed as a promising indicator of positive animal welfare. We aimed to study play in rats across contexts (conspecific/heterospecific) and types (social: pinning, being pinned; solitary: scampering), and we investigated its structure using behavioral sequence analysis. Group-housed (three per cage) adolescent male Lister Hooded rats (n = 21) were subjected to a Play-In-Pairs test: after a 3 hour isolation period, a pair of cage-mates was returned to the home cage and both social and solitary play were scored for 20 min. This procedure was repeated for each pair combination across three consecutive days, and individual play scores were calculated. Heterospecific play was measured using a Tickling test: rats were individually tickled by the experimenter through bouts of gentle, rapid finger movements on their underside, and the number of positive 50 kHz frequency modulated vocalizations and experimenter-directed approach behaviors were recorded. Both of the above tests were compared with social play in the home cage. While conspecific play in both the Play-In-Pairs test and home cage were correlated, both seemed to be unrelated to heterospecific play in the Tickling test. During the Play-In-Pairs test, although both solitary and social play types occurred, they were unrelated, and solitary locomotor play of one rat did not predict the subsequent play behavior of its cage mate. Analysis of play structure revealed that social play occurred more often in bouts of repeated behaviors while solitary play sequences did not follow a specific pattern. If play is to be used as an indicator of positive welfare in rats, context, type and structure differences should be taken into account.

Assessing the Potential of Multiple Imputation in Sequence Analysis

Sequence analysis and optimal matching are useful heuristic tools for the descriptive analysis of heterogeneous individual pathways such as educational careers, job sequences or patterns of family formation. However, to date it remains unclear how to handle the inevitable problems caused by missing values with regard to such analysis. Multiple Imputation (MI) offers a possible solution for this problem but it has not been tested in the context of sequence analysis. Against this background, we contribute to the literature by assessing the potential of MI in the context of sequence analyses using an empirical example. Methodologically, we draw upon the work of Brendan Halpin and extend it to additional types of missing value patterns. Our empirical case is a sequence analysis of panel data with substantial attrition that examines the typical patterns and the persistence of sex segregation in school-to-work transitions in Switzerland. The preliminary results indicate that MI is a valuable methodology for handling missing values due to panel mortality in the context of sequence analysis. MI is especially useful in facilitating a sound interpretation of the resulting sequence types.