Population genetics of CAPN10 and GPR35: implications for the evolution of type 2 diabetes variants

A positional cloning study of type 2 diabetes in Mexican Americans identified a region, termed "NIDDM1," on chromosome 2q37 with significant linkage evidence. Haplotype combinations at the calpain-10 gene (CAPN10) within this region were shown to increase diabetes risk in several populations. On the basis of the thrifty genotype hypothesis, variants that increase susceptibility to type 2 diabetes under modern lifestyle conditions provided a survival advantage in past environments by increasing the efficiency of energy use and storage. Here, our goal is to make inferences about the evolutionary forces shaping variation in genes in the NIDDM1 region and to investigate the population genetics models that may underlie the thrifty genotype hypothesis. To this end, we surveyed sequence variation in CAPN10 and in an adjacent gene, G-protein-coupled receptor 35 (GPR35), in four population samples from different ethnic groups. These data revealed two distinct deviations from the standard neutral model in CAPN10, whereas GPR35 variation was largely consistent with neutrality. CAPN10 showed a significant deficit of variation in the haplotype class defined by the derived allele at SNP44, a polymorphism that is significantly associated with diabetes in meta-analysis studies. This suggests that this haplotype class was quickly driven to high frequency by positive natural selection. Interestingly, the derived allele at SNP44 is protective against diabetes. CAPN10 also showed a local excess of polymorphism and linkage disequilibrium decay in intron 13. Simulations show that this pattern may be explained by long-standing balancing selection that maintains multiple selected alleles. Alternatively, it is possible that the local mutation and recombination rates changed since the divergence of human and chimpanzee; this scenario does not require the action of natural selection on intron 13 variation.     

Combining sperm typing and linkage disequilibrium analyses reveals differences in selective pressures or recombination rates across human populations

A previous polymorphism survey of the type 2 diabetes gene CAPN10 identified a segment showing an excess of polymorphism levels in all population samples, coinciding with localized breakdown of linkage disequilibrium (LD) in a sample of Hausa from Cameroon, but not in non-African samples. This raised the possibility that a recombination hotspot is present in all populations and we had insufficient power to detect it in the non-African data. To test this possibility, we estimated the crossover rate by sperm typing in five non-African men; these estimates were consistent with the LD decay in the non-African, but not in the Hausa data. Moreover, resequencing the orthologous region in a sample of Western chimpanzees did not show either an excess of polymorphism level or rapid LD decay, suggesting that the processes underlying the patterns observed in humans operated only on the human lineage. These results suggest that a hotspot of recombination has recently arisen in humans and has reached higher frequency in the Hausa than in non-Africans, or that there is no elevation in crossover rate in any human population, and the observed variation results from long-standing balancing selection.     

Disentangling the effects of recombination,selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen‐mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d_n/d_s ratio and the presence of trans‐species polymorphisms suggest that a strong long‐term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D‐loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000–30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.     

Polymorphism and divergence at the prune locus in Drosophila melanogaster and D. simulans

The prune locus of Drosophila melanogaster lies at the tip of the X chromosome, in a region of reduced recombination in which nearby loci show reduced variation relative to evolutionary divergence from D. simulans. DNA sequencing of prune alleles from D. melanogaster and D. simulans reveals extremely low variation in D. melanogaster but greater variation in D. simulans. Divergence between the two species is not reduced. This pattern may be explained by either positive selection leading to hitchhiking of neutral variation or background selection against deleterious mutations. The pattern of silent versus replacement polymorphism and divergence at prune is consistent with either a model of weakly deleterious selection against amino acid substitutions or balancing selection.      

Are variants in the CAPN10 gene related to risk of type 2 diabetes? A quantitative assessment of population and family-based association studies

