Recent positive selection of a human androgen receptor/ectodysplasin A2 receptor haplotype and its relationship to male pattern baldness

Genetic variants in the human androgen receptor gene (AR) are associated with male pattern baldness (androgenetic alopecia, AGA) in Europeans. Previous observations of long-range linkage disequilibrium at the AR locus are consistent with the hypothesis of recent positive selection. Here, we further investigate this signature and its relationship to the AGA risk haplotype. The haplotype homozygosity suggests that the AGA risk haplotype was driven to high frequency by positive selection in Europeans although a low meiotic recombination rate contributed to the high haplotype homozygosity. Further, we find high levels of population differentiation as measured by F _ST and a series of fixed derived alleles along an extended region centromeric to AR in the Asian HapMap sample. The predominant AGA risk haplotype also carries the putatively functional variant 57K in the flanking ectodysplasin A2 receptor gene (EDA2R). It is therefore probable that the AGA risk haplotype rose to high frequency in combination with this EDA2R variant, possibly by hitchhiking on a positively selected 57K haplotype. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Polymorphism in the neurofibromin gene,Nf1, is associated with antagonistic selection on wing size and development time in Drosophila melanogaster

In many invertebrates, body size shows genetically based clines, with size increasing in colder climates. Large body size is typically associated with prolonged development times. We consider variation in the CNS‐specific gene neurofibromin 1 (Nf1) and its association with body size and development time. We identified two major Nf1 haplotypes in natural populations, Nf1‐insertion‐A and Nf1‐deletion‐G. These haplotypes are characterized by a 45‐base insertion/deletion (INDEL) in Nf1 intron 2 and an A/G synonymous substitution (locus L17277). Linkage disequilibrium (LD) between the INDEL and adjacent sites is high but appears to be restricted within the Nf1 gene interval. In Australia, the frequency of the Nf1‐insertion‐A haplotype increases with latitude where wing size is larger, independent of the chromosomal inversion In(3R)Payne. Unexpectedly, the Nf1‐insertion‐A haplotype is negatively associated with wing size. We found that the Nf1‐insertion‐A haplotype is enriched in females with shorter development time. This suggests that the Nf1 haplotype cline may be driven by selection for development time rather than size; females from southern (higher latitude) D. melanogaster populations maintain a rapid development time despite being relatively larger, and the higher incidence of Nf1‐insertion‐A in Southern Australia may contribute to this pattern, whereas the effects of the Nf1 haplotypes on size may be countered by other loci with antagonistic effects on size and development time. Our results point to the potential complexity involved in identifying selection on genetic variants exhibiting pleiotropic effects when studies are based on spatial patterns or association studies.     

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life‐history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness‐related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O₃₊₄ arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Individual inversions or their combinations: which is the main selective target in a natural population of Drosophila subobscura?

It is generally accepted that chromosomal inversions have been key elements in adaptation and speciation processes. In this context, Drosophila subobscura has been, and still is, an excellent model species due to its rich chromosomal polymorphism. In this species, many analyses from natural populations have demonstrated the adaptive potential of individual inversions (and their overlapped combinations, the so‐called arrangements). However, little information is available on the evolutionary role of combinations generated by inversions located in homologous and nonhomologous chromosomes. The aim of this research was to ascertain whether these combinations are also a target for natural selection. For this objective, we have studied the inversion composition of homologous and nonhomologous chromosomes from a D. subobscura sample collected in a well‐studied population, Mount Avala (Serbia). No significant deviation from H‐W expectations was detected, and when comparing particular karyotypic combinations, likelihood ratios close to 1 were obtained. Thus, it seems that for each pair of homologous chromosomes inversions no deviation from randomness was detected. Finally, no linkage disequilibrium was observed between inversions located in different chromosomes of the karyotype. For all these reasons, it can be assumed that, at the cytological level, the individual inversions rather than their combinations in different chromosomes are the main target of selection.     

Integrative taxonomy of the Rhinolophus macrotis complex (Chiroptera,Rhinolophidae) in Vietnam and nearby regions

