Panel Of Recombinant Inbred Lines Research Articles

Towards mouse genetic-specific RNA-sequencing read mapping.

Genetic variations affect behavior and cause disease but understanding how these variants drive complex traits is still an open question. A common approach is to link the genetic variants to intermediate molecular phenotypes such as the transcriptome using RNA-sequencing (RNA-seq). Paradoxically, these variants between the samples are usually ignored at the beginning of RNA-seq analyses of many model organisms. This can skew the transcriptome estimates that are used later for downstream analyses, such as expression quantitative trait locus (eQTL) detection. Here, we assessed the impact of reference-based analysis on the transcriptome and eQTLs in a widely-used mouse genetic population: the BXD panel of recombinant inbred lines. We highlight existing reference bias in the transcriptome data analysis and propose practical solutions which combine available genetic variants, genotypes, and genome reference sequence. The use of custom BXD line references improved downstream analysis compared to classical genome reference. These insights would likely benefit genetic studies with a transcriptomic component and demonstrate that genome references need to be reassessed and improved.

Genetic mapping of novel modifiers for ApcMin induced intestinal polyps\u2019 development using the genetic architecture power of the collaborative cross mice

BackgroundFamilial adenomatous polyposis is an inherited genetic disease, characterized by colorectal polyps. It is caused by inactivating mutations in the Adenomatous polyposis coli (Apc) gene. Mice carrying a nonsense mutation in the Apc gene at R850, which is designated ApcMin/+ (Multiple intestinal neoplasia), develop intestinal adenomas. Several genetic modifier loci of Min (Mom) were previously mapped, but so far, most of the underlying genes have not been identified. To identify novel modifier loci associated with ApcMin/+, we performed quantitative trait loci (QTL) analysis for polyp development using 49 F1 crosses between different Collaborative Cross (CC) lines and C57BL/6 J-ApcMin/+mice. The CC population is a genetic reference panel of recombinant inbred lines, each line independently descended from eight genetically diverse founder strains. C57BL/6 J-ApcMin/+ males were mated with females from 49 CC lines. F1 offspring were terminated at 23 weeks and polyp counts from three sub-regions (SB1–3) of small intestinal and colon were recorded.ResultsThe number of polyps in all these sub-regions and colon varied significantly between the different CC lines. At 95% genome-wide significance, we mapped nine novel QTL for variation in polyp number, with distinct QTL associated with each intestinal sub-region. QTL confidence intervals varied in width between 2.63–17.79 Mb. We extracted all genes in the mapped QTL at 90 and 95% CI levels using the BioInfoMiner online platform to extract, significantly enriched pathways and key linker genes, that act as regulatory and orchestrators of the phenotypic landscape associated with the ApcMin/+ mutation.ConclusionsGenomic structure of the CC lines has allowed us to identify novel modifiers and confirmed some of the previously mapped modifiers. Key genes involved mainly in metabolic and immunological processes were identified. Future steps in this analysis will be to identify regulatory elements – and possible epistatic effects – located in the mapped QTL.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses.

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans.

Zinc is an essential trace element that acts as a co-factor for many enzymes and transcription factors required for cellular growth and development. Altering intracellular zinc levels can produce dramatic effects ranging from cell proliferation to cell death. To avoid such fates, cells have evolved mechanisms to handle both an excess and a deficiency of zinc. Zinc homeostasis is largely maintained via zinc transporters, permeable channels, and other zinc-binding proteins. Variation in these proteins might affect their ability to interact with zinc, leading to either increased sensitivity or resistance to natural zinc fluctuations in the environment. We can leverage the power of the roundworm nematode Caenorhabditis elegans as a tractable metazoan model for quantitative genetics to identify genes that could underlie variation in responses to zinc. We found that the laboratory-adapted strain (N2) is resistant and a natural isolate from Hawaii (CB4856) is sensitive to micromolar amounts of exogenous zinc supplementation. Using a panel of recombinant inbred lines, we identified two large-effect quantitative trait loci (QTL) on the left arm of chromosome III and the center of chromosome V that are associated with zinc responses. We validated and refined both QTL using near-isogenic lines (NILs) and identified a naturally occurring deletion in sqst-5, a sequestosome-related gene, that is associated with resistance to high exogenous zinc. We found that this deletion is relatively common across strains within the species and that variation in sqst-5 is associated with zinc resistance. Our results offer a possible mechanism for how organisms can respond to naturally high levels of zinc in the environment and how zinc homeostasis varies among individuals.

