Infrared nanoimaging of neuronal ultrastructure and nanoparticle interaction with cells

Here we introduce scattering-type scanning near-field optical microscopy (s-SNOM) as a novel tool for nanoscale chemical-imaging of sub-cellular organelles, nanomaterials and of the interactions between them. Our setup uses a tuneable mid-infrared laser and a sharp scanning probe to image at a resolution substantially surpassing the diffraction limit. The laser can be tuned to excite vibrational modes of functional groups in biomolecules, (e.g. amide moieties), in a way that enables direct chemical mapping without the need for labelling. We, for the first time, chemically image neuronal ultrastructure, identify neuronal organelles and sub-organelle structures as small as 10 nm and validate our findings using transmission electron microscopy (TEM). We produce chemical and morphological maps of neurons treated with gold nanospheres and characterize nanoparticle size and intracellular location, and their interaction with the plasma membrane. Our results show that the label-free nature of s-SNOM means it has a 'true' chemical resolution of up to 20 nm which can be further improved. We argue that it offers significant potential in nanomedicine for nanoscale chemical imaging of cell ultrastructure and the subcellular distribution of nanomaterials within tissues.

Autocrine IL-6 drives cell and extracellular matrix anisotropy in scar fibroblasts

Fibrosis is associated with dramatic changes in extracellular matrix (ECM) architecture of unknown etiology. Here we exploit keloid scars as a paradigm to understand fibrotic ECM organization. We reveal that keloid patient fibroblasts uniquely produce a globally aligned ECM network in 2-D culture as observed in scar tissue. ECM anisotropy develops after rapid initiation of a fibroblast supracellular actin network, suggesting that cell alignment initiates ECM patterning. Keloid fibroblasts produce elevated levels of IL-6, and autocrine IL-6 production is both necessary and sufficient to induce cell and ECM alignment, as evidenced by ligand stimulation of normal dermal fibroblasts and treatment of keloid fibroblasts with the function blocking IL-6 receptor monoclonal antibody, tocilizumab. Downstream of IL-6, supracellular organization of keloid fibroblasts is controlled by activation of cell-cell adhesion. Adhesion formation inhibits contact-induced cellular overlap leading to nematic organization of cells and an alignment of focal adhesions. Keloid fibroblasts placed on isotropic ECM align the pre-existing matrix, suggesting that focal adhesion alignment leads to active anisotropic remodeling. These results show that IL-6-induced fibroblast cooperativity can control the development of a nematic ECM, highlighting both IL-6 signaling and cell-cell adhesions as potential therapeutic targets to inhibit this common feature of fibrosis.

Does hyperthermic povidone-iodine lavage increase the apoptotic rate of residual cancer cell in patients with malignant pleural mesothelioma?-a prospective pilot study

Background

Malignant pleural mesothelioma (MPM) is an incurable, late presenting primary cancer, conferring a survival of 8-14 months. Different intrapleural treatments have been tested as part of a multimodality approach to treat a select group of patients with limited disease, increasing survival. Recently, povidone-iodine has been shown to induce apoptosis in microscopic tumour cells in vitro, with no reported complications. This is the first in vivo study assessing the apoptotic rate caused by intraoperative hyperthermic betadine lavage using routine immunohistochemistry combined with transmission electron microscopy (TEM).

Methods

We included surgically fit patients aged >18, undergoing minimally invasive video-assisted thoracoscopic surgery (VATS) pleural biopsy between December 2016 and February 2018, for confirmed or presumed pleural malignancy. Parietal pleural biopsies were obtained at 7.5, 15 and 30 minutes after hyperthermic betadine lavage, and compared to pre-lavage biopsy samples, for apoptotic changes. Viable tumour samples underwent histological, immunohistochemical and ultrastructural analysis as well as TEM for features of apoptosis.

Results

N=6. Median age was 76 years. Median overall survival was 26.7 months. There was no statistical impact on survival of side of disease (left vs. right). There was no significant difference in expressions of markers of apoptotic index pre and post betadine treatment upon immunohistochemical analysis. There was no discernible effect on morphological features of apoptosis seen with betadine treatment, on TEM analysis. No side effects were identified post betadine lavage.

Conclusions

Although hyperthermic betadine lavage is a safe antiseptic solution with no toxicity when performed intraoperatively, it confers no effect on apoptotic rate or necrosis. It is therefore unlikely that hyperthermic betadine lavage will have an impact on reducing the microscopic residual disease after pleurectomy decortication and enhancing survival.

Cortical wiring by synapse type-specific control of local protein synthesis

