Gene and genotypic diversity of Phytophthora cinnamomi in South Africa and Australia revealed by DNA polymorphisms

Phytophthora cinnamomi isolates from South Africa and Australia were compared to assess genetic differentiation between the two populations. These two populations were analysed for levels of phenotypic diversity using random amplified polymorphic DNAs (RAPDs) and gene and genotypic diversity using restriction fragment length polymorphisms (RFLPs). Sixteen RAPD markers from four decanucleotide Operon primers and 34 RFLP alleles from 15 putative loci were used. A few isolates from Papua New Guinea known to posses alleles different from Australian isolates were also included for comparative purposes. South African and Australian P. cinnamomi populations were almost identical with an extremely low level of genetic distance between them (D-m = 0.003). Common features for the two populations include shared alleles, low levels of phenotypic/genotypic diversity, high clonality, and low observed and expected levels of heterozygosity. Furthermore, relatively high levels of genetic differentiation between mating type populations (D-m South Africa = 0.020 and D-m Australia = 0.025 respectively), negative fixation indices, and significant deviations from Hardy-Weinberg equilibrium, all provided evidence for the lack of frequent sexual reproduction in both populations. The data strongly suggest that both the South African and Australian P. cinnamomi populations are introduced.

Population structure of Phytophthora cinnamomi in South Africa

Phytophthora cinnamomi isolates collected from 1977 to 1986 and 1991 to 1993 in two regions in South Africa were analyzed using isozymes. A total of 135 isolates was analyzed for 14 enzymes representing 20 putative loci, of which four were polymorphic. This led to the identification of nine different multilocus isozyme genotypes. Both mating types of P. cinnamomi occurred commonly in the Cape region, whereas, predominantly, the A2 mating type occurred in the Mpumalanga region of South Africa. A2 mating type isolates could be resolved into seven multilocus isozyme genotypes, compared with only two multilocus isozyme genotypes for the A1 mating type isolates. Low levels of gene (0.115) and genotypic (2.4%) diversity and a low number of alleles per locus (1.43) were observed for the South African P. cinnamomi population. The genetic distance between the Cape and Mpumalanga P. cinnamomi populations was relatively low (D-m = 0.165), and no specific pattern in regional distribution of multilocus isozyme genotypes could be observed. The genetic distance between the ''old'' (isolated between 1977 and 1986) and ''new'' (isolated between 1991 and 1993) P. cinnamomi populations from the Cape was low (D-m = 0.164), indicating a stable population over time. Three of the nine multilocus isozyme genotypes were specific to the ''old'' population, and only one multilocus isozyme genotype was specific to the ''new'' population. Significant differences in allele frequencies, a high genetic distance (D-m = 0.581) between the Cape A1 and A2 mating type isolates, significant deviations from Hardy-Weinberg equilibrium, a low overall level of heterozygosity, and a high fixation index (0.71) all indicate that sexual reproduction occurs rarely, if at all, in the South African P. cinnamomi population.

Panmixia in Pocillopora verrucosa from South Africa

The genetic structure of six local collections of Pocillopora verrrucosa from six coral reefs in KwaZulu-Natal, South Africa, was examined using allozyme electrophoresis. The six separate reefs lie within two different reef complexes. Twenty-two enzymes were screened on five buffer systems, but only five polymorphic loci (Gpi-1, Gdh-1, Lgg-2, Lpp-1, Est-1) could be consistently resolved. No significant differences in allelic frequencies were detected among the six sites. All local collections were genotypically diverse, with evidence of only very limited clonal replication at each site. Indeed, the ratio of observed to expected genotypic diversity (mean Go:Ge=0.64 +/-0.05 SD), the ratio of observed number of genotypes to the number of individuals (mean Ng:N = 0.65 +/-0.04 SE), and deviations from the Hardy-Weinberg equilibrium indicate that sexual reproduction plays a major role in the maintenance of the populations. No genetic differentiation was found either within (FSR = 0.026 +/-0.003 SE) or between (FRT = 0.000 +/-0.001 SE) reef complexes. The homogeneity of the gene frequencies across the six reefs strongly supports the assumption that the KwaZulu-Natal reef complexes are highly connected by gene flow (Nem=44). The reefs in the southern and central reef complexes along the northern Maputaland coastline can therefore be considered part of a single population.

When vicars meet: A narrow contact zone between morphologically cryptic phylogeographic lineages of the rainforest skink, Carlia rubrigularis

