Short report: historical analysis of a near disaster: Anopheles gambiae in Brazil

Attributed to human-mediated dispersal, a species of the Anopheles gambiae complex invaded northeastern Brazil in 1930. This event is considered unique among the intercontinental introductions of disease vectors and the most serious one: "Few threats to the future health of the Americas have equalled that inherent in the invasion of Brazil, in 1930, by Anopheles gambiae." Because it was only in the 1960s that An. gambiae was recognized as a species complex now including seven species, the precise species identity of the Brazilian invader remains a mystery. Here we used historical DNA analysis of museum specimens, collected at the time of invasion from Brazil, and aimed at the identification of the Brazilian invader. Our results identify the arid-adapted Anopheles arabiensis as being the actual invading species. Establishing the identity of the species, in addition to being of intrinsic historical interest, can inform future threats of this sort especially in a changing environment. Furthermore, these results highlight the potential danger of human-mediated range expansions of insect disease vectors and the importance of museum collections in retrieving historical information

Phylogeography of the common vampire bat (Desmodus rotundus): Marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers

Background: The common vampire bat Desmodus rotundus is an excellent model organism for studying ecological vicariance in the Neotropics due to its broad geographic range and its preference for forested areas as roosting sites. With the objective of testing for Pleistocene ecological vicariance, we sequenced a mitocondrial DNA (mtDNA) marker and two nuclear markers (RAG2 and DRB) to try to understand how Pleistocene glaciations affected the distribution of intraspecific lineages in this bat. Results: Five reciprocally monophyletic clades were evident in the mitochondrial gene tree, and in most cases with high bootstrap support: Central America (CA), Amazon and Cerrado (AMC), Pantanal (PAN), Northern Atlantic Forest (NAF) and Southern Atlantic Forest (SAF). The Atlantic forest clades formed a monophyletic clade with high bootstrap support, creating an east/west division for this species in South America. On the one hand, all coalescent and non-coalescent estimates point to a Pleistocene time of divergence between the clades. On the other hand, the nuclear markers showed extensive sharing of haplotypes between distant localities, a result compatible with male-biased gene flow. In order to test if the disparity between the mitochondrial and nuclear markers was due to the difference in mutation rate and effective size, we performed a coalescent simulation to examine the feasibility that, given the time of separation between the observed lineages, even with a gene flow rate close to zero, there would not be reciprocal monophyly for a neutral nuclear marker. We used the observed values of theta and an estimated mutation rate for the nuclear marker gene to perform 1000 iterations of the simulation. The results of this simulation were inconclusive: the number of iterations with and without reciprocal monophyly of one or more clades are similar. Conclusions: We therefore conclude that the pattern exhibited by the common vampire bat, with marked geographical structure for a mitochondrial marker and no phylogeographic structure for nuclear markers is compatible with a historical scenario of complete isolation of refuge-like populations during the Pleistocene. The results on demographic history on this species is compatible with the Carnaval-Moritz model of Pleistocene vicariance, with demographic expansions in the southern Atlantic forest.

Industrial apiculture in the Jordan valley during Biblical times with Anatolian honeybees

Although texts and wall paintings suggest that bees were kept in the Ancient Near East for the production of precious wax and honey, archaeological evidence for beekeeping has never been found. The Biblical term ""honey"" commonly was interpreted as the sweet product of fruits, such as dates and figs. The recent discovery of unfired clay cylinders similar to traditional hives still used in the Near East at the site of Tel Rehov in the Jordan valley in northern Israel suggests that a large-scale apiary was located inside the town, dating to the 10th-early 9th centuries B.C.E. This paper reports the discovery of remains of honeybee workers, drones, pupae, and larvae inside these hives. The exceptional preservation of these remains provides unequivocal identification of the clay cylinders as the most ancient beehives yet found. Morphometric analyses indicate that these bees differ from the local subspecies Apis mellifera syriaca and from all subspecies other than A. m. anatoliaca, which presently resides in parts of Turkey. This finding suggests either that the Western honeybee subspecies distribution has undergone rapid change during the last 3,000 years or that the ancient inhabitants of Tel Rehov imported bees superior to the local bees in terms of their milder temper and improved honey yield.