Neurons use local protein synthesis to support their morphological complexity, which requires independent control across multiple subcellular compartments up to the level of individual synapses. We identify a signaling pathway that regulates the local synthesis of proteins required to form excitatory synapses on parvalbumin-expressing (PV+) interneurons in the mouse cerebral cortex. This process involves regulation of the TSC subunit 2 (Tsc2) by the Erb-B2 receptor tyrosine kinase 4 (ErbB4), which enables local control of messenger RNA {mRNA} translation in a cell type-specific and synapse type-specific manner. Ribosome-associated mRNA profiling reveals a molecular program of synaptic proteins downstream of ErbB4 signaling required to form excitatory inputs on PV+ interneurons. Thus, specific connections use local protein synthesis to control synapse formation in the nervous system.

Bacteriological Survey of Fresh Minced Beef on Sale at Retail Outlets in Scotland in 2019: Three Foodborne Pathogens, Hygiene Process Indicators, and Phenotypic Antimicrobial Resistance

ABSTRACT

The health and economic burden of foodborne illness is high, with approximately 2.4 million cases occurring annually in the United Kingdom. A survey to understand the baseline microbial quality and prevalence of food-related hazards of fresh beef mince on retail sale could inform risk assessment, management, and communication to ensure the safety of this commodity. In such a survey, a two-stage sampling design was used to reflect variations in population density and the market share of five categories of retail outlets in Scotland. From January to December 2019, 1,009 fresh minced beef samples were collected from 15 geographic areas. The microbial quality of each sample was assessed using aerobic colony count and Escherichia coli count. Samples were cultured for Campylobacter and Salmonella, and PCR was used to detect target genes (stx1 all variants, stx2 a to g, and rfbO157) for Shiga toxin-producing E. coli (STEC). The presence of viable E. coli O157 and STEC in samples with a positive PCR signal was confirmed via culture and isolation. Phenotypic antimicrobial sensitivity patterns of cultured pathogens and 100 E. coli isolates were determined, mostly via disk diffusion. The median aerobic colony count and E. coli counts were 6.4 × 105 (interquartile range, 6.9 × 104 to 9.6 × 106) and <10 CFU/g (interquartile range, <10 to 10) of minced beef, respectively. The prevalence was 0.1% (95% confidence interval [CI], 0 to 0.7%) for Campylobacter, 0.3% (95% CI, 0 to 1%) for Salmonella, 22% (95% CI, 20 to 25%) for PCR-positive STEC, and 4% (95% CI, 2 to 5%) for culture-positive STEC. The evidence for phenotypic antimicrobial resistance detected did not give cause for concern, mainly occurring in a few E. coli isolates as single nonsusceptibilities to first-line active substances. The low prevalence of pathogens and phenotypic antimicrobial resistance is encouraging, but ongoing consumer food safety education is necessary to mitigate the residual public health risk.

HIGHLIGHTS

Pseudomonas aeruginosa adapts to octenidine via a combination of efflux and membrane remodelling

Pseudomonas aeruginosa is an opportunistic pathogen capable of stably adapting to the antiseptic octenidine by an unknown mechanism. Here we characterise this adaptation, both in the laboratory and a simulated clinical setting, and identify a novel antiseptic resistance mechanism. In both settings, 2 to 4-fold increase in octenidine tolerance was associated with stable mutations and a specific 12 base pair deletion in a putative Tet-repressor family gene (smvR), associated with a constitutive increase in expression of the Major Facilitator Superfamily (MFS) efflux pump SmvA. Adaptation to higher octenidine concentrations led to additional stable mutations, most frequently in phosphatidylserine synthase pssA and occasionally in phosphatidylglycerophosphate synthase pgsA genes, resulting in octenidine tolerance 16- to 256-fold higher than parental strains. Metabolic changes were consistent with mitigation of oxidative stress and altered plasma membrane composition and order. Mutations in SmvAR and phospholipid synthases enable higher level, synergistic tolerance of octenidine.

Impacts of Metabolism and Organic Acids on Cell Wall Composition and Pseudomonas aeruginosa Susceptibility to Membrane Active Antimicrobials

Reliable antimicrobial susceptibility testing is essential in informing both clinical antibiotic therapy decisions and the development of new antibiotics. Mammalian cell culture media have been proposed as an alternative to bacteriological media, potentially representing some critical aspects of the infection environment more accurately. Here, we use a combination of NMR metabolomics and electron microscopy to investigate the response of Escherichia coli and Pseudomonas aeruginosa to growth in differing rich media to determine whether and how this determines metabolic strategies, the composition of the cell wall, and consequently susceptibility to membrane active antimicrobials including colistin and tobramycin. The NMR metabolomic approach is first validated by characterizing the expected E. coli acid stress response to fermentation and the accompanying changes in the cell wall composition, when cultured in glucose rich mammalian cell culture media. Glucose is not a major carbon source for P. aeruginosa but is associated with a response to osmotic stress and a modest increase in colistin tolerance. Growth of P. aeruginosa in a range of bacteriological media is supported by consumption of formate, an important electron donor in anaerobic respiration. In mammalian cell culture media, however, the overall metabolic strategy of P. aeruginosa is instead dependent on consumption of glutamine and lactate. Formate doping of mammalian cell culture media does not alter the overall metabolic strategy but is associated with polyamine catabolism, remodelling of both inner and outer membranes, and a modest sensitization of P. aeruginosa PAO1 to colistin. Further, in a panel of P. aeruginosa isolates an increase between 2- and 3-fold in sensitivity to tobramycin is achieved through doping with other organic acids, notably propionate which also similarly enhances the activity of colistin. Organic acids are therefore capable of nonspecifically influencing the potency of membrane active antimicrobials.

