Rainbow trout resistance to bacterial cold water disease: two new quantitative trait loci identified after a natural disease outbreak on a French farm

In rainbow trout farming, Flavobacterium psychrophilum, the causative agent of bacterial cold water disease, is responsible for important economic losses. Resistance to F. psychrophilum is heritable, and several quantitative trait loci (QTL) with moderate effects have been detected, opening up promising perspectives for the genetic improvement of resistance. In most studies however, resistance to F. psychrophilum was assessed in experimental infectious challenges using injection as the infection route, which is not representative of natural infection. Indeed, injection bypasses external barriers, such as mucus and skin, that likely play a protective role against the infection. In this study, we aimed at describing the genetic architecture of the resistance to F. psychrophilum after a natural disease outbreak. In a 2000‐fish cohort, reared on a French farm, 720 fish were sampled and genotyped using the medium‐throughput Axiom? Trout Genotyping Array. Overall mortality at the end of the outbreak was 25%. Genome‐wide association studies were performed under two different models for time to death measured on 706 fish with validated genotypes for 30 060 SNPs. This study confirms the polygenic inheritance of resistance to F. psychrophilum with a few QTL with moderate effects and a large polygenic background, the heritability of the trait being estimated at 0.34. Two new chromosome‐wide significant QTL and three suggestive QTL were detected, each of them explaining between 1% and 4% of genetic variance.     

Localization of quantitative trait loci for diapause and other photoperiodically regulated life history traits important in adaptation to seasonally varying environments

Seasonally changing environments at high latitudes present great challenges for the reproduction and survival of insects, and photoperiodic cues play an important role in helping them to synchronize their life cycle with prevalent and forthcoming conditions. We have mapped quantitative trait loci (QTL) responsible for the photoperiodic regulation of four life history traits, female reproductive diapause, cold tolerance, egg‐to‐eclosion development time and juvenile body weight in Drosophila montana strains from different latitudes in Canada and Finland. The F2 progeny of the cross was reared under a single photoperiod (LD cycle 16:8), which the flies from the Canadian population interpret as early summer and the flies from the Finnish population as late summer. The analysis revealed a unique QTL for diapause induction on the X chromosome and several QTL for this and the other measured traits on the 4th chromosome. Flies’ cold tolerance, egg‐to‐eclosion development time and juvenile body weight had several QTL also on the 2nd, 3rd and 5th chromosome, some of the peaks overlapping with each other. These results suggest that while the downstream output of females’ photoperiodic diapause response is partly under a different genetic control from that of the other traits in the given day length, all traits also share some QTL, possibly involving genes with pleiotropic effects and/or multiple tightly linked genes. Nonoverlapping QTL detected for some of the traits also suggest that the traits are potentially capable of independent evolution, even though this may be restricted by epistatic interactions and/or correlations and trade‐offs between the traits.     

QTL affecting stress response to crowding in a rainbow trout broodstock population

Background

Genomic analyses have the potential to impact selective breeding programs by identifying markers that serve as proxies for traits which are expensive or difficult to measure. Also, identifying genes affecting traits of interest enhances our understanding of their underlying biochemical pathways. To this end we conducted genome scans of seven rainbow trout families from a single broodstock population to identify quantitative trait loci (QTL) having an effect on stress response to crowding as measured by plasma cortisol concentration. Our goal was to estimate the number of major genes having large effects on this trait in our broodstock population through the identification of QTL.

Results

A genome scan including 380 microsatellite markers representing 29 chromosomes resulted in the de novo construction of genetic maps which were in good agreement with the NCCCWA genetic map. Unique sets of QTL were detected for two traits which were defined after observing a low correlation between repeated measurements of plasma cortisol concentration in response to stress. A highly significant QTL was detected in three independent analyses on Omy16, many additional suggestive and significant QTL were also identified. With linkage-based methods of QTL analysis such as half-sib regression interval mapping and a variance component method, we determined that the significant and suggestive QTL explain about 40-43% and 13-27% of the phenotypic trait variation, respectively.