Land planarians (Platyhelminthes) as a model organism for fine-scale phylogeographic studies: understanding patterns of biodiversity in the Brazilian Atlantic Forest hotspot

The Brazilian Atlantic Forest is one of the richest biodiversity hotspots of the world. Paleoclimatic models have predicted two large stability regions in its northern and central parts, whereas southern regions might have suffered strong instability during Pleistocene glaciations. Molecular phylogeographic and endemism studies show, nevertheless, contradictory results: although some results validate these predictions, other data suggest that paleoclimatic models fail to predict stable rainforest areas in the south. Most studies, however, have surveyed species with relatively high dispersal rates whereas taxa with lower dispersion capabilities should be better predictors of habitat stability. Here, we have used two land planarian species as model organisms to analyse the patterns and levels of nucleotide diversity on a locality within the Southern Atlantic Forest. We find that both species harbour high levels of genetic variability without exhibiting the molecular footprint of recent colonization or population expansions, suggesting a long-term stability scenario. The results reflect, therefore, that paleoclimatic models may fail to detect refugia in the Southern Atlantic Forest, and that model organisms with low dispersal capability can improve the resolution of these models.

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Immunogenicity of a Whole-Cell Pertussis Vaccine with Low Lipopolysaccharide Content in Infants

The lack of a clear correlation between the levels of antibody to pertussis antigens and protection against disease lends credence to the possibility that cell-mediated immunity provides primary protection against disease. This phase I comparative trial had the aim of comparing the in vitro cellular immune response and anti-pertussis toxin (anti-PT) immunoglobulin G (IgG) titers induced by a cellular pertussis vaccine with low lipopolysaccharide (LPS) content (wP(low) vaccine) with those induced by the conventional whole-cell pertussis (wP) vaccine. A total of 234 infants were vaccinated at 2, 4, and 6 months with the conventional wP vaccine or the wP(low) vaccine. Proliferation of CD3(+) T cells was evaluated by flow cytometry after 6 days of peripheral blood mononuclear cell culture with stimulation with heat-killed Bordetella pertussis or phytohemagglutinin (PHA). CD3(+), CD4(+), CD8(+), and T-cell receptor gamma delta-positive (gamma delta(+)) cells were identified in the gate of blast lymphocytes. Gamma interferon, tumor necrosis factor alpha, interleukin-4 (IL-4), and IL-10 levels in super-natants and serum anti-PT IgG levels were determined using enzyme-linked immunosorbent assay (ELISA). The net percentage of CD3(+) blasts in cultures with B. pertussis in the group vaccinated with wP was higher than that in the group vaccinated with the wP(low) vaccine (medians of 6.2% for the wP vaccine and 3.9% for the wP(low) vaccine; P = 0.029). The frequencies of proliferating CD4(+), CD8(+), and gamma delta(+) cells, cytokine concentrations in supernatants, and the geometric mean titers of anti-PT IgG were similar for the two vaccination groups. There was a significant difference between the T-cell subpopulations for B. pertussis and PHA cultures, with a higher percentage of gamma delta(+) cells in the B. pertussis cultures (P < 0.001). The overall data did suggest that wP vaccination resulted in modestly better specific CD3(+) cell proliferation, and gamma delta(+) cell expansions were similar with the two vaccines.

Temperature dependence of the magnetic susceptibility for triangular-lattice antiferromagnets with spatially anisotropic exchange constants