Effects of Germline VHL Deficiency on Growth, Metabolism, and Mitochondria

Mutations in VHL, which encodes von Hippel-Lindau tumor suppressor (VHL), are associated with divergent diseases. We describe a patient with marked erythrocytosis and prominent mitochondrial alterations associated with a severe germline VHL deficiency due to homozygosity for a novel synonymous mutation (c.222C→A, p.V74V). The condition is characterized by early systemic onset and differs from Chuvash polycythemia (c.598C→T) in that it is associated with a strongly reduced growth rate, persistent hypoglycemia, and limited exercise capacity. We report changes in gene expression that reprogram carbohydrate and lipid metabolism, impair muscle mitochondrial respiratory function, and uncouple oxygen consumption from ATP production. Moreover, we identified unusual intermitochondrial connecting ducts. Our findings add unexpected information on the importance of the VHL-hypoxia-inducible factor (HIF) axis to human phenotypes. (Funded by Associazione Italiana Ricerca sul Cancro and others.).

Mitochondrial abnormalities and disruption of the neuromuscular junction precede the clinical phenotype and motor neuron loss in hFUSWT transgenic mice

FUS (fused in sarcoma) mislocalization and cytoplasmic aggregation are hallmark pathologies in FUS-related amyotrophic lateral sclerosis and frontotemporal dementia. Many of the mechanistic hypotheses have focused on a loss of nuclear function in the FUS-opathies, implicating dysregulated RNA transcription and splicing in driving neurodegeneration. Recent studies describe an additional somato-dendritic localization for FUS in the cerebral cortex implying a regulatory role in mRNA transport and local translation at the synapse. Here, we report that FUS is also abundant at the pre-synaptic terminal of the neuromuscular junction (NMJ), suggesting an important function for this protein at peripheral synapses. We have previously reported dose and age-dependent motor neuron degeneration in transgenic mice overexpressing human wild-type FUS, resulting in a motor phenotype detected by ∼28 days and death by ∼100 days. Now, we report the earliest structural events using electron microscopy and quantitative immunohistochemistry. Mitochondrial abnormalities in the pre-synaptic motor nerve terminals are detected at postnatal day 6, which are more pronounced at P15 and accompanied by a loss of synaptic vesicles and synaptophysin protein coupled with NMJs of a smaller size at a time when there is no detectable motor neuron loss. These changes occur in the presence of abundant FUS and support a peripheral toxic gain of function. This appearance is typical of a ‘dying-back’ axonopathy, with the earliest manifestation being mitochondrial disruption. These findings support our hypothesis that FUS has an important function at the NMJ, and challenge the ‘loss of nuclear function’ hypothesis for disease pathogenesis in the FUS-opathies.

VTEC infections and livestock-related exposures in Scotland, 2004

Scotland generally reports higher rates of infection with verotoxin-producing Escherichia coli O157 (VTEC O157) than other parts of the United Kingdom

Evaluation of the Escherichia coli threonine deaminase gene as a selectable marker for plant transformation

The initial step in the synthesis of isoleucine (Ile) is the conversion of threonine to alpha-ketobutyrate. This reaction is carried out by threonine deaminase (TD), which is feedback-regulated by Ile. Mutations in TD that manifest insensitivity to Ile feedback inhibition result in intracellular accumulation of Ile. Previous reports have shown that in planta expression of the wild-type Escherichia coli TD, ilvA, or an Ile-insensitive mutant designated ilvA-466, increased cellular concentrations of Ile. A structural analog of Ile, l-O-methylthreonine (OMT), is able to compete effectively with Ile during translation and induce cell death. It has been postulated that OMT could therefore be utilized as an effective selective agent in plant engineering studies. To test this concept, we designed two binary plasmids that harbored an nptII cassette and either the wild-type ilvA or mutant ilvA-466. The ilvA coding sequences were fused to a plastid transit peptide down stream of a modified 35S CaMV promoter. Tobacco transformations were set up implementing a selection protocol based on either kanamycin or OMT. The ilvA gene was effectively utilized as a selectable marker gene to identify tobacco transformants when coupled with OMT as the selection agent. However, the transformation efficiency was substantially lower than that observed with nptII using kanamycin as the selection agent. Moreover, in a subset of the ilvA transformants and in a majority of the ilvA-466 transgenic lines, a severe off-type was observed under greenhouse conditions that correlated with increased levels of expression of the ilvA transgene.