Conclusions

The cortisol response to crowding stress is a complex trait controlled in a sub-sample of our broodstock population by multiple QTL on at least 8 chromosomes. These QTL are largely different from others previously identified for a similar trait, documenting that population specific genetic variants independently affect cortisol response in ways that may result in different impacts on growth. Also, mapping QTL for multiple traits associated with stress response detected trait specific QTL which indicate the significance of the first plasma cortisol measurement in defining the trait. Fine mapping these QTL can lead towards the identification of genes affecting stress response and may influence approaches to selection for this economically important stress response trait.

     

Quantitative trait loci for growth and body size in the nine‐spined stickleback Pungitius pungitius L.

Body size is an ecologically important trait shown to be genetically variable both within and among different animal populations as revealed by quantitative genetic studies. However, few studies have looked into underlying genetic architecture of body size variability in the wild using genetic mapping methods. With the aid of quantitative trait loci (QTL) analyses based on 226 microsatellite markers, we mapped body size and growth rate traits in the nine‐spined stickleback (Pungitius pungitius) using an F₂‐intercross (n = 283 offspring) between size‐divergent populations. In total, 17 QTL locations were detected. The proportion of phenotypic variation explained by individual body size‐related QTL ranged from 3% to 12% and those related to growth parameters and increments from 3% to 10%. Several of the detected QTL affected either early or late growth. These results provide a solid starting point for more in depth investigations of structure and function of genomic regions involved in determination of body size in this popular model of ecological and evolutionary research.     

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.     

More than one way to evolve a weed: parallel evolution of US weedy rice through independent genetic mechanisms

Many different crop species were selected for a common suite of ‘domestication traits’, which facilitates their use for studies of parallel evolution. Within domesticated rice (Oryza sativa), there has also been independent evolution of weedy strains from different cultivated varieties. This makes it possible to examine the genetic basis of parallel weed evolution and the extent to which this process occurs through shared genetic mechanisms. We performed comparative QTL mapping of weediness traits using two recombinant inbred line populations derived from crosses between an indica crop variety and representatives of each of the two independently evolved weed strains found in US rice fields, strawhull (S) and blackhull awned (B). Genotyping‐by‐sequencing provided dense marker coverage for linkage map construction (average marker interval <0.25 cM), with 6016 and 13 730 SNPs mapped in F₅ lines of the S and B populations, respectively. For some weediness traits (awn length, hull pigmentation and pericarp pigmentation), QTL mapping and sequencing of underlying candidate genes confirmed that trait variation was largely attributable to individual loci. However, for more complex quantitative traits (including heading date, panicle length and seed shattering), we found multiple QTL, with little evidence of shared genetic bases between the S and B populations or across previous studies of weedy rice. Candidate gene sequencing revealed causal genetic bases for 8 of 27 total mapped QTL. Together these findings suggest that despite the genetic bottleneck that occurred during rice domestication, there is ample genetic variation in this crop to allow agricultural weed evolution through multiple genetic mechanisms.     

QTL for body weight,morphometric traits and stress response in European sea bass Dicentrarchus labrax

Natural mating and mass spawning in the European sea bass (Dicentrarchus labrax L., Moronidae, Teleostei) complicate genetic studies and the implementation of selective breeding schemes. We utilized a two‐step experimental design for detecting QTL in mass‐spawning species: 2122 offspring from natural mating between 57 parents (22 males, 34 females and one missing) phenotyped for body weight, eight morphometric traits and cortisol levels, had been previously assigned to parents based on genotypes of 31 DNA microsatellite markers. Five large full‐sib families (five sires and two dams) were selected from the offspring (570 animals), which were genotyped with 67 additional markers. A new genetic map was compiled, specific to our population, but based on the previously published map. QTL mapping was performed with two methods: half‐sib regression analysis (paternal and maternal) and variance component analysis accounting for all family relationships. Two significant QTL were found for body weight on linkage group 4 and 6, six significant QTL for morphometric traits on linkage groups 1B, 4, 6, 7, 15 and 23 and three suggestive QTL for stress response on linkage groups 3, 14 and 23. The QTL explained between 8% and 38% of phenotypic variance. The results are the first step towards identifying genes involved in economically important traits like body weight and stress response in European sea bass.     