QTL mapping of natural variation reveals that the developmental regulator bruno reduces tolerance to P-element transposition in the Drosophila female germline.

Transposable elements (TEs) are obligate genetic parasites that propagate in host genomes by replicating in germline nuclei, thereby ensuring transmission to offspring. This selfish replication not only produces deleterious mutations—in extreme cases, TE mobilization induces genotoxic stress that prohibits the production of viable gametes. Host genomes could reduce these fitness effects in two ways: resistance and tolerance. Resistance to TE propagation is enacted by germline-specific small-RNA-mediated silencing pathways, such as the Piwi-interacting RNA (piRNA) pathway, and is studied extensively. However, it remains entirely unknown whether host genomes may also evolve tolerance by desensitizing gametogenesis to the harmful effects of TEs. In part, the absence of research on tolerance reflects a lack of opportunity, as small-RNA-mediated silencing evolves rapidly after a new TE invades, thereby masking existing variation in tolerance. We have exploited the recent historical invasion of the Drosophila melanogaster genome by P-element DNA transposons in order to study tolerance of TE activity. In the absence of piRNA-mediated silencing, the genotoxic stress imposed by P-elements disrupts oogenesis and, in extreme cases, leads to atrophied ovaries that completely lack germline cells. By performing quantitative trait locus (QTL) mapping on a panel of recombinant inbred lines (RILs) that lack piRNA-mediated silencing of P-elements, we uncovered multiple QTL that are associated with differences in tolerance of oogenesis to P-element transposition. We localized the most significant QTL to a small 230-kb euchromatic region, with the logarithm of the odds (LOD) peak occurring in the bruno locus, which codes for a critical and well-studied developmental regulator of oogenesis. Genetic, cytological, and expression analyses suggest that bruno dosage modulates germline stem cell (GSC) loss in the presence of P-element activity. Our observations reveal segregating variation in TE tolerance for the first time, and implicate gametogenic regulators as a source of tolerant variants in natural populations.

Circadian rhythms are associated with variation in photosystem II function and photoprotective mechanisms.

The circadian clock regulates many aspects of leaf gas supply and biochemical demand for CO2 , and is hypothesized to improve plant performance. Yet the extent to which the clock may regulate the efficiency of photosystem II (PSII) and photoprotective mechanisms such as heat dissipation is less explored. Based on measurements of chlorophyll a fluorescence, we estimated the maximum efficiency of PSII in light (Fv'/Fm') and heat dissipation by nonphotochemical quenching (NPQ). We further dissected total NPQ into its main components, qE (pH-dependent quenching), qT (state-transition quenching), and qI (quenching related to photoinhibition), in clock mutant genotypes of Arabidopsis thaliana, the cognate wild-type genotypes, and a panel of recombinant inbred lines expressing quantitative variation in clock period. Compared with mutants with altered clock function, we observed that wild-type genotypes with clock period lengths of approximately 24hr had both higher levels of Fv'/Fm', indicative of improved PSII function, and reduced NPQ, suggestive of lower stress on PSII light harvesting complexes. In the recombinant inbred lines, genetic variances were significant for Fv'/Fm' and all 3 components of NPQ, with qE explaining the greatest proportion of NPQ. Bivariate tests of association and structural equation models of hierarchical trait relationships showed that quantitative clock variation was empirically associated with Fv'/Fm' and NPQ, with qE mediating the relationship with gas exchange. The results demonstrate significant segregating variation for all photoprotective components, and suggest the adaptive significance of the clock may partly derive from its regulation of the light reactions of photosynthesis and of photoprotective mechanisms.

Determining Phenological Patterns Associated with the Onset of Senescence in a Wheat MAGIC Mapping Population.

