Evolution of the 28S ribosomal RNA gene in anurans: regions of variability and their phylogenetic implications

Fifteen restriction sites were mapped to the 28S ribosomal RNA gene of individuals representing 54 species of frogs, two species of salamanders, a caecilian, and a lungfish. Eight of these sites were present in all species examined, and two were found in all but one species. Alignment of these conserved restriction sites revealed, among anuran 28S rRNA genes, five regions of major length variation that correspond to four of 12 previously identified divergent domains of this gene. One of the divergent domains (DD8) consists of two regions of length variation separated by a short segment that is conserved at least throughout tetrapods. Most of the insertions, deletions, and restriction-site variations identified in the 28S gene will require sequence-level analysis for a detailed reconstruction of their history. However, an insertion in DD9 that is coextensive with frogs in the suborder Neobatrachia, a BstEII site that is limited to representatives of two leptodactylid subfamilies, and a deletion in DD10 that is found only in three ranoid genera are probably synapomorphies.      

A Generalized Approach for Estimating Effective Population Size from Temporal Changes in Allele Frequency

The temporal method for estimating effective population size (Ne) from the standardized variance in allele frequency change (F) is presented in a generalized form. Whereas previous treatments of this method have adopted rather limiting assumptions, the present analysis shows that the temporal method is generally applicable to a wide variety of organisms. Use of a revised model of gene sampling permits a more generalized interpretation of Ne than that used by some other authors studying this method. It is shown that two sampling plans (individuals for genetic analysis taken before or after reproduction) whose differences have been stressed by previous authors can be treated in a uniform way. Computer simulations using a wide variety of initial conditions show that different formulas for computing F have much less effect on Ne than do sample size (S), number of generations between samples (t), or the number of loci studied (L). Simulation results also indicate that (1) bias of F is small unless alleles with very low frequency are used; (2) precision is typically increased by about the same amount with a doubling of S, t, or L; (3) confidence intervals for Ne computed using a chi 2 approximation are accurate and unbiased under most conditions; (4) the temporal method is best suited for use with organisms having high juvenile mortality and, perhaps, a limited effective population size.     

Characterization of cross-bridge elasticity and kinetics of cross-bridge cycling during force development in single smooth muscle cells

Force development in smooth muscle, as in skeletal muscle, is believed to reflect recruitment of force-generating myosin cross-bridges. However, little is known about the events underlying cross-bridge recruitment as the muscle cell approaches peak isometric force and then enters a period of tension maintenance. In the present studies on single smooth muscle cells isolated from the toad (Bufo marinus) stomach muscularis, active muscle stiffness, calculated from the force response to small sinusoidal length changes (0.5% cell length, 250 Hz), was utilized to estimate the relative number of attached cross-bridges. By comparing stiffness during initial force development to stiffness during force redevelopment immediately after a quick release imposed at peak force, we propose that the instantaneous active stiffness of the cell reflects both a linearly elastic cross-bridge element having 1.5 times the compliance of the cross-bridge in frog skeletal muscle and a series elastic component having an exponential length-force relationship. At the onset of force development, the ratio of stiffness to force was 2.5 times greater than at peak isometric force. These data suggest that, upon activation, cross-bridges attach in at least two states (i.e., low-force-producing and high-force-producing) and redistribute to a steady state distribution at peak isometric force. The possibility that the cross-bridge cycling rate was modulated with time was also investigated by analyzing the time course of tension recovery to small, rapid step length changes (0.5% cell length in 2.5 ms) imposed during initial force development, at peak force, and after 15 s of tension maintenance. The rate of tension recovery slowed continuously throughout force development following activation and slowed further as force was maintained. Our results suggest that the kinetics of force production in smooth muscle may involve a redistribution of cross-bridge populations between two attached states and that the average cycling rate of these cross-bridges becomes slower with time during contraction.     

Estimation of Allele Frequencies at Isoloci