Discovering patterns in Plasmodium falciparum genomic DNA

A method has been developed for discovering patterns in DNA sequences. Loosely based on the well-known Lempel Ziv model for text compression, the model detects repeated sequences in DNA. The repeats can be forward or inverted, and they need not be exact. The method is particularly useful for detecting distantly related sequences, and for finding patterns in sequences of biased nucleotide composition, where spurious patterns are often observed because the bias leads to coincidental nucleotide matches. We show here the utility of the method by applying it to genomic sequences of Plasmodium falciparum. A single scan of chromosomes 2 and 3 of P. falciparum, using our method and no other a priori information about the sequences, reveals regions of low complexity in both telomeric and central regions, long repeats in the subtelomeric regions, and shorter repeat areas in dense coding regions. Application of the method to a recently sequenced contig of chromosome 10 that has a particularly biased base composition detects a long internal repeat more readily than does the conventional dot matrix plot. Space requirements are linear, so the method can be used on large sequences. The observed repeat patterns may be related to large-scale chromosomal organization and control of gene expression. The method has general application in detecting patterns of potential interest in newly sequenced genomic material.

Fast, optimal alignment of three sequences using linear gap costs

Alignment algorithms can be used to infer a relationship between sequences when the true relationship is unknown. Simple alignment algorithms use a cost function that gives a fixed cost to each possible point mutation-mismatch, deletion, insertion. These algorithms tend to find optimal alignments that have many small gaps. It is more biologically plausible to have fewer longer gaps rather than many small gaps in an alignment. To address this issue, linear gap cost algorithms are in common use for aligning biological sequence data. More reliable inferences are obtained by aligning more than two sequences at a time. The obvious dynamic programming algorithm for optimally aligning k sequences of length n runs in O(n(k)) time. This is impractical if k>/=3 and n is of any reasonable length. Thus, for this problem there are many heuristics for aligning k sequences, however, they are not guaranteed to find an optimal alignment. In this paper, we present a new algorithm guaranteed to find the optimal alignment for three sequences using linear gap costs. This gives the same results as the dynamic programming algorithm for three sequences, but typically does so much more quickly. It is particularly fast when the (three-way) edit distance is small. Our algorithm uses a speed-up technique based on Ukkonen's greedy algorithm (Ukkonen, 1983) which he presented for two sequences and simple costs.

Clinical effectiveness: testicular self-examination

Evidence must be used to support a change in practice. Incorporating evidence into practice is essential for clinical effectiveness. The findings should be disseminated to others involved. The effect on patients should be evaluated.

Sequence complexity for biological sequence analysis

A new statistical model for DNA considers a sequence to be a mixture of regions with little structure and regions that are approximate repeats of other subsequences, i.e. instances of repeats do not need to match each other exactly. Both forward- and reverse-complementary repeats are allowed. The model has a small number of parameters which are fitted to the data. In general there are many explanations for a given sequence and how to compute the total probability of the data given the model is shown. Computer algorithms are described for these tasks. The model can be used to compute the information content of a sequence, either in total or base by base. This amounts to looking at sequences from a data-compression point of view and it is argued that this is a good way to tackle intelligent sequence analysis in general.

Genetic heterogeneity of Escherichia coli O157:H7 in Scotland and its utility in strain subtyping

From April 1994 to March 1995, seven outbreaks of Escherichia coli O157:H7 infection occurred throughout Scotland, including the largest milk-borne outbreak to date worldwide. Various vehicles of infection were identified, and there were 144 confirmed cases in total. All isolates associated with the outbreaks were subjected to detailed subtyping: phage typing, testing for carriage of verotoxin genes (VT), and pulsed-field gel electrophoresis. The outbreak strains were of three different phage types (2, 4, and 28). Those of phage type 2 and 28 were VT1-/VT2+, those of phage type 4 were VT1+/VT2+. To discriminate outbreak-associate isolates from the high sporadic background, real-time pulsed-field gel electrophoresis analyses were performed. The results demonstrated that, within each of the seven outbreak groups, the macrorestriction profiles observed were indistinguishable, whereas profiles for sporadic isolates were not. The consistent genetic heterogeneity observed within the Scottish Escherichia coli O157 population can be exploited in epidemiological investigations.

rbcL Transcript levels in tobacco plastids are independent of light: reduced dark transcription rate is compensated by increased mRNA stability

