The X chromosome in quantitative trait locus mapping

QTLViewer: an interactive webtool for genetic analysis in the Collaborative Cross and Diversity Outbred mouse populations

The Collaborative Cross and the Diversity Outbred mouse populations are related multiparental populations, derived from the same 8 isogenic founder strains. They carry >50 M known genetic variants, which makes them ideal tools for mapping genetic loci that regulate phenotypes, including physiological and molecular traits. Mapping quantitative trait loci requires statistical and computational training, which can present a barrier to access for some researchers. The QTLViewer software allows users to graphically explore Collaborative Cross and Diversity Outbred quantitative trait locus mapping and related analyses performed through the R/qtl2 package. Additionally, the QTLViewer website serves as a repository for published Collaborative Cross and Diversity Outbred studies, increasing the accessibility of these genetic resources to the broader scientific community.

Genetic Evidence for a Causal Relationship between Hyperlipidemia and Type 2 Diabetes in Mice

Dyslipidemia is considered a risk factor for type 2 diabetes (T2D), yet studies with statins and candidate genes suggest that circulating lipids may protect against T2D development. Apoe-null (Apoe^-/-) mouse strains develop spontaneous dyslipidemia and exhibit a wide variation in susceptibility to diet-induced T2D. We thus used Apoe^-/- mice to elucidate phenotypic and genetic relationships of circulating lipids with T2D. A male F2 cohort was generated from an intercross between LP/J and BALB/cJ Apoe^-/- mice and fed 12 weeks of a Western diet. Fasting, non-fasting plasma glucose, and lipid levels were measured and genotyping was performed using miniMUGA arrays. We uncovered a major QTL near 60 Mb on chromosome 15, Nhdlq18, which affected non-HDL cholesterol and triglyceride levels under both fasting and non-fasting states. This QTL was coincident with Bglu20, a QTL that modulates fasting and non-fasting glucose levels. The plasma levels of non-HDL cholesterol and triglycerides were closely correlated with the plasma glucose levels in F2 mice. Bglu20 disappeared after adjustment for non-HDL cholesterol or triglycerides. These results demonstrate a causative role for dyslipidemia in T2D development in mice.

Quantitative trait loci for growth-related traits in Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study aimed to identify quantitative trait loci (QTLs) for growth-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers in Japanese quail. A QTL mapping population of 277 F₂ birds was obtained from an intercross between a male of a large-sized strain and three females of a normal-sized strain. Body weight (BW) was measured weekly from hatching to 16 weeks of age. Non-linear regression growth models of Weibull, Logistic, Gompertz, Richards, and Brody were analyzed, and growth curve parameters of Richards was selected as the best model to describe the quail growth curve of the F₂ birds. Restriction-site associated DNA sequencing developed 125 SNP markers that were informative between their parental strains. The SNP markers were distributed on 16 linkage groups that spanned 795.9 centiMorgan (cM) with an average marker interval of 7.3 cM. QTL analysis of phenotypic traits revealed four main-effect QTLs. Detected QTLs were located on chromosomes 1 and 3 and were associated with BW from 4 to 16 weeks of age and asymptotic weight of Richards model at genome-wide significant at 1% or 5% level. No QTL was detected for BW from 0 to 3 weeks of age. This is the first report identified QTLs for asymptotic weight of the Richards parameter in Japanese quail. These results highlight that the combination of QTL studies and the RAD-seq method will aid future breeding programs identify genes underlying the QTL and the application of marker-assisted selection in the poultry industry, particularly the Japanese quail.

Mapping of Quantitative Trait Loci Controlling Egg-Quality and -Production Traits in Japanese Quail (Coturnix japonica) Using Restriction-Site Associated DNA Sequencing

This research was conducted to identify quantitative trait loci (QTL) associated with egg-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers using restriction-site associated DNA sequencing (RAD-seq) in Japanese quail. A total of 138 F₂ females were produced by full-sib mating of F₁ birds derived from an intercross between a male of the large-sized strain with three females of the normal-sized strain. Eggs were investigated at two different stages: the beginning stage of egg-laying and at 12 weeks of age (second stage). Five eggs were analyzed for egg weight, lengths of the long and short axes, egg shell strength and weight, yolk weight and diameter, albumen weight, egg equator thickness, and yolk color (L*, a*, and b* values) at each stage. Moreover, the age at first egg, the cumulative number of eggs laid, and egg production rate were recorded. RAD-seq developed 118 SNP markers and mapped them to 13 linkage groups using the Map Manager QTX b20 software. Markers were spanned on 776.1 cM with an average spacing of 7.4 cM. Nine QTL were identified on chromosomes 2, 4, 6, 10, 12, and Z using the simple interval mapping method in the R/qtl package. The QTL detected affected 10 egg traits of egg weight, lengths of the long and short axes of egg, egg shell strength, yolk diameter and weight, albumen weight, and egg shell weight at the beginning stage, yellowness of the yolk color at the second stage, and age at first egg. This is the first report to perform a quail QTL analysis of egg-related traits using RAD-seq. These results highlight the effectiveness of RAD-seq associated with targeted QTL and the application of marker-assisted selection in the poultry industry, particularly in the Japanese quail.

Treatment of the X chromosome in mapping multiple quantitative trait loci

Statistical methods to map quantitative trait loci (QTL) often neglect the X chromosome and may focus exclusively on autosomal loci. But the X chromosome often requires special treatment: sex and cross-direction covariates may need to be included to avoid spurious evidence of linkage, and the X chromosome may require a separate significance threshold. In multiple-QTL analyses, including the consideration of epistatic interactions, the X chromosome also requires special care and consideration. We extend a penalized likelihood method for multiple-QTL model selection, to appropriately handle the X chromosome. We examine its performance in simulation and by application to a large eQTL data set. The method has been implemented in the package R/qtl.

Neighbor QTL: an interval mapping method for quantitative trait loci underlying plant neighborhood effects