The calpain-10 gene (CAPN10) on chromosome 2q37.3 was the first candidate gene for type 2 diabetes (T2D) identified through a genomewide screen and positional cloning. One polymorphism (UCSNP-43: G-->A) and a specific haplotype combination defined by three polymorphisms (UCSNP-43, -19, and -63) were linked to an increased risk of T2D in several populations. To quantitatively assess the collective evidence for the effects of CAPN10 on risk of T2D, we conducted a meta-analysis of both population-based and family-based association studies. We retrieved data from the MEDLINE, PubMed, and Online Mendelian Inheritance in Man databases, as well as from other relevant reports and abstracts published up to July 2003. From a total of 26 studies with primary data (21 population-based studies: 5,013 cases and 5,876 controls; 5 family-based studies: 487 parent-offspring trios), we developed a summary database that contains variables of study design, study population/ethnicity, specific polymorphisms and haplotype combinations in CAPN10, and diabetes-related metabolic phenotypes. For population-based studies, we used both fixed-effects and random-effects models to calculate the pooled odds ratio (OR) and 95% confidence interval (CI) for the associations of CAPN10 genotypes with the risk of T2D. We also calculated weighted mean differences for the associations between CAPN10 and diabetes-related quantitative traits. Under either an additive or a dominant effect model, we found no statistically significant relation between CAPN10 genotypes in the UCSNP-43 locus and T2D risk. However, under a recessive model, individuals homozygous for the common G allele had a statistically significant 19% higher risk of T2D than carriers of the A allele (OR 1.19; 95% CI 1.07-1.33). The association between the 112/121 haplotype combination and T2D risk appeared to be overestimated by several initial small studies with positive findings (OR 1.38; 95% CI 1.04-1.84). After we removed these initial studies, this association became nonsignificant (OR 1.11; 95% CI 0.91-1.35). Moreover, we found no evidence for the associations between the UCSNP-43 G/G genotype and the 112/121 haplotype combination and metabolic phenotypes. Our meta-analysis of family-based studies showed only an overtransmission of the rare allele C in UCSNP-44 from heterozygous parents to their affected offspring with T2D. Our analysis indicates that inadequate statistical power, racial/ethnic differences in frequencies of alleles, haplotypes and haplotype combinations, potential gene-gene or gene-environment interactions, publication bias, and multiple hypothesis testing may contribute to the significant heterogeneity in previous studies of CAPN10 and T2D. Our findings also suggest that both large-scale, well-designed association studies and functional studies are warranted to either reliably confirm or conclusively refute the initial hypothesis regarding the role of CAPN10 in T2D risk.     

<Emphasis Type="Italic">CAPN1</Emphasis> markers in three Argentinean cattle breeds: report of a new InDel polymorphism within intron 17

In this study, the genotype distribution and allelic frequencies of CAPN1 (Calcium activated neutral protease) single nucleotide polymorphisms (SNPs) were analyzed taking advantage of the different genetic backgrounds provided by Hereford, Brahman and Braford cattle. We report a new insertion/deletion (InDel) polymorphism, consisting of a change of seven nucleotides for only one nucleotide (TCTGGGT → C) within intron 17 of the CAPN1 gene. The segregation pattern of this polymorphism was analyzed together with the markers CAPN316, CAPN530 and CAPN4751 already described. The allele distribution of CAPN1 markers in the Braford crossbreed (3/8 Brahman 5/8 Hereford) is described for the first time. Four assays of allelic discrimination were designed: the tetra primer ARMS-PCR technique for genotyping the new InDel and the CAPN4751 marker, and a PCR-RFLP method for genotyping the markers CAPN316 and CAPN530. The genotypic and minor allele frequencies (MAFs) obtained showed that the InDel polymorphism does not provide redundant information to that already provided by the other CAPN1 markers and segregates differently between breeds, being a common SNP (MAF ≥ 0.05) in the herds with a high percentage of Bos indicus background. The high percentage of heterozygous individuals found in the Braford crossbreed for the markers assessed reveals enough genetic variation that could help to solve the tenderness problem of tropical-adapted cattle.     

Genetic Variation of Microsatellite Loci in the Major Histocompatibility Complex (MHC) Region in the Southeast Asian House Mouse (<Emphasis Type="Italic">Mus musculus castaneus</Emphasis>)