The pulse of the tree is under genetic control: eucalyptus as a case study

The pulse of the tree (diurnal cycle of stem radius fluctuations) has been widely studied as a way of analyzing tree responses to the environment, including the phenotypic plasticity of tree–water relationships in particular. However, the genetic basis of this daily phenotype and its interplay with the environment remain largely unexplored. We characterized the genetic and environmental determinants of this response, by monitoring daily stem radius fluctuation (dSRF) on 210 trees from a Eucalyptus urophylla × E. grandis full‐sib family over 2 years. The dSRF signal was broken down into hydraulic capacitance, assessed as the daily amplitude of shrinkage (DA), and net growth, estimated as the change in maximum radius between two consecutive days (ΔR). The environmental determinants of these two traits were clearly different: DA was positively correlated with atmospheric variables relating to water demand, while ΔR was associated with soil water content. The heritability for these two traits ranged from low to moderate over time, revealing a time‐dependent or environment‐dependent complex genetic determinism. We identified 686 and 384 daily quantitative trait loci (QTL) representing 32 and 31 QTL regions for DA and ΔR, respectively. The identification of gene networks underlying the 27 major genomics regions for both traits generated additional hypotheses concerning the biological mechanisms involved in response to water demand and supply. This study highlights that environmentally induced changes in daily stem radius fluctuation are genetically controlled in trees and suggests that these daily responses integrated over time shape the genetic architecture of mature traits.     

Selection on QTL and complex traits in complex environments

Understanding genetic variation for complex traits in heterogeneous environments is a fundamental problem in biology. In this issue of Molecular Ecology, Fournier‐Level et al. ( 2013 ) analyse quantitative trait loci (QTL) influencing ecologically important phenotypes in mapping populations of Arabidopsis thaliana grown in four habitats across its native European range. They used causal modelling to quantify the selective consequences of life history and morphological traits and QTL on components of fitness. They found phenology QTL colocalizing with known flowering time genes as well as novel loci. Most QTL influenced fitness via life history and size traits, rather than QTL having direct effects on fitness. Comparison of phenotypes among environments found no evidence for genetic trade‐offs for phenology or growth traits, but genetic trade‐offs for fitness resulted because flowering time had opposite fitness effects in different environments. These changes in QTL effects and selective consequences may maintain genetic variation among populations.     

Genome‐wide association study for body weight in cattle populations from Siberia

Body weight is a complex trait in cattle associated with commonly used commercial breeding measurements related to growth. Although many quantitative trait loci (QTL) for body weight have been identified in cattle so far, searching for genetic determinants in different breeds or environments is promising. Therefore, we carried out a genome‐wide association study (GWAS) in two cattle populations from the Russian Federation (Siberian region) using the GGP HD150K array containing 139 376 single nucleotide polymorphism (SNP) markers. Association tests for 107 550 SNPs left after filtering revealed five statistically significant SNPs on BTA5, considering a false discovery rate of less than 0.05. The chromosomal region containing these five SNPs contains the CCND2 gene, which was previously associated with average daily weight gain and body mass index in US beef cattle populations and in humans respectively. Our study is the first GWAS for body weight in beef cattle populations from the Russian Federation. The results provided here suggest that, despite the existence of breed‐ and species‐specific QTL, the genetic architecture of body weight could be evolutionarily conserved in mammals.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