Phenotypes of sessile organisms, such as plants, rely not only on their own genotypes but also on those of neighboring individuals. Previously, we incorporated such neighbor effects into a single-marker regression using the Ising model of ferromagnetism. However, little is known regarding how neighbor effects should be incorporated in quantitative trait locus (QTL) mapping. In this study, we propose a new method for interval QTL mapping of neighbor effects, designated "neighbor QTL," the algorithm of which includes: (1) obtaining conditional self-genotype probabilities with recombination fraction between flanking markers; (2) calculating conditional neighbor genotypic identity using the self-genotype probabilities; and (3) estimating additive and dominance deviations for neighbor effects. Our simulation using F2 and backcross lines showed that the power to detect neighbor effects increased as the effective range decreased. The neighbor QTL was applied to insect herbivory on Col × Kas recombinant inbred lines of Arabidopsis thaliana. Consistent with previous results, the pilot experiment detected a self-QTL effect on the herbivory at the GLABRA1 locus. Regarding neighbor QTL effects on herbivory, we observed a weak QTL on the top of chromosome 4, at which a weak self-bolting QTL was also identified. The neighbor QTL method is available as an R package (https://cran.r-project.org/package=rNeighborQTL), providing a novel tool to investigate neighbor effects in QTL studies.

Genomic landscape of reproductive isolation in Lucania killifish: The role of sex loci and salinity

Adaptation to different environments can directly and indirectly generate reproductive isolation between species. Bluefin killifish (Lucania goodei) and rainwater killifish (L. parva) are sister species that have diverged across a salinity gradient and are reproductively isolated by habitat, behavioural, extrinsic and intrinsic post-zygotic isolation. We asked if salinity adaptation contributes indirectly to other forms of reproductive isolation via linked selection and hypothesized that low recombination regions, such as sex chromosomes or chromosomal rearrangements, might facilitate this process. We conducted QTL mapping in backcrosses between L. parva and L. goodei to explore the genetic architecture of salinity tolerance, behavioural isolation and intrinsic isolation. We mapped traits relative to a chromosome that has undergone a centric fusion in L. parva (relative to L. goodei). We found that the sex locus appears to be male determining (XX-XY), was located on the fused chromosome and was implicated in intrinsic isolation. QTL associated with salinity tolerance were spread across the genome and did not overly co-localize with regions associated with behavioural or intrinsic isolation. This preliminary analysis of the genetic architecture of reproductive isolation between Lucania species does not support the hypothesis that divergent natural selection for salinity tolerance led to behavioural and intrinsic isolation as a by-product. Combined with previous studies in this system, our work suggests that adaptation as a function of salinity contributes to habitat isolation and that reinforcement may have contributed to the evolution of behavioural isolation instead, possibly facilitated by linkage between behavioural isolation and intrinsic isolation loci on the fused chromosome.

Naive Pluripotent Stem Cells Exhibit Phenotypic Variability that Is Driven by Genetic Variation

Variability among pluripotent stem cell (PSC) lines is a prevailing issue that hampers not only experimental reproducibility but also large-scale applications and personalized cell-based therapy. This variability could result from epigenetic and genetic factors that influence stem cell behavior. Naive culture conditions minimize epigenetic fluctuation, potentially overcoming differences in PSC line differentiation potential. Here we derived PSCs from distinct mouse strains under naive conditions and show that lines from distinct genetic backgrounds have divergent differentiation capacity, confirming a major role for genetics in PSC phenotypic variability. This is explained in part through inconsistent activity of extra-cellular signaling, including the Wnt pathway, which is modulated by specific genetic variants. Overall, this study shows that genetic background plays a dominant role in driving phenotypic variability of PSCs.

Synchronization of senescence and desynchronization of flowering in Arabidopsis thaliana

In a recent publication, we proposed that adjusting lifespan in order to synchronize senescence is important for timing of reproduction, and we quantified the synchrony of reproductive timing relative to germination timing. Here, in a second sequential seeding experiment (SSE), the germination timing of Arabidopsis thaliana accessions was manipulated and plants were then grown under two different temperature regimes. Life stage traits of plants in each temperature regime were analysed and it was evaluated whether the cohorts were grouped according to age and/or environmental conditions. While flowering-related traits showed desynchrony among cohorts, striking synchrony in the timing of senescence among cohorts for each group was found. A quantitative trait locus (QTL) analysis using a genotyped population of 'Cvi/Ler' recombinant inbred lines (RILs) was then conducted. Novel and known loci were assigned to flowering and senescence timing. However, senescence synchrony resulted in low variation in senescence time and weak QTL detection for flowering termination. Overlapping flowering and senescence genes with loci affecting either of those traits were found and suggest a potential interdependency of reproductive traits.

Quantitative trait loci mapping for the shear force value in breast muscle of F2 chickens

The shear force value is one of the major traits that determine meat quality. In the present study, we performed QTL analysis for chicken breast muscle shear force value at 7 wk of age using 545 single nucleotide polymorphism (SNP) markers developed via restriction-site associated DNA sequencing (RAD-seq). An F2 resource family was generated by mating Oh-Shamo, a native Japanese chicken breed, and the White Plymouth Rock chicken breed. A total of 215 F2 birds were produced. Simple interval mapping revealed one significant main-effect QTL between 6.28 and 8.10 Mb SNPs on the chromosome Z with a logarithm of odds score of 5.53 at the genome-wide 5% level. At this QTL, the confidence interval, phenotypic variance explained, and additive effect were 26 cM, 12.24%, and -0.31 in males and -0.34 in females, respectively. No QTL with epistatic interaction effects were detected. To our knowledge, this is the first report on a QTL affecting the shear force value in the chicken breast muscle, using SNP markers derived from RAD-seq.

Replication and discovery of musculoskeletal QTLs in LG/J and SM/J advanced intercross lines

The genetics underlying variation in health‐related musculoskeletal phenotypes can be investigated in a mouse model. Quantitative trait loci (QTLs) affecting musculoskeletal traits in the LG/J and SM/J strain lineage remain to be refined and corroborated. The aim of this study was to map muscle and bone traits in males (n = 506) of the 50th filial generation of advanced intercross lines (LG/SM AIL) derived from the two strains. Genetic contribution to variation in all musculoskeletal traits was confirmed; the SNP heritability of muscle mass ranged between 0.46 and 0.56; and the SNP heritability of tibia length was 0.40. We used two analytical software, GEMMA and QTLRel, to map the underlying QTLs. GEMMA required substantially less computation and recovered all the QTLs identified by QTLRel. Seven significant QTLs were identified for muscle weight (Chr 1, 7, 11, 12, 13, 15, and 16), and two for tibia length, (Chr 1 and 13). Each QTL explained 4–5% of phenotypic variation. One muscle and both bone loci replicated previous findings; the remaining six were novel. Positional candidates for the replicated QTLs were prioritized based on in silico analyses and gene expression in muscle tissue. In summary, we replicated existing QTLs and identified novel QTLs affecting muscle weight, and replicated bone length QTLs in LG/SM AIL males. Heritability estimates substantially exceed the cumulative effect of the QTLs, hence a richer genetic architecture contributing to muscle and bone variability could be uncovered with a larger sample size.

Mapping of Quantitative Trait Loci for Growth and Carcass-Related Traits in Chickens Using a Restriction-Site Associated DNA Sequencing Method

In the present study, quantitative trait loci (QTLs) analysis was performed to identify the chromosomal positions of growth and carcass-related trait QTLs using 319 F₂ chickens obtained from intercrosses of an Oh-Shamo male and four White Plymouth Rock females. Body weight was measured weekly until the birds were 7 weeks old. Carcass-related traits were also measured at this timepoint. A genetic linkage map was constructed using 545 single nucleotide polymorphism (SNP) markers that were developed using a restriction-site associated DNA sequencing method. The linkage map included the 23 autosomes and the Z chromosome. Using simple interval QTL mapping, we were able to identify 10 significant and suggestive main-effect QTLs for growth and carcass-related traits present on chromosomes 1, 2, 3, 5, 8, 19, 24, and Z. These loci explained 5.60–16.52% of the phenotypic variances. The chromosomal positions of the 10 QTLs overlapped with those of previously reported QTLs, whereas the targeted traits varied. Our QTLs will aid future breeding programs in improving growth and meat yield of chickens (e.g., via marker-assisted selection), particularly in the Japanese brand chicken industry.

R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely-used R/qtl software package to include multiparental populations, better handles modern high-dimensional data....

R/qtl2 is an interactive software environment for mapping quantitative trait loci (QTL) in experimental populations. The R/qtl2 software expands the scope of the widely used R/qtl software package to include multiparent populations derived from more than two founder strains, such as the Collaborative Cross and Diversity Outbred mice, heterogeneous stocks, and MAGIC plant populations. R/qtl2 is designed to handle modern high-density genotyping data and high-dimensional molecular phenotypes, including gene expression and proteomics. R/qtl2 includes the ability to perform genome scans using a linear mixed model to account for population structure, and also includes features to impute SNPs based on founder strain genomes and to carry out association mapping. The R/qtl2 software provides all of the basic features needed for QTL mapping, including graphical displays and summary reports, and it can be extended through the creation of add-on packages. R/qtl2, which is free and open source software written in the R and C++ programming languages, comes with a test framework.

Genetic Dissection of Hybrid Male Sterility Across Stages of Spermatogenesis

Hybrid sterility is a common form of reproductive isolation between nascent species. Although hybrid sterility is routinely documented and genetically dissected in speciation studies, its developmental basis is rarely examined, especially in generations beyond the F1 generation. To identify phenotypic and genetic determinants of hybrid male sterility from a developmental perspective, we characterized testis histology in 312 F2 hybrids generated by intercrossing inbred strains of Mus musculus domesticus and M. m. musculus, two subspecies of house mice. Hybrids display a range of histologic abnormalities that indicate defective spermatogenesis. Among these abnormalities, we quantified decreased testis size, reductions in spermatocyte and spermatid number, increased apoptosis of meiosis I spermatocytes, and more multinucleated syncytia. Collectively, our phenotypic data point to defects in meiosis I as a primary barrier to reproduction. We identified seven quantitative trait loci (QTL) controlling five histologic traits. A region of chromosome 17 that contains Prdm9, a gene known to confer F1 hybrid male sterility, affects multinucleated syncytia and round spermatids, potentially extending the phenotypic outcomes of this incompatibility. The X chromosome also plays a key role, with loci affecting multinucleated syncytia, apoptosis of round spermatids, and round spermatid numbers. We detected an epistatic interaction between QTL on chromosomes 17 and X for multinucleated syncytia. Our results refine the developmental basis of a key reproductive barrier in a classic model system for speciation genetics.

Genetic Mapping of Quantitative Trait Loci for Egg Production and Egg Quality Traits in Chickens: a Review

Chickens display a wide spectrum of phenotypic variations in quantitative traits such as egg-related traits. Quantitative trait locus (QTL) analysis is a statistical method used to understand the relationship between phenotypic (trait measurements) and genotypic data (molecular markers). We have performed QTL analyses for egg-related traits using an original resource population based on the Japanese Large Game (Oh-Shamo) and the White Leghorn breeds of chickens. In this article, we summarize the results of our extensive QTL analyses for 11 and 66 traits for egg production and egg quality, respectively. We reveal that at least 30 QTL regions on 17 different chromosomes affect phenotypic variation in egg-related traits. Each locus had an age-specific effect on traits, and a variety in effects was also apparent, such as additive, dominance, and epistatic-interaction effects. Although genome-wide association study (GWAS) is suitable for gene-level resolution mapping of GWAS loci with additive effects, QTL mapping studies enable us to comprehensively understand genetic control, such as chromosomal regions, genetic contribution to phenotypic variance, mode of inheritance, and age-specificity of both common and rare alleles. QTL analyses also describe the relationship between genotypes and phenotypes in experimental populations. Accumulation of QTL information, including GWAS loci, is also useful for studies of population genomics approached without phenotypic data in order to validate the identified genomic signatures of positive selection. The combination of QTL studies and next-generation sequencing techniques with uncharacterized genetic resources will enhance current understanding of the relationship between genotypes and phenotypes in livestock animals.

Natural Variation in Brachypodium Links Vernalization and Flowering Time Loci as Major Flowering Determinants

The domestication of plants is underscored by the selection of agriculturally favorable developmental traits, including flowering time, which resulted in the creation of varieties with altered growth habits. Research into the pathways underlying these growth habits in cereals has highlighted the role of three main flowering regulators: VERNALIZATION1 (VRN1), VRN2, and FLOWERING LOCUS T (FT). Previous reverse genetic studies suggested that the roles of VRN1 and FT are conserved in Brachypodium distachyon yet identified considerable ambiguity surrounding the role of VRN2 To investigate the natural diversity governing flowering time pathways in a nondomesticated grass, the reference B. distachyon accession Bd21 was crossed with the vernalization-dependent accession ABR6. Resequencing of ABR6 allowed the creation of a single-nucleotide polymorphism-based genetic map at the F4 stage of the mapping population. Flowering time was evaluated in F4:5 families in five environmental conditions, and three major loci were found to govern flowering time. Interestingly, two of these loci colocalize with the B. distachyon homologs of the major flowering pathway genes VRN2 and FT, whereas no linkage was observed at VRN1 Characterization of these candidates identified sequence and expression variation between the two parental genotypes, which may explain the contrasting growth habits. However, the identification of additional quantitative trait loci suggests that greater complexity underlies flowering time in this nondomesticated system. Studying the interaction of these regulators in B. distachyon provides insights into the evolutionary context of flowering time regulation in the Poaceae as well as elucidates the way humans have utilized the natural variation present in grasses to create modern temperate cereals.

Genetics of Skeletal Evolution in Unusually Large Mice from Gough Island

Organisms on islands often undergo rapid morphological evolution, providing a platform for understanding mechanisms of phenotypic change. Many examples of evolution on islands involve the vertebrate skeleton. Although the genetic basis of skeletal variation has been studied in laboratory strains, especially in the house mouse Mus musculus domesticus, the genetic determinants of skeletal evolution in natural populations remain poorly understood. We used house mice living on the remote Gough Island-the largest wild house mice on record-to understand the genetics of rapid skeletal evolution in nature. Compared to a mainland reference strain from the same subspecies (WSB/EiJ), the skeleton of Gough Island mice is considerably larger, with notable expansions of the pelvis and limbs. The Gough Island mouse skeleton also displays changes in shape, including elongations of the skull and the proximal vs. distal elements in the limbs. Quantitative trait locus (QTL) mapping in a large F₂ intercross between Gough Island mice and WSB/EiJ reveals hundreds of QTL that control skeletal dimensions measured at 5, 10, and/or 16 weeks of age. QTL exhibit modest, mostly additive effects, and Gough Island alleles are associated with larger skeletal size at most QTL. The QTL with the largest effects are found on a few chromosomes and affect suites of skeletal traits. Many of these loci also colocalize with QTL for body weight. The high degree of QTL colocalization is consistent with an important contribution of pleiotropy to skeletal evolution. Our results provide a rare portrait of the genetic basis of skeletal evolution in an island population and position the Gough Island mouse as a model system for understanding mechanisms of rapid evolution in nature.

Mapping QTL Contributing to Variation in Posterior Lobe Morphology between Strains of Drosophila melanogaster

Closely-related, and otherwise morphologically similar insect species frequently show striking divergence in the shape and/or size of male genital structures, a phenomenon thought to be driven by sexual selection. Comparative interspecific studies can help elucidate the evolutionary forces acting on genital structures to drive this rapid differentiation. However, genetic dissection of sexual trait divergence between species is frequently hampered by the difficulty generating interspecific recombinants. Intraspecific variation can be leveraged to investigate the genetics of rapidly-evolving sexual traits, and here we carry out a genetic analysis of variation in the posterior lobe within D. melanogaster. The lobe is a male-specific process emerging from the genital arch of D. melanogaster and three closely-related species, is essential for copulation, and shows radical divergence in form across species. There is also abundant variation within species in the shape and size of the lobe, and while this variation is considerably more subtle than that seen among species, it nonetheless provides the raw material for QTL mapping. We created an advanced intercross population from a pair of phenotypically-different inbred strains, and after phenotyping and genotyping-by-sequencing the recombinants, mapped several QTL contributing to various measures of lobe morphology. The additional generations of crossing over in our mapping population led to QTL intervals that are smaller than is typical for an F₂ mapping design. The intervals we map overlap with a pair of lobe QTL we previously identified in an independent mapping cross, potentially suggesting a level of shared genetic control of trait variation. Our QTL additionally implicate a suite of genes that have been shown to contribute to the development of the posterior lobe. These loci are strong candidates to harbor naturally-segregating sites contributing to phenotypic variation within D. melanogaster, and may also be those contributing to divergence in lobe morphology between species.

The importance of having two X chromosomes

Historically, it was thought that the number of X chromosomes plays little role in causing sex differences in traits. Recently, selected mouse models have been used increasingly to compare mice with the same type of gonad but with one versus two copies of the X chromosome. Study of these models demonstrates that mice with one X chromosome can be strikingly different from those with two X chromosomes, when the differences are not attributable to confounding group differences in gonadal hormones. The number of X chromosomes affects adiposity and metabolic disease, cardiovascular ischaemia/reperfusion injury and behaviour. The effects of X chromosome number are likely the result of inherent differences in expression of X genes that escape inactivation, and are therefore expressed from both X chromosomes in XX mice, resulting in a higher level of expression when two X chromosomes are present. The effects of X chromosome number contribute to sex differences in disease phenotypes, and may explain some features of X chromosome aneuploidies such as in Turner and Klinefelter syndromes.

Integration of genome-wide association and extant brain expression QTL identifies candidate genes influencing prepulse inhibition in inbred F1 mice

Genetic association mapping in structured populations of model organisms can offer a fruitful complement to human genetic studies by generating new biological hypotheses about complex traits. Here we investigated prepulse inhibition (PPI), a measure of sensorimotor gating that is disrupted in a number of psychiatric disorders. To identify genes that influence PPI, we constructed a panel of half-sibs by crossing 30 females from common inbred mouse strains with inbred C57BL/6J males to create male and female F1 offspring. We used publicly available single nucleotide polymorphism (SNP) genotype data from these inbred strains to perform a genome-wide association scan using a dense panel of over 150,000 SNPs in a combined sample of 604 mice representing 30 distinct F1 genotypes. We identified two independent PPI-associated loci on Chromosomes 2 and 7, each of which explained 12-14% of the variance in PPI. Searches of available databases did not identify any plausible causative coding polymorphisms within these loci. However, previously collected expression quantitative trait locus (eQTL) data from hippocampus and striatum indicated that the SNPs on Chromosomes 2 and 7 that showed the strongest association with PPI were also strongly associated with expression of several transcripts, some of which have been implicated in human psychiatric disorders. This integrative approach successfully identified a focused set of genes which can be prioritized for follow-up studies. More broadly, our results show that F1 crosses among common inbred strains can be used in combination with other informatics and expression datasets to identify candidate genes for complex behavioral traits.

Genetics of Rapid and Extreme Size Evolution in Island Mice

Organisms on islands provide a revealing window into the process of adaptation. Populations that colonize islands often evolve substantial differences in body size from their mainland relatives. Although the ecological drivers of this phenomenon have received considerable attention, its genetic basis remains poorly understood. We use house mice (subspecies: Mus musculus domesticus) from remote Gough Island to provide a genetic portrait of rapid and extreme size evolution. In just a few hundred generations, Gough Island mice evolved the largest body size among wild house mice from around the world. Through comparisons with a smaller-bodied wild-derived strain from the same subspecies (WSB/EiJ), we demonstrate that Gough Island mice achieve their exceptional body weight primarily by growing faster during the 6 weeks after birth. We use genetic mapping in large F(2) intercrosses between Gough Island mice and WSB/EiJ to identify 19 quantitative trait loci (QTL) responsible for the evolution of 16-week weight trajectories: 8 QTL for body weight and 11 QTL for growth rate. QTL exhibit modest effects that are mostly additive. We conclude that body size evolution on islands can be genetically complex, even when substantial size changes occur rapidly. In comparisons to published studies of laboratory strains of mice that were artificially selected for divergent body sizes, we discover that the overall genetic profile of size evolution in nature and in the laboratory is similar, but many contributing loci are distinct. Our results underscore the power of genetically characterizing the entire growth trajectory in wild populations and lay the foundation necessary for identifying the mutations responsible for extreme body size evolution in nature.

Ontogenetic stage-specific quantitative trait loci contribute to divergence in developmental trajectories of sexually dimorphic fins between medaka populations

Sexual dimorphism can evolve when males and females differ in phenotypic optima. Genetic constraints can, however, limit the evolution of sexual dimorphism. One possible constraint is derived from alleles expressed in both sexes. Because males and females share most of their genome, shared alleles with different fitness effects between sexes are faced with intralocus sexual conflict. Another potential constraint is derived from genetic correlations between developmental stages. Sexually dimorphic traits are often favoured at adult stages, but selected against as juvenile, so developmental decoupling of traits between ontogenetic stages may be necessary for the evolution of sexual dimorphism in adults. Resolving intralocus conflicts between sexes and ages is therefore a key to the evolution of age-specific expression of sexual dimorphism. We investigated the genetic architecture of divergence in the ontogeny of sexual dimorphism between two populations of the Japanese medaka (Oryzias latipes) that differ in the magnitude of dimorphism in anal and dorsal fin length. Quantitative trait loci (QTL) mapping revealed that few QTL had consistent effects throughout ontogenetic stages and the majority of QTL change the sizes and directions of effects on fin growth rates during ontogeny. We also found that most QTL were sex-specific, suggesting that intralocus sexual conflict is almost resolved. Our results indicate that sex- and age-specific QTL enable the populations to achieve optimal developmental trajectories of sexually dimorphic traits in response to complex natural and sexual selection.

Quantitative trait locus mapping methods for diversity outbred mice

Genetic mapping studies in the mouse and other model organisms are used to search for genes underlying complex phenotypes. Traditional genetic mapping studies that employ single-generation crosses have poor mapping resolution and limit discovery to loci that are polymorphic between the two parental strains. Multiparent outbreeding populations address these shortcomings by increasing the density of recombination events and introducing allelic variants from multiple founder strains. However, multiparent crosses present new analytical challenges and require specialized software to take full advantage of these benefits. Each animal in an outbreeding population is genetically unique and must be genotyped using a high-density marker set; regression models for mapping must accommodate multiple founder alleles, and complex breeding designs give rise to polygenic covariance among related animals that must be accounted for in mapping analysis. The Diversity Outbred (DO) mice combine the genetic diversity of eight founder strains in a multigenerational breeding design that has been maintained for >16 generations. The large population size and randomized mating ensure the long-term genetic stability of this population. We present a complete analytical pipeline for genetic mapping in DO mice, including algorithms for probabilistic reconstruction of founder haplotypes from genotyping array intensity data, and mapping methods that accommodate multiple founder haplotypes and account for relatedness among animals. Power analysis suggests that studies with as few as 200 DO mice can detect loci with large effects, but loci that account for <5% of trait variance may require a sample size of up to 1000 animals. The methods described here are implemented in the freely available R package DOQTL.

Genomic networks of hybrid sterility

Hybrid dysfunction, a common feature of reproductive barriers between species, is often caused by negative epistasis between loci ("Dobzhansky-Muller incompatibilities"). The nature and complexity of hybrid incompatibilities remain poorly understood because identifying interacting loci that affect complex phenotypes is difficult. With subspecies in the early stages of speciation, an array of genetic tools, and detailed knowledge of reproductive biology, house mice (Mus musculus) provide a model system for dissecting hybrid incompatibilities. Male hybrids between M. musculus subspecies often show reduced fertility. Previous studies identified loci and several X chromosome-autosome interactions that contribute to sterility. To characterize the genetic basis of hybrid sterility in detail, we used a systems genetics approach, integrating mapping of gene expression traits with sterility phenotypes and QTL. We measured genome-wide testis expression in 305 male F2s from a cross between wild-derived inbred strains of M. musculus musculus and M. m. domesticus. We identified several thousand cis- and trans-acting QTL contributing to expression variation (eQTL). Many trans eQTL cluster into eleven 'hotspots,' seven of which co-localize with QTL for sterility phenotypes identified in the cross. The number and clustering of trans eQTL-but not cis eQTL-were substantially lower when mapping was restricted to a 'fertile' subset of mice, providing evidence that trans eQTL hotspots are related to sterility. Functional annotation of transcripts with eQTL provides insights into the biological processes disrupted by sterility loci and guides prioritization of candidate genes. Using a conditional mapping approach, we identified eQTL dependent on interactions between loci, revealing a complex system of epistasis. Our results illuminate established patterns, including the role of the X chromosome in hybrid sterility. The integrated mapping approach we employed is applicable in a broad range of organisms and we advocate for widespread adoption of a network-centered approach in speciation genetics.

Genetic mapping of two components of reproductive isolation between two sibling species of moths, Ostrinia nubilalis and O. scapulalis

We report the quantitative trait loci (QTL) mapping of reproductive isolation traits between Ostrinia nubilalis (the European corn borer) and its sibling species O. scapulalis (the Adzuki bean borer), focusing on two traits: mating isolation (mi) and pheromone production (Pher). Four genetic maps were generated from two backcross families, with two maps (one chromosomal map and one linkage map) per backcross. We located 165-323 AFLP markers on these four maps, resulting in the identification of 27-31 linkage groups, depending on the map considered. No-choice mating experiments with the offspring of each backcross led to the detection of at least two QTLs for mi in different linkage groups. QTLs underlying Pher were located in a third linkage group. The Z heterochromosome was identified by a specific marker (Tpi) and did not carry any of these QTLs. Finally, we considered the global divergence between the two sibling species, distortions of segregation throughout the genome, and the location and effect of mi and Pher QTLs in light of the known candidate genes for reproductive isolation within the genus Ostrinia and, more broadly, in phytophagous insects.

Identification of a novel polymorphism in X-linked sterol-4-alpha-carboxylate 3-dehydrogenase (Nsdhl) associated with reduced high-density lipoprotein cholesterol levels in I/LnJ mice

Loci controlling plasma lipid concentrations were identified by performing a quantitative trait locus analysis on genotypes from 233 mice from a F2 cross between KK/HlJ and I/LnJ, two strains known to differ in their high-density lipoprotein (HDL) cholesterol levels. When fed a standard diet, HDL cholesterol concentration was affected by two significant loci, the Apoa2 locus on Chromosome (Chr) 1 and a novel locus on Chr X, along with one suggestive locus on Chr 6. Non-HDL concentration also was affected by loci on Chr 1 and X along with a suggestive locus on Chr 3. Additional loci that may be sex-specific were identified for HDL cholesterol on Chr 2, 3, and 4 and for non-HDL cholesterol on Chr 5, 7, and 14. Further investigation into the potential causative gene on Chr X for reduced HDL cholesterol levels revealed a novel, I/LnJ-specific nonsynonymous polymorphism in Nsdhl, which codes for sterol-4-alpha-carboxylate 3-dehydrogenase in the cholesterol synthesis pathway. Although many lipid quantitative trait locus have been reported previously, these data suggest there are additional genes left to be identified that control lipid levels and that can provide new pharmaceutical targets.

A major X-linked locus affects kidney function in mice

Chronic kidney disease is a common disease with increasing prevalence in the western population. One common reason for chronic kidney failure is diabetic nephropathy. Diabetic nephropathy and hyperglycemia are characteristics of the mouse inbred strain KK/HlJ, which is predominantly used as a model for metabolic syndrome due to its inherited glucose intolerance and insulin resistance. We used KK/HlJ, an albuminuria-sensitive strain, and C57BL/6J, an albuminuria-resistant strain, to perform a quantitative trait locus (QTL) cross to identify the genetic basis for chronic kidney failure. Albumin-creatinine ratio (ACR) was measured in 130 F2 male offspring. One significant QTL was identified on chromosome (Chr) X and four suggestive QTL were found on Chrs 6, 7, 12, and 13. Narrowing of the QTL region was focused on the X-linked QTL and performed by incorporating genotype and expression analyses for genes located in the region. From the 485 genes identified in the X-linked QTL region, a few candidate genes were identified using a combination of bioinformatic evidence based on genomic comparison of the parental strains and known function in urine homeostasis. Finally, this study demonstrates the significance of the X chromosome in the genetic determination of albuminuria.

Mapping an X-linked locus that influences heat-induced febrile seizures in mice

PURPOSE

Febrile seizures (FS) are the most common seizure type in children between the age of 6 months and 5 years. Although FS are largely benign, recurrent FS are a major risk factor for developing temporal lobe epilepsy (TLE) later in life. The mechanisms underlying FS are largely unknown; however, family and twin studies indicate that FS susceptibility is under complex genetic control. We have recently developed a phenotypic screen to study the genetics of FS susceptibility in mice. Using this screen in a phenotype-driven genetic strategy we analyzed the C57BL/6J-Chr #(A)/NaJ chromosome substitution strain (CSS) panel. In each CSS line one chromosome of the A/J strain is substituted in a genetically homogeneous C57BL/6J background. The analysis of the CSS panel revealed that A/J chromosomes 1, 2, 6, 10, 13, and X carry at least one quantitative trait locus (QTL) for heat-induced FS susceptibility. The fact that many X-linked genes are highly expressed in the brain and have been implicated in human developmental disorders often presenting with seizures (like fragile X mental retardation) prompted us to map the chromosome X QTL.

METHODS

C57BL/6J mice were mated with C57BL/6J-Chr X(A) /NaJ (CSSX) to generate F(2)-generations-CXBL6 and BL6CX-originating from CSSX or C57BL/6J mothers, respectively. Heat-induced FS were elicited on postnatal day 14 by exposure to a controlled warm airstream of 50°C. The latency to heat-induced FS is our phenotype. This phenotype has previously been validated by video-electroencephalography (EEG) monitoring. After phenotyping and genotyping the F(2)-population, QTL analysis was performed using R/QTL software.

KEY FINDINGS

QTL analysis revealed a significant peak with an LOD-score of 3.25. The 1-LOD confidence interval (149,886,866-158,836,462 bp) comprises 52 protein coding genes, of which 34 are known to be brain expressed. Two of these brain-expressed genes have previously been linked to X-linked epilepsies, namely Cdkl5 and Pdha1.

SIGNIFICANCE

Our results show that the mouse genetics of X-linked FS susceptibility is complex, and that our heat-induced FS-driven genetic approach is a powerful tool for use in unraveling the complexities of this trait in mice. Fine-mapping and functional studies will be required to further identify the X-linked FS susceptibility genes.

Using bioinformatics and systems genetics to dissect HDL-cholesterol genetics in an MRL/MpJ x SM/J intercross

A higher incidence of coronary artery disease is associated with a lower level of HDL-cholesterol. We searched for genetic loci influencing HDL-cholesterol in F2 mice from a cross between MRL/MpJ and SM/J mice. Quantitative trait loci (QTL) mapping revealed one significant HDL QTL (Apoa2 locus), four suggestive QTL on chromosomes 10, 11, 13, and 18 and four additional QTL on chromosomes 1 proximal, 3, 4, and 7 after adjusting HDL for the strong Apoa2 locus. A novel nonsynonymous polymorphism supports Lipg as the QTL gene for the chromosome 18 QTL, and a difference in Abca1 expression in liver tissue supports it as the QTL gene for the chromosome 4 QTL. Using weighted gene co-expression network analysis, we identified a module that after adjustment for Apoa2, correlated with HDL, was genetically determined by a QTL on chromosome 11, and overlapped with the HDL QTL. A combination of bioinformatics tools and systems genetics helped identify several candidate genes for both the chromosome 11 HDL and module QTL based on differential expression between the parental strains, cis regulation of expression, and causality modeling. We conclude that integrating systems genetics to a more-traditional genetics approach improves the power of complex trait gene identification.

Multiple Quantitative Trait Loci Influence the Shape of a Male-Specific Genital Structure in Drosophila melanogaster

The observation that male genitalia diverge more rapidly than other morphological traits during evolution is taxonomically widespread and likely due to some form of sexual selection. One way to elucidate the evolutionary forces acting on these traits is to detail the genetic architecture of variation both within and between species, a program of research that is considerably more tractable in a model system. Drosophila melanogaster and its sibling species, D. simulans, D. mauritiana, and D. sechellia, are morphologically distinguishable only by the shape of the posterior lobe, a male-specific elaboration of the genital arch. We extend earlier studies identifying quantitative trait loci (QTL) responsible for lobe divergence across species and report the first genetic dissection of lobe shape variation within a species. Using an advanced intercross mapping design, we identify three autosomal QTL contributing to the difference in lobe shape between a pair of D. melanogaster inbred lines. The QTL each contribute 4.6–10.7% to shape variation, and two show a significant epistatic interaction. Interestingly, these intraspecific QTL map to the same locations as interspecific lobe QTL, implying some shared genetic control of the trait within and between species. As a first step toward a mechanistic understanding of natural lobe shape variation, we find an association between our QTL data and a set of genes that show sex-biased expression in the developing genital imaginal disc (the precursor of the adult genitalia). These genes are good candidates to harbor naturally segregating polymorphisms contributing to posterior lobe shape.

Identification of candidate gene regions in the rat by co-localization of QTLs for bone density, size, structure and strength

Susceptibility to osteoporotic fracture is influenced by genetic factors that can be dissected by whole-genome linkage analysis in experimental animal crosses. The aim of this study was to characterize quantitative trait loci (QTLs) for biomechanical and two-dimensional dual-energy X-ray absorptiometry (DXA) phenotypes in reciprocal F2 crosses between diabetic GK and normo-glycemic F344 rat strains and to identify possible co-localization with previously reported QTLs for bone size and structure. The biomechanical measurements of rat tibia included ultimate force, stiffness and work to failure while DXA was used to characterize tibial area, bone mineral content (BMC) and areal bone mineral density (aBMD). F2 progeny (108 males, 98 females) were genotyped with 192 genome-wide markers followed by sex- and reciprocal cross-separated whole-genome QTL analyses. Significant QTLs were identified on chromosome 8 (tibial area; logarithm of odds (LOD) = 4.7 and BMC; LOD = 4.1) in males and on chromosome 1 (stiffness; LOD = 5.5) in females. No QTLs showed significant sex-specific interactions. In contrast, significant cross-specific interactions were identified on chromosome 2 (aBMD; LOD = 4.7) and chromosome 6 (BMC; LOD = 4.8) for males carrying F344mtDNA, and on chromosome 15 (ultimate force; LOD = 3.9) for males carrying GKmtDNA, confirming the effect of reciprocal cross on osteoporosis-related phenotypes. By combining identified QTLs for biomechanical-, size- and qualitative phenotypes (pQCT and 3D CT) from the same population, overlapping regions were detected on chromosomes 1, 3, 4, 6, 8 and 10. These are strong candidate regions in the search for genetic risk factors for osteoporosis.

Genetic dissection of a key reproductive barrier between nascent species of house mice

Reproductive isolation between species is often caused by deleterious interactions among loci in hybrids. Finding the genes involved in these incompatibilities provides insight into the mechanisms of speciation. With recently diverged subspecies, house mice provide a powerful system for understanding the genetics of reproductive isolation early in the speciation process. Although previous studies have yielded important clues about the genetics of hybrid male sterility in house mice, they have been restricted to F1 sterility or incompatibilities involving the X chromosome. To provide a more complete characterization of this key reproductive barrier, we conducted an F2 intercross between wild-derived inbred strains from two subspecies of house mice, Mus musculus musculus and Mus musculus domesticus. We identified a suite of autosomal and X-linked QTL that underlie measures of hybrid male sterility, including testis weight, sperm density, and sperm morphology. In many cases, the autosomal loci were unique to a specific sterility trait and exhibited an effect only when homozygous, underscoring the importance of examining reproductive barriers beyond the F1 generation. We also found novel two-locus incompatibilities between the M. m. musculus X chromosome and M. m. domesticus autosomal alleles. Our results reveal a complex genetic architecture for hybrid male sterility and suggest a prominent role for reproductive barriers in advanced generations in maintaining subspecies integrity in house mice.

Integration of QTL and bioinformatic tools to identify candidate genes for triglycerides in mice

To identify genetic loci influencing lipid levels, we performed quantitative trait loci (QTL) analysis between inbred mouse strains MRL/MpJ and SM/J, measuring triglyceride levels at 8 weeks of age in F2 mice fed a chow diet. We identified one significant QTL on chromosome (Chr) 15 and three suggestive QTL on Chrs 2, 7, and 17. We also carried out microarray analysis on the livers of parental strains of 282 F2 mice and used these data to find cis-regulated expression QTL. We then narrowed the list of candidate genes under significant QTL using a "toolbox" of bioinformatic resources, including haplotype analysis; parental strain comparison for gene expression differences and nonsynonymous coding single nucleotide polymorphisms (SNP); cis-regulated eQTL in livers of F2 mice; correlation between gene expression and phenotype; and conditioning of expression on the phenotype. We suggest Slc25a7 as a candidate gene for the Chr 7 QTL and, based on expression differences, five genes (Polr3 h, Cyp2d22, Cyp2d26, Tspo, and Ttll12) as candidate genes for Chr 15 QTL. This study shows how bioinformatics can be used effectively to reduce candidate gene lists for QTL related to complex traits.

A new mapping method for quantitative trait loci of silkworm

BACKGROUND

Silkworm is the basis of sericultural industry and the model organism in insect genetics study. Mapping quantitative trait loci (QTLs) underlying economically important traits of silkworm is of high significance for promoting the silkworm molecular breeding and advancing our knowledge on genetic architecture of the Lepidoptera. Yet, the currently used mapping methods are not well suitable for silkworm, because of ignoring the recombination difference in meiosis between two sexes.

RESULTS

A mixed linear model including QTL main effects, epistatic effects, and QTL × sex interaction effects was proposed for mapping QTLs in an F2 population of silkworm. The number and positions of QTLs were determined by F-test and model selection. The Markov chain Monte Carlo (MCMC) algorithm was employed to estimate and test genetic effects of QTLs and QTL × sex interaction effects. The effectiveness of the model and statistical method was validated by a series of simulations. The results indicate that when markers are distributed sparsely on chromosomes, our method will substantially improve estimation accuracy as compared to the normal chiasmate F2 model. We also found that a sample size of hundreds was sufficiently large to unbiasedly estimate all the four types of epistases (i.e., additive-additive, additive-dominance, dominance-additive, and dominance-dominance) when the paired QTLs reside on different chromosomes in silkworm.

CONCLUSION

The proposed method could accurately estimate not only the additive, dominance and digenic epistatic effects but also their interaction effects with sex, correcting the potential bias and precision loss in the current QTL mapping practice of silkworm and thus representing an important addition to the arsenal of QTL mapping tools.

Genetic loci for bone architecture determined by three-dimensional CT in crosses with the diabetic GK rat

The F344 rat carries alleles contributing to bone fragility while the GK rat spontaneously develops type-2 diabetes. These characteristics make F344×GK crosses well suited for the identification of genes related to bone size and allow for future investigation on the association with type-2 diabetes. The aim of this study was to identify quantitative trait loci (QTLs) for bone size phenotypes measured by a new application of three-dimensional computed tomography (3DCT) and to investigate the effects of sex- and reciprocal cross. Tibia from male and female GK and F344 rats, representing the parental, F1 and F2 generations, were examined with 3DCT and analyzed for: total and cortical volumetric BMD, straight and curved length, peri- and endosteal area at mid-shaft. F2 progeny (108 male and 98 female) were genotyped with 192 genome-wide microsatellite markers (average distance 10 cM). Sex- and reciprocal cross-separated QTL analyses were performed for the identification of QTLs linked to 3DCT phenotypes and true interactions were confirmed by likelihood ratio analysis in all F2 animals. Several genome-wide significant QTLs were found in the sex- and reciprocal cross-separated progeny on chromosomes (chr) 1, 3, 4, 9, 10, 14, and 17. Overlapping QTLs for both males and females in the (GK×F344)F2 progeny were located on chr 1 (39-67 cM). This region confirms previously reported pQCT QTLs and overlaps loci for fasting glucose. Sex separated linkage analysis confirmed a male specific QTL on chr 9 (67-82 cM) for endosteal area at the fibula site. Analyses separating the F2 population both by sex and reciprocal cross identified cross specific QTLs on chr 14 (males) and chr 3 and 4 (females). Two loci, chr 4 and 6, are unique to 3DCT and separate from pQCT generated loci. The 3DCT method was highly reproducible and provided high precision measurements of bone size in the rat enabling identification of new sex- and cross-specific loci. The QTLs on chr 1 indicate potential genetic association between bone-related phenotypes and traits affecting type-2 diabetes. The results illustrate the complexity of the genetic architecture of bone size phenotypes and demonstrate the importance of complementary methods for bone analysis.

Multiple loci contribute to genome-wide recombination levels in male mice

Recent linkage-based studies in humans suggest the presence of loci that affect either genome-wide recombination rates, utilization of recombination hotspots, or both. We have been interested in utilizing cytological methodology to directly assess recombination in mammalian meiocytes and to identify recombination-associated loci. In the present report we summarize studies in which we combined a cytological assay of recombination in mouse pachytene spermatocytes with QTL analyses to identify loci that contribute to genome-wide levels of recombination in male meiosis. Specifically, we analyzed MLH1 foci, a marker of crossovers, in 194 F2 male mice derived from a subspecific cross between CAST/EiJ and C57BL/6J parental strains. We then used these data to uncover loci associated with individual variation in mean MLH1 values. We identified seven recombination-associated loci across the genome (on chromosomes 2, 3, 4, 14, 15, 17, and X), indicating that there are multiple recombination “setting” loci in mammalian male meiosis.

Uncovering genes and regulatory pathways related to urinary albumin excretion

Identifying the genes underlying quantitative trait loci (QTL) for disease is difficult, mainly because of the low resolution of the approach and the complex genetics involved. However, recent advances in bioinformatics and the availability of genetic resources now make it possible to narrow the genetic intervals, test candidate genes, and define pathways affected by these QTL. In this study, we mapped three significant QTL and one suggestive QTL for an increased albumin-to-creatinine ratio on chromosomes (Chrs) 1, 4, 15, and 17, respectively, in a cross between the inbred MRL/MpJ and SM/J strains of mice. By combining data from several sources and by utilizing gene expression data, we identified Tlr12 as a likely candidate for the Chr 4 QTL. Through the mapping of 33,881 transcripts measured by microarray on kidney RNA from each of the 173 male F2 animals, we identified several downstream pathways associated with these QTL, including the glycan degradation, leukocyte migration, and antigen-presenting pathways. We demonstrate that by combining data from multiple sources, we can identify not only genes that are likely to be causal candidates for QTL but also the pathways through which these genes act to alter phenotypes. This combined approach provides valuable insights into the causes and consequences of renal disease.

A modifier locus on chromosome 5 contributes to L1 cell adhesion molecule X-linked hydrocephalus in mice

Humans with L1 cell adhesion molecule (L1CAM) mutations exhibit X-linked hydrocephalus, as well as other severe neurological disorders. L1-6D mutant mice, which are homozygous for a deletion that removes the sixth immunoglobulin-like domain of L1cam, seldom display hydrocephalus on the 129/Sv background. However, the same L1-6D mutation produces severe hydrocephalus on the C57BL/6J background. To begin to understand how L1cam deficiencies result in hydrocephalus and to identify modifier loci that contribute to X-linked hydrocephalus by genetically interacting with L1cam, we conducted a genome-wide scan on F2 L1-6D mice, bred from L1-6D 129S2/SvPasCrlf and C57BL/6J mice. Linkage studies, utilizing chi-square tests and quantitative trait loci mapping techniques, were performed. Candidate modifier loci were further investigated in an extension study. Linkage was confirmed for a locus on chromosome 5, which we named L1cam hydrocephalus modifier 1 (L1hydro1), p = 4.04 X 10(-11).

eQTL analysis in mice and rats

Since the introduction of genetical genomics in 2001, many studies have been published on various organisms, including mouse and rat. Genetical genomics makes use of the latest microarray profiling technologies and combines vast amounts of genotype and gene expression information, a strategy that has proven very successful in inbred line crosses. The data are analyzed using standard tools for linkage analysis to map the genetic determinants of gene expression variation. Typically, studies have singled out hundreds of genomic loci regulating the expression of nearby and distant genes (called local and distant expression quantitative trait loci, respectively; eQTLs). In this chapter, we provide a step-by-step guide to performing genome-wide linkage analysis in an eQTL mapping experiment by using the R statistical software framework.

A method for mapping intralocus interactions influencing excessive alcohol drinking

Excessive alcohol (ethanol) consumption is the hallmark of alcohol use disorders. The F1 hybrid cross between the C57BL/6J (B6) and FVB/NJ (FVB) inbred mouse strains consumes more ethanol than either progenitor strain. The purpose of this study was to utilize ethanol-drinking data and genetic information to map genes that result in overdominant (or heterotic) ethanol drinking. About 600 B6 x FVB F2 mice, half of each sex, were tested for ethanol intake and preference in a 24-h, two-bottle water versus ethanol choice procedure, with ascending ethanol concentrations. They were then tested for ethanol intake in a Drinking in the Dark (DID) procedure, first when there was no water choice and then when ethanol was offered versus water. DNA samples were obtained and genome-wide QTL analyses were performed to search for single QTLs (both additive and dominance effects) and interactions between pairs of QTLs, or epistasis. On average, F2 mice consumed excessive amounts of ethanol in the 24-h choice procedure, consistent with high levels of consumption seen in the F1 cross. Consumption in the DID procedure was similar or higher than amounts reported previously for the B6 progenitor. QTLs resulting in heightened consumption in heterozygous compared to homozygous animals were found on Chrs 11, 15, and 16 for 24-h choice 30% ethanol consumption, and on Chr 11 for DID. No evidence was found for epistasis between any pair of significant or suggestive QTLs. This indicates that the hybrid overdominance is due to intralocus interactions at the level of individual QTL.

Modifiers of von Willebrand factor identified by natural variation in inbred strains of mice

Type 1 von Willebrand disease (VWD) is the most common inherited human bleeding disorder. However, diagnosis is complicated by incomplete penetrance and variable expressivity, as well as wide variation in von Willebrand factor (VWF) levels among the normal population. Previous work has exploited the highly variable plasma VWF levels among inbred strains of mice to identify 2 major regulators, Mvwf1 and Mvwf2 (modifier of VWF). Mvwf1 is a glycosyltransferase and Mvwf2 is a natural variant in Vwf that alters biosynthesis. We report the identification of an additional alteration at the Vwf locus (Mvwf5), as well as 2 loci unlinked to Vwf (Mvwf6-7) using a backcross approach with the inbred mouse strains WSB/EiJ and C57BL/6J. Through positional cloning, we show that Mvwf5 is a cis-regulatory variant that alters Vwf mRNA expression. A similar mechanism could potentially explain a significant percentage of human VWD cases, especially those with no detectable mutation in the VWF coding sequence. Mvwf6 displays conservation of synteny with potential VWF modifier loci identified in human pedigrees, suggesting that its ortholog may modify VWF in human populations.

Sex-specific quantitative trait loci linked to autoresuscitation failure in SWR/J mice

Autoresuscitation (AR) is a highly conserved response among mammals, which allows survival from transient extreme hypoxia. During hypoxia, bradycardia, and hypoxic gasping develop after a brief period of hyperactivity. Normally, AR occurs if oxygen is restored during the gasping period where an initial heart rate increase is rapidly followed resumption or eupneic breathing. Humans and other mammals can survive multiple immediately repeated AR. A defective AR capacity has been implicated in Sudden Infant Death Syndrome. We had reported earlier that inbred strains of mice such as BALB/cJ could survive a characteristic number of immediately repeated AR trials, but that SWR/J mice failed to AR from a single hypoxic episode. We now report that strains closely related to SWR/J, FVB/N and SJL/J exhibit partial resuscitation defects relative to BALB/cJ or other mouse strains, establishing a genetic basis for variation in AR failure. The AR trial phenotype of BALB/cJ x SWR/J intercross F(1) and F(2) mice was consistent with BALB/cJ dominance and a discrete number of loci. Genome-wide mapping conducted with 60 intercross F(2) animals linked two loci to the number of AR trials survived, including one sex-specific locus with male expression, consistent with the observed 50% male bias for Sudden Infant Death Syndrome in humans. A locus carried on SWR/J chromosome 10 seems to be particularly important in AR failure and was confirmed in a partial consomic line. These results establish a genetic basis for AR failure phenotype in mice, with relevance to Sudden Infant Death Syndrome.

A Cmv2 QTL on chromosome X affects MCMV resistance in New Zealand male mice

NK cell-mediated resistance to viruses is subject to genetic control in humans and mice. Here we used classical and quantitative genetic strategies to examine NK-mediated murine cytomegalovirus (MCMV) control in genealogically related New Zealand white (NZW) and black (NZB) mice. NZW mice display NK cell-dependent MCMV resistance while NZB NK cells fail to limit viral replication after infection. Unlike Ly49H(+) NK resistance in C57BL/6 mice, NZW NK-mediated MCMV control was Ly49H-independent. Instead, MCMV resistance in NZW (Cmv2) involves multiple genetic factors. To establish the genetic basis of Cmv2 resistance, we further characterized a major chromosome X-linked resistance locus (DXMit216) responsible for innate MCMV control in NZW x NZB crosses. We found that the DXMit216 locus affects early MCMV control in New Zealand F(2) crosses and demonstrate that the NZB-derived DXMit216 allele enhances viral resistance in F(2) males. The evolutionary conservation of the DXMit216 region in mice and humans suggests that a Cmv2-related mechanism may affect human antiviral responses.

yellow0, a marker for low body weight in Drosophila melanogaster

Marker-assisted selection (MAS) is an important modern breeding technique, but it has been found that the effect of the markers for quantitative trait loci (QTL) is inconsistent, leading in some cases to MAS failure and raising doubts about its effectiveness. Here the model organism Drosophila melanogaster was employed to study whether an effective marker could be found and applied to MAS. We crossed the stock carrying the y (0) marker (a recessive mutation allele of the yellow gene on the X chromosome) with three other stocks carrying corresponding wild-type markers in an F2 design, and found that the y (0) marker was in significant association with low body weight (P<0.001). This association was consistent across different backgrounds and the marker effects in female and male were approximately 0.95 sigma (P) (phenotypic standard deviation) and 0.68 sigma (P), respectively. We next introgressed a fragment via the y (0) marker into a wild stock background over 20 generations of marker-assisted introgression (MAI), and constructed the introgression stock y (0)(OR)20 in which body weight decreased by 13% and 7%, in female and male, respectively, compared to the wild stock (P<0.0001). This indicated that there must be a single QTL for low body weight that is tightly linked to the y (0) marker. We then shortened the introgressed fragment to less than 1.5 cM by a deeper MAI using the y (0) marker and the white marker. This narrower fragment also resulted in a similar decrease in body weight to that induced by y (0)(OR)20, indicating that the QTL for low body weight is located within this less-than-1.5 cM interval. Molecular characteristics of the y (0) marker by PCR amplification and Southern blotting revealed that yellow gene was deficient in the y (0) stock, leading to disappearance of melanin from the cuticle and probably influencing the developmental process. The above results confirmed the existence of effective QTL markers applicable to MAS breeding schemes, and their potential application in breeding new stocks.

Reciprocal backcross mice confirm major loci linked to hyperoxic acute lung injury survival time

Morbidity and mortality associated with acute lung injury (ALI) and acute respiratory distress syndrome remain substantial. Although many candidate genes have been tested, a clear understanding of the pathogenesis is lacking, as is our ability to predict individual outcome. Because ALI is a complex disease, single gene approaches cannot easily identify effectors that must be treated concurrently. We employed a strategy to help identify critical genes and gene combinations involved in ALI mortality. Using hyperoxia to induce ALI, a mouse model for genetic analyses of ALI survival time was identified: C57BL/6J (B) mice are sensitive (i.e., die early), whereas 129X1/SvJ (S) mice are significantly more resistant, but with low penetrance. Segregation analysis of reciprocal F(2) mice generated from B and S strains revealed significant sex, cross, and parent of origin effects. Quantitative trait locus (QTL) analysis identified five chromosomal regions significantly linked to hyperoxic ALI survival time (named Shali1-Shali5). Further analyses demonstrated that both parental strains contribute resistance alleles to their offspring and that the phenotype demonstrated parent of origin effects. To validate earlier findings, we generated and tested mice from all eight possible B-S-derived backcrosses. Results from segregation and QTL analyses of 935 backcrosses, alone and combined with the previous 840 B-S-derived F(2) population, further supported the highly significant QTLs on chromosomes 1 (Shali1) and 4 (Shali2) and confirmed that the sex, cross, and parent of origin all contribute to survival time with hyperoxic ALI.

A model selection approach for the identification of quantitative trait loci in experimental crosses, allowing epistasis

The identification of quantitative trait loci (QTL) and their interactions is a crucial step toward the discovery of genes responsible for variation in experimental crosses. The problem is best viewed as one of model selection, and the most important aspect of the problem is the comparison of models of different sizes. We present a penalized likelihood approach, with penalties on QTL and pairwise interactions chosen to control false positive rates. This extends the work of Broman and Speed to allow for pairwise interactions among QTL. A conservative version of our penalized LOD score provides strict control over the rate of extraneous QTL and interactions; a more liberal criterion is more lenient on interactions but seeks to maintain control over the rate of inclusion of false loci. The key advance is that one needs only to specify a target false positive rate rather than a prior on the number of QTL and interactions. We illustrate the use of our model selection criteria as exploratory tools; simulation studies demonstrate reasonable power to detect QTL. Our liberal criterion is comparable in power to two Bayesian approaches.

A dominant X-linked QTL regulating pubertal timing in mice found by whole genome scanning and modified interval-specific congenic strain analysis

BACKGROUND

Pubertal timing in mammals is triggered by reactivation of the hypothalamic-pituitary-gonadal (HPG) axis and modulated by both genetic and environmental factors. Strain-dependent differences in vaginal opening among inbred mouse strains suggest that genetic background contribute significantly to the puberty timing, although the exact mechanism remains unknown.

METHODOLOGY/PRINCIPAL FINDINGS

We performed a genome-wide scanning for linkage in reciprocal crosses between two strains, C3H/HeJ (C3H) and C57BL6/J (B6), which differed significantly in the pubertal timing. Vaginal opening (VO) was used to characterize pubertal timing in female mice, and the age at VO of all female mice (two parental strains, F1 and F2 progeny) was recorded. A genome-wide search was performed in 260 phenotypically extreme F2 mice out of 464 female progeny of the F1 intercrosses to identify quantitative trait loci (QTLs) controlling this trait. A QTL significantly associated was mapped to the DXMit166 marker (15.5 cM, LOD = 3.86, p<0.01) in the reciprocal cross population (C3HB6F2). This QTL contributed 2.1 days to the timing of VO, which accounted for 32.31% of the difference between the original strains. Further study showed that the QTL was B6-dominant and explained 10.5% of variation to this trait with a power of 99.4% at an alpha level of 0.05.The location of the significant ChrX QTL found by genome scanning was then fine-mapped to a region of approximately 2.5 cM between marker DXMit68 and rs29053133 by generating and phenotyping a panel of 10 modified interval-specific congenic strains (mISCSs).

CONCLUSIONS/SIGNIFICANCE

Such findings in our study lay a foundation for positional cloning of genes regulating the timing of puberty, and also reveal the fact that chromosome X (the sex chromosome) does carry gene(s) which take part in the regulative pathway of the pubertal timing in mice.