Single Linkage Group Research Articles

Population genomic analyses reveal hybridization and marked differences in genetic structure of Scurria limpet sister species with parapatric distributions across the South Eastern Pacific.

The study of sister species that occur in parapatry around biogeographic transition zones can help understand the evolutionary processes that underlie the changes in species composition across biogeographic transition zones. The South Eastern Pacific (SEP) coast is a highly productive coastal system that exhibits a broad biogeographic transition zone around 30–35°S. Here, we present a comparative genome‐wide analysis of the sister species Scurria viridula and Scurria zebrina, that occur in parapatry and whose poleward and equatorward range edges intersect in the 30–35°S SEP biogeographic transition zone. We sampled 118 specimens sourced from nine sites from Tocopilla (22°S) to Chiloé (41°S) including one site where both species overlap and analyzed over 8000 biallelic single nucleotide polymorphisms. We found evidence of hybridization between these species in the contact zone and found significant but contrasting population structures for both species. Our results indicate that the genetic structure in S. viridula, which is currently expanding its range poleward, follows a simple isolation by distance model with no traces of natural selection (no evidence of outlier loci). In contrast, S. zebrina, which has its equatorward range edge at the transition zone, displayed a pronounced genetic break approximately at 32–34°S, along a region of marked environmental heterogeneity in association with a semi‐permanent coastal upwelling regime. For S. zebrina we also found 43 outlier loci associated with this genetic break, with a significant proportion of them clustering in a single linkage group. This marked difference in the presence of outlier loci between species suggests that they could be responding differently to local environmental challenges found at their overlapping geographic range edges, thus providing important new insights about genomic changes around biogeographic transition zones in sister species and the forces that shape genetic diversity in intertidal marine species.

QTL mapping and marker identification for sex determination in the ridgetail white prawn, Exopalaemon carinicauda

Sex determination is an important and intriguing research topic in the field of evolutionary and developmental biology. Quantitative trait locus (QTL) mapping for sex is helpful in clarifying the sex determination system of species. In this study, a second high-resolution genetic linkage map was constructed for the ridgetail white prawn, Exopalaemon carinicauda, which included 9280 markers, covering 99.98% of the complete genome. Based on the linkage map, a highly significant sex-related QTL was first mapped to a single linkage group (LG3, LOD > 55.6). Fifty-two markers in the QTL region were significantly associated with sex (p ≤ 10−40), of which heterogametic genotypes in females supported the ZW sex determination mechanism. Six markers were verified to be significantly associated with sex in the wild population. Some sex-related genes were identified, including phospholipase D, protein kinase shaggy, and longitudinals lacking protein. These results inform our understanding of the mechanisms of sex determination in E. carinicauda.

Pod and Seed Trait QTL Identification To Assist Breeding for Peanut Market Preferences.

Although seed and pod traits are important for peanut breeding, little is known about the inheritance of these traits. A recombinant inbred line (RIL) population of 156 lines from a cross of Tifrunner x NC 3033 was genotyped with the Axiom_Arachis1 SNP array and SSRs to generate a genetic map composed of 1524 markers in 29 linkage groups (LG). The genetic positions of markers were compared with their physical positions on the peanut genome to confirm the validity of the linkage map and explore the distribution of recombination and potential chromosomal rearrangements. This linkage map was then used to identify Quantitative Trait Loci (QTL) for seed and pod traits that were phenotyped over three consecutive years for the purpose of developing trait-associated markers for breeding. Forty-nine QTL were identified in 14 LG for seed size index, kernel percentage, seed weight, pod weight, single-kernel, double-kernel, pod area and pod density. Twenty QTL demonstrated phenotypic variance explained (PVE) greater than 10% and eight more than 20%. Of note, seven of the eight major QTL for pod area, pod weight and seed weight (PVE >20% variance) were attributed to NC 3033 and located in a single linkage group, LG B06_1. In contrast, the most consistent QTL for kernel percentage were located on A07/B07 and derived from Tifrunner.

Genome Synteny Has Been Conserved Among the Octoploid Progenitors of Cultivated Strawberry Over Millions of Years of Evolution.

