Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

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.

Local adaptation and species segregation in two mussel (Mytilus edulis X M. trossulus) hybrid zones

Few marine hybrid zones have been studied extensively, the major exception being the hybrid zone between the mussels Mytilus edulis and M. galloprovincialis in southwestern Europe. Here, we focus on two less studied hybrid zones that also involve Mytilus spp.; M. edulis and M. trossulus are sympatric and hybridize on both western and eastern coasts of the Atlantic Ocean. We review the dynamics of hybridization in these two hybrid zones and evaluate the role of local adaptation for maintaining species boundaries. In Scandinavia, hybridization and gene introgression is so extensive that no individuals with pure M. trossulus genotypes have been found. However, M. trossulus alleles are maintained at high frequencies in the extremely low salinity Baltic Sea for some allozyme genes. A synthesis of reciprocal transplantation experiments between different salinity regimes shows that unlinked Gpi and Pgm alleles change frequency following transplantation, such that post-transplantation allelic composition resembles native populations found in the same salinity. These experiments provide strong evidence for salinity adaptation at Gpi and Pgm (or genes linked to them). In the Canadian Maritimes, pure M. edulis and M. trossulus individuals are abundant, and limited data suggest that M. edulis predominates in low salinity and sheltered conditions, whereas M. trossulus are more abundant on the wave-exposed open coasts. We suggest that these conflicting patterns of species segregation are, in part, caused by local adaptation of Scandinavian M. trossulus to the extremely low salinity Baltic Sea environment.

Differential patterns of Male and Female mtDNA exchange across the Atlantic Ocean in the blue mussel, Mytilus edulis

Comparisons among loci with differing modes of inheritance can reveal unexpected aspects of population history. We employ a multilocus approach to ask whether two types of independently assorting mitochondrial DNAs (maternally and paternally inherited: F- and M-mtDNA) and a nuclear locus (ITS) yield concordant estimates of gene flow and population divergence. The blue mussel, Mytilus edulis, is distributed on both North American and European coastlines and these populations are separated by the waters of the Atlantic Ocean. Gene flow across the Atlantic Ocean differs among loci, with F-mtDNA and ITS showing an imprint of some genetic interchange and M-mtDNA showing no evidence for gene flow. Gene flow of F-mtDNA and ITS causes trans-Atlantic population divergence times to be greatly underestimated for these loci, although a single trans-Atlantic population divergence time (1.2 MYA) can be accommodated by considering all three loci in combination in a coalescent framework. The apparent lack of gene flow for M-mtDNA is not readily explained by different dispersal capacities of male and female mussels. A genetic barrier to M-mtDNA exchange between North American and European mussel populations is likely to explain the observed pattern, perhaps associated with the double uniparental system of mitochondrial DNA inheritance.

Ecological speciation

Ecological processes are central to the formation of new species when barriers to gene flow (reproductive isolation) evolve between populations as a result of ecologically-based divergent selection. Although laboratory and field studies provide evidence that 'ecological speciation' can occur, our understanding of the details of the process is incomplete. Here we review ecological speciation by considering its constituent components: an ecological source of divergent selection, a form of reproductive isolation, and a genetic mechanism linking the two. Sources of divergent selection include differences in environment or niche, certain forms of sexual selection, and the ecological interaction of populations. We explore the evidence for the contribution of each to ecological speciation. Forms of reproductive isolation are diverse and we discuss the likelihood that each may be involved in ecological speciation. Divergent selection on genes affecting ecological traits can be transmitted directly (via pleiotropy) or indirectly (via linkage disequilibrium) to genes causing reproductive isolation and we explore the consequences of both. Along with these components, we also discuss the geography and the genetic basis of ecological speciation. Throughout, we provide examples from nature, critically evaluate their quality, and highlight areas where more work is required.

Long-term flow rates and biomat zone hydrology in soil columns receiving septic tank effluent

Soil absorption systems (SAS) are used commonly to treat and disperse septic tank effluent (STE). SAS can hydraulically fail as a result of the low permeable biomat zone that develops on the infiltrative surface. The objectives of this experiment were to compare the hydraulic properties of biomats grown in soils of different textures, to investigate the long-term acceptance rates (LTAR) from prolonged application of STE, and to assess if soils were of major importance in determining LTAR. The STE was applied to repacked sand, Oxisol and Vertisol soil columns over a period of 16 months, at equivalent hydraulic loading rates of 50, 35 and 8 L/m(2)/d, respectively Infiltration rates, soil matric potentials, and biomat hydraulic properties were measured either directly from the soil columns or calculated using established soil physics theory. Biomats 1 to 2 cm thick developed in all soils columns with hydraulic resistances of 27 to 39 d. These biomats reduced a 4 order of magnitude variation in saturated hydraulic conductivity (K.) between the soils to a one order of magnitude variation in LTAR. A relationship between biomat resistance and organic loading rate was observed in all soils. Saturated hydraulic conductivity influenced the rate and extent of biomat development. However, once the biomat was established, the LTAR was governed by the resistance of the biomat and the sub-biomat soil unsaturated flow regime induced by the biomat. Results show that whilst initial soil K. is likely to be important in the establishment of the biomat zone in a trench, LTAR is determined by the biomat resistance and the unsaturated soil hydraulic conductivity, not the K, of a soil. The results call into question the commonly used approach of basing the LTAR, and ultimately trench length in SAS, on the initial K, of soils. (c) 2006 Elsevier Ltd. All rights reserved.

