BioOne.org will be down briefly for maintenance on 14 May 2025 between 18:00-22:00 Pacific Time US. We apologize for any inconvenience.
Registered users receive a variety of benefits including the ability to customize email alerts, create favorite journals list, and save searches.
Please note that a BioOne web account does not automatically grant access to full-text content. An institutional or society member subscription is required to view non-Open Access content.
Contact helpdesk@bioone.org with any questions.
Common principal components (CPC) analysis is a new tool for the comparison of phenotypic and genetic variance-covariance matrices. CPC was developed as a method of data summarization, but frequently biologists would like to use the method to detect analogous patterns of trait correlation in multiple populations or species. To investigate the properties of CPC, we simulated data that reflect a set of causal factors. The CPC method performs as expected from a statistical point of view, but often gives results that are contrary to biological intuition. In general, CPC tends to underestimate the degree of structure that matrices share. Differences of trait variances and covariances due to a difference in a single causal factor in two otherwise identically structured datasets often cause CPC to declare the two datasets unrelated. Conversely, CPC could identify datasets as having the same structure when causal factors are different. Reordering of vectors before analysis can aid in the detection of patterns. We urge caution in the biological interpretation of CPC analysis results.
A hypothesis has been advanced recently predicting that, in evolution, as higher-level entities arise from associations of lower-level organisms, and as these entities acquire the ability to feed, reproduce, defend themselves, and so on, the lower-level organisms will tend to lose much of their internal complexity (McShea 2001a). In other words, in hierarchical transitions, there is a drain on numbers of part types at the lower level. One possible rationale is that the transfer of functional demands to the higher level renders many part types at the lower level useless, and thus their loss in evolution is favored by selection for economy. Here, a test is conducted at the cell level, comparing numbers of part types in free-living eukaryotic cells (protists) and the cells of metazoans and land plants. Differences are significant and consistent with the hypothesis, suggesting that tests at other hierarchical levels may be worthwhile.
As the ultimate source of genetic variation, spontaneous mutation is essential to evolutionary change. Theoretical studies over several decades have revealed the dependence of evolutionary consequences of mutation on specific mutational properties, including genomic mutation rates, U, and the effects of newly arising mutations on individual fitness, s. The recent resurgence of empirical effort to infer these properties for diverse organisms has not achieved consensus. Estimates, which have been obtained by methods that assume mutations are unidirectional in their effects on fitness, are imprecise. Both because a general approach must allow for occurrence of fitness-enhancing mutations, even if these are rare, and because recent evidence demands it, we present a new method for inferring mutational parameters. For the distribution of mutational effects, we retain Keightley's assumption of the gamma distribution, to take advantage of the flexibility of its shape. Because the conventional gamma is one sided, restricting it to unidirectional effects, we include an additional parameter, ρ, as an amount it is displaced from zero. Estimation is accomplished by Markov chain Monte Carlo maximum likelihood. Through a limited set of simulations, we verify the accuracy of this approach. We apply it to analyze data on two reproductive fitness components from a 17-generation mutation-accumulation study of a Columbia accession of Arabidopsis thaliana in which 40 lines sampled in three generations were assayed simultaneously. For these traits, U ≃ 0.1–0.2, with distributions of mutational effects broadly spanning zero, such that roughly half the mutations reduce reproductive fitness. One evolutionary consequence of these results is lower extinction risks of small populations of A. thaliana than expected from the process of mutational meltdown. A comprehensive view of the evolutionary consequences of mutation will depend on quantitatively accounting for fitness-enhancing, as well as fitness-reducing, mutations.
The Winteraceae are traditionally regarded as the least-specialized descendents of the first flowering plants, based largely on their lack of xylem vessels. Since vessels have been viewed as a key innovation for angiosperm diversification, Winteraceae have been portrayed as declining relicts, limited to wet forest habitats where their tracheid-based wood does not impose a significant hydraulic constraints. In contrast, phylogenetic analyses place Winteraceae among angiosperm clades with vessels, indicating that their vesselless wood is derived rather than primitive, whereas extension of the Winteraceae fossil record into the Early Cretaceous suggests a more complex ecological history than has been deduced from their current distribution. However, the selective regime and ecological events underlying the possible loss of vessels in Winteraceae have remained enigmatic. Here we examine the hypothesis that vessels were lost as an adaptation to freezing-prone environments in Winteraceae by measuring the responses of xylem water transport to freezing for a diverse group of Winteraceae taxa as compared to Canella winterana (Canellaceae, a close relative with vessels) and sympatric conifer taxa. We found that mean percent loss of xylem water transport capacity following freeze-thaw varied from 0% to 6% for Winteraceae species from freezing-prone temperate climates and approximately 20% in those taxa from tropical (nonfreezing) climates. Similarly, conifers exhibit almost no decrease in xylem hydraulic conductivity following freezing. In contrast, water transport in Canella stems is nearly 85% blocked after freeze-thaw. Although vessel-bearing wood of Canella possesses considerably greaterhydraulic capacity than Winteraceae, nearly 20% of xylem hydraulic conductance remains, a value that is comparable to the hydraulic capacity of vesselless Winteraceae xylem, if the proportion of hydraulic flow through vessels (modeled as ideal capillaries) is removed. Thus, the evolutionary removal of vessels may not necessarily require a deleterious shift to an ineffective vascular system. By integrating Winteraceae's phylogenetic relationships and fossil history with physiological and ecological observations, we suggest that, as ancestors of modern Winteraceae passed through temperate conditions present in Southern Gondwana during the Early Cretaceous, they were exposed to selective pressures against vessel-possession and returned to a vascular system relying on tracheids. These results suggest that the vesselless condition is advantageous in freezing-prone areas, which is supported by the strong bias in the ecological abundance of Winteraceae to wet temperate and tropical alpine habitats, rather than a retained feature from the first vesselless angiosperms. We believe that vesselless wood plays an important role in the ecological abundance of Winteraceae in Southern Hemisphere temperate environments by enabling the retention of leaves and photosynthesis in the face of frequent freeze-thaw events.
The Pleistocene extinction of the widespread organ-pipe Montastraea coral had measurable morphological and ecological effects on surviving lineages of the Montastraea “annularis” species complex. Extinction of the organ-pipe Montastraea occurred after more than 500,000 years of dominance in the shallow-water reef habitat of Barbados. Extinction resulted in a morphological shift of the columnar Montastraea lineage from thick to thin columns in modern reef environments. Pleistocene colonies of the columnar morphotype sympatric with organ-pipe Montastraea showed greater column widths than those in allopatry. We subjected our data to a number of criteria for interpreting the morphological shift as character release following lifting of competitive pressure after extinction. The morphological differences do not appear to be due either to chance or to physical properties of the marine environment. Differential local extinction and recolonization of four members of the species complex did not occur on Barbados, so that the species coexisted and appear to have coevolved between more than 600,000 and 82,000 years ago. The morphological shift is related to coral growth form and growth rate, and thus reflects the acquisition of a primary resource in corals—light. Character release occurred at the same oceanic Caribbean island (Barbados) where environments have fluctuated with similar variance throughout the period of coexistence. Not only has competition among living members of the Montastraea “annularis” species complex been convincingly demonstrated, but trends in relative abundance among fossil members of the species complex strongly suggest that a competitive hierarchy was operating during their Pleistocene coexistence on Barbados. We also observed an ecological analogue to character release on another Caribbean island, Curaçao. The distribution and abundance of living columnar M. annularis s.s. and massive M. faveolata from the leeward reef crest in Curaçao is greater now than in the Pleistocene, when organ-pipe Montastraea dominated this shallow-water reef habitat. Extinction of the faster growing, shallow-water organ-pipe Montastraea resulted in higher abundance of the columnar Montastraea lineage in shallow-water habitats, where it shifted its morphology to one adapted to high light levels. The species extinction released surviving lineages from a competitive network that had resulted in lower rank abundance in the Pleistocene community and enhanced abundance of both columnar M. annularis s.s. and M. faveolata in modern communities. Full validation of our interpretation of character release must await experiments that demonstrate whether phenotypic differences between populations have a genetic basis. However, we believe the results of this study point to the important, yet heretofore neglected, role that biological interactions have played in the evolution of closely related reef coral species.
Many members of the cnidarian subclass Zoantharia (sea anemones, corals, and their allies) pass through a larval stage with eight complete mesenteries and without posterior musculature. This larva is usually transient, developing into an adult with 12 or more mesenteries. The adults of one family of sea anemones, the Edwardsiidae, bear the larval number and arrangement of mesenteries and lack the pedal disc seen in other sea anemones. The morphology of the Edwardsiidae has been interpreted in a number of ways: (1) the Edwardsiidae are the most basal extant zoantharian, having diverged before the evolution of additional mesenteries and basal musculature; (2) they are relatively advanced sea anemones that have secondarily simplified because they burrow in sand or mud rather than attaching to a hard substrate; or (3) edwardsiids are derived anemones that have retained a juvenile morphology through paedomorphosis. Phylogenetic analyses of small subunit ribosomal gene sequences reveal that the Edwardsiidae are derived zoantharians, nested within sea anemones. None of the proposed explanations fully explain the edwardsiid's body plan; edwardsiid anatomy is a mosaic of retained primitive and derived features. The results of the present study provide insight into zoantharian phylogeny and illustrate how phylogenetic tests can be used to study the evolution of cnidarian body plans.
Although the sex ratios of many groups conform to Fisher's (1930) prediction that parents should invest equally in daughters and sons, a number of taxa are characterized by excesses of one gender. A variety of mechanisms may lead to sex ratio biases, but in organisms that reproduce clonally as well as sexually, gender differences in the rate of cloning could drive the development of sex-ratio biases. In this study, I demonstrate that males of the clonal brittle star Ophiactis savignyi were significantly more likely to divide than females and that the magnitude of this difference was sufficient to explain the consistent and significant excess of males in natural populations. Females were significantly more likely to lose sexual reproductive capabilities following division, and this greater cost associated with division may explain why females are less likely to divide. Gender differences in mortality rates are unlikely to explain the excess of males in this species. Because of their potential influence on the operational sex ratio, gender differences in division rates may have important ecological and evolutionary implications including effects on the direction and strength of selection.
The deleterious effects of inbreeding have long been known, and inbreeding can increase the risk of extinction for local populations in metapopulations. However, other consequences of inbreeding in metapopulations are still not well understood. Here we show the presence of strong inbreeding depression in a rockpool metapopulation of the planktonic freshwater crustacean Daphnia magna, which reproduces by cyclical parthenogenesis. We conducted three experiments in real and artificial rockpools to quantify components of inbreeding depression in the presence and the absence of competition between clonal lines of selfed and outcrossed genotypes. In replicated asexual populations, we recorded strong selection against clones produced by selfing in competition with clones produced by outcrossing. In contrast, inbreeding depression was much weaker in single-clone populations, that is, in the absence of competition between inbred and outbred clones. The finding of a competitive advantage of the outbred genotypes in this metapopulation suggests that if rockpool populations are inbred, hybrid offspring resulting from crosses between immigrants and local genotypes might have a strong selective advantage. This would increase the effective gene flow in the metapopulation. However, the finding of low inbreeding depression in the monoclonal populations suggests that inbred and outbred genotypes might have about equal chances of establishing new populations.
We ask whether the observed mitochondrial DNA (mtDNA) population subdivision of Drosophila simulans is indicative of organismal structure or of specific processes acting on the mitochondrial genome. Factors either intrinsic or extrinsic to the host genome may influence the evolutionary dynamics of mtDNA. Potential intrinsic factors include adaptation of the mitochondrial genome and of nucleomitochondrial gene complexes specific to the local environment. An extrinsic force that has been shown to influence mtDNA evolution in invertebrates is the bacterial endosymbiont Wolbachia. Evidence presented in this study suggests that mtDNA is not a good indicator of organismal subdivision in D. simulans. Furthermore, there is no evidence to suggest that Wolbachia causes any reduction in nuclear gene flow in this species. The observed differentiation in mtDNA is not corroborated by data from NADH: ubiquinone reductase 75kD subunit precursor or the Alcohol dehydrogenase-related loci, from the shape or size of the male genital arch, or from assortative premating behavior. We discuss these results in relation to a mitochondrial genetic species concept and the potential for Wolbachia-induced incompatibility to be a mechanism of speciation in insects. We conclude with an iterated appeal to include phylogenetic and statistical tests of neutrality as a supplement to phylogenetic and population genetic analyses when using mtDNA as an evolutionary marker.
Between sister species of Drosophila, both pre- and postzygotic reproductive isolation commonly appear by the time a Nei's genetic distance of 0.5 is observed. The degree of genetic differentiation present when allopatric populations of the same Drosophila species exhibit incipient reproductive isolation has not been systematically investigated. Here we compare the relationship between genetic differentiation and pre- and postzygotic isolation among allopatric populations of three cactophilic desert Drosophila: D. mettleri, D. nigrospiracula, and D. mojavensis. The range of all three is interrupted by the Gulf of California, while two species, D. mettleri and D. mojavensis, have additional allopatric populations residing on distant Santa Catalina Island, off the coast of southern California. Significant population structure exists within all three species, but only for allopatric populations of D. mojavensis is significant isolation at the prezygotic level observed. The genetic distances for the relevant populations of D. mojavensis were in the range of 0.12, similar to that for D. mettleri whose greatest D = 0.11 was unassociated with any form of isolation. These observations suggest further investigations of Drosophila populations with genetic distances in this range be undertaken to identify any potential patterns in the relationship between degree of genetic differentiation and the appearance of pre- and/or postzygotic isolation.
Although multiple mating most likely increases mortality risk for social insect queens and lowers the kin benefits for nonreproductive workers, a significant proportion of hymenopteran queens mate with several males. It has been suggested that queens may mate multiply as a means to manipulate sex ratios to their advantage. Multiple paternity reduces the extreme relatedness value of females for workers, selecting for workers to invest more in males. In populations with female-biased sex ratios, queens heading such male-producing colonies would achieve a higher fitness. We tested this hypothesis in a Swiss and a Swedish population of the ant Lasius niger. There was substantial and consistent variation in queen mating frequency and colony sex allocation within and among populations, but no evidence that workers regulated sex allocation in response to queen mating frequency; the investment in females did not differ among paternity classes. Moreover, population-mean sex ratios were consistently less female biased than expected under worker control and were close to the queen optimum. Queens therefore had no incentive to manipulate sex ratios because their fitness did not depend on the sex ratio of their colony. Thus, we found no evidence that the sex-ratio manipulation theory can explain the evolution and maintenance of multiple mating in L. niger.
The recently developed geometric morphometrics methods represent an important contribution of statistics and geometry to the study of biological shapes. We propose simple protocols using shape distances that incorporate geometric techniques into linear quantitative genetic models that should provide insights into the contribution of genetics to shape variation in organisms. The geometric approaches use Procrustes distances in a curved shape space and distances in tangent spaces within and among families to estimate shape heritability. We illustrate the protocols with an example of wing shape variation in the honeybee, Apis mellifera. The heritability of overall shape variation was small, but some localized components depicting shape changes on distal wing regions showed medium to large heritabilities. The genetic variance-covariance matrix of the geometric shape variables was significantly correlated with the phenotypic shape variance-covariance matrix. A comparison of the results of geometric methods with the traditional multivariate analysis of interlandmark distances indicated that even with a larger dimensionality, the interlandmark distances were not as rich in shape information as the landmark coordinates. Quantitative genetics studies of shape should greatly benefit from the application of geometric methods.
Mitochondrial DNA and allozyme variation was examined in populations of two Neotropical butterflies, Heliconius charithonia and Dryas iulia. On the mainland, both species showed evidence of considerable gene flow over huge distances. The island populations, however, revealed significant genetic divergence across some, but not all, ocean passages. Despite the phylogenetic relatedness and broadly similar ecologies of these two butterflies, their intraspecific biogeography clearly differed. Phylogenetic analyses of mitochondrial DNA sequences revealed that populations of D. iulia north of St. Vincent are monophyletic and were probably derived from South America. By contrast, the Jamaican subspecies of H. charithonia rendered West Indian H. charithonia polyphyletic with respect to the mainland populations; thus, H. charithonia seems to have colonized the Greater Antilles on at least two separate occasions from Central America. Colonization velocity does not correlate with subsequent levels of gene flow in either species. Even where range expansion seems to have been instantaneous on a geological timescale, significant allele frequency differences at allozyme loci demonstrate that gene flow is severely curtailed across narrow ocean passages. Stochastic extinction, rapid (re)colonization, but low gene flow probably explain why, in the same species, some islands support genetically distinct and nonexpanding populations, while nearby a single lineage is distributed across several islands. Despite the differences, some common biogeographic patterns were evident between these butterflies and other West Indian taxa; such congruence suggests that intraspecific evolution in the West Indies has been somewhat constrained by earth history events, such as changes in sea level.
In many insect systems, males donate nuptial gifts to insure an effective copulation or as a form of paternal investment. However, if gift magnitude is both body size-limited and positively related to fitness, then the opportunity exists for the gift to promote the evolution of large male size. In the striped ground cricket, Allonemobius socius, males transfer a body size-limited, somatic nuptial gift that is comprised primarily of hemolymph. To address the implications of this gift on male size evolution, we quantified the intensity and direction of natural (fecundity) and sexual (mating success) selection over multiple generations. We found that male size was under strong positive sexual selection throughout the breeding season. This pattern of selection was similar in successive generations spanning multiple years. Male size was also under strong natural selection, with the largest males siring the most offspring. However, multivariate selection gradients indicated that gift size, and not male size, was the best predictor of female fecundity. In other words, direct fecundity selection for larger gifts placed indirect positive selection on male body size, supporting the hypothesis that nuptial gifts can influence the evolution of male body size in this system. Although female size was also under strong selection due to a size related fecundity advantage, it did not exceed selection on male size. The implications of these results with regard to the maintenance of the female-biased size dimorphic system are discussed.
The evolution of species-specific mate recognition signals is of particular interest within speciose monophyletic groups with restricted distributions (known as “species flocks”). However, the explosive nature of speciation in these clades makes difficult the reconstruction of their phylogenetic history. Here we describe a species flock of riverine mormyrid fishes from west-central Africa in which electric signals may play a role in the reproductive isolation of sympatric species. In our recent field collections, totaling more than 1400 specimens from many localities, we recognize 38 forms that are distinct in their morphologies and electric organ discharge (EOD) characteristics. Of these 38, only four clearly correspond to described species. Here we treat these forms as operational taxonomic units (OTUs) in a phylogenetic analysis of cytochrome b sequence data from a sample of 86 specimens. We examined support in the molecular data for the monophyly of these 38 OTUs considered together, the monophyly of each phenotypically delimited OTU considered individually, and for relationships among OTUs congruent with those inferred from the distribution of morphological and EOD character states. Trees obtained by both maximum-parsimony and maximum-likelihood analyses, rooted with sequence data from outgroup taxa, provide evidence for the monophyly of these 38 OTUs with respect to other mormyrid fishes. The small genetic distances between many distinct forms suggest their recent divergence. However, in many instances the cytochrome b tree topology fails to support the monophyly of individual OTUs and close relationships between OTUs that are similar in morphology and EOD characteristics. In other cases, individuals from distinct OTUs share identical or nearly identical haplotypes. Close examination of these cases suggests that unnatural OTU definition is not the sole cause of this pattern, and we infer an incongruence between the mitochondrial gene tree and the organismal phylogeny caused by incomplete mitochondrial lineage sorting and/or introgression across forms. The apparently rapid diversification in this clade of riverine electric fishes and the problems associated with recovering a meaningful species-level phylogeny from mitochondrial data parallel findings in other species flocks. Selection on EOD waveforms as mate recognition signals may be involved in the radiation of these fishes. This is the first description of a freshwater fish species flock from a riverine, as opposed to a lacustrine, environment.
Phenotypic plasticity provides means for adapting to environmental unpredictability. In terms of accelerated development in the face of pond-drying risk, phenotypic plasticity has been demonstrated in many amphibian species, but two issues of evolutionary interest remain unexplored. First, the heritable basis of plastic responses is poorly established. Second, it is not known whether interpopulational differences in capacity to respond to pond-drying risk exist, although such differences, when matched with differences in desiccation risk would provide strong evidence for local adaptation. We investigated sources of within- and among-population variation in plastic responses to simulated pond-drying risk (three desiccation treatments) in two Rana temporaria populations originating from contrasting environments: (1) high desiccation risk with weak seasonal time constraint (southern population); and (2) low desiccation risk with severe seasonal time constraint (northern population). The larvae originating from the environment with high desiccation risk responded adaptively to the fast decreasing water treatment by accelerating their development and metamorphosing earlier, but this was not the case in the larvae originating from the environment with low desiccation risk. In both populations, metamorphic size was smaller in the high-desiccation-risk treatment, but the effect was larger in the southern population. Significant additive genetic variation in development rate was found in the northern and was nearly significant in the southern population, but there was no evidence for genetic variation in plasticity for development rates in either of the populations. No genetic variation for plasticity was found either in size at metamorphosis or growth rate. All metamorphic traits were heritable, and additive genetic variances were generally somewhat higher in the southern population, although significantly so in only one trait. Dominance variances were also significant in three of four traits, but the populations did not differ. Maternal effects in metamorphic traits were generally weak in both populations. Within-environment phenotypic correlations between larval period and metamorphic size were positive and genetic correlations negative in both populations. These results suggest that adaptive phenotypic plasticity is not a species-specific fixed trait, but evolution of interpopulational differences in plastic responses are possible, although heritability of plasticity appears to be low. The lack of adaptive response to desiccation risk in northern larvae is consistent with the interpretation that selection imposed by shorter growing season has favored rapid development in north (∼8% faster development in north as compared to south) or a minimum metamorphic size at the expense of phenotypic plasticity.
Miniaturization has evolved numerous times and reached impressive extremes in the Anura. I compared the skeletons of miniature frog species to those of closely related larger species to assess patterns of morphological change, sampling 129 species from 12 families. Two types of morphological data were examined: (1) qualitative data on bone presence and absence; and (2) thin-plate spline morphometric descriptions of skull structure and bone shape. Phylogenetic comparative methods were used to address the shared history of species. Miniature anurans were more likely to lose skull bones and phalangeal elements of the limbs. Their skulls also showed consistent differences compared to those of their larger relatives, including relatively larger braincases and sensory capsules, verticalization of lateral elements, rostral displacement of the jaw joint, and reduction of some skull elements. These features are explained by functional constraints and by paedomorphosis. Variation among lineages in the morphological response to miniaturization was also explored. Certain lineages appear to be unusually resistant to the morphological trends that characterize miniature frogs as a whole. This study represents the first large-scale examination of morphology and miniaturization across a major, diverse group of organisms conducted in a phylogenetic framework and with statistical rigor.
The genetic basis of fluctuating asymmetry (FA), or nondirectional variation in the subtle differences between left and right sides of bilateral characters, continues to be of considerable theoretical interest. FA generally has been thought to arise from random noise during development and therefore to have a largely or entirely environmental origin. Whereas additive genetic variation for FA generally has been small and often insignificant, a number of investigators have hypothesized that interactions between loci, or epistasis, significantly influence FA. We tested this hypothesis by conducting a whole-genome scan to detect any epistasis in FA of centroid size in the mandibles of more than 400 mice from an F2 intercross population formed from crossing the Large (LG/J) and Small (SM/J) inbred strains. Genotypic deviations were imputed at each site 2 cM apart on all 19 autosomes, and these and centroid size asymmetry values were used in canonical correlation analyses for each of the 171 possible pairs of 19 autosomes to identify the most probable sites for epistasis. Epistasis for centroid size asymmetry was abundant, occurring far more often than was expected by chance alone (there were 30 separate instances of epistasis at the 0.001 significance level, when only two were expected by chance alone). The contributions of epistasis from 30 pairwise combinations of loci tended to suppress the additive and dominance genetic variance, but greatly increased the epistatic genetic variance for FA in centroid size given the intermediate allele frequencies of an F2 intercross population.
It has been assumed, based on theoretical studies, that lethals with the level of dominance estimated from experimental studies would have an allele frequency that is virtually independent of effective population size. However, here it is shown numerically that the expected frequency of lethals with low levels of dominance is also dependent on finite population size, although not as much as completely recessive lethals. This finding is significant in determining the standing level of inbreeding depression and the consequent potential for the evolution of self-fertilization. In addition, the architecture of genetic variation influencing inbreeding depression in populations with a history of small size may be of important consequence in endangered species. Finally, it is shown that the loss of lethal genetic variation often occurs much more quickly than the regeneration of lethal variation by mutation. This asymmetry may result in a lower standing genetic variation for inbreeding depression than expected from mutation rates and contemporary population size data.
This article is only available to subscribers. It is not available for individual sale.
Access to the requested content is limited to institutions that have
purchased or subscribe to this BioOne eBook Collection. You are receiving
this notice because your organization may not have this eBook access.*
*Shibboleth/Open Athens users-please
sign in
to access your institution's subscriptions.
Additional information about institution subscriptions can be foundhere