We present the temperature dependence of the uniform susceptibility of spin-half quantum antiferromagnets on spatially anisotropic triangular lattices, using high-temperature series expansions. We consider a model with two exchange constants J1 and J2 on a lattice that interpolates between the limits of a square lattice (J1=0), a triangular lattice (J2=J1), and decoupled linear chains (J2=0). In all cases, the susceptibility, which has a Curie-Weiss behavior at high temperatures, rolls over and begins to decrease below a peak temperature Tp. Scaling the exchange constants to get the same peak temperature shows that the susceptibilities for the square lattice and linear chain limits have similar magnitudes near the peak. Maximum deviation arises near the triangular-lattice limit, where frustration leads to much smaller susceptibility and with a flatter temperature dependence. We compare our results to the inorganic materials Cs2CuCl4 and Cs2CuBr4 and to a number of organic molecular crystals. We find that the former (Cs2CuCl4 and Cs2CuBr4) are weakly frustrated and their exchange parameters determined through the temperature dependence of the susceptibility are in agreement with neutron-scattering measurements. In contrast, the organic materials considered are strongly frustrated with exchange parameters near the isotropic triangular-lattice limit.

Fixed bed reactors with catalyst poisoning - first order kinetics.

The long performance of an isothermal fixed bed reactor undergoing catalyst poisoning is theoretically analyzed using the dispersion model. First order reaction with dth order deactivation is assumed and the model equations are solved by matched asymptotic expansions for large Peclet number. Simple closed-form solutions, uniformly valid in time, are obtained.

Substrate-inhibited kinetics with catalyst deactivation in an isothermal CSTR

Analytical expressions are developed for the time-dependent reactant concentration and catalyst activity in an isothermal CSTR with Langmuir-Hinshelwood kinetics of deactivation and reaction. Several parallel and series posioning mechanisms are considered for a reactor which, without poisoning, would operate at a unique steady state. The use of matched asymptotic expansions and abandonment of the usual initial-steady-state assumption give results, valid from startup to final loss of activity, whose accuracy can be improved systematically.

Asymmetry and duration dependence in Australian GDP and unemployment

The per iodic structure of business cycles suggests that significant asymmetries are present over different phases of the cycle. This paper uses markov regime-switching models with fixed and duration dependent transition probabilities to directly model expansions, contractions and durations in Australian GDP growth and unemployment growth. Evidence is found of significant asymmetry in growth rates across expansions and contractions for both series. GDP contractions exhibit duration dependence implying that as output recessions age the likelihood of switching into an expansion phase increases. Unemployment growth does not exhibit duration dependence in either phase. Evidence is also presented that non-linearities in unemployment growth are well explained by the asymmetries in the GDP growth cycle. The analysis suggests that recessions are periods of rapid and intense job destruction, that Australian unemployment tends to ratchet up in recessionary periods and, in contrast to US and UK studies, that shocks to Australian unemployment growth are more persistent in recessions than expansions. [E37 C5 C41].

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Is there an association between normal cognitive functioning and trinucleotide repeat expansion at loci involved in neurodegenerative disorders?

Recent reports have shown neurodegenerative disorders to be associated with abnormal expansions of a CAG trinucleotide repeat allele at various autosomal loci. While normal chromosomes have 14 to 44 repeats, disease chromosomes may have 60 to 84 repeats. The number of CAG repeats on mutant chromosomes correlates with increasing severity of disease or decreasing age at onset of symptoms. Since we are interested in identifying the many quantitative trait loci (QTL) influencing brain functioning, we examined the possibility that the number of CAG repeats in the normal size range at these loci are relevant to "normal" neural functioning. We have used 150 pairs of adolescent (aged 16 years) twins and their parents to examine allele size at the MJD, SCA1, and DRPLA loci in heterozygous normal individuals. These are part of a large ongoing project using cognitive and physiological measures to investigate the genetie influences on cognition, and an extensive protocol of tests is employed to assess some of the key components of intellectual functioning. This study selected to examine full-scale psychometric IQ (FSIQ) and a measure of information processing (choice reaction time) and working memory (slow wave amplitude). CAG repeat size was determined on an ABI Genescan system following multiplex PCR amplification. Quantitative genetic analyses were performed to determine QTL effects of MJD, SCA1, and DRPLA on cognitive functioning. Analyses are in progress and will be discussed.