The role of mobile belts for the longevity of deep cratonic lithosphere: The crumple zone model

Many Archean cratons are surrounded by Proterozoic mobile belts that have experienced episodes of tectonic re-activation over their lifetimes. This suggests that mobile belt lithosphere may be associated with long lived, inherited weakness. It is proposed that the proximity of this weakness can increase the longevity of deep Archean lithosphere by buffering Archean cratons from mantle derived stresses. The physical plausibility of this idea is explored through numerical simulations of mantle convection that include continents and allow for material rheologies that model the combined brittle and ductile behavior of the lithosphere. Within the simulations, the longevity of deep cratonic lithosphere does increase if it is buffered by mobile belts that can fail at relatively low stress levels.

In familial hyperaldosteronism type I, hybrid gene-induced aldosterone production dominates that induced by wild-type genes

We compared the aldosterone-producing potency of the angiotensin II-sensitive wild-type aldosterone synthase genes and the ACTH-sensitive hybrid 11 beta-hydroxylase/aldosterone synthase gene by examining aldosterone, PRA, and cortisol day-curves (2-hourly levels over 24 h) in patients with familial hyperaldosteronism type I, before and during long-term (0.8-13.5 yr) glucocorticoid treatment. In 8 untreated patients, PRA levels were usually suppressed, and aldosterone correlated strongly with cortisol (r = 0.69-0.99). Fourteen studies were performed on 10 patients receiving glucocorticoid treatment that corrected hypertension, hypokalemia, and PRA suppression in all. ACTH was markedly and continuously suppressed in 6 studies, 3 of which demonstrated strong correlations between aldosterone and PRA (r = 0.77-0.92), ACTH was only partially suppressed in the remaining 8 studies; aldosterone correlated strongly: 1) with cortisol alone in 5 (r = 0.71-0.98); 2) with cortisol (r = 0.90) and PRA (r = 0.74) in one; 3) with PRA only in one (r = 0.80); and 4) with neither PRA nor cortisol in one. Unless ACTH is markedly and continuously suppressed, aldosterone is more responsive to ACTH than to renin/angiotensin II, despite the latter being unsuppressed. This is consistent with the hybrid gene being more powerfully expressed than the wild-type aldosterone synthase genes in familial hyperaldosteronism type I.

Geographic parthenogenesis in the Australian arid zone: I. A climatic analysis of the Heteronotia binoei complex (Gekkonidae)

Patterns of geographic parthenogenesis can provide insight into the ecological implications of the transition from sexual to parthenogenetic reproduction. We analysed quantitatively the environmental niches occupied by sexual and parthenogenetic geckos of the Heteronotia binoei complex in the Australian and zone. This complex consists of two independently derived maternal lineages of hybrid parthenogens, which, in turn, include two different triploid races that resulted from reciprocal backcrossing with the parental sexual taxa. The sexual progenitors are still extant and occupy very distinct environmental niches. The triploid parthenogenetic races are biased in their environmental niche towards those of the sexual races for which their genomes are biased and this dosage effect is apparent in both maternal lineages. Thus triploidy may have benefited the parthenogens through partial recovery of the parental niches. Although the parthenogens have a broader geographic distribution than their sexual progenitors, their environmental niche is narrower and biased towards one of the sexual races. In keeping with general patterns of geographic parthenogenesis. parthenogenetic H. binoei occupy a harsher environment than the sexual forms. occurring in regions of persistently low rainfall. Bioclimatic modelling suggests patterns of rainfall are important in limiting the distribution of sexual and parthenogenetic taxa. and extrapolation from the current bioclimatic profiles indicates potential for further eastward range expansion by the parthenogens.

Reproductive isolation of a new hybrid species, Senecio eboracensis Abbott & Lowe (Asteraceae)

The nature and extent of reproductive isolation was examined between a new self-compatible hybrid species Senecio eboracensis (2n = 40) and its parents, self-incompatible S. squalidus (2n = 20) and self-compatible S. vulgaris (2n = 40). The triploid F-1 of S. eboracensis x S. squalidus exhibited very low seed set ((x) over bar = 0.63%), and F-2 and F-3 progeny were able to recover nominal levels of fertility ((x) over bar = 23.9 and 9.7%), while F-1 and F-2 offspring of S. eboracensis x S. vulgaris showed reduced seed set ((x) over bar = 63.8 and 58.8%). In both cases, evidence from previous work indicates that reduced fertility is associated with meiotic chromosome mispairing, and is a likely consequence of recombining both parental genomes within this new taxon. No hybrid offspring between S. eboracensis and S. squalidus were found in the wild, and only one such hybrid was recorded among 769 progeny produced by S. eboracensis surrounded by S. squalidus on an experimental plot. Natural crossing between S. eboracensis and S. vulgaris was recorded to be very low (between 0 and 1.46%) in the wild, but rose to 18.3% when individuals of S. eboracensis were surrounded by plants of S. vulgaris. It was concluded that strong breeding barriers exist between the new hybrid species and its two parents. Prezygotic isolation between S. eboracensis and S. vulgaris is likely to be largely due to both species reproducing by predominant self-fertilisation. However, differences recorded for germination, seedling survival, time of flowering and characters associated with pollinator attraction, plus significant clumping of juvenile and adult conspecifics in the wild, probably also contribute to reproductive isolation and ecological differentiation.