The taxonomic status of Rhinolophus macrotis sensu lato (s.l.) in Vietnam and adjacent territories remains problematic. To address this issue, we performed an integrated study of morphological, acoustic, and genetic characters of R. macrotis s.l. specimens and compared these with sympatric species within the philippinensis group (R. marshalli, R. paradoxolophus, and R. rex). Our results reveal that in addition to a cryptic species of R. macrotis previously found in Jiangxi and Jingmen, China, R. macrotis s.l. in continental Asia includes three further species, namely R. cf. siamensis, R. cf. macrotis, and R. cf. macrotis “Phia Oac.” These four taxa are distinguished from genuine R. macrotis in Nepal and R. siamensis in Thailand by their morphological and/or genetic features. Further taxonomic evaluation of the subspecies of R. macrotis s.l. is needed to determine their affinities with recently recognized cryptic species and to possibly describe new taxa. Our results also show that interspecific divergences in mitochondrial DNA sequences (Cytb and COI genes) among taxa within the philippinensis group (particularly between R. cf. siamensis/R. cf. macrotis and R. rex/R. paradoxolophus) are significantly lower than those of other morphological groups in the genus. These phylogenetic patterns might be explained by recent allopatric speciation or ancient introgression events among ancestors of the taxa during the Pleistocene. However, further investigations including genetic analyses of nuclear genes are needed to test the latter hypothesis.     

Molecular,genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana

Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17‐Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F‐box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col‐FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South‐Eastern Europe approximately 5000 years ago and the FRI‐Col allele to North‐West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The ‘American haplotype’ was identified from several European localities, potentially due to return migration.     

Invasion genetics of American cherry fruit fly in Europe and signals of hybridization with the European cherry fruit fly

The American cherry fruit fly is an invasive pest species in Europe, of serious concern in tart cherry production as well as for the potential to hybridize with the European cherry fruit fly, Rhagoletis cerasi L. (Diptera: Tephritidae), which might induce new pest dynamics. In the first European reports, the question arose whether only the eastern American cherry fruit fly, Rhagoletis cingulata (Loew) (Diptera: Tephritidae), is present, or also the closely related western American cherry fruit fly, Rhagoletis indifferens Curran. In this study, we investigate the species status of European populations by comparing these with populations of both American species from their native ranges, the invasion dynamics in German (first report in 1993) and Hungarian (first report in 2006) populations, and we test for signals of hybridization with the European cherry fruit fly. Although mtDNA sequence genealogy could not separate the two American species, cross‐species amplification of 14 microsatellite loci separated them with high probabilities (0.99–1.0) and provided evidence for R. cingulata in Europe. German and Hungarian R. cingulata populations differed significantly in microsatellite allele frequencies, mtDNA haplotype and wing pattern distributions, and both were genetically depauperate relative to North American populations. The diversity suggests independent founding events in Germany and Hungary. Within each country, R. cingulata displayed little or no structure in any trait, which agrees with rapid local range expansions. In cross‐species amplifications, signals of hybridization between R. cerasi and R. cingulata were found in 2% of R. cingulata individuals and in 3% of R. cerasi. All putative hybrids had R. cerasi mtDNA indicating that the original between‐species mating involved R. cerasi females and R. cingulata males.     

Beta‐globin gene evolution in the ruminants: evidence for an ancient origin of sheep haplotype B

Domestic sheep (Ovis aries) can be divided into two groups with significantly different responses to hypoxic environments, determined by two allelic beta‐globin haplotypes. Haplotype A is very similar to the goat beta‐globin locus, whereas haplotype B has a deletion spanning four globin genes, including beta‐C globin, which encodes a globin with high oxygen affinity. We surveyed the beta‐globin locus using resequencing data from 70 domestic sheep from 42 worldwide breeds and three Ovis canadensis and two Ovis dalli individuals. Haplotype B has an allele frequency of 71.4% in O. aries and was homozygous (BB) in all five wild sheep. This shared ancestry indicates haplotype B is at least 2–3 million years old. Approximately 40 kb of the sequence flanking the ~37‐kb haplotype B deletion had unexpectedly low identity between haplotypes A and B. Phylogenetic analysis showed that the divergent region of sheep haplotype B is remarkably distinct from the beta‐globin loci in goat and cattle but still groups with the Ruminantia. We hypothesize that this divergent ~40‐kb region in haplotype B may be from an unknown ancestral ruminant and was maintained in the lineage to O. aries, but not other Bovidae, evolving independently of haplotype A. Alternatively, the ~40‐kb sequence in haplotype B was more recently acquired by an ancestor of sheep from an unknown non‐Bovidae ruminant, replacing part of haplotype A. Haplotype B has a lower nucleotide diversity than does haplotype A, suggesting a recent bottleneck, whereas the higher frequency of haplotype B suggests a subsequent spread through the global population of O. aries.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