Major histocompatibility complex (MHC) genes are the most polymorphic loci known for vertebrates. Here we employed five microsatellite loci closely linked to the MHC region in an attempt to study the amount of genetic variation in 19 populations of the southeast Asian house mouse (Mus musculus castaneus) in Taiwan. The overall polymorphism at the five loci was high (H_e = 0.713), and the level of polymorphism varied from locus to locus. Furthermore, in order to investigate if selection is operating on MHC genes in natural mouse populations, we compared the extent and pattern of genetic variation for the MHC-linked microsatellite loci (the MHC loci) with those for the microsatellite loci located outside the MHC region (the non-MHC loci). The number of alleles and the logarithm of variance in repeat number were significantly higher for the MHC loci than for the non-MHC loci, presumably reflecting linkage to a locus under balancing selection. Although three statistical tests used do not provide support for selection, their lack of support may be due to low statistical power of the tests, to weakness of selection, or to a profound effect of genetic drift reducing the signature of balancing selection. Our results also suggested that the populations in the central and the southwestern regions of Taiwan might be one part of a metapopulation structure.     

A Comparison of Intraspecific Patterns of DNA Sequence Variation in Mitochondrial DNA, α-Enolase,and MHC Class II B Loci in Auklets (Charadriiformes: Alcidae)

Patterns of DNA sequence variation can be used to learn about mechanisms of organismal evolution, but only if mechanisms of sequence evolution are well understood. Although theories of molecular evolution are well developed, few empirical studies have addressed patterns and mechanisms of sequence evolution in nuclear genes within species. In the present study, we compared DNA sequences among three loci with different evolutionary constraints to determine the influences of effective population size, balancing selection, and linkage on intraspecific patterns of sequence variation. Specifically, we assessed the degree and nature of polymorphism in a 307-base pair (bp) fragment of the mitochondrial cytochrome b gene, intron VIII of the gene for -enolase (a presumably neutral nuclear gene), and an ~600-bp fragment of an MHC class II B gene, including 155 bp of the hypervariable peptide binding region (a nuclear locus thought to be under balancing selection) for least and crested auklets (Aethia pusilla and A. cristatella; Charadriiformes: Alcidae). Transspecies polymorphism was found in both -enolase and the MHC but not cytochrome b and, given estimates of effective population size, probably represents retained ancestral variation. Biases in nucleotide composition suggested that mutational bias, tRNA availability, and the secondary structure of mRNA and/or DNA may influence base usage. Several lines of evidence indicated that balancing selection may be acting on the MHC II B exon 2. However, no evidence of balancing selection was observed in the intron and exon sequences immediately downstream of MHC II B exon 2. Current address (Hollie E. Walsh): Department of Zoology, University of Washington, Box 351800, Seattle, WA 98195-1800, USA     

Long‐term balancing selection on chromosomal variants associated with crypsis in a stick insect

How polymorphisms are maintained within populations over long periods of time remains debated, because genetic drift and various forms of selection are expected to reduce variation. Here, we study the genetic architecture and maintenance of phenotypic morphs that confer crypsis in Timema cristinae stick insects, combining phenotypic information and genotyping‐by‐sequencing data from 1,360 samples across 21 populations. We find two highly divergent chromosomal variants that span megabases of sequence and are associated with colour polymorphism. We show that these variants exhibit strongly reduced effective recombination, are geographically widespread and probably diverged millions of generations ago. We detect heterokaryotype excess and signs of balancing selection acting on these variants through the species’ history. A third chromosomal variant in the same genomic region likely evolved more recently from one of the two colour variants and is associated with dorsal pattern polymorphism. Our results suggest that large‐scale genetic variation associated with crypsis has been maintained for long periods of time by potentially complex processes of balancing selection.     

The evolution of TEP1, an exceptionally polymorphic immunity gene in <Emphasis Type="Italic">Anopheles gambiae</Emphasis>

Background  

Host-parasite coevolution can result in balancing selection, which maintains genetic variation in the susceptibility of hosts to parasites. It has been suggested that variation in a thioester-containing protein called TEP1 (AGAP010815) may alter the ability of Anopheles mosquitoes to transmit Plasmodium parasites, and high divergence between alleles of this gene suggests the possible action of long-term balancing selection. We studied whether TEP1 is a case of an ancient balanced polymorphism in an animal immune system.     

