A genomic island linked to ecotype divergence in Atlantic cod

The genomic architecture underlying ecological divergence and ecological speciation with gene flow is still largely unknown for most organisms. One central question is whether divergence is genome-wide or localized in 'genomic mosaics' during early stages when gene flow is still pronounced. Empirical work has so far been limited, and the relative impacts of gene flow and natural selection on genomic patterns have not been fully explored. Here, we use ecotypes of Atlantic cod to investigate genomic patterns of diversity and population differentiation in a natural system characterized by high gene flow and large effective population sizes, properties which theoretically could restrict divergence in local genomic regions. We identify a genomic region of strong population differentiation, extending over approximately 20 cM, between pairs of migratory and stationary ecotypes examined at two different localities. Furthermore, the region is characterized by markedly reduced levels of genetic diversity in migratory ecotype samples. The results highlight the genomic region, or 'genomic island', as potentially associated with ecological divergence and suggest the involvement of a selective sweep. Finally, we also confirm earlier findings of localized genomic differentiation in three other linkage groups associated with divergence among eastern Atlantic populations. Thus, although the underlying mechanisms are still unknown, the results suggest that 'genomic mosaics' of differentiation may even be found under high levels of gene flow and that marine fishes may provide insightful model systems for studying and identifying initial targets of selection during ecological divergence.

Inbreeding uncovers fundamental differences in the genetic load affecting male and female fertility in a butterfly

Inbreeding depression is most pronounced for traits closely associated with fitness. The traditional explanation is that natural selection eliminates deleterious mutations with additive or dominant effects more effectively than recessive mutations, leading to directional dominance for traits subject to strong directional selection. Here we report the unexpected finding that, in the butterfly Bicyclus anynana, male sterility contributes disproportionately to inbreeding depression for fitness (complete sterility in about half the sons from brother-sister matings), while female fertility is insensitive to inbreeding. The contrast between the sexes for functionally equivalent traits is inconsistent with standard selection arguments, and suggests that trait-specific developmental properties and cryptic selection play crucial roles in shaping genetic architecture. There is evidence that spermatogenesis is less developmentally stable than oogenesis, though the unusually high male fertility load in B. anynana additionally suggests the operation of complex selection maintaining male sterility recessives. Analysis of the precise causes of inbreeding depression will be needed to generate a model that reliably explains variation in directional dominance and reconciles the gap between observed and expected genetic loads carried by populations. This challenging evolutionary puzzle should stimulate work on the occurrence and causes of sex differences in fertility load.

Effects of selection and drift on G matrix evolution in a heterogeneous environment: a multivariate Qst-Fst Test with the freshwater snail Galba truncatula.

Unraveling the effect of selection vs. drift on the evolution of quantitative traits is commonly achieved by one of two methods. Either one contrasts population differentiation estimates for genetic markers and quantitative traits (the Q(st)-F(st) contrast) or multivariate methods are used to study the covariance between sets of traits. In particular, many studies have focused on the genetic variance-covariance matrix (the G matrix). However, both drift and selection can cause changes in G. To understand their joint effects, we recently combined the two methods into a single test (accompanying article by Martin et al.), which we apply here to a network of 16 natural populations of the freshwater snail Galba truncatula. Using this new neutrality test, extended to hierarchical population structures, we studied the multivariate equivalent of the Q(st)-F(st) contrast for several life-history traits of G. truncatula. We found strong evidence of selection acting on multivariate phenotypes. Selection was homogeneous among populations within each habitat and heterogeneous between habitats. We found that the G matrices were relatively stable within each habitat, with proportionality between the among-populations (D) and the within-populations (G) covariance matrices. The effect of habitat heterogeneity is to break this proportionality because of selection for habitat-dependent optima. Individual-based simulations mimicking our empirical system confirmed that these patterns are expected under the selective regime inferred. We show that homogenizing selection can mimic some effect of drift on the G matrix (G and D almost proportional), but that incorporating information from molecular markers (multivariate Q(st)-F(st)) allows disentangling the two effects.