X-linked myoclonic epilepsy with spasticity and intellectual disability - Mutation in the homeobox gene ARX

Objective: To describe a new syndrome of X-linked myoclonic epilepsy with generalized spasticity and intellectual disability (XMESID) and identify the gene defect underlying this disorder. Methods: The authors studied a family in which six boys over two generations had intractable seizures using a validated seizure questionnaire, clinical examination, and EEG studies. Previous records and investigations were obtained. Information on seizure disorders was obtained on 271 members of the extended family. Molecular genetic analysis included linkage studies and mutational analysis using a positional candidate gene approach. Results: All six affected boys had myoclonic seizures and TCS; two had infantile spasms, but only one had hypsarrhythmia. EEG studies show diffuse background slowing with slow generalized spike wave activity. All affected boys had moderate to profound intellectual disability. Hyperreflexia was observed in obligate carrier women. A late-onset progressive spastic ataxia in the matriarch raises the possibility of late clinical manifestations in obligate carriers. The disorder was mapped to Xp11.2-22.2 with a maximum lod score of 1.8. As recently reported, a missense mutation (1058C>T/P353L) was identified within the homeodomain of the novel human Aristaless related homeobox gene (ARX). Conclusions: XMESID is a rare X-linked recessive myoclonic epilepsy with spasticity and intellectual disability in boys. Hyperreflexia is found in carrier women. XMESID is associated with a missense mutation in ARX. This disorder is allelic with X-linked infantile spasms (ISSX; MIM 308350) where polyalanine tract expansions are the commonly observed molecular defect. Mutations of ARX are associated with a wide range of phenotypes; functional studies in the future may lend insights to the neurobiology of myoclonic seizures and infantile spasms.

Folate-sensitive fragile site FRA10A is due to an expansion of a CGG repeat in a novel gene, FRA10AC1, encoding a nuclear protein

Fragile sites appear visually as nonstaining gaps on chromosomes that are inducible by specific cell culture conditions. Expansion of CGG/ CCG repeats has been shown to be the molecular basis of all five folate-sensitive fragile sites characterized molecularly so far, i.e., FRAXA, FRAXE, FRAXF, FRA11B, and FRA16A. In the present study we have refined the localization of the FRA10A folate-sensitive fragile site by fluorescence in situ hybridization. Sequence analysis of a BAC clone spanning FRA10A identified a single, imperfect, but polymorphic CGG repeat that is part of a CpG island in the 5'UTR of a novel gene named FRA10ACl. The number of CGG repeats varied in the population from 8 to 13. Expansions exceeding 200 repeat units were methylated in all FRA10A fragile site carriers tested. The FRA10ACl gene consists of 19 exons and is transcribed in the centromeric direction from the FRA10A repeat. The major transcript of similar to 1450 nt is ubiquitously expressed and codes for a highly conserved protein, FRA10ACl, of unknown function. Several splice variants leading to alternative 3' ends were identified (particularly in testis). These give rise to FRA10ACl proteins with altered COOH-termini. Immunofluorescence analysis of full-length, recombinant EGFP-tagged FRA10ACl protein showed that it was present exclusively in the nucleoplasm. We show that the expression of FRA10A, in parallel to the other cloned folate-sensitive fragile sites, is caused by an expansion and subsequent methylation of an unstable CGG trinucleotide repeat. Taking advantage of three cSNPs within the FRA10ACl gene we demonstrate that one allele of the gene is not transcribed in a FRA10A carrier. Our data also suggest that in the heterozygous state FRA10A is likely a benign folate-sensitive fragile site. (C) 2004 Elsevier Inc. All rights reserved.