Detecting selected haplotype blocks in evolve and resequence experiments

Shifting from the analysis of single nucleotide polymorphisms to the reconstruction of selected haplotypes greatly facilitates the interpretation of evolve and resequence (E&R) experiments. Merging highly correlated hitchhiker SNPs into haplotype blocks reduces thousands of candidates to few selected regions. Current methods of haplotype reconstruction from Pool‐seq data need a variety of data‐specific parameters that are typically defined ad hoc and require haplotype sequences for validation. Here, we introduce haplovalidate, a tool which detects selected haplotypes in Pool‐seq time series data without the need for sequenced haplotypes. Haplovalidate makes data‐driven choices of two key parameters for the clustering procedure, the minimum correlation between SNPs constituting a cluster and the window size. Applying haplovalidate to simulated E&R data reliably detects selected haplotype blocks with low false discovery rates. Importantly, our analyses identified a restriction of the haplotype block‐based approach to describe the genomic architecture of adaptation. We detected a substantial fraction of haplotypes containing multiple selection targets. These blocks were considered as one region of selection and therefore led to underestimation of the number of selection targets. We demonstrate that the separate analysis of earlier time points can significantly increase the separation of selection targets into individual haplotype blocks. We conclude that the analysis of selected haplotype blocks has great potential for the characterization of the adaptive architecture with E&R experiments.     

Polymorphic chromosomal inversions in Anopheles moucheti,a major malaria vector in Central Africa

Anopheles moucheti Evans (Diptera: Culicidae) is a major vector of malaria in forested areas of Central Africa. However, few genetic tools are available for this species. The present study represents the first attempt to characterize chromosomes in An. moucheti females collected in Cameroon. Ovarian nurse cells contained polytene chromosomes, which were suitable for standard cytogenetic applications. The presence of three polymorphic chromosomal inversions in An. moucheti was revealed. Two of these inversions were located on the 2R chromosome arm. The homology between the 2R chromosome arms of An. moucheti and Anopheles gambiae Giles was established by fluorescent in situ hybridization of six An. gambiae genic sequences. Mapping of the probes on chromosomes of An. moucheti detected substantial gene order reshuffling between the two species. The presence of polytene chromosomes and polymorphic inversions in An. moucheti provides a new basis for further population genetic, taxonomic and ecological studies of this neglected malaria vector.     

A targeted genotyping approach enhances identification of variants in taste receptor and appetite/reward genes of potential functional importance for obesity‐related porcine traits

Taste receptors (TASRs) and appetite and reward (AR) mechanisms influence eating behaviour, which in turn affects food intake and risk of obesity. In a previous study, we used next generation sequencing to identify potentially functional mutations in TASR and AR genes and found indications for genetic associations between identified variants and growth and fat deposition in a subgroup of animals (n = 38) from the UNIK resource pig population. This population was created for studying obesity and obesity‐related diseases. In the present study we validated results from our previous study by investigating genetic associations between 24 selected single nucleotide variants in TASR and AR gene variants and 35 phenotypes describing obesity and metabolism in the entire UNIK population (n = 564). Fifteen variants showed significant association with specific obesity‐related phenotypes after Bonferroni correction. Six of the 15 genes, namely SIM1, FOS, TAS2R4, TAS2R9, MCHR2 and LEPR, showed good correlation between known biological function and associated phenotype. We verified a genetic association between potentially functional variants in TASR/AR genes and growth/obesity and conclude that the combination of identification of potentially functional variants by next generation sequencing followed by targeted genotyping and association studies is a powerful and cost‐effective approach for increasing the power of genetic association studies.     