Genetic diversity and differentiation processes in the ploidy series of Olea europaea L.: a multiscale approach from subspecies to insular populations.

Geographical isolation and polyploidization are central concepts in plant evolution. The hierarchical organization of archipelagos in this study provides a framework for testing the evolutionary consequences for polyploid taxa and populations occurring in isolation. Using amplified fragment length polymorphism and simple sequence repeat markers, we determined the genetic diversity and differentiation patterns at three levels of geographical isolation in Olea europaea: mainland-archipelagos, islands within an archipelago, and populations within an island. At the subspecies scale, the hexaploid ssp. maroccana (southwest Morocco) exhibited higher genetic diversity than the insular counterparts. In contrast, the tetraploid ssp. cerasiformis (Madeira) displayed values similar to those obtained for the diploid ssp. guanchica (Canary Islands). Geographical isolation was associated with a high genetic differentiation at this scale. In the Canarian archipelago, the stepping-stone model of differentiation suggested in a previous study was partially supported. Within the western lineage, an east-to-west differentiation pattern was confirmed. Conversely, the easternmost populations were more related to the mainland ssp. europaea than to the western guanchica lineage. Genetic diversity across the Canarian archipelago was significantly correlated with the date of the last volcanic activity in the area/island where each population occurs. At the island scale, this pattern was not confirmed in older islands (Tenerife and Madeira), where populations were genetically homogeneous. In contrast, founder effects resulted in low genetic diversity and marked genetic differentiation among populations of the youngest island, La Palma.

Genes mirror geography within Europe.

Understanding the genetic structure of human populations is of fundamental interest to medical, forensic and anthropological sciences. Advances in high-throughput genotyping technology have markedly improved our understanding of global patterns of human genetic variation and suggest the potential to use large samples to uncover variation among closely spaced populations. Here we characterize genetic variation in a sample of 3,000 European individuals genotyped at over half a million variable DNA sites in the human genome. Despite low average levels of genetic differentiation among Europeans, we find a close correspondence between genetic and geographic distances; indeed, a geographical map of Europe arises naturally as an efficient two-dimensional summary of genetic variation in Europeans. The results emphasize that when mapping the genetic basis of a disease phenotype, spurious associations can arise if genetic structure is not properly accounted for. In addition, the results are relevant to the prospects of genetic ancestry testing; an individual's DNA can be used to infer their geographic origin with surprising accuracy-often to within a few hundred kilometres.

Male mutation bias and possible long-term effects of human activities.

The ability of a population to adapt to changing environments depends critically on the amount and kind of genetic variability it possesses. Mutations are an important source of new genetic variability and may lead to new adaptations, especially if the population size is large. Mutation rates are extremely variable between and within species, and males usually have higher mutation rates as a result of elevated rates of male germ cell division. This male bias affects the overall mutation rate. We examined the factors that influence male mutation bias, and focused on the effects of classical life-history parameters, such as the average age at reproduction and elevated rates of sperm production in response to sexual selection and sperm competition. We argue that human-induced changes in age at reproduction or in sexual selection will affect male mutation biases and hence overall mutation rates. Depending on the effective population size, these changes are likely to influence the long-term persistence of a population.

Heterochromatic banding pattern in two Brazilian populations of Aedes aegypti

We analysed samples of Aedes aegypti from Sao Jose do Rio Preto and Franca (Brazil) by C-banding and Ag-banding staining techniques. C-banding pattern of Ae. aegypti from Sao Jose do Rio Preto examined in metaphase cells differed from Franca. The chromosomes 2, 3 and X showed centromeric C-bands in both populations, but a slightly stained centromeric band in the Y chromosome was observed only in Sao Jose do Rio Preto. In addition, the X chromosome in both populations and the Y chromosome of all individuals from Sao Jose do Rio Preto showed an intercalary band on one of the arms that was absent in Franca. An intercalary, new band, lying on the secondary constriction of chromosome 3 was also present in mosquitoes of both populations. The comparison of the present data with data in the literature for Ae. aegypti from other regions of the world showed that they differ as to the banding pattern of sex chromosomes and the now described intercalary band in chromosome 3. The observations suggested that the heterochromatic regions of all chromosomes are associated to constitute a single C-banded body in interphase cells. Ag-banding technique stained the centromeric regions of all chromosomes (including the Y) and the intercalary C-band region of the X chromosome in both populations. As Ae. aegypti populations are widespread in a great part of the world, the banding pattern variations indicate environmental interactions and may reveal both the chromosome evolutionary patterns in this species and the variations that may interfere with its vector activity.