The plastid rbcL gene, encoding the large subunit of ribulose-1, 5-bisphosphate carboxylase, in higher plants is transcribed from a sigma70 promoter by the eubacterial-type RNA polymerase. To identify regulatory elements outside of the rbcL -10/-35 promoter core, we constructed transplastomic tobacco plants with uidA reporter genes expressed from rbcL promoter derivatives. Promoter activity was characterized by measuring steady state levels of uidA mRNA on RNA gel blots and by measuring promoter strength in run-on transcription assays. We report here that the rbcL core promoter is sufficient to obtain wild-type rates of transcription. Furthermore, the rates of transcription were up to 10-fold higher in light-grown leaves than in dark-adapted plants. Although the rates of transcription were lower in the dark, rbcL mRNA accumulated to similar levels in light-grown and dark-adapted leaves. Accumulation of uidA mRNA from most rbcL promoter deletion derivatives directly reflected the relative rates of transcription: high in the light-grown and low in the dark-adapted leaves. However, uidA mRNA accumulated to high levels in a light-independent fashion as long as a segment encoding a stem-loop structure in the 5' untranslated region was included in the promoter construct. This finding indicates that lower rates of rbcL transcription in the dark are compensated by increased mRNA stability.

Discovering simple DNA sequences by compression

An information-theoretic DNA compression scheme devised by Milosavljevic and Jurka (1993) has been used in many places in the literature for both the discovery of new genes and the compression of DNA. Their compression method applies an encoding of previously occurring runs. They use 5 different code-words: four being the DNA bases, A, C, G and T, and the other being a pointer to a previously occurring run. They advocate a code-word of length log2 5 for each of these and then encoding a run by a code-word of length 2 x log2 n, where n is the length of the sequence. This scheme encodes the start of the sequence with a code-word of length log2 n and likewise encodes the end of the sequence with a code-word of length log2 n. In this paper, we show the above coding scheme to be inefficient in various ways and improve upon it so that it can compress DNA. We discuss our implementation of various schemes some of which run in linear time.

An MML classification of protein structure that knows about angles and sequence

The MML classification program, Snob, deals with mixture modelling (or clustering) of circular data. It has recently been extended to do Markov modelling of the serial correlation between clusters such as modelling the fact that a Helix cluster favours being followed by another Helix cluster. Such a model is better known as a Hidden Markov Model. The search for the most appropriate secondary structure classification of protein data is of significant importance and was addressed by Hunter and States (1992) using the Bayesian classifier, AutoClass, on Cartesian co-ordinate data of protein residues. Dowe et al. (1996) improved upon this earlier work by using Snob to cluster dihedral angle data, with the advantage that 3 x 3 = 9 Cartesian co-ordinates can be represented by the 2 orientation-invariant angles, phi and psi. The Hidden Markov Model used here is shown to be a more appropriate way again of modelling protein data and results in the selection of a simpler class model with 17 structure classes. We report on this classification, including the class transition matrix, and relate it back to the amino-acid sequence and the simple Helix, Beta, Turn classification. We find 3 types of Helix, 2 types of Beta and many types of Turn. The msot numerous Turn class defines a continuous flexible structure that is negatively correlated to all the other classes.

Compression of strings with approximate repeats

We describe a model for strings of characters that is loosely based on the Lempel Ziv model with the addition that a repeated substring can be an approximate match to the original substring; this is close to the situation of DNA, for example. Typically there are many explanations for a given string under the model, some optimal and many suboptimal. Rather than commit to one optimal explanation, we sum the probabilities over all explanations under the model because this gives the probability of the data under the model. The model has a small number of parameters and these can be estimated from the given string by an expectation-maximization (EM) algorithm. Each iteration of the EM algorithm takes O(n2) time and a few iterations are typically sufficient. O(n2) complexity is impractical for strings of more than a few tens of thousands of characters and a faster approximation algorithm is also given. The model is further extended to include approximate reverse complementary repeats when analyzing DNA strings. Tests include the recovery of parameter estimates from known sources and applications to real DNA strings.

Circular clustering of protein dihedral angles by Minimum Message Length

Early work on proteins identified the existence of helices and extended sheets in protein secondary structures, a high-level classification which remains popular today. Using the Snob program for information-theoretic Minimum Message Length (MML) classification, we are able to take the protein dihedral angles as determined by X-ray crystallography, and cluster sets of dihedral angles into groups. Previous work by Hunter and States has applied a similar Bayesian classification method, AutoClass, to protein data with site position represented by 3 Cartesian co-ordinates for each of the alpha-Carbon, beta-Carbon and Nitrogen, totalling 9 co-ordinates. By using the von Mises circular distribution in the Snob program, we are instead able to represent local site properties by the two dihedral angles, phi and psi. Since each site can be modelled as having 2 degrees of freedom, this orientation-invariant dihedral angle representation of the data is more compact than that of nine highly-correlated Cartesian co-ordinates. Using the information-theoretic message length concepts discussed in the paper, such a more concise model is more likely to represent the underlying generating process from which the data came. We report on the results of our classification, plotting the classes in (phi, psi) space; and introducing a symmetric information-theoretic distance measure to build a minimum spanning tree between the classes. We also give a transition matrix between the classes and note the existence of three classes in the region phi approximately -1.09 rad and psi approximately -0.75 rad which are close on the spanning tree and have high inter-transition probabilities. This gives rise to a tight, abundant and self-perpetuating structure.