In some polyploid animals and plants, pairs of duplicated loci occur that share alleles encoding proteins with identical electrophoretic mobilities. Except in cases where these ``isoloci' are known to be inherited tetrasomically, individual genotypes cannot be determined unambiguously, and there is no direct way to assign observed variation to a particular locus of the pair. For a pair of diallelic isoloci, nine genotypes are possible but only five phenotypes can be identified, corresponding to individuals with 0-4 doses of the variant allele. A maximum likelihood (ML) approach is used here to identify the set of allele frequencies (p, q) at the individual gene loci with the highest probability of producing the observed phenotypic distribution. A likelihood ratio test is used to generate the asymmetrical confidence intervals around ML estimates. Simulations indicate that the standard error of p is typically about twice the binomial sampling error associated with single locus allele frequency estimates. ML estimates can be used in standard indices of genetic diversity and differentiation and in goodness-of-fit tests of genetic hypotheses. The noncentral χ(2) distribution is used to evaluate the power of a test of apparent heterozygote deficiency that results from attributing all variation to one locus when both loci are polymorphic.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Evidence for "twisted plane" undulations in golden hamster sperm tails

Motile spermatozoa from the golden hamster have been arrested by rapid freezing and then fixed with glutaraldehyde at low temperature after substitution with ethylene glycol. As far as can be judged, the flagellar waveforms thus stabilized are similar to those seen in living sperm; in contrast, fixation in glutaraldehyde, without prior freezing, induces agonal changes in flagellar conformation. The characteristics waveform after freeze substitution contains three bends. Approx. half of these flagella are entirely planar. The rest are three dimensional, with the third bend displaced in a regular way from the plane containing the proximal two bends. From the geometry of these flagella, it is concluded that the plane of action of a given bending cycle undergoes a clockwise twist (from a forward viewpoint) as the cycle is succeeded by new bending cycles. This "twisted plane" undulation is quite different from helical movement. The twisting seems to occur abruptly, between cycles, as if each bending cycle has a preferred plane of action. The mechanism underlying the twisting is uncertain. However, on the basis of the angular displacements between the preferred planes, and the findings from electron microscopy, the following idea is presented as a working hypothesis: that, if the most proximal plane of bending is topographically determined by peripheral doublet 1, then successive distal planes of action are influenced predominantly by doublets 2, 3, etc., in clockwise sequence. The merits and weaknesses of this hypothesis are discussed.     

The evolution of microendemism in a reef fish (Hypoplectrus maya)

Marine species tend to have extensive distributions, which are commonly attributed to the dispersal potential provided by planktonic larvae and the rarity of absolute barriers to dispersal in the ocean. Under this paradigm, the occurrence of marine microendemism without geographic isolation in species with planktonic larvae poses a dilemma. The recently described Maya hamlet (Hypoplectrus maya, Serranidae) is exactly such a case, being endemic to a 50‐km segment of the Mesoamerican Barrier Reef System (MBRS). We use whole‐genome analysis to infer the demographic history of the Maya hamlet and contrast it with the sympatric and pan‐Caribbean black (H. nigricans), barred (H. puella) and butter (H. unicolor) hamlets, as well as the allopatric but phenotypically similar blue hamlet (H. gemma). We show that H. maya is indeed a distinct evolutionary lineage, with genomic signatures of inbreeding and a unique demographic history of continuous decrease in effective population size since it diverged from congeners just ~3,000 generations ago. We suggest that this case of microendemism may be driven by the combination of a narrow ecological niche and restrictive oceanographic conditions in the southern MBRS, which is consistent with the occurrence of an unusually high number of marine microendemics in this region. The restricted distribution of the Maya hamlet, its decline in both census and effective population sizes, and the degradation of its habitat place it at risk of extinction. We conclude that the evolution of marine microendemism can be a fast and dynamic process, with extinction possibly occurring before speciation is complete.     

Impact of human activities on the reproduction of Hooded Vultures Necrosyrtes monachus in Burkina Faso