Allo-octoploid cultivated strawberry (Fragaria × ananassa) originated through a combination of polyploid and homoploid hybridization, domestication of an interspecific hybrid lineage, and continued admixture of wild species over the last 300 years. While genes appear to flow freely between the octoploid progenitors, the genome structures and diversity of the octoploid species remain poorly understood. The complexity and absence of an octoploid genome frustrated early efforts to study chromosome evolution, resolve subgenomic structure, and develop a single coherent linkage group nomenclature. Here, we show that octoploid Fragaria species harbor millions of subgenome-specific DNA variants. Their diversity was sufficient to distinguish duplicated (homoeologous and paralogous) DNA sequences and develop 50K and 850K SNP genotyping arrays populated with co-dominant, disomic SNP markers distributed throughout the octoploid genome. Whole-genome shotgun genotyping of an interspecific segregating population yielded 1.9M genetically mapped subgenome variants in 5,521 haploblocks spanning 3,394 cM in F. chiloensis subsp. lucida, and 1.6M genetically mapped subgenome variants in 3,179 haploblocks spanning 2,017 cM in F. × ananassa. These studies provide a dense genomic framework of subgenome-specific DNA markers for seamlessly cross-referencing genetic and physical mapping information and unifying existing chromosome nomenclatures. Using comparative genomics, we show that geographically diverse wild octoploids are effectively diploidized, nearly completely collinear, and retain strong macro-synteny with diploid progenitor species. The preservation of genome structure among allo-octoploid taxa is a critical factor in the unique history of garden strawberry, where unimpeded gene flow supported its origin and domestication through repeated cycles of interspecific hybridization.

Genomic analysis of powdery mildew resistance in a hop (Humulus lupulus L.) bi-parental population segregating for \u201cR6-locus\u201d

Genetic response in hop to fungal pathogen infection has been evaluated at the chromosomal level through QTL analyses but very little information exists on the expression of genes during infection periods. Raw GBS reads and phenotypic data from a previously published QTL analysis along with a newly assembled PacBio-derived hop reference genome were used to re-evaluate resistance to races v4/v6 of powdery mildew (PM; Podosphaera humuli). QTL analyses revealed two tightly linked regions of association on a single linkage group. The three SNP markers most tightly linked to PM resistance (found on contig 000559F) were observed downstream from a putative R-gene locus for powdery mildew resistance. This 230 kb region contained a series of seven putative R-genes surrounded by seven putative peroxidase-3 genes downstream and seven putative glucan endo-1,3-beta-glucosidase upstream and an expressed F-box domain protein. RNAseq data showed all putative R-genes along with all putative glucan endo-1,3-beta-glucosidase genes were expressed under diseased conditions, while none of the peroxidase genes were expressed. The second region contained three SNPs found on contig 002916F next to two putative R-genes. RNAseq data showed complex expression of exons contained in putative isoforms of R-genes. This preliminary information will prove valuable information for development of precise markers located either within or next to genes responsible for race v4/v6 PM resistance.

QTL Mapping and Marker Identification for Sex-Determining: Indicating XY Sex Determination System in the Swimming Crab (Portunus trituberculatus).

Sex determination is an important area of research, which has always had an intriguing aspect in evolutionary and developmental biology. Quantitative trait locus (QTL) mapping for sex will be helpful in clarifying the sex determination system. In this study, the sex QTL mapping of the swimming crab (Portunus trituberculatus) was performed based on a high-density linkage map, and a highly significant QTL specifically mapped on a single linkage group (LG) was firstly identified (LG24, LOD > 14). Twenty markers in the QTL region showed significant associations with sex by association analysis, of which heterogametic genotypes in males supported the XY sex determination mechanism. Two sex-specific markers at the family level were identified via segregation distortion analysis, which were known to be the most closely linked to the sex of P. trituberculatus. Based on sex marker sequences (Marker3840, Marker20320, and Marker10494), three potential sex-related genes were identified, and the quantitative real-time PCR results suggested that these genes were important in spermatogenesis or sex characteristics in males. Our results will contribute to the fine-mapping of sex-determining genes and clarify the sex determination mechanism of P. trituberculatus.