Two-trait selection response with marker-based assortative mating

 Marker-based assortative mating (MAM) – the mating of individuals that have similar genotypes at random marker loci – can increase selection response for a single trait by 3–8% over random mating (RM). Genetic gain is usually desired for multiple traits rather than for a single trait. My objectives in this study were to (1) compare MAM, phenotypic assortative mating (PAM), and RM of selected individuals for improving two traits and (2) determine when MAM will be most useful for improving two traits. I simulated 20 generations of selecting 32 out of 200 individuals in an F₂ population. The individuals were selected based on an index (SI) of two traits and were intermated by MAM, PAM, or RM. I studied eight genetic models that differed in three contrasts: (1) weight, number of quantitative trait loci (QTL), and heritability (h ²) for each trait; (2) linkage of QTL for each trait; and (3) trait means of the inbred parents of the F₂. For SI and the two component traits, MAM increased short-term selection response by 5–8% in six out of the eight genetic models. The MAM procedure was least effective in two genetic models, wherein the QTL for one trait were unlinked to the QTL for the other trait and the parents of the F₂ had divergent means for each trait. The loss of QTL heterozygosity was much greater with MAM than with PAM or RM. Consequently, the advantage of MAM over RM dissipated after 5–7 generations. Differences were small between selection responses with PAM and RM. The MAM procedure can enhance short-term selection response for two traits when selection is not stringent, h ² is low, and the means of the parents of the F₂ are equal for each trait. Received: 10 June 1998 / Accepted: 5 August 1998     

Mapping of post-flowering drought resistance traits in grain sorghum: association between QTLs influencing premature senescence and maturity

The identification of genetic factors underlying the complex responses of plants to drought stress provides a solid basis for improving drought resistance. The stay-green character in sorghum (Sorghum bicolor L. Moench) is a post-flowering drought resistance trait, which makes plants resistant to premature senescence under drought stress during the grainfilling stage. The objective of this study was to identify quantitative trait loci (QTLs) that control premature senescence and maturity traits, and to investigate their association under post-flowering drought stress in grain sorghum. A genetic linkage map was developed using a set of recombinant inbred lines (RILs) obtained from the cross B35 × Tx430, which were scored for 142 restriction fragment length polymorphism (RFLP) markers. The RILs and their parental lines were evaluated for post-flowering drought resistance and maturity in four environments. Simple interval mapping identified seven stay-green QTLs and two maturity QTLs. Three major stay-green QTLs (SGA, SGD and SGG) contributed to 42% of the phenotypic variability (LOD 9.0) and four minor QTLs (SGB, SGI.1, SGI.2, and SGJ) significantly contributed to an additional 25% of the phenotypic variability in stay-green ratings. One maturity QTL (DFB) alone contributed to 40% of the phenotypic variability (LOD 10.0), while the second QTL (DFG) significantly contributed to an additional 17% of the phenotypic variability (LOD 4.9). Composite interval mapping confirmed the above results with an additional analysis of the QTL × Environment interaction. With heritability estimates of 0.72 for stay-green and 0.90 for maturity, the identified QTLs explained about 90% and 63% of genetic variability for stay-green and maturity traits, respectively. Although stay-green ratings were significantly correlated (r=0.22, P ≤ 0.05) with maturity, six of the seven stay-green QTLs were independent of the QTLs influencing maturity. Similarly, one maturity QTL (DFB) was independent of the stay-green QTLs. One stay-green QTL (SGG), however, mapped in the vicinity of a maturity QTL (DFG), and all markers in the vicinity of the independent maturity QTL (DFB) were significantly (P ≤ 0.1) correlated with stay-green ratings, confounding the phenotyping of stay-green. The molecular genetic analysis of the QTLs influencing stay-green and maturity, together with the association between these two inversely related traits, provides a basis for further study of the underlying physiological mechanisms and demonstrates the possibility of improving drought resistance in plants by pyramiding the favorable QTLs. Received: 10 October 1998 / Accepted: 12 July 1999     

Mapping QTL for production traits in segregating Piétrain pig populations using genome‐wide association study results of F2 crosses