The appropriate timing of developmental transitions is critical for adapting many crops to their local climatic conditions. Therefore, understanding the genetic basis of different aspects of phenology could be useful in highlighting mechanisms underpinning adaptation, with implications in breeding for climate change. For bread wheat (Triticum aestivum), the transition from vegetative to reproductive growth, the start and rate of leaf senescence and the relative timing of different stages of flowering and grain filling all contribute to plant performance. In this study we screened under Smart house conditions a large, multi-founder “NIAB elite MAGIC” wheat population, to evaluate the genetic elements that influence the timing of developmental stages in European elite varieties. This panel of recombinant inbred lines was derived from eight parents that are or recently have been grown commercially in the UK and Northern Europe. We undertook a detailed temporal phenotypic analysis under Smart house conditions of the population and its parents, to try to identify known or novel Quantitative Trait Loci associated with variation in the timing of key phenological stages in senescence. This analysis resulted in the detection of QTL interactions with novel traits such the time between “half of ear emergence above flag leaf ligule” and the onset of senescence at the flag leaf as well as traits associated with plant morphology such as stem height. In addition, strong correlations between several traits and the onset of senescence of the flag leaf were identified. This work establishes the value of systematically phenotyping genetically unstructured populations to reveal the genetic architecture underlying morphological variation in commercial wheat.

Identifying Loci Contributing to Natural Variation in Xenobiotic Resistance in Drosophila.

Natural populations exhibit a great deal of interindividual genetic variation in the response to toxins, exemplified by the variable clinical efficacy of pharmaceutical drugs in humans, and the evolution of pesticide resistant insects. Such variation can result from several phenomena, including variable metabolic detoxification of the xenobiotic, and differential sensitivity of the molecular target of the toxin. Our goal is to genetically dissect variation in the response to xenobiotics, and characterize naturally-segregating polymorphisms that modulate toxicity. Here, we use the Drosophila Synthetic Population Resource (DSPR), a multiparent advanced intercross panel of recombinant inbred lines, to identify QTL (Quantitative Trait Loci) underlying xenobiotic resistance, and employ caffeine as a model toxic compound. Phenotyping over 1,700 genotypes led to the identification of ten QTL, each explaining 4.5–14.4% of the broad-sense heritability for caffeine resistance. Four QTL harbor members of the cytochrome P450 family of detoxification enzymes, which represent strong a priori candidate genes. The case is especially strong for Cyp12d1, with multiple lines of evidence indicating the gene causally impacts caffeine resistance. Cyp12d1 is implicated by QTL mapped in both panels of DSPR RILs, is significantly upregulated in the presence of caffeine, and RNAi knockdown robustly decreases caffeine tolerance. Furthermore, copy number variation at Cyp12d1 is strongly associated with phenotype in the DSPR, with a trend in the same direction observed in the DGRP (Drosophila Genetic Reference Panel). No additional plausible causative polymorphisms were observed in a full genomewide association study in the DGRP, or in analyses restricted to QTL regions mapped in the DSPR. Just as in human populations, replicating modest-effect, naturally-segregating causative variants in an association study framework in flies will likely require very large sample sizes.

Using Drosophila melanogaster to identify chemotherapy toxicity genes.

The severity of the toxic side effects of chemotherapy shows a great deal of interindividual variability, and much of this variation is likely genetically based. Simple DNA tests predictive of toxic side effects could revolutionize the way chemotherapy is carried out. Due to the challenges in identifying polymorphisms that affect toxicity in humans, we use Drosophila fecundity following oral exposure to carboplatin, gemcitabine and mitomycin C as a model system to identify naturally occurring DNA variants predictive of toxicity. We use the Drosophila Synthetic Population Resource (DSPR), a panel of recombinant inbred lines derived from a multiparent advanced intercross, to map quantitative trait loci affecting chemotoxicity. We identify two QTL each for carboplatin and gemcitabine toxicity and none for mitomycin. One QTL is associated with fly orthologs of a priori human carboplatin candidate genes ABCC2 and MSH2, and a second QTL is associated with fly orthologs of human gemcitabine candidate genes RRM2 and RRM2B. The third, a carboplatin QTL, is associated with a posteriori human orthologs from solute carrier family 7A, INPP4A&B, and NALCN. The fourth, a gemcitabine QTL that also affects methotrexate toxicity, is associated with human ortholog GPx4. Mapped QTL each explain a significant fraction of variation in toxicity, yet individual SNPs and transposable elements in the candidate gene regions fail to singly explain QTL peaks. Furthermore, estimates of founder haplotype effects are consistent with genes harboring several segregating functional alleles. We find little evidence for nonsynonymous SNPs explaining mapped QTL; thus it seems likely that standing variation in toxicity is due to regulatory alleles.