Target sequence capture of nuclear-encoded genes for phylogenetic analysis in ferns.

Premise of the StudyUntil recently, most phylogenetic studies of ferns were based on chloroplast genes. Evolutionary inferences based on these data can be incomplete because the characters are from a single linkage group and are uniparentally inherited. These limitations are particularly acute in studies of hybridization, which is prevalent in ferns; fern hybrids are common and ferns are able to hybridize across highly diverged lineages, up to 60 million years since divergence in one documented case. However, it not yet clear what effect such hybridization has on fern evolution, in part due to a paucity of available biparentally inherited (nuclear‐encoded) markers.MethodsWe designed oligonucleotide baits to capture 25 targeted, low‐copy nuclear markers from a sample of 24 species spanning extant fern diversity.ResultsMost loci were successfully sequenced from most accessions. Although the baits were designed from exon (transcript) data, we successfully captured intron sequences that should be useful for more focused phylogenetic studies. We present phylogenetic analyses of the new target sequence capture data and integrate these into a previous transcript‐based data set.DiscussionWe make our bait sequences available to the community as a resource for further studies of fern phylogeny.

Factors affecting the accuracy of genomic selection for growth and wood quality traits in an advanced-breeding population of black spruce (Picea mariana)

BackgroundGenomic selection (GS) uses information from genomic signatures consisting of thousands of genetic markers to predict complex traits. As such, GS represents a promising approach to accelerate tree breeding, which is especially relevant for the genetic improvement of boreal conifers characterized by long breeding cycles. In the present study, we tested GS in an advanced-breeding population of the boreal black spruce (Picea mariana [Mill.] BSP) for growth and wood quality traits, and concurrently examined factors affecting GS model accuracy.ResultsThe study relied on 734 25-year-old trees belonging to 34 full-sib families derived from 27 parents and that were established on two contrasting sites. Genomic profiles were obtained from 4993 Single Nucleotide Polymorphisms (SNPs) representative of as many gene loci distributed among the 12 linkage groups common to spruce. GS models were obtained for four growth and wood traits. Validation using independent sets of trees showed that GS model accuracy was high, related to trait heritability and equivalent to that of conventional pedigree-based models. In forward selection, gains per unit of time were three times higher with the GS approach than with conventional selection. In addition, models were also accurate across sites, indicating little genotype-by-environment interaction in the area investigated. Using information from half-sibs instead of full-sibs led to a significant reduction in model accuracy, indicating that the inclusion of relatedness in the model contributed to its higher accuracies. About 500 to 1000 markers were sufficient to obtain GS model accuracy almost equivalent to that obtained with all markers, whether they were well spread across the genome or from a single linkage group, further confirming the implication of relatedness and potential long-range linkage disequilibrium (LD) in the high accuracy estimates obtained. Only slightly higher model accuracy was obtained when using marker subsets that were identified to carry large effects, indicating a minor role for short-range LD in this population.ConclusionsThis study supports the integration of GS models in advanced-generation tree breeding programs, given that high genomic prediction accuracy was obtained with a relatively small number of markers due to high relatedness and family structure in the population. In boreal spruce breeding programs and similar ones with long breeding cycles, much larger gain per unit of time can be obtained from genomic selection at an early age than by the conventional approach. GS thus appears highly profitable, especially in the context of forward selection in species which are amenable to mass vegetative propagation of selected stock, such as spruces.

Interspecific hybridization causes long‐term phylogenetic discordance between nuclear and mitochondrial genomes in freshwater fishes

Classification, phylogeography and the testing of evolutionary hypotheses rely on correct estimation of species phylogeny. Early molecular phylogenies often relied on mtDNA alone, which acts as a single linkage group with one history. Over the last decade, the use of multiple nuclear sequences has often revealed conflict among gene trees. This observation can be attributed to hybridization, lineage sorting, paralogy or selection. Here, we use 54 groups of fishes from 48 studies to estimate the degree of concordance between mitochondrial and nuclear gene trees in two ecological grades of fishes: marine and freshwater. We test the hypothesis that freshwater fish phylogenies should, on average, show more discordance because of their higher propensity for hybridization in the past. In keeping with this idea, concordance between mitochondrial and nuclear gene trees (as measured by proportion of components shared) is on average 50% higher in marine fishes. We discuss why this difference almost certainly results from introgression caused by greater historical hybridization among lineages in freshwater groups, and further emphasize the need to use multiple nuclear genes, and identify conflict among them, in estimation of species phylogeny.