Relationship between head growth and neurodevelopmental outcome of Malaysian very low birthweight infants during the 1st year of life

A prospective study was carried out to (i) compare head growth patterns of 103 very low birthweight (VLBW, < 1500 g) Malaysian infants and 98 normal birthweight (NBW, 2500- < 4500 g) controls during the 1st year of life; and (ii) examine the relationship between neurodevelopmental outcome at 1 year of age and occipito-frontal head circumferences (OFC) at birth and at 1 year of age in VLBW babies. When compared with those of NBW infants at birth, mid-infancy and 1 year of age, the mean OFC ratios (observed/expected OFC at 50th percentile) of VLBW infants were significantly lower (p < 0.001). Small-for-gestational-age (SGA) VLBW babies had significantly lower mean OFC ratios than their appropriate-for-gestational-age (AGA) VLBW counterparts at birth (p < 0.001), but this difference was no longer seen at mid-infancy or at 1 year of age. Logistic regression analysis showed that abnormal late neonatal cranial ultrasound findings (odds ratio 8.5, 95% confidence interval 4.12-22.07; p < 0.001) and each additional day of oxygen therapy (odds ratio 1.15, 95% confidence interval 1.00-4.45; p = 0.045) were significant risk factors associated with neurodevelopmental disability at 1 year of age, while mean OFC ratios at birth or at 1 year of age were not. Poor postnatal head growth per se did not predict disability, but probably reflected the consequences of "brain injury" as evidenced by abnormal brain scans.

Neurodevelopmental outcome of Malaysian very low birth weight infants: predictive value of cranial ultrasound appearances

The aim of the study was to determine the predictive value of cranial ultrasound scans done in the neonatal period for neurodevelopmental outcome of the Malaysian very low birthweight (VLBW, < 1500 grams) infants assessed at 12 months of corrected age. Of the 101 infants studied, 68 (67.3%) were neurodevelopmentally normal at one year of age, 18 (17.8%) had major and 15 (14.9%) had minor neurodevelopmental impairment. Neurodevelopmental outcome was normal in 66/88 (75.0%) infants who did not have severe intraventricular haemorrhage (IVH) or periventricular intraparenchymal echo densities (PVE) in the first week of life, and in 57/73 (78.1%) with uncomplicated scans at discharge. In contrast, 11/13 (84.6%) with parenchymal echo densities or severe intraventricular bleed in the early neonatal period and 17/28 (60.7%) with complicated scans at discharge had adverse sequelae. There was a significant association between lesions seen on cranial ultrasound in the neonatal period and subsequent neurodevelopmental impairment. Late neonatal ultrasound scans appear to be a better predictor of short-term neurodevelopmental outcome than early scans.

Comparison of morbidities in very low birthweight and normal birthweight infants during the first year of life in a developing country

OBJECTIVE

To compare the morbidities in the very low birthweight (VLBW; < 1500 g) and normal birthweight (NBW; > or = 2500 g) Malaysian infants during the first year of life.

METHODOLOGY

Prospective observational cohort study of consecutive surviving VLBW infants and randomly sampled NBW infants born in the Kuala Lumpur Maternity Hospital between 1 December 1989 and 31 December 1992. Infants were followed up regularly during the first year of life, after correction for prematurity.

RESULTS

Compared with NBW infants (n = 106), VLBW infants (n = 127) had significantly higher risk of failure to thrive (odds ratio [OR] = 8.0, 95% confidence intervals [CI]: 1.1 to 354.3), wheezing (OR = 3.7, 95% CI: 1.6 to 9.3), rehospitalization (OR = 2.3, 95% CI: 1.1 to 5.0), cerebral palsy (OR = 8.6, 95% CI: 2.0 to 77.6), neurosensory hearing loss (OR = 12.0, 95% CI: 1.7 to 513.6) and visual loss (7.9 vs 0%, P = 0.002). The mean mental developmental index (MDI) and mean psychomotor developmental index (PDI) at 1 year of age were significantly lower among VLBW infants (MDI 99 [SD = 28], PDI 89 [SD = 25]) than NBW infants (MDI 106 [SD = 18], PDI 101 [SD = 18]) (95% CI for difference between means being MDI: -14.1 to -1.7; and PDI: -17.6 to -6.0). Logistic regression analysis showed that among VLBW infants: (i) male sex, Malay ethnicity and bronchopulmonary dysplasia were significant risk factors associated with wheezing; (ii) longer duration of oxygen therapy during the neonatal period, seizures after the post-neonatal period and wheezing were significant risk factors associated with rehospitalization; and (iii) longer duration of oxygen therapy during the neonatal period was a significant risk factor associated with adverse neurodevelopmental outcome during the first year of life.