In this study, genome‐wide association study (GWAS) results of porcine F2 crosses were used to map QTL in outcross Piétrain populations. For this purpose, two F2 crosses (Piétrain × Meishan, n = 304; Piétrain × Wild Boar, n = 291) were genotyped with the PorcineSNP60v2 BeadChip and phenotyped for the dressing yield, carcass length, daily gain and drip loss traits. GWASs were conducted in the pooled F2 cross applying single marker mixed linear models. For the investigated traits, between two and five (in total 15) QTL core regions, spanning 250 segregating SNPs around a significant trait‐associated peak SNP, were identified. The SNPs within the QTL core regions were subsequently tested for trait association in two outcross Piétrain populations consisting of 771 progeny‐tested boars and 210 sows with their own performance records. In the sow (boar) dataset, five (eight) of the 15 mapped QTL were validated. Hence, many QTL mapped in the F2 crosses (with Piétrain as a common founder breed) are still segregating in the current Piétrain breed. This confirms the usefulness of existing F2 crosses for mapping QTL that are still segregating in the recent founder breed generation. The approach utilizes the high power of an F2 cross to map QTL in a breeding population for which it is not guaranteed that they would be found using a GWAS in this population.     

Stress responsiveness affects dominant-subordinate relationships in rainbow trout.

The magnitude by which plasma cortisol levels increase following exposure to a stressor is a heritable trait in rainbow trout. The relative growth in coculture of F1 lines selected for high responsiveness (HR) and low responsiveness (LR) to a confinement stressor suggested that behavioral characteristics related to food acquisition, aggression, or competitive ability might differ between the two lines. This hypothesis was tested using the F2 generation of the selected lines. The F2 lines clearly exhibited the characteristics of the F1 parents, displaying significantly divergent plasma cortisol responses to a 1-h confinement stressor and a high heritability for the trait. Behavioral differences between the lines were assessed by observing the outcome of staged fights for dominance in size-matched pairs of HR and LR fish. The identification of dominant and subordinate fish within each pair on the basis of their behavior was supported by the levels of blood cortisol in the fish attributed to each group (dominant < subordinate). Fish from the LR line were identified as dominant in significantly more trials than were HR individuals. The results suggest that behavioral attributes that affect the outcome of rank-order fights are closely linked to the magnitude of the plasma cortisol response to stress in rainbow trout. Whether the link is causal or circumstantial is not yet evident.     

Quantitative trait loci associated with constitutive traits control water use in pearl millet [Pennisetum glaucum (L.) R. Br.]

There is substantial genetic variation for drought adaption in pearl millet in terms of traits controlling plant water use. It is important to understand genomic regions responsible for these traits. Here, F₇ recombinant inbred lines were used to identify quantitative trait loci (QTL) and allelic interactions for traits affecting plant water use, and their relevance is discussed for crop productivity in water‐limited environments. Four QTL contributed to increased transpiration rate under high vapour pressure deficit (VPD) conditions, all with alleles from drought‐sensitive parent ICMB 841. Of these four QTL, a major QTL (35.7%) was mapped on linkage group (LG) 6. The alleles for 863B at this QTL decreased transpiration rate and this QTL co‐mapped to a previously detected LG 6 QTL, with alleles from 863B for grain weight and panicle harvest index across severe terminal drought stress environments. This provided additional support for a link between water saving from a lower transpiration rate under high VPD and drought tolerance. 863B alleles in this same genomic region also increased shoot weight, leaf area and total transpiration under well‐watered conditions. One unexpected outcome was reduced transpiration under high VPD (15%) from the interaction of two alleles for high VPD transpiration (LG 6 (B), 40.7) and specific leaf mass and biomass (LG 7 (A), 35.3), (A, allele from ICMB 841, B, allele from 863B, marker position). The LG 6 QTL appears to combine alleles for growth potential, beneficial for non‐stress conditions, and for saving water under high evaporative demand, beneficial under stressful conditions. Mapping QTL for water‐use traits, and assessing their interactions offers considerable potential for improving pearl millet adaptation to specific stress conditions through physiology‐informed marker‐assisted selection.     