Development of a Method to Produce Chromosome Lacking Lines (CLLs) in <i>Nicotiana tabacum</i> L. “Red Russian”

Monosomic lines of Nicotiana tabacum are helpful to confirm the location of genes on specific chromosomes. In the cross N. nudicaulis and N. tabacum, hybrid seedlings express lethal symptoms, which are controlled by the S subgenome of N. tabacum. To identify the responsible chromosome, we needed to produce chromosome lacking lines (CLLs) of N. tabacum L. “Red Russian” and use them to cross with N. nudicaulis. From a cross of (N. tabacum × N. tomentosiformis) × N. tabacum, 380 BC1 individuals were obtained. Using a Haplo-Q line (a monosomic line lacking the single linkage group 11) and N. tabacum, we found that qPCR is a simple and reliable screening method for CLLs of N. tabacum. The marker PT30342 is located on linkage group 11, and the -Ct value (Ct Actin - Ct PT30342) was 2.0 for a disomic line and was 1.097 for a Haplo-Q line. By the use of flow cytometry, qPCR and chromosome counting together as a screening method, we identified 6 CLLs lacking 2 to 6 chromosomes. Compared with conventional methods, our method is a rapid technique for making and screening CLLs ofthe S or S/T subgenome of N. tabacum. Further, these CLLs will be useful to identify the location of two or more factors on chromosomes controlling a variety of genetic problems affecting breeding. Here, we only made CLLs of the S or S/T subgenome of N. tabacum. We will use the method established in this study to produce CLLs of the T subgenome of N. tabacum, and gather a full set of CLLs of N. tabacum. qPCR could also be applied to the identification of chromosome aberrations in other plants.

Second-Generation Linkage Maps for the Pacific Oyster Crassostrea gigas Reveal Errors in Assembly of Genome Scaffolds.

The Pacific oyster Crassostrea gigas, a widely cultivated marine bivalve mollusc, is becoming a genetically and genomically enabled model for highly fecund marine metazoans with complex life-histories. A genome sequence is available for the Pacific oyster, as are first-generation, low-density, linkage and gene-centromere maps mostly constructed from microsatellite DNA markers. Here, higher density, second-generation, linkage maps are constructed from more than 1100 coding (exonic) single-nucleotide polymorphisms (SNPs), as well as 66 previously mapped microsatellite DNA markers, all typed in five families of Pacific oysters (nearly 172,000 genotypes). The map comprises 10 linkage groups, as expected, has an average total length of 588 cM, an average marker-spacing of 1.0 cM, and covers 86% of a genome estimated to be 616 cM. All but seven of the mapped SNPs map to 618 genome scaffolds; 260 scaffolds contain two or more mapped SNPs, but for 100 of these scaffolds (38.5%), the contained SNPs map to different linkage groups, suggesting widespread errors in scaffold assemblies. The 100 misassembled scaffolds are significantly longer than those that map to a single linkage group. On the genetic maps, marker orders and intermarker distances vary across families and mapping methods, owing to an abundance of markers segregating from only one parent, to widespread distortions of segregation ratios caused by early mortality, as previously observed for oysters, and to genotyping errors. Maps made from framework markers provide stronger support for marker orders and reasonable map lengths and are used to produce a consensus high-density linkage map containing 656 markers.

Identification of Single-Locus PCR-Based Markers Linked to Shell Background Color in the Pacific Oyster (Crassostrea gigas).