CONCLUSIONS

Compared with NBW infants, VLBW Malaysian infants had significantly higher risks of physical and neuro-developmental morbidities.

Long-term efficacy of continuing hepatitis B vaccination in infancy in two Gambian villages

In 1984, all non-immune children under the age of 5 years in the Gambian villages of Keneba and Manduar were vaccinated against hepatitis B virus (HBV). All children born in these villages since 1984 have been vaccinated in infancy. Despite a rapid fall in antibody concentrations, vaccine efficacy against HBV infection and chronic carriage of HBsAg has increased with time. Overall, vaccine efficacies in 1993 against HBV infection and chronic HBsAg carriage were 94.7% (95% Cl 93.0-96.0) and 95.3% (91.0-97.5), respectively. Breakthrough infections in vaccinated children largely originate from chronic HBsAg carriers. Thus, we tested 261 chronic carriers for HBV DNA and e antigen. The prevalence of these markers of infectivity, and the amount of HBV DNA, decreased greatly with age. Detailed studies of breakthrough infections over two 4-year periods revealed that in the second period there were fewer than half the expected numbers of infections. Our findings suggest that in Keneba and Manduar long-term vaccination is progressively decreasing HBV transmission by chronic carriers, since their infectivity diminishes with time.

The posterior probability distribution of alignments and its application to parameter estimation of evolutionary trees and to optimization of multiple alignments

How to sample alignments from their posterior probability distribution given two strings is shown. This is extended to sampling alignments of more than two strings. The result is first applied to the estimation of the edges of a given evolutionary tree over several strings. Second, when used in conjunction with simulated annealing, it gives a stochastic search method for an optimal multiple alignment.

Breakthrough infections and identification of a viral variant in Gambian children immunized with hepatitis B vaccine

Hepatitis B (HB) breakthrough infections, identified by the presence of HB core (c) antibody, were found in 32 of 358 Gambian children vaccinated with plasma-derived HB vaccine. Over 2 years, 15 of these children lost their HBc antibodies. These children had significantly higher HB surface antibody levels before infection than those who retained HBc antibodies. One child, who responded well to the vaccine, had HB viral DNA detected in the presence of HBs antibodies. The S gene sequence of this DNA showed nucleotide changes that resulted in an amino acid substitution at residue 141 (lysine to glutamic acid) of the surface antigen. This finding suggests the child was infected with a variant virus that was not neutralized by antibodies resulting from HB vaccination.

A fast algorithm for the optimal alignment of three strings

Ukkonen's (pair-wise) string alignment technique is extended to the problem of finding an optimal alignment for three strings. The resulting algorithm has worst-case time-complexity O(nd2) and space-complexity O(d3), where the string lengths are ñ and d is the three-way edit-distance based on tree-costs. In practice, the algorithm usually runs in O(n + d3) time. The algorithm is particularly fast when the strings are similar, in which case, d << n. Three-way alignment is an important special case in string alignment. Each internal node in an unrooted, binary evolutionary-tree has three neighbours. The algorithm presented can be used as an iterative step in a heuristic multiple-alignment program for more than three strings.

Normalization of affine gap costs used in optimal sequence alignment

It is shown how to normalize the costs of an alignment algorithm that employs affine or linear gap costs. The normalized costs are interpreted as the -log probabilities of the instructions of a finite-state edit-machine. This gives an explicit model relating sequences that can be linked to processes of mutation and evolution.

Reconstruction of strings past

A major use of string-alignment algorithms is to compare macromolecules that are thought to have evolved from a common ancestor to estimate the duration of, or the amount of mutation in, their separate evolution and to infer as much as possible about their most recent common ancestor. Minimum message length encoding, a method of inductive inference, is applied to the string-alignment problem. It leads to an alignment method that averages over all alignments in a weighted fashion. Experiments indicates that this method can recover the actual parameters of evolution with high accuracy and over a wide range of values, whereas the use of a single optimal alignment gives biased results.

Finite-state models in the alignment of macromolecules

Minimum message length encoding is a technique of inductive inference with theoretical and practical advantages. It allows the posterior odds-ratio of two theories or hypotheses to be calculated. Here it is applied to problems of aligning or relating two strings, in particular two biological macromolecules. We compare the r-theory, that the strings are related, with the null-theory, that they are not related. If they are related, the probabilities of the various alignments can be calculated. This is done for one-, three-, and five-state models of relation or mutation. These correspond to linear and piecewise linear cost functions on runs of insertions and deletions. We describe how to estimate parameters of a model. The validity of a model is itself an hypothesis and can be objectively tested. This is done on real DNA strings and on artificial data. The tests on artificial data indicate limits on what can be inferred in various situations. The tests on real DNA support either the three- or five-state models over the one-state model. Finally, a fast, approximate minimum message length string comparison algorithm is described.