High nucleotide sequence variation in a region of low recombination in Drosophila simulans is consistent with the background selection model

We surveyed nucleotide sequence variation at glucose dehydrogenase (Gld), in a region of low recombination on chromosome 3R, from a population sample of Drosophila simulans. The levels of nucleotide variation were surprisingly high. There was no departure from the expectation of a neutral model for the level of polymorphism, indicating no evidence of a selective sweep in this region. There was a significant deficiency of singleton polymorphisms according to the Fu and Li test, although Tajima and Hudson, Kreitman, and Aguade (HKA) tests do not provide evidence of a significant elevation of variation due to balancing selection. Genetic map data for the D. simulans third chromosome were used to calculate expected values of pi for Gld under a current model of background selection, varying the values for the parameter sh (selection coefficient against deleterious mutations). We show that the recombinational landscape of D. simulans is sufficiently different from that of D. melanogaster that we expect higher variation under the background selection model, even when effective population sizes are assumed to be equal. The data for Gld were tested against the predictions using computer simulations of the distribution of the number of segregating sites conditioned on pi. Background selection alone can explain our observations as long as sh is larger than 0.005 and species-level effective population size is assumed to be several- fold larger than in D. melanogaster. Alternatively, the deleterious mutation rate may be smaller in D. simulans, or balancing selection may be acting nearby, thereby reducing the effect of background selection.      

Genome‐wide association analysis for quantitative trait loci influencing Warner–Bratzler shear force in five taurine cattle breeds

We performed a genome‐wide association study for Warner–Bratzler shear force (WBSF), a measure of meat tenderness, by genotyping 3360 animals from five breeds with 54 790 BovineSNP50 and 96 putative single‐nucleotide polymorphisms (SNPs) within μ‐calpain [HUGO nomenclature calpain 1, (mu/I) large subunit; CAPN1] and calpastatin (CAST). Within‐ and across‐breed analyses estimated SNP allele substitution effects (ASEs) by genomic best linear unbiased prediction (GBLUP) and variance components by restricted maximum likelihood under an animal model incorporating a genomic relationship matrix. GBLUP estimates of ASEs from the across‐breed analysis were moderately correlated (0.31–0.66) with those from the individual within‐breed analyses, indicating that prediction equations for molecular estimates of breeding value developed from across‐breed analyses should be effective for genomic selection within breeds. We identified 79 genomic regions associated with WBSF in at least three breeds, but only eight were detected in all five breeds, suggesting that the within‐breed analyses were underpowered, that different quantitative trait loci (QTL) underlie variation between breeds or that the BovineSNP50 SNP density is insufficient to detect common QTL among breeds. In the across‐breed analysis, CAPN1 was followed by CAST as the most strongly associated WBSF QTL genome‐wide, and associations with both were detected in all five breeds. We show that none of the four commercialized CAST and CAPN1 SNP diagnostics are causal for associations with WBSF, and we putatively fine‐map the CAPN1 causal mutation to a 4581‐bp region. We estimate that variation in CAST and CAPN1 explains 1.02 and 1.85% of the phenotypic variation in WBSF respectively.     

Sequence diversity patterns suggesting balancing selection in partially sex‐linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex‐linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within‐species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex‐determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation‐based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex‐linked genes, which requires closer linkage with the fully sex‐linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.     

Introgression from Domestic Goat Generated Variation at the Major Histocompatibility Complex of Alpine Ibex