A number of Pacific oyster (Crassostrea gigas) with golden shell background color were obtained which show great potential to develop a niche market. To improve the selective breeding progress of true-breeding strains with complete golden oysters, research was conducted to identify genetic markers linked to the shell color locus. An F1-segregating population was obtained by crossing two oysters with golden shell and white shell. Genomic DNA from eight progenies with golden shell and eight progenies with white shell were equally pooled for amplified fragment length polymorphism (AFLP) screening. In bulked segregant analysis, six out of 225 selective primer pair combinations produced seven polymorphic fragments tightly linked to shell color across the segregating population. The seven AFLP markers were all derived from the golden dam and mapped onto a single linkage group flanking the shell color locus. In conversion of the AFLPs into single-locus PCR-based markers, a sequence-characterized amplified region (SCAR) marker, named SCARJ8-2, a single nucleotide polymorphism (SNP) marker, named SNPL2-4, and a simple sequence repeat (SSR) marker, named SSRO11-2, were obtained. These markers obtained in this study will be useful for marker-assisted selection of the Pacific oyster.

SimpleMap: A Pipeline to Streamline High-Density Linkage Map Construction.

The recent development in high-throughput genotyping techniques requires new statistical methods to analyze large datasets. The current available linkage mapping software are time consuming and limited in terms of the maximum number of markers that can be mapped on a single linkage group. In this paper, we propose the Perl pipeline, SimpleMap. This tool can significantly improve the speed of currently available linkage mapping software with minimal impact on marker order and map length by limiting the consideration of duplicated and tightly linked molecular markers during linkage group development. SimpleMap works with the following three main steps: (i) generating a subset of markers for which each pair has a number of recombinants higher than a threshold determined by the user (the repulsion threshold), (ii) mapping this subset with any external mapping tool, and (iii) intersecting the remaining unmapped markers to the constructed map. The script was tested on 15 wheat (Triticum aestivum L.) linkage groups derived from two different crosses. In 13 genetic groups, the computational time was reduced from ∼8 h to ∼8 min, while it was impossible to map the remaining two linkage groups without applying SimpleMap first. SimpleMap is a very time-efficient tool, and considering a repulsion threshold equivalent to 1 cM results in a number of markers similar to map lengths that can be analyzed on a simple personal computer. SimpleMap can be downloaded from http://simplemap-aj.sourceforge.net/.

High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis.

The sex determination system in crabs is believed to be XY-XX from karyotypy, but centromeres could not be identified in some chromosomes and their morphology is not completely clear. Using quantitative trait locus mapping of the gender phenotype, we revealed a ZW-ZZ sex determination system in Eriocheir sinensis and presented a high-density linkage map covering ~98.5% of the genome, with 73 linkage groups corresponding to the haploid chromosome number. All sex-linked markers in the family we used were located on a single linkage group, LG60, and sex linkage was confirmed by genome-wide association studies (GWAS). Forty-six markers detected by GWAS were heterozygous and segregated only in the female parent. The female LG60 was thus the putative W chromosome, with the homologous male LG60 as the Z chromosome. The putative Z and W sex chromosomes were identical in size and carried many homologous loci. Sex ratio (5:1) skewing towards females in induced triploids using unrelated animals also supported a ZW-ZZ system. Transcriptome data were used to search for candidate sex-determining loci, but only one LG60 gene was identified as an ankyrin-2 gene. Double sex- and mab3-related transcription factor 1 (Dmrt1), a Z-linked gene in birds, was located on a putative autosome. With complete genome sequencing and transcriptomic data, more genes on putative sex chromosomes will be characterised, thus leading towards a comprehensive understanding of the sex determination and differentiation mechanisms of E. sinensis, and decapod crustaceans in general.

Genetic markers for western corn rootworm resistance to Bt toxin.