Restriction site mapping for three or more enzymes

Restriction site mapping requires a generator to put forward possible maps and a constraint checker to reject false maps. Ideally these combine to give an algorithm which calculates a sound and complete solution set. Three algorithms for generation are presented and compared. Two decompose a multi-enzyme problem (greater than or equal to 3) into subproblems. The constraint checker is based on separation theory. Some insights into the extent of constraint checking involved in and feasibility of more checking for three or more enzymes are discussed. The trade-off between computation time and the soundness of the solution set is examined.

Minimum message length encoding and the comparison of macromolecules

A comparison of inductive inference known as minimum message length encoding is applied to string comparison in molecular biology. The question of whether or not two strings are related and, if so, of how they are related and the problem of finding a good theory of string mutation are treated as inductive inference problems. The method allows the posterior odds-ratio of two string alignments or of two models of string mutation to be computed. The connection between models of mutation and existing string alignment algorithms is made explicit. A fast minimum message length alignment algorithm is also described.

Restriction site mapping is in separation theory

A computer algorithm for restriction-site mapping consists of a generator of partial maps and a consistency checker. This paper examines consistency checking and argues that a method based on separation theory extracts the maximum amount of information from fragment lengths in digest data. It results in the minimum number of false maps being generated.

Mapping of the pattern of DNA replication in polytene chromosome from Chironomus thummi using monoclonal anti-bromodeoxyuridine antibodies

We present results from a nonautoradiographic study of DNA replication in polytene chromosomes from dipteran larvae. Monoclonal antibodies with specificity for 5-bromodeoxyuridine (BrdUrd) were used to localize by indirect immunofluorescence the sites of BrdUrd incorporation and to follow the dynamics of DNA synthesis in salivary gland cells of 4th instar Chironomus thummi larvae. This technique presents numerous advantages over autoradiographic procedures and allows mapping of DNA synthesis patterns at the level of resolution of one chromosomal band. Several replication patterns were observed, classified according to characteristic features, and tentatively assigned to specific periods of the S-phase. In early S-phase, DNA synthesis is first detectable in puffs and interbands, later in bands. Most chromosomal bands appear to initiate DNA synthesis synchronously; however, in bands within centromeric and heterochromatic regions the start of synthesis is delayed. At mid S-phase, all the bands show uniform staining. Subsequent staining patterns are increasingly differential with the bands displaying characteristic fluorescence intensities. As replication progresses through the late S-phase period, the chromosomes show a decreasing number of fluorescent bands. The last bands to terminate replication are located in centromeric and heterochromatic DNA-rich regions and a few bands of low DNA content in region IIAa-c.

Properties and characterization of monoclonal antibodies to Tacaribe virus

Monoclonal antibodies prepared against Tacaribe and Junin viruses have been used to define further the serological relationships between arenaviruses of the Tacaribe complex. A close relationship was found between these two viruses and the heterologous Amapari and Machupo viruses, with Pichinde virus and Parana virus being more distantly related. Among the antibodies specific for Tacaribe virus, five were found to react with viral antigens at the surface of infected cells and to neutralize virus infectivity in vitro. These five antibodies could be differentiated by competitive immunoassay as recognizing at least two antigenically distinct epitopes. The kinetics of reaction between antibody and virus were examined for all five neutralizing antibodies. One antibody (2.25.4) effectively neutralized all infectious virus. The remaining four directed against a second epitope gave significant persistent fractions which could be reduced by addition of complement, anti-mouse immunoglobulin, or antibody 2.25.4. Variants of Tacaribe virus resistant to neutralization by antibody 2.25.4 were obtained by growth in the presence of this antibody and neutralization kinetics were reexamined using the heterologous monoclonal neutralizing antibodies. Several different neutralization profiles were obtained, suggesting that point mutations resulted in conformational changes at topographically selected distinct epitopes recognized by the remaining antibodies.

Thermal inactivation of Pichinde virus

Detailed information regarding the kinetics of thermal inactivation of Pichinde, an arenavirus, is presented. Inactivation of virus infectivity proceeded as a first order reaction over the temperature range 22-53 degrees C. The determined inactivation rates analysed as a function of absolute temperature revealed that two different reactions were involved. Below 37 degrees C, the energy of activation was determined to be compatible with RNA degradation, whereas at higher temperatures a correspondingly greater value suggests that protein inactivation contributes significantly to loss of infectivity. Both inactivation reactions were retarded in the presence of foetal calf serum to a final concentration of 1%. The relatively short half-life of 12-24 h at 22 degrees C suggests transmission in nature via contaminated foodstuffs and soil may be inefficient.