During the last decades, the critically endangered Hooded Vulture Necrosyrtes monachus has strongly declined across its African range. Although direct persecution has been suggested as a major cause of this decline, little is known about the impact of humans on reproductive output in West Africa. We studied the impact of human activities on the reproductive output of Hooded Vultures in the Garango area of Burkina Faso. Twenty and 56 nesting attempts were monitored, respectively, during the breeding season in 2013/14 and 2014/15, to determine reproductive success and identify causes of nest failure. Annual breeding success varied between 0.68 and 0.71 chicks fledged per breeding pair per year and productivity was assessed at 0.57 chicks fledged per territorial pair in 2014/15. The main threats imposed by humans were poaching of eggs, chicks and collection of nest materials, leading to 20% (13 out of 64 breeding attempts) of nest failures over the two years. An additional important reason for nest failure was the pruning and (partial) cutting of nest trees. Despite this high level of human interference, we found that Hooded Vulture nest success increased with proximity to human settlements, probably because breeding vultures benefit from protection by people against persecution and disturbance.     

Production, recovery and immunogenicity of the protective antigen from a recombinant strain of Bacillus anthracis

The protective antigen (PA) is one of the three components of the anthrax toxin. It is a secreted nontoxic protein with a molecular weight of 83 kDa and is the major component of the currently licensed human vaccine for anthrax. Due to limitations found in the existing vaccine formulation, it has been proposed that genetically modified PA may be more effective as a vaccine. The expression and the stability of two recombinant PA (rPA) variants, PA-SNKE-ΔFF-E308D and PA-N657A, were studied. These proteins were expressed in the nonsporogenic avirulent strain BH445. Initial results indicated that PA-SNKE-ΔFF-E308D, which lacks two proteolysis-sensitive sites, is more stable than PA-N657A. Process development was conducted to establish an efficient production and purification process for PA-SNKE-ΔFF-E308D. pH, media composition, growth strategy and protease inhibitors composition were analyzed. The production process chosen was based on batch growth of B. anthracis using tryptone and yeast extract as the only source of carbon, pH control at 7.5, and antifoam 289. Optimal harvest time was 14–18 h after inoculation, and EDTA (5 mM) was added upon harvest for proteolysis control. Recovery of the rPA was performed by expanded-bed adsorption (EBA) on a hydrophobic interaction chromatography (HIC) resin, eliminating the need for centrifugation, microfiltration and diafiltration. The EBA step was followed by ion exchange and gel filtration. rPA yields before and after purification were 130 and 90 mg/l, respectively. The purified rPA, without further treatment, treated with small amounts of formalin or adsorbed on alum, induced, high levels of IgG anti-PA with neutralization activities. Journal of Industrial Microbiology & Biotechnology (2002) 28, 232–238 DOI: 10.1038/sj/jim/7000239 Received 28 August 2001/ Accepted in revised form 20 December 2001     

Effects of Overlapping Generations on Linkage Disequilibrium Estimates of Effective Population Size

Use of single-sample genetic methods to estimate effective population size has skyrocketed in recent years. Although the underlying models assume discrete generations, they are widely applied to age-structured species. We simulated genetic data for 21 iteroparous animal and plant species to evaluate two untested hypotheses regarding performance of the single-sample method based on linkage disequilibrium (LD): (1) estimates based on single-cohort samples reflect the effective number of breeders in one reproductive cycle (N_b), and (2) mixed-age samples reflect the effective size per generation (N_e). We calculated true N_e and N_b, using the model species’ vital rates, and verified these with individual-based simulations. We show that single-cohort samples should be equally influenced by N_b and N_e and confirm this with simulated results: N^b was a linear (r² = 0.98) function of the harmonic mean of N_e and N_b. We provide a quantitative bias correction for raw N^b based on the ratio N_b/N_e, which can be estimated from two or three simple life history traits. Bias-adjusted estimates were within 5% of true N_b for all 21 study species and proved robust when challenged with new data. Mixed-age adult samples produced downwardly biased estimates in all species, which we attribute to a two-locus Wahlund effect (mixture LD) caused by combining parents from different cohorts in a single sample. Results from this study will facilitate interpretation of rapidly accumulating genetic estimates in terms of both N_e (which influences long-term evolutionary processes) and N_b (which is more important for understanding eco-evolutionary dynamics and mating systems).