Western corn rootworm (WCR) is a major maize (Zea mays L.) pest leading to annual economic losses of more than 1 billion dollars in the United States. Transgenic maize expressing insecticidal toxins derived from the bacterium Bacillus thuringiensis (Bt) are widely used for the management of WCR. However, cultivation of Bt-expressing maize places intense selection pressure on pest populations to evolve resistance. Instances of resistance to Bt toxins have been reported in WCR. Developing genetic markers for resistance will help in characterizing the extent of existing issues, predicting where future field failures may occur, improving insect resistance management strategies, and in designing and sustainably implementing forthcoming WCR control products. Here, we discover and validate genetic markers in WCR that are associated with resistance to the Cry3Bb1 Bt toxin. A field-derived WCR population known to be resistant to the Cry3Bb1 Bt toxin was used to generate a genetic map and to identify a genomic region associated with Cry3Bb1 resistance. Our results indicate that resistance is inherited in a nearly recessive manner and associated with a single autosomal linkage group. Markers tightly linked with resistance were validated using WCR populations collected from Cry3Bb1 maize fields showing significant WCR damage from across the US Corn Belt. Two markers were found to be correlated with both diet (R2 = 0.14) and plant (R2 = 0.23) bioassays for resistance. These results will assist in assessing resistance risk for different WCR populations, and can be used to improve insect resistance management strategies.

Computationally efficient map construction in the presence of segregation distortion.

We present a novel estimator for map construction in the presence of segregation distortion which is highly computationally efficient. For multi-parental designs this estimator outperforms methods that do not account for segregation distortion, at no extra computational cost. Inclusion of genetic markers exhibiting segregation distortion in a linkage map can result in biased estimates of genetic distance and distortion of map positions. Removal of distorted markers is hence a typical filtering criterion; however, this may result in exclusion of biologically interesting regions of the genome such as introgressions and translocations. Estimation of additional parameters characterizing the distortion is computationally slow, as it relies on estimation via the Expectation Maximization algorithm or a higher dimensional numerical optimisation. We propose a robust M-estimator (RM) capable of handling tens of thousands of distorted markers from a single linkage group. We show via simulation that for multi-parental designs the RM estimator can perform much better than uncorrected estimation, at no extra computational cost. We then apply the RM estimator to chromosome 2B in wheat in a multi-parent population segregating for the Sr36 introgression, a known transmission distorter. The resulting map contains over 700 markers, and is consistent with maps constructed from crosses which do not exhibit segregation distortion.

Identification and mapping of a SCAR marker linked to a locus involved in shell pigmentation of the Pacific oyster (Crassostrea gigas)

The Pacific oyster (Crassostrea gigas) is one of the most important farmed oysters worldwide. To improve breeding progress, research was done to identify molecular markers linked to shell pigmentation, which has been viewed as a key factor for product value. An F1-segregating population exhibiting a bimodal distribution of shell pigmentation intensity was obtained by crossing two wild oysters with opposite shell pigmentation. Genomic DNA from nine individuals with lightest shell pigmentation and nine individuals with darkest shell pigmentation were equally pooled for amplified fragment length polymorphism (AFLP) screening. In bulked segregant analysis, six out of 225 selective primer pair combinations produced seven polymorphic fragments tightly associated with shell pigmentation across the segregating population. The seven AFLP markers were mapped onto a single linkage group, and 80% of phenotypic variance could be explained by this locus. In conversion of the seven fragments into sequence-characterized amplified region (SCAR) markers, only one was successfully converted into a co-dominant SCAR marker, named SP-170. The SCAR marker was integrated into a previously constructed linkage map. The SCAR marker obtained in this study will be useful for marker-assisted selection of the Pacific oyster.

QTL mapping reveals the genetic architecture of loci affecting pre- and post-zygotic isolating barriers in Louisiana Iris