HRAS1 and LASS1 with APOE are associated with human longevity and healthy aging

The search for longevity‐determining genes in human has largely neglected the operation of genetic interactions. We have identified a novel combination of common variants of three genes that has a marked association with human lifespan and healthy aging. Subjects were recruited and stratified according to their genetically inferred ethnic affiliation to account for population structure. Haplotype analysis was performed in three candidate genes, and the haplotype combinations were tested for association with exceptional longevity. An HRAS1 haplotype enhanced the effect of an APOE haplotype on exceptional survival, and a LASS1 haplotype further augmented its magnitude. These results were replicated in a second population. A profile of healthy aging was developed using a deficit accumulation index, which showed that this combination of gene variants is associated with healthy aging. The variation in LASS1 is functional, causing enhanced expression of the gene, and it contributes to healthy aging and greater survival in the tenth decade of life. Thus, rare gene variants need not be invoked to explain complex traits such as aging; instead rare congruence of common gene variants readily fulfills this role. The interaction between the three genes described here suggests new models for cellular and molecular mechanisms underlying exceptional survival and healthy aging that involve lipotoxicity.     

Haplotype block and superblock structures of the alpha1-adrenergic receptor genes reveal echoes from the chromosomal past

A significant proportion of the human genome is contained within haplotype blocks across which pairwise linkage disequilibrium (LD) is very high. However, LD is also often high between markers at more remote distances, and within different haplotype blocks. Here, we evaluate the origins of haplotype block structure in the three genes for alpha1 adrenergic receptors (alpha1-AR) in the human genome ( ADRA1A, ADRA1B and ADRA1D) by genotyping dense single-nucleotide polymorphism (SNP) marker maps, and show that LD signals between distant markers are due to the presence of extended haplotype superblocks in individuals with ancient chromosomes which have escaped historic recombination. ARs mediate the physiological effects of epinephrine and norepinephrine, and are targets of many therapeutic drugs. This work has identified haplotype backgrounds of alpha1-AR missense variants, haplotype block structures in US Caucasians and African Americans, and haplotype tag SNPs for each block, and we present strong evidence for ancient haplotype block superstructure at these genes which has been partially disrupted by recombination, and evidence for reinstatement of linkage disequilibrium by subsequent recombination events. ADRA1A is comprised of four haplotype blocks in US Caucasians, while in African Americans Block 1 is split. ADRA1B has four blocks in US Caucasians, but in African Americans only the first two blocks are present. ADRA1D has two blocks in US Caucasians, and the first block is replaced by two smaller blocks in African Americans. For both ADRA1A and ADRA1B, haplotype superstructures may represent a novel, higher-level hierarchy in the human genome, which may reduce redundancy of testing by further aggregation of genotype data.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by W. R. McCombie     

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site‐specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole‐genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree‐based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.     

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement‐specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.     

Contrasting population genetic structures using allozymes and the inversion polymorphism in Drosophila buzzatii

Second chromosome inversion and genotypic frequencies at seven allozyme loci, differentially associated with inversions, were determined in seven natural populations of Drosophila buzzatii. The patterns of variation of allozymes and the inversion polymorphisms were significantly different, indicating the role of adaptive differentiation for the latter. Moreover, the patterns of population structure varied among allozyme loci, suggesting the operation of diversifying selection for certain loci. Differentiation was negligible for Leucyl‐amino peptidase (Lap) and Peptidase‐2 (Pep‐2), low to moderate for Aldehyde oxidase (Aldox), Peptidase‐1 (Pep‐1) and Esterase‐1 (Est‐1) and high for Esterase‐2 (Est‐2) and Xanthine dehydrogenase (Xdh). Significant linkage disequilibria were detected between inversions and Aldox, Est‐1, Est‐2 and Xdh. Multiple regression analyses of inversion and allele frequencies on environmental variables revealed the existence of clines for inversions, Est‐1, Est‐2, Xdh and Aldox along altitudinal, latitudinal and/or climatic gradients. Tests using conditional allele frequencies showed that Est‐1 and Aldox clines could be accounted for by hitchhiking with inversions, whereas natural selection should be invoked to explain the clines observed for Est‐2 and Xdh.