Phylogeographic analyses of the fauna of the Australian wet tropics rainforest have provided strong evidence for long-term isolation of populations among allopatric refugia, yet typically there is no corresponding divergence in morphology. This system provides an opportunity to examine the consequences of geographic isolation, independent of morphological divergence, and thus to assess the broader significance of historical subdivisions revealed through mitochondrial DNA phylogeography. We have located and characterized a zone of secondary contact between two long isolated (mtDNA divergence > 15%) lineages of the skink Carlia rubrigularis using one mitochondrial and eight nuclear (two intron, six microsatellite) markers. This revealed a remarkably narrow (width < 3 km) hybrid zone with substantial linkage disequilibrium and strong deficits of heterozygotes at two of three nuclear loci with diagnostic alleles. Cline centers were coincident across loci. Using a novel form of likelihood analysis, we were unable to distinguish between sigmoidal and stepped cline shapes except at one nuclear locus for which the latter was inferred. Given estimated dispersal rates of 90-133 m x gen(-1/2) and assuming equilibrium, the observed cline widths suggest effective selection against heterozygotes of at least 22-49% and possibly as high as 70%. These observations reveal substantial postmating isolation, although the absence of consistent deviations from Hardy-Weinberg equilibrium at diagnostic loci suggests that there is little accompanying premating isolation. The tight geographic correspondence between transitions in mtDNA and those for nuclear genes and corresponding evidence for selection against hybrids indicates that these morphologically cryptic phylogroups could be considered as incipient species. Nonetheless, we caution against the use of mtDNA phylogeography as a sole criterion for defining species boundaries.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Analysis and application of an equilibrium model for in vitro bioassay systems with three components: Receptor, hormone and hormone-binding-protein

A simple theoretical framework is presented for bioassay studies using three component in vitro systems. An equilibrium model is used to derive equations useful for predicting changes in biological response after addition of hormone-binding-protein or as a consequence of increased hormone affinity. Sets of possible solutions for receptor occupancy and binding protein occupancy are found for typical values of receptor and binding protein affinity constants. Unique equilibrium solutions are dictated by the initial condition of total hormone concentration. According to the occupancy theory of drug action, increasing the affinity of a hormone for its receptor will result in a proportional increase in biological potency. However, the three component model predicts that the magnitude of increase in biological potency will be a small fraction of the proportional increase in affinity. With typical initial conditions a two-fold increase in hormone affinity for its receptor is predicted to result in only a 33% increase in biological response. Under the same conditions an Ii-fold increase in hormone affinity for receptor would be needed to produce a two-fold increase in biological potency. Some currently used bioassay systems may be unrecognized three component systems and gross errors in biopotency estimates will result if the effect of binding protein is not calculated. An algorithm derived from the three component model is used to predict changes in biological response after addition of binding protein to in vitro systems. The algorithm is tested by application to a published data set from an experimental study in an in vitro system (Lim et al., 1990, Endocrinology 127, 1287-1291). Predicted changes show good agreement (within 8%) with experimental observations. (C) 1998 Academic Press Limited.

Experimental study of phase equilibria in the "FeO"-ZnO-(CaO+SiO2) system with CaO/SiO2 weight ratios of 0.33, 0.93, and 1.2 in equilibrium with metallic iron

The pseudoternary sections FeO-ZnO-(CaO + SiO2) with CaO/SiO2 weight ratios of 0.33, 0.93, and 1.2 in equilibrium with metallic iron have been experimentally investigated in the temperature range from 1000 degreesC to 1300 degreesC (1273 to 1573 K). The liquidus surfaces in these pseudoternary sections have been experimentally determined in the composition range from 0 to 33 wt pct ZnO and 30 to 70 wt pct (CaO + SiO2). The sections contain primary-phase fields of wustite (FexZn1-xO1+y), zincite (ZnzFe1-zO), fayalite (Fu(w)Zn(2-w)SiO(4)), melilite (Ca2ZnuFe1-uSi2O7), willemite (ZnvFe2-vSiO4), dicalcium silicate (Ca2SiO4), pseudowollastonite and wollastonite (CaSiO3), and tridymite (SiO2). The phase equilibria involving the liquid phase and the solid solutions-have also been measured.

Non-equilibrium energy and momentum accommodation coefficients of Ar atoms scattered from Ni(001) in the thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study energy and momentum transfer of low-energy Ar atoms scattered from the Ni(001) surface. The investigation concentrates on the dependence of these processes on incident energy, angles of incidence and surface temperature. Energy transfer exhibits a strong dependence on the surface temperature, at incident energies below 500 meV, and incident angles close to specular incidence. Above 500 meV, the surface temperature dependence vanishes, and a limiting value in the amount of energy transferred to the surface is attained. Momentum exchange is investigated in terms of tangential and normal components. Both components exhibit a weak surface temperature dependence, but they have opposite behaviours at all incidence angles. In each component, momentum can be lost or gained following the interaction with the surface. (C) 1997 Elsevier Science B.V.