BackgroundHybridization among Louisiana Irises has been well established and the genetic architecture of reproductive isolation is known to affect the potential for and the directionality of introgression between taxa. Here we use co-dominant markers to identify regions where QTL are located both within and between backcross maps to compare the genetic architecture of reproductive isolation and fitness traits across treatments and years.ResultsQTL mapping was used to elucidate the genetic architecture of reproductive isolation between Iris fulva and Iris brevicaulis. Homologous co-dominant EST-SSR markers scored in two backcross populations between I. fulva and I. brevicaulis were used to generate genetic linkage maps. These were used as the framework for mapping QTL associated with variation in 11 phenotypic traits likely responsible for reproductive isolation and fitness. QTL were dispersed throughout the genome, with the exception of one region of a single linkage group (LG) where QTL for flowering time, sterility, and fruit production clustered. In most cases, homologous QTL were not identified in both backcross populations, however, homologous QTL for flowering time, number of growth points per rhizome, number of nodes per inflorescence, and number of flowers per node were identified on several linkage groups.ConclusionsTwo different traits affecting reproductive isolation, flowering time and sterility, exhibit different genetic architectures, with numerous QTL across the Iris genome controlling flowering time and fewer, less distributed QTL affecting sterility. QTL for traits affecting fitness are largely distributed across the genome with occasional overlap, especially on LG 4, where several QTL increasing fitness and decreasing sterility cluster. Given the distribution and effect direction of QTL affecting reproductive isolation and fitness, we have predicted genomic regions where introgression may be more likely to occur (those regions associated with an increase in fitness and unlinked to loci controlling reproductive isolation) and those that are less likely to exhibit introgression (those regions linked to traits decreasing fitness and reproductive isolation).

A high-resolution map of the Nile tilapia genome: a resource for studying cichlids and other percomorphs

BackgroundThe Nile tilapia (Oreochromis niloticus) is the second most farmed fish species worldwide. It is also an important model for studies of fish physiology, particularly because of its broad tolerance to an array of environments. It is a good model to study evolutionary mechanisms in vertebrates, because of its close relationship to haplochromine cichlids, which have undergone rapid speciation in East Africa. The existing genomic resources for Nile tilapia include a genetic map, BAC end sequences and ESTs, but comparative genome analysis and maps of quantitative trait loci (QTL) are still limited.ResultsWe have constructed a high-resolution radiation hybrid (RH) panel for the Nile tilapia and genotyped 1358 markers consisting of 850 genes, 82 markers corresponding to BAC end sequences, 154 microsatellites and 272 single nucleotide polymorphisms (SNPs). From these, 1296 markers could be associated in 81 RH groups, while 62 were not linked. The total size of the RH map is 34,084 cR3500 and 937,310 kb. It covers 88% of the entire genome with an estimated inter-marker distance of 742 Kb. Mapping of microsatellites enabled integration to the genetic map. We have merged LG8 and LG24 into a single linkage group, and confirmed that LG16-LG21 are also merged. The orientation and association of RH groups to each chromosome and LG was confirmed by chromosomal in situ hybridizations (FISH) of 55 BACs. Fifty RH groups were localized on the 22 chromosomes while 31 remained small orphan groups. Synteny relationships were determined between Nile tilapia, stickleback, medaka and pufferfish.ConclusionThe RH map and associated FISH map provide a valuable gene-ordered resource for gene mapping and QTL studies. All genetic linkage groups with their corresponding RH groups now have a corresponding chromosome which can be identified in the karyotype. Placement of conserved segments indicated that multiple inter-chromosomal rearrangements have occurred between Nile tilapia and the other model fishes. These maps represent a valuable resource for organizing the forthcoming genome sequence of Nile tilapia, and provide a foundation for evolutionary studies of East African cichlid fishes.

QTL linkage analysis of connected populations using ancestral marker and pedigree information

The common assumption in quantitative trait locus (QTL) linkage mapping studies that parents of multiple connected populations are unrelated is unrealistic for many plant breeding programs. We remove this assumption and propose a Bayesian approach that clusters the alleles of the parents of the current mapping populations from locus-specific identity by descent (IBD) matrices that capture ancestral marker and pedigree information. Moreover, we demonstrate how the parental IBD data can be incorporated into a QTL linkage analysis framework by using two approaches: a Threshold IBD model (TIBD) and a Latent Ancestral Allele Model (LAAM). The TIBD and LAAM models are empirically tested via numerical simulation based on the structure of a commercial maize breeding program. The simulations included a pilot dataset with closely linked QTL on a single linkage group and 100 replicated datasets with five linkage groups harboring four unlinked QTL. The simulation results show that including parental IBD data (similarly for TIBD and LAAM) significantly improves the power and particularly accuracy of QTL mapping, e.g., position, effect size and individuals’ genotype probability without significantly increasing computational demand.