Evolution of the regular zone of histone H1 in Fabaceae plants

An analysis of the historic H1 subtype, H1-1, in eight legumes belonging to four genera of the tribe Vicieae (Pisum, Lathyrus, Lens, and Vicia), revealed an extended region consisting of the tandemly repeated AKPAAK motifs. We named this region the Regular zone (RZ). The AKPAAK motifs are organized into two blocks separated by a short (two or six amino acids) intervening sequence (IS). The distal block contains six AKPAAK motifs, while the number of repeats in the proximal block varies from six in V. faba to seven in the other species. In V. hirsuta, the first two repeated units of the proximal block are octapeptides AKAKPAAK. The apparent rate of synonymous substitutions in the blocks of RZ is much higher than in the rest of the gene. This can be explained by repeat shuffling within each block. In the C-domain of the orthologous H1 subtype froth Medicago truncatula (tribe Trifolieae), a region corresponding to the RZ of Vicieae species was found. It also consists of two blocks of AKPAAK motifs (four and three repeats in the proximal and distal blocks, respectively). These blocks are separated by a 20-amino acid IS. The first 20 amino acids of the Medicago RZ are not part of AKPAAK repeats. We hypothesise that the RZ has most probably evolved as a result of an expansion of AKPAAK repeats from two separate sites in the C-domain. This process started tens of millions of years ago and was most likely directed by positive selection.

The estrogen receptor 1 G594A polymorphism is associated with migraine susceptibility in two independent case/control groups

Migraine is a painful and debilitating disorder with a significant genetic component. Steroid hormones, in particular estrogen, have long been considered to play a role in migraine, as variations in hormone levels are associated with migraine onset in many sufferers of the disorder. Steroid hormones mediate their activity via hormone receptors, which have a wide tissue distribution. Estrogen receptors have been localized to the brain in regions considered to be involved in migraine pathogenesis. Hence it is possible that genetic variation in the estrogen receptor gene may play a role in migraine susceptibility. This study thus examined the estrogen receptor 1 (ESRalpha) gene for a potential role in migraine pathogenesis and susceptibility. A population-based cohort of 224 migraine sufferers and 224 matched controls were genotyped for the G594A polymorphism located in exon 8 of the ESR1 gene. Statistical analysis indicated a significant difference between migraineurs and non-migraineurs in both the allele frequencies (P=0.003) and genotype distributions (P=0.008) in this sample. An independent follow-up study was then undertaken using this marker in an additional population-based cohort of 260 migraine sufferers and 260 matched controls. This resulted in a significant association between the two groups with regard to allele frequencies (P=8x10(-6)) and genotype distributions (P=4x10(-5)). Our findings support the hypothesis that genetic variation in hormone receptors, in particular the ESR1 gene, may play a role in migraine.

Field observations of instantaneous water slopes and horizontal pressure gradients in the swash-zone

Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within +/- 0.001, or 1/20th -1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05-0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly Suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., delta u/delta t > 0 similar to delta p/delta x < 0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.

Long wave generation by the shoaling and breaking of transient wave groups on a beach

This paper presents new laboratory data on the generation of long waves by the shoaling and breaking of transient-focused short-wave groups. Direct offshore radiation of long waves from the breakpoint is shown experimentally for the first time. High spatial resolution enables identification of the relationship between the spatial gradients of the short-wave envelope and the long-wave surface. This relationship is consistent with radiation stress theory even well inside the surf zone and appears as a result of the strong nonlinear forcing associated with the transient group. In shallow water, the change in depth across the group leads to asymmetry in the forcing which generates significant dynamic setup in front of the group during shoaling. Strong amplification of the incident dynamic setup occurs after short-wave breaking. The data show the radiation of a transient long wave dominated by a pulse of positive elevation, preceded and followed by weaker trailing waves with negative elevation. The instantaneous cross-shore structure of the long wave shows the mechanics of the reflection process and the formation of a transient node in the inner surf zone. The wave run-up and relative amplitude of the radiated and incident long waves suggests significant modification of the incident bound wave in the inner surf zone and, the dominance of long waves generated by the breaking process. It is proposed that these conditions occur when the primary short waves and bound wave are not shallow water waves at the breakpoint. A simple criterion is given to determine these conditions, which generally occur for the important case of storm waves.