The major histocompatibility complex (MHC) is a crucial component of the vertebrate immune system and shows extremely high levels of genetic polymorphism. The extraordinary genetic variation is thought to be ancient polymorphisms maintained by balancing selection. However, introgression from related species was recently proposed as an additional mechanism. Here we provide evidence for introgression at the MHC in Alpine ibex (Capra ibex ibex). At a usually very polymorphic MHC exon involved in pathogen recognition (DRB exon 2), Alpine ibex carried only two alleles. We found that one of these DRB alleles is identical to a DRB allele of domestic goats (Capra aegagrus hircus). We sequenced 2489 bp of the coding and non-coding regions of the DRB gene and found that Alpine ibex homozygous for the goat-type DRB exon 2 allele showed nearly identical sequences (99.8%) to a breed of domestic goats. Using Sanger and RAD sequencing, microsatellite and SNP chip data, we show that the chromosomal region containing the goat-type DRB allele has a signature of recent introgression in Alpine ibex. A region of approximately 750 kb including the DRB locus showed high rates of heterozygosity in individuals carrying one copy of the goat-type DRB allele. These individuals shared SNP alleles both with domestic goats and other Alpine ibex. In a survey of four Alpine ibex populations, we found that the region surrounding the DRB allele shows strong linkage disequilibria, strong sequence clustering and low diversity among haplotypes carrying the goat-type allele. Introgression at the MHC is likely adaptive and introgression critically increased MHC DRB diversity in the genetically impoverished Alpine ibex. Our finding contradicts the long-standing view that genetic variability at the MHC is solely a consequence of ancient trans-species polymorphism. Introgression is likely an underappreciated source of genetic diversity at the MHC and other loci under balancing selection.     

A signature of balancing selection in the region upstream to the human <Emphasis Type="Italic">UGT2B4</Emphasis> gene and implications for breast cancer risk

UDP-glucuronosyltransferase 2 family, polypeptide B4 (UGT2B4) is an important metabolizing enzyme involved in the clearance of many xenobiotics and endogenous substrates, especially steroid hormones and bile acids. The HapMap data show that numerous SNPs upstream of UGT2B4 are in near-perfect linkage disequilibrium with each other and occur at intermediate frequency, indicating that this region might contain a target of natural selection. To investigate this possibility, we chose three regions (4.8 kb in total) for resequencing and observed a striking excess of intermediate-frequency alleles that define two major haplotypes separated by many mutation events and with little differentiation across populations, thus suggesting that the variation pattern upstream UGT2B4 is highly unusual and may be the result of balancing selection. We propose that this pattern is due to the maintenance of a regulatory polymorphism involved in the fine tuning of UGT2B4 expression so that heterozygous genotypes result in optimal enzyme levels. Considering the important role of steroid hormones in breast cancer susceptibility, we hypothesized that variation in this region could predispose to breast cancer. To test this hypothesis, we genotyped tag SNP rs13129471 in 1,261 patients and 825 normal women of African ancestry from three populations. The frequency comparison indicated that rs13129471 was significantly associated with breast cancer after adjusting for ethnicity [P = 0.003; heterozygous odds ratio (OR) 1.02, 95% confidence interval (CI) 0.81–1.28; homozygous OR 1.50, 95% CI 1.15–1.95]. Our results provide new insights into UGT2B4 sequence variation and indicate that a signal of natural selection may lead to the identification of disease susceptibility variants.     

Lineage pattern, trans-species polymorphism, and selection pressure among the major lineages of feline Mhc-DRB peptide-binding region

The long-term evolution of major histocompatibility complex (MHC) involves the birth-and-death process and independent divergence of loci during episodes punctuated by natural selection. Here, we investigated the molecular signatures of natural selection at exon-2 of MHC class II DRB gene which includes a part of the peptide-binding region (PBR) in seven of eight putative extant Felidae lineages. The DRB alleles in felids can be mainly divided into five lineages. Signatures of trans-species polymorphism among major allelic lineages indicate that balancing selection has maintained the MHC polymorphism for a long evolutionary time. Analysis based on maximum likelihood models of codon substitution revealed overall purifying selection acting on the feline DRB. Sites that have undergone positive selection and those that are under divergent selective pressure among lineages were detected and found to fall within the putative PBR. This study increased our understanding of the nature of selective forces acting on DRB during feline radiation.     

SELECTION ON THE COLOR POLYMORPHISM IN HAWAIIAN HAPPY-FACE SPIDERS: EVIDENCE FROM GENETIC STRUCTURE AND TEMPORAL FLUCTUATIONS