Genomic architecture of habitat‐related divergence and signature of directional selection in the body shapes of Gnathopogon fishes

Evolution of ecomorphologically relevant traits such as body shapes is important to colonize and persist in a novel environment. Habitat‐related adaptive divergence of these traits is therefore common among animals. We studied the genomic architecture of habitat‐related divergence in the body shape of Gnathopogon fishes, a novel example of lake–stream ecomorphological divergence, and tested for the action of directional selection on body shape differentiation. Compared to stream‐dwelling Gnathopogon elongatus, the sister species Gnathopogon caerulescens, exclusively inhabiting a large ancient lake, had an elongated body, increased proportion of the caudal region and small head, which would be advantageous in the limnetic environment. Using an F₂ interspecific cross between the two Gnathopogon species (195 individuals), quantitative trait locus (QTL) analysis with geometric morphometric quantification of body shape and restriction‐site associated DNA sequencing‐derived markers (1622 loci) identified 26 significant QTLs associated with the interspecific differences of body shape‐related traits. These QTLs had small to moderate effects, supporting polygenic inheritance of the body shape‐related traits. Each QTL was mostly located on different genomic regions, while colocalized QTLs were detected for some ecomorphologically relevant traits that are proxy of body and caudal peduncle depths, suggesting different degree of modularity among traits. The directions of the body shape QTLs were mostly consistent with the interspecific difference, and QTL sign test suggested a genetic signature of directional selection in the body shape divergence. Thus, we successfully elucidated the genomic architecture underlying the adaptive changes of the quantitative and complex morphological trait in a novel system.     

Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L.

Gilthead sea bream (Sparus aurata L.) is an important marine fish in Mediterranean aquaculture. Sex determination by age and/or body weight is a critical life‐history trait, the genetic basis for which is largely unknown in this sequential hermaphrodite species. Herein, we performed a partial genome scan to map quantitative trait loci (QTL) affecting body weight and sex using 74 informative microsatellite markers from 10 paternal half‐sib families to construct nine linkage groups (LG). In total, four growth‐related QTL (two chromosome‐wide and two genome‐wide) and six QTL related to sex determination (three pairs in three different LGs) were detected (two chromosome‐wide and one genome‐wide). The proportion of phenotypic variation explained by the body‐weight QTL ranged from 9.3% to 17.2%, showing their potential for use in marker‐assisted selection. The results obtained offer solid ground to investigate the structure and function of the genomic regions involved in the mechanisms of sex reversal.     

QTL for the species-specific male and female genital morphologies in Ohomopterus ground beetles

Animals with internal fertilization often exhibit marked diversification in genital morphology among closely related species. However, our knowledge of the genetic architecture underlying genital evolution is still limited. We constructed genetic linkage maps and analysed quantitative trait loci (QTL) for F₂ hybrids of two closely related species of the carabid beetles Carabus (Ohomopterus) iwawakianus and C. (O.) maiyasanus, which show matching male and female genital shapes within species, but marked differences in genital morphologies between species. The linkage maps comprised both amplified fragment length polymorphism and microsatellite markers. Composite interval mapping to detect QTL for three traits of male copulatory piece (length, width, weight) and two traits for female vaginal appendix (length, width) resulted in the detection of one to five significant QTL for each trait. The QTL explained large proportions of phenotypic variance. Thus, the interspecific difference in the genital morphologies appeared to be determined by relatively small numbers of genes with large genetic effects. QTL of different traits for the same or different sexes co‐occurred on five of eight linkage groups with significant QTL; in particular, three QTL for different male and female genital traits occurred almost at the same position. Each of the male genital traits showed uniform signs of additive genetic effects, suggesting that directional selection has led to species‐specific morphologies. However, the signs of additive genetic effects in each female genital trait were not uniform, suggesting that coevolution between sexes is not necessarily concerted. This result requires further assessment because the sample size of F₂ females was small.