A High-Resolution Genetic Map of Yellow Monkeyflower Identifies Chemical Defense QTLs and Recombination Rate Variation

Genotyping-by-sequencing methods have vastly improved the resolution and accuracy of genetic linkage maps by increasing both the number of marker loci as well as the number of individuals genotyped at these loci. Using restriction-associated DNA sequencing, we construct a dense linkage map for a panel of recombinant inbred lines derived from a cross between divergent ecotypes of Mimulus guttatus. We used this map to estimate recombination rate across the genome and to identify quantitative trait loci for the production of several secondary compounds (PPGs) of the phenylpropanoid pathway implicated in defense against herbivores. Levels of different PPGs are correlated across recombinant inbred lines suggesting joint regulation of the phenylpropanoid pathway. However, the three quantitative trait loci identified in this study each act on a distinct PPG. Finally, we map three putative genomic inversions differentiating the two parental populations, including a previously characterized inversion that contributes to life-history differences between the annual/perennial ecotypes.

The Genetic Architecture of Methotrexate Toxicity Is Similar in Drosophila melanogaster and Humans

The severity of the toxic side effects of chemotherapy varies among patients, and much of this variation is likely genetically based. Here, we use the model system Drosophila melanogaster to genetically dissect the toxicity of methotrexate (MTX), a drug used primarily to treat childhood acute lymphoblastic leukemia and rheumatoid arthritis. We use the Drosophila Synthetic Population Resource, a panel of recombinant inbred lines derived from a multiparent advanced intercross, and quantify MTX toxicity as a reduction in female fecundity. We identify three quantitative trait loci (QTL) affecting MTX toxicity; two colocalize with the fly orthologs of human genes believed to mediate MTX toxicity and one is a novel MTX toxicity gene with a human ortholog. A fourth suggestive QTL spans a centromere. Local single-marker association scans of candidate gene exons fail to implicate amino acid variants as the causative single-nucleotide polymorphisms, and we therefore hypothesize the causative variation is regulatory. In addition, the effects at our mapped QTL do not conform to a simple biallelic pattern, suggesting multiple causative factors underlie the QTL mapping results. Consistent with this observation, no single single-nucleotide polymorphism located in or near a candidate gene can explain the QTL mapping signal. Overall, our results validate D. melanogaster as a model for uncovering the genetic basis of chemotoxicity and suggest the genetic basis of MTX toxicity is due to a handful of genes each harboring multiple segregating regulatory factors.

Status and access to the Collaborative Cross population

The Collaborative Cross (CC) is a panel of recombinant inbred lines derived from eight genetically diverse laboratory inbred strains. Recently, the genetic architecture of the CC population was reported based on the genotype of a single male per line, and other publications reported incompletely inbred CC mice that have been used to map a variety of traits. The three breeding sites, in the US, Israel, and Australia, are actively collaborating to accelerate the inbreeding process through marker-assisted inbreeding and to expedite community access of CC lines deemed to have reached defined thresholds of inbreeding. Plans are now being developed to provide access to this novel genetic reference population through distribution centers. Here we provide a description of the distribution efforts by the University of North Carolina Systems Genetics Core, Tel Aviv University, Israel and the University of Western Australia.

HIGH-DIMENSIONAL VARIANCE PARTITIONING REVEALS THE MODULAR GENETIC BASIS OF ADAPTIVE DIVERGENCE IN GENE EXPRESSION DURING REPRODUCTIVE CHARACTER DISPLACEMENT