Variation of genetic expression during development, revealed by esterase patterns in Aedes aegypti (Diptera, Culicidae)

Polyacrylamide gel electrophoresis was used to analyze esterase patterns during development of Aedes aegypti from the cities of Marília and São José do Rio Preto (SJRP), Brazil. The zymograms showed a total of 23 esterase bands, 22 of which were in the specimens from Marília and 19 in those from SJRP. These esterase bands were considered to be the product of 23 alleles distributed tentatively in eight genetic loci. Most of the alleles were developmentally regulated. The larval stage expressed the greatest number of them (19 alleles, from the eight loci, in Marília; and 17 alleles, from seven loci, in SJRP). The pupal stage expressed 10 alleles from seven loci, in both populations, and the adult stage expressed 8 alleles from five and six loci in SJRP and Marília, respectively. Some alleles that were active in every stage were developmentally controlled at the level of expression (amount of product). A single allele was constitutively and highly expressed, in larvae, pupae, and adults, in both populations. Differences in esterase synthesis among stages are probably due to regulatory mechanisms acting in agreement with the requirements of a variable number of processes in which esterases are involved. The larval stage is the most active in developmental processes and shows very intense intake of food and very high mobility. These features may demand increased esterase production at that stage. Comparison of the two populations examined showed (besides the existence of alleles that they do not share) that they exhibit differences in the control of expression of other alleles. Such findings may reflect genetic differences between founders in each population, but the possibility of involvement of the intensive use of insecticides in SJRP is also discussed.

Strong selection against hybrids at a hybrid zone in the Ensatina ring species complex and its evolutionary implications

The analysis of interactions between lineages at varying levels of genetic divergence can provide insights into the process of speciation through the accumulation of incompatible mutations. Ring species, and especially the Ensatina eschscholtzii system exemplify this approach. The plethodontid salamanders E. eschscholtzii xanthoptica and E. eschscholtzii platensis hybridize in the central Sierran foothills of California. We compared the genetic structure across two transects (southern and northern Calaveras Co.), one of which was resampled over 20 years, and examined diagnostic molecular markers (eight allozyme loci and mitochondrial DNA) and a diagnostic quantitative trait (color pattern). Key results across all studies were: (1) cline centers for all markers were coincident and the zones were narrow, with width estimates of 730 m to 2000 m; (2) cline centers at the northern Calaveras transect were coincident between 1981 and 2001, demonstrating repeatability over five generations; (3) there were very few if any putative F1s, but a relatively high number of backcrossed individuals in the central portion of transects: and (4) we found substantial linkage disequilibrium in all three studies and strong heterozygote deficit both in northern Calaveras, in 2001, and southern Calaveras. Both linkage disequilibrium and heterozygote deficit showed maximum values near the center of the zones. Using estimates of cline width and dispersal, we infer strong selection against hybrids. This is sufficient to promote accumulation of differences at loci that are neutral or under divergent selection, but would still allow for introgression of adaptive alleles. The evidence for strong but incomplete isolation across this centrally located contact is consistent with theory suggesting a gradual increase in postzygotic incompatibility between allopatric populations subject to divergent selection and reinforces the value of Ensatina as a system for the study of divergence and speciation at multiple stages. © 2005 The Society for the Study of Evolution. All rights reserved.

Estimation of variance components and prediction of breeding values in rubber tree breeding using the REML/BLUP procedure