Throughout this century genetic polymorphisms for color have been widely used as a research tool to allow insights into key evolutionary processes. Although color variants can often be diverse within populations, frequencies of different morphs may be similar across populations, either as a result of balancing selection or gene flow. Under these circumstances selection can be extremely difficult to demonstrate. Here we test for balancing selection on the naturally occurring color forms of the Hawaiian happy-face spider, Theridion grallator with two approaches. First, allozyme loci are used to generate a null model against which to test selection. Frequencies of alleles involved in the color polymorphism of T. grallator are used to generate another estimate for comparison. The results suggest that statistically similar frequencies of color morphs among populations of T. grallator may be maintained by some form of balancing selection. Second, we make use of an unusual event in which the normally stable frequencies of unpatterned and patterned morphs within a population were found to have shifted toward an excess of unpatterned morphs. We scored offspring of all fertilized, unpatterned (bottom-recessive) females found during this period of skewed morph frequencies and also in a year when morph frequencies were normal to deduce paternal color phenotypes. Mating was found to be random in the normal year, but in the perturbed year females had mated with rare (patterned) males twice as frequently as expected on the basis of the frequency of this morph type in the population. Both of these results are consistent with selection operating on the color polymorphism, and we speculate that apostatic selection, perhaps mediated by bird predators, may provide the mechanism.     

Balancing selection maintains cryptic colour morphs

Animals display incredibly diverse colour patterns, a testament to evolution's endless innovation in shaping life. In many species, the interplay between males and females in the pursuit of mates has driven the evolution of a myriad of colour forms, from the flashy peacock tail feathers to the tiniest colour markings in damselflies. In others, colour provides crypsis by allowing to blend into the background and to escape the eyes of predators. While the obvious benefits of this dazzling diversity for reproduction and survival seem straightforward, its maintenance is not. Theory predicts that genetic drift and various forms of selection reduce variation over time, making the persistence of colour variants over generations a puzzle. In this issue of Molecular Ecology, Lindtke et al. ( 2017 ) study the cryptic colour morphs of Timema cristinae walking sticks to shed light on the genetic architecture and mechanisms that allow colour polymorphism maintenance over long timescales. By combining genome‐wide data with phenotyping information from natural populations, they were able to map the green and melanistic colour to one genomic region with highly reduced effective recombination rate between two main chromosomal variants, consistent with an inversion polymorphism. These two main chromosomal variants showed geographically widespread heterozygote excess, and genomic signatures consistent with long‐term balancing selection. A younger chromosomal variant was detected for the third morph, the green‐striped colour morphs, in the same genomic regions as the melanistic and the green‐unstriped morphs. Together, these results suggest that the genetic architecture of cryptic T. cristinae morphs is caused by nonrecombining genomic blocks that have been maintained over extended time periods by balancing selection making this study one of the few available empirical examples documenting that balancing selection of various forms may play an important role in maintaining adaptive genetic variation in nature.     

Unconstrained evolution in short introns? – An analysis of genome‐wide polymorphism and divergence data from Drosophila

An unconstrained reference sequence facilitates the detection of selection. In Drosophila, sequence variation in short introns seems to be least influenced by selection and dominated by mutation and drift. Here, we test this with genome‐wide sequences using an African population (Malawi) of D. melanogaster and data from the related outgroup species D. simulans, D. sechellia, D. erecta and D. yakuba. The distribution of mutations deviates from equilibrium, and the content of A and T (AT) nucleotides shows an excess of variance among introns. We explain this by a complex mutational pattern: a shift in mutational bias towards AT, leading to a slight nonequilibrium in base composition and context‐dependent mutation rates, with G or C (GC) sites mutating most frequently in AT‐rich introns. By comparing the corresponding allele frequency spectra of AT‐rich vs. GC‐rich introns, we can rule out the influence of directional selection or biased gene conversion on the mutational pattern. Compared with neutral equilibrium expectations, polymorphism spectra show an excess of low frequency and a paucity of intermediate frequency variants, irrespective of the direction of mutation. Combining the information from different outgroups with the polymorphism data and using a generalized linear model, we find evidence for shared ancestral polymorphism between D. melanogaster and D. simulans, D. sechellia, arguing against a bottleneck in D. melanogaster. Generally, we find that short introns can be used as a neutral reference on a genome‐wide level, if the spatially and temporally varying mutational pattern is accounted for.