Although adaptive change is usually associated with complex changes in phenotype, few genetic investigations have been conducted on adaptations that involve sets of high-dimensional traits. Microarrays have supplied high-dimensional descriptions of gene expression, and phenotypic change resulting from adaptation often results in large-scale changes in gene expression. We demonstrate how genetic analysis of large-scale changes in gene expression generated during adaptation can be accomplished by determining high-dimensional variance partitioning within classical genetic experimental designs. A microarray experiment conducted on a panel of recombinant inbred lines (RILs) generated from two populations of Drosophila serrata that have diverged in response to natural selection, revealed genetic divergence in 10.6% of 3762 gene products examined. Over 97% of the genetic divergence in transcript abundance was explained by only 12 genetic modules. The two most important modules, explaining 50% of the genetic variance in transcript abundance, were genetically correlated with the morphological traits that are known to be under selection. The expression of three candidate genes from these two important genetic modules was assessed in an independent experiment using qRT-PCR on 430 individuals from the panel of RILs, and confirmed the genetic association between transcript abundance and morphological traits under selection.

Collaborative Cross mice and their power to map host susceptibility to Aspergillus fumigatus infection.

The Collaborative Cross (CC) is a genetic reference panel of recombinant inbred lines of mice, designed for the dissection of complex traits and gene networks. Each line is independently descended from eight genetically diverse founder strains such that the genomes of the CC lines, once fully inbred, are fine-grained homozygous mosaics of the founder haplotypes. We present an analysis of 120 CC lines, from a cohort of the CC bred at Tel Aviv University in collaboration with the University of Oxford, which at the time of this study were between the sixth and 12th generations of inbreeding and substantially homozygous at 170,000 SNPs. We show how CC genomes decompose into mosaics, and we identify loci that carry a deficiency or excess of a founder, many being deficient for the wild-derived strains WSB/EiJ and PWK/PhJ. We phenotyped 371 mice from 66 CC lines for a susceptibility to Aspergillus fumigatus infection. The survival time after infection varied significantly between CC lines. Quantitative trait locus (QTL) mapping identified genome-wide significant QTLs on chromosomes 2, 3, 8, 10 (two QTLs), 15, and 18. Simulations show that QTL mapping resolution (the median distance between the QTL peak and true location) varied between 0.47 and 1.18 Mb. Most of the QTLs involved contrasts between wild-derived founder strains and therefore would not segregate between classical inbred strains. Use of variation data from the genomes of the CC founder strains refined these QTLs further and suggested several candidate genes. These results support the use of the CC for dissecting complex traits.

Systems genetics of liver fibrogenesis in a murine reference population: Association mapping and pair scans in BXD recombinant inbred lines

Aims: Complex phenotypes are not only due to single locus variation, but are influenced by environmental covariates as well as epistatic interactions. Association mapping is based on phenotypic variations between inbred strains to identify genetic determinants of complex traits. Within the framework of the German Network of Systems Genetics (GeNeSys), we are investigating the BXD panel of recombinant inbred lines, which were generated by intercrossing C57BL/6J and DBA/2J mice, in a systems genetics approach to identify genetic loci and loci interactions that confer susceptibility to liver fibrosis.

Quantitative genetic analysis suggests causal association between cuticular hydrocarbon composition and desiccation survival in Drosophila melanogaster.

Survival to low relative humidity is a complex adaptation, and many repeated instances of evolution to desiccation have been observed among Drosophila populations and species. One general mechanism for desiccation resistance is Cuticular Hydrocarbon (CHC) melting point. We performed the first Quantitative Trait Locus (QTL) map of population level genetic variation in desiccation resistance in D. melanogaster. Using a panel of Recombinant Inbred Lines (RILs) derived from a single natural population, we mapped QTL in both sexes throughout the genome. We found that in both sexes, CHCs correlated strongly with desiccation resistance. At most desiccation resistance loci there was a significant association between CHCs and desiccation resistance of the sort predicted from clinal patterns of CHC variation and biochemical properties of lipids. This association was much stronger in females than males, perhaps because of greater overall abundance of CHCs in females, or due to correlations between CHCs used for waterproofing and sexual signalling in males. CHC evolution may be a common mechanism for desiccation resistance in D. melanogaster. It will be interesting to compare patterns of CHC variation and desiccation resistance in species which adapt to desiccation, and rainforest restricted species which cannot.

The Genetics of Phenotypic Plasticity in Plant Defense: Trichome Production inMimulus guttatus