The present paper deals with estimation of variance components, prediction of breeding values and selection in a population of rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Müell.-Arg.] from Rio Branco, State of Acre, Brazil. The REML/BLUP (restricted maximum likelihood/best linear unbiased prediction) procedure was applied. For this purpose, 37 rubber tree families were obtained and assessed in a randomized complete block design, with three unbalanced replications. The field trial was carried out at the Experimental Station of UNESP, located in Selvíria, State of Mato Grosso do Sul, Brazil. The quantitative traits evaluated were: girth (G), bark thickness (BT), number of latex vessel rings (NR), and plant height (PH). Given the unbalanced condition of the progeny test, the REML/BLUP procedure was used for estimation. The narrow-sense individual heritability estimates were 0.43 for G, 0.18 for BT, 0.01 for NR, and 0.51 for PH. Two selection strategies were adopted: one short-term (ST - selection intensity of 8.85%) and the other long-term (LT - selection intensity of 26.56%). For G, the estimated genetic gains in relation to the population average were 26.80% and 17.94%, respectively, according to the ST and LT strategies. The effective population sizes were 22.35 and 46.03, respectively. The LT and ST strategies maintained 45.80% and 28.24%, respectively, of the original genetic diversity represented in the progeny test. So, it can be inferred that this population has potential for both breeding and ex situ genetic conservation as a supplier of genetic material for advanced rubber tree breeding programs. Copyright by the Brazilian Society of Genetics.

Cytogenetic and random amplified polymorphic DNA analysis of Leptodactylus species from rural and urban environments (Anura, Amphibia)

Cytogenetic and random amplified polymorphic DNA analyses carried out in the species Leptodactylus podicipinus, L. ocellatus, L. labyrinthicus, and L. fuscus from rural and urban habitats of the northwest region of São Paulo State, Brazil, showed that the karyotypes (2n = 22), constitutive heterochromatin distribution and nucleolus organizer region (NOR) location did not differ between the populations from the two environments. The in situ hybridization with an rDNA probe confirmed the location of the NORs on chromosome 8 revealing an in tandem duplication of that region in one of the chromosomes of L. fuscus. DAPI showed that part of the C-band-positive heterochromatin is rich in AT, including that in the proximity the NORs in L. podicipinus and L. ocellatus. The molecular analyses showed that the two populations (urban and rural) of L. podicipinus and L. fuscus are similar from a genetic point of view. The urban and rural populations of species L. ocellatus and L. labyrinthicus showed differences in genetic structures, probably due to urbanization which interferes with the dispersion of those frogs. The marked differences observed between the two populations of L. ocellatus can be representing the cryptic condition of the species. Unweighted pair-group method of analysis and genetic distance analysis detected the genetic proximity between L. ocellatus and L. fuscus. The results indicate that there was no reduction in the genetic diversity in the populations from the urban environment; however, the survival of these frogs would not be guaranteed in the case of an increase in human impact especially for populations of L. labyrinthicus and L. ocellatus. ©FUNPEC-RP.

Descrição cariotípica de possíveis novas espécies do complexo Hypostomus ancistroides e outros Hypostominae

Cytogenetic analyses were performed in four species of the Hypostominae subfamily, three from Hypostomus (Hypostomini) genus and Rhinelepis aspera (Rhinelepini). Three populations of Hypostomus ancistroides were analyzed, which had 2n=68 chromosomes, but presented different karyotype formulas. Hypostomus regani and H. strigaticeps, both from Ivaí river, showed 2n=72 chromosomes with two distinct cytotypes. In turn, R. aspera of the upper Paraná river basin presented 2n=54 chromosome. Multiple Nucleolar Organizer Regions (NORs) have been evidenced by silver nitrate staining in species of Hypostomus and single NOR in R. aspera. The observed variation in the chromosome number and the marked variability in karyotype formulas and NORs reveal a certain amount of karyotype variation in the genus Hypostomus suggesting the probable existence of cryptic species with independent chromosome traits. Therefore, our data can be of great value in discriminating species and understanding their chromosomal evolution.