Insect herbivory is a major driving force of plant evolution. Phenotypic plasticity and developmental variation provide a means for plants to cope with variable herbivory. We characterized the genetics of developmental variation and phenotypic plasticity in trichome density, a putative defensive trait of Mimulus guttatus (yellow monkeyflower). Our results are evaluated in relation to the optimal defense theory, which provides testable predictions for plastic and developmental patterns in defense traits. We found that both developmental stage and simulated insect damage affected trichome production, but in different ways. Plants were more likely to produce at least some trichomes on later leaves than on earlier leaves, regardless of damage. Damage did not affect the average probability of producing trichomes, but it did increase the density of hairs on trichome-positive plants. We mapped trichome quantitative trait loci (QTL) by selectively genotyping a large panel of recombinant inbred lines derived from two highly divergent populations. Several highly pleiotropic QTL influenced multiple aspects of the trichome phenotype (constitutive, developmental, and/or plastic responses). Only one of the QTL influenced trichome induction following damage. In a result that is consistent with a central prediction of optimal defense theory, the high allele at this location was from the ancestral population with low constitutive trichome production.

A Forward Genetics Approach to Identify Genetic Regulators of Liver Mineral Accumulation in Mice

We integrated ionomic profiling of mouse liver with quantitative trait loci (QTL) analysis using a panel of recombinant inbred lines (RIL) resulting from a cross of C57BL/6J (B6) and DBA/2J (DBA) inbred mouse lines (BXD). Following acid digestion, the levels of 13 elements (Na, Mg, P, K, Ca, Mn, Fe, Co, Cu, Zn, Se, Mo, and Cd) were determined by ICP‐MS in 400 samples from 27 RIL. Bayesian mapping using gender and analytical run as covariates and incorporating two stages of marker selection based on probabilities amongst prior distributions was used to identify QTL. 37 QTL were located and 12 strong peaks with a Bayes Factor (BF) >10 affected liver Fe (3 QTL), Mg (1), P (1), Mn (2), Cu (2), Mo (1), and Cd (2). Few of these loci contained genes known to regulate the mineral metabolism, demonstrating the value of this approach for gene discovery. The peaks controlling liver iron level on Chr 16, 17 and 18 each contain >100 candidate genes. The Chr 16 QTL peak (BF>4.5e+15) included the transferrin receptor (Tfrc) gene and polymorphisms between B6 and DBA mice were found within the 3′ UTR that normally regulates mRNA stability. However, transcript‐level QTL analysis revealed no difference in Tfrc mRNA levels between the BXD lines and suggests a Tfrc‐independent mechanism. Our data highlight the potential of coupling ionomics and forward genetics to identify natural genetic variants controlling mineral metabolism.

Bayesian Quantitative Trait Locus Mapping Using Inferred Haplotypes

We describe a fast hierarchical Bayesian method for mapping quantitative trait loci by haplotype-based association, applicable when haplotypes are not observed directly but are inferred from multiple marker genotypes. The method avoids the use of a Monte Carlo Markov chain by employing priors for which the likelihood factorizes completely. It is parameterized by a single hyperparameter, the fraction of variance explained by the quantitative trait locus, compared to the frequentist fixed-effects model, which requires a parameter for the phenotypic effect of each combination of haplotypes; nevertheless it still provides estimates of haplotype effects. We use simulation to show that the method matches the power of the frequentist regression model and, when the haplotypes are inferred, exceeds it for small QTL effect sizes. The Bayesian estimates of the haplotype effects are more accurate than the frequentist estimates, for both known and inferred haplotypes, which indicates that this advantage is independent of the effect of uncertainty in haplotype inference and will hold in comparison with frequentist methods in general. We apply the method to data from a panel of recombinant inbred lines of Arabidopsis thaliana, descended from 19 inbred founders.

Does segregating variation in sexual or microhabitat preferences lead to non-random mating within a population of Drosophila melanogaster ?

Variation in female choice for mates has implications for the maintenance of genetic variation and the evolution of male traits. Yet, estimates of population-level variation in male mating success owing to female genotype are rare. Here, we used a panel of recombinant inbred lines to estimate the strength of selection at many genetic loci in a single generation and attempt to assess differences between females with respect to the males they mated with. We performed selection assays in a complex environment to allow differences in habitat or social group preference to be expressed. We detected directional selection at loci across the genome, but are unable to provide support for differential male success because of variation in female genotype.