Detection of nonrandom association of alleles from the distribution of the number of heterozygous loci in a sample.

The distribution of the number of heterozygous loci in two randomly chosen gametes or in a random diploid zygote provides information regarding the nonrandom association of alleles among different genetic loci. Two alternative statistics may be employed for detection of nonrandom association of genes of different loci when observations are made on these distributions: observed variance of the number of heterozygous loci (s2k) and a goodness-of-fit criterion (X2) to contrast the observed distribution with that expected under the hypothesis of random association of genes. It is shown, by simulation, that s2k is statistically more efficient than X2 to detect a given extent of nonrandom association. Asymptotic normality of s2k is justified, and X2 is shown to follow a chi-square (chi 2) distribution with partial loss of degrees of freedom arising because of estimation of parameters from the marginal gene frequency data. Whenever direct evaluations of linkage disequilibrium values are possible, tests based on maximum likelihood estimators of linkage disequilibria require a smaller sample size (number of zygotes or gametes) to detect a given level of nonrandom association in comparison with that required if such tests are conducted on the basis of s2k. Summarization of multilocus genotype (or haplotype) data, into the different number of heterozygous loci classes, thus, amounts to appreciable loss of information.

Parentage analysis with genetic markers in natural populations. I. The expected proportion of offspring with unambiguous paternity.

Recent studies indicate that polymorphic genetic markers are potentially helpful in resolving genealogical relationships among individuals in a natural population. Genetic data provide opportunities for paternity exclusion when genotypic incompatibilities are observed among individuals, and the present investigation examines the resolving power of genetic markers in unambiguous positive determination of paternity. Under the assumption that the mother for each offspring in a population is unambiguously known, an analytical expression for the fraction of males excluded from paternity is derived for the case where males and females may be derived from two different gene pools. This theoretical formulation can also be used to predict the fraction of births for each of which all but one male can be excluded from paternity. We show that even when the average probability of exclusion approaches unity, a substantial fraction of births yield equivocal mother-father-offspring determinations. The number of loci needed to increase the frequency of unambiguous determinations to a high level is beyond the scope of current electrophoretic studies in most species. Applications of this theory to electrophoretic data on Chamaelirium luteum (L.) shows that in 2255 offspring derived from 273 males and 70 females, only 57 triplets could be unequivocally determined with eight polymorphic protein loci, even though the average combined exclusionary power of these loci was 73%. The distribution of potentially compatible male parents, based on multilocus genotypes, was reasonably well predicted from the allele frequency data available for these loci. We demonstrate that genetic paternity analysis in natural populations cannot be reliably based on exclusionary principles alone. In order to measure the reproductive contributions of individuals in natural populations, more elaborate likelihood principles must be deployed.

Estimation of the incidence of a rare genetic disease through a two-tier mutation survey.

Recent attempts to detect mutations involving single base changes or small deletions that are specific to genetic diseases provide an opportunity to develop a two-tier mutation-screening program through which incidence of rare genetic disorders and gene carriers may be precisely estimated. A two-tier survey consists of mutation screening in a sample of patients with specific genetic disorders and in a second sample of newborns from the same population in which mutation frequency is evaluated. We provide the statistical basis for evaluating the incidence of affected and gene carriers in such two-tier mutation-screening surveys, from which the precision of the estimates is derived. Sample-size requirements of such two-tier mutation-screening surveys are evaluated. Considering examples of cystic fibrosis (CF) and medium-chain acyl-CoA dehydrogenase deficiency (MCAD), the two most frequent autosomal recessive disease in Caucasian populations and the two most frequent mutations (delta F508 and G985) that occur on these disease allele-bearing chromosomes, we show that, with 50-100 patients and a 20-fold larger sample of newborns screened for these mutations, the incidence of such diseases and their gene carriers in a population may be quite reliably estimated. The theory developed here is also applicable to rare autosomal dominant diseases for which disease-specific mutations are found.