The genetic map and comparative analysis with the physical map of Trypanosoma brucei

Trypanosoma brucei is the causative agent of African sleeping sickness in humans and contributes to the debilitating disease ‘Nagana’ in cattle. To date we know little about the genes that determine drug resistance, host specificity, pathogenesis and virulence in these parasites. The availability of the complete genome sequence and the ability of the parasite to undergo genetic exchange have allowed genetic investigations into this parasite and here we report the first genetic map of T.brucei for the genome reference stock TREU 927, comprising of 182 markers and 11 major linkage groups, that correspond to the 11 previously identified chromosomes. The genetic map provides 90% probability of a marker being 11 cM from any given locus. Its comparison to the available physical map has revealed the average physical size of a recombination unit to be 15.6 Kb/cM. The genetic map coupled with the genome sequence and the ability to undertake crosses presents a new approach to identifying genes relevant to the disease and its prevention in this important pathogen through forward genetic analysis and positional cloning.     

Whole-Genome Sequencing of Theileria parva Strains Provides Insight into Parasite Migration and Diversification in the African Continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.     

A re-assigned American mink (Neovison vison) map optimal for genome-wide studies

Our previously published second generation genetic map for the American mink (Neovison vison) has been used and redesigned in its best for genome-wide studies with maximum of efficiency. A number of 114 selected markers, including 33 newly developed microsatellite markers from the CHORI-231 mink Bacterial Artificial Chromosome (BAC) library, have been genotyped in a two generation population composed of 1200 individuals. The outcome reassigns the position of some markers on the chromosomes and it produces a more reliable map with a convenient distance between markers. A total of 104 markers mapped to 14 linkage groups corresponding to the mink autosomes. Six markers are unlinked and four markers are allocated to the X chromosome by homology but no linkage was detected. The sex-average linkage map spans 1192 centiMorgans (cM) with an average intermarker distance of 11.4 cM and 1648 cM when the ends of the linkage groups and the autosomal unlinked markers are added. Sex-specific genetic linkage maps were also generated. The male sex-specific map had a total length of 1014.6 cM between the linked markers and an average inter-marker interval of 9.7 cM. The female map has a corresponding length of 1378.6 cM and an average inter-marker interval of 13.3 cM. The study is complemented with additional anchorage for most of the chromosomes of the map by BAC in situ hybridization with clones containing microsatellites strategically selected from the various parts of the genome. This map provides an improved tool for genetic mapping and comparative genomics in mink, also useful for the future assembly of the mink genome sequence when this will be taken forward.     

A genetic linkage map for the apicomplexan protozoan parasite Eimeria maxima and comparison with Eimeria tenella

Eimeria maxima is one of the seven Eimeria spp. that infect the chicken and cause the disease coccidiosis. The well characterised immunogenicity and genetic diversity associated with E. maxima promote its use in genetics-led studies on avian coccidiosis. The development of a genetic map for E. maxima, presented here based upon 647 amplified fragment length polymorphism markers typed from 22 clonal hybrid lines and assembled into 13 major linkage groups, is a major new resource for work with this parasite. Comparison with genetic maps produced for other coccidial parasites indicates relatively high levels of genetic recombination. Conversion of ∼14% of the markers representing the major linkage groups to sequence characterised amplified region markers can provide a scaffold for the assembly of future genomic sequences as well as providing a foundation for more detailed genetic maps. Comparison with the Eimeria tenella genetic map produced 10 years ago has revealed a less biased marker distribution, with no more than nine markers mapped within any unresolved heritable unit. Nonetheless, preliminary bioinformatic characterisation of the three largest publicly available genomic E. maxima sequences suggest that the feature-poor/feature-rich structure which has previously been found to define the first sequenced E. tenella chromosome also defines the E. maxima genome. The significance of such a segmented genome and the apparent potential for variation in genetic recombination will be relevant to haplotype stability and the longevity of future anticoccidial strategies based upon multiple loci targeted by novel chemotherapeutic drugs or recombinant subunit vaccines.     

Recombination and Its Impact on the Genome of the Haplodiploid Parasitoid Wasp Nasonia

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4–1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.     

A New Genetic Linkage Map of the Zygomycete Fungus Phycomyces blakesleeanus

Phycomyces blakesleeanus is a member of the subphylum Mucoromycotina. A genetic map was constructed from 121 progeny of a cross between two wild type isolates of P. blakesleeanus with 134 markers. The markers were mostly PCR-RFLPs. Markers were located on 46 scaffolds of the genome sequence, covering more than 97% of the genome. Analysis of the alleles in the progeny revealed nine or 12 linkage groups, depending on the log of the odds (LOD) score, across 1583.4 cM at LOD 5. The linkage groups were overlaid on previous mapping data from crosses between mutants, aided by new identification of the mutations in primary metabolism mutant strains. The molecular marker map, the phenotype map and the genome sequence are overall congruent, with some exceptions. The new genetic map provides a genome-wide estimate for recombination, with the average of 33.2 kb per cM. This frequency is one piece of evidence for meiosis during zygospore development in Mucoromycotina species. At the same time as meiosis, transmission of non-recombinant chromosomes is also evident in the mating process in Phycomyces. The new map provides scaffold ordering for the genome sequence and a platform upon which to identify the genes in mutants that are affected in traits of interest, such as carotene biosynthesis, phototropism or gravitropism, using positional cloning.     

Development of forward genetics in Toxoplasma gondii

The development of forward genetics as a functional system in Toxoplasma gondii spanned more than three decades from the mid-1970s until now. The initial demonstration of experimental genetics relied on chemically induced drug-resistant mutants that were crossed by co-infecting cats, collecting oocysts, sporulating and hatching progeny in vitro. To capitalise on this, genetic markers were employed to develop linkage maps by tracking inheritance through experimental crosses. In all, three generations of genetic maps were developed to define the chromosomes, estimate recombination rates and provide a system for linkage analysis. Ultimately this genetic map would become the foundation for the assembly of the T. gondii genome, which was derived from whole genome shotgun sequencing, into a chromosome-centric view. Finally, application of forward genetics to multigenic biological traits showed the potential to map and identify specific genes that control complex phenotypes including virulence.     

A genetic map of Peromyscus with chromosomal assignment of linkage groups (a Peromyscus genetic map)

The rodent genus Peromyscus is the most numerous and species-rich mammalian group in North America. The naturally occurring diversity within this genus allows opportunities to investigate the genetic basis of adaptation, monogamy, behavioral and physiological phenotypes, growth control, genomic imprinting, and disease processes. Increased genomic resources including a high quality genetic map are needed to capitalize on these opportunities. We produced interspecific hybrids between the prairie deer mouse (P. maniculatus bairdii) and the oldfield mouse (P. polionotus) and scored meiotic recombination events in backcross progeny. A genetic map was constructed by genotyping of backcross progeny at 185 gene-based and 155 microsatellite markers representing all autosomes and the X-chromosome. Comparison of the constructed genetic map with the molecular maps of Mus and Rattus and consideration of previous results from interspecific reciprocal whole chromosome painting allowed most linkage groups to be unambiguously assigned to specific Peromyscus chromosomes. Based on genomic comparisons, this Peromyscus genetic map covers ~83 % of the Rattus genome and 79 % of the Mus genome. This map supports previous results that the Peromyscus genome is more similar to Rattus than Mus. For example, coverage of the 20 Rattus autosomes and the X-chromosome is accomplished with only 28 segments of the Peromyscus map, but coverage of the 19 Mus autosomes and the X-chromosome requires 40 chromosomal segments of the Peromyscus map. Furthermore, a single Peromyscus linkage group corresponds to about 91 % of the rat and only 76 % of the mouse X-chromosomes.     

An expanded genetic linkage map of an intervarietal Agaricus bisporus var. bisporus × A. bisporus var. burnettii hybrid based on AFLP,SSR and CAPS markers sheds light on the recombination behaviour of the species

A genetic linkage map for the edible basidiomycete Agaricus bisporus was constructed from 118 haploid homokaryons derived from an intervarietal A. bisporus var. bisporus × A. bisporus var. burnettii hybrid. Two hundred and thirty-one AFLP, 21 SSR, 68 CAPS markers together with the MAT, BSN, PPC1 loci and one allozyme locus (ADH) were evenly spread over 13 linkage groups corresponding to the chromosomes of A. bisporus. The map covers 1156 cM, with an average marker spacing of 3.9 cM and encompasses nearly the whole genome. The average number of crossovers per chromosome per individual is 0.86. Normal recombination over the entire genome occurs in the heterothallic variety, burnettii, contrary to the homothallic variety, bisporus, which showed adaptive genome-wide suppressed recombination. This first comprehensive genetic linkage map for A. bisporus provides foundations for quantitative trait analyses and breeding programme monitoring, as well as genome organisation studies.     

A Consensus Linkage Map Provides Insights on Genome Character and Evolution in Common Carp (Cyprinus carpio L.)

Common carp (Cyprinus carpio L.) is cultured worldwide and is a major contributor to the world’s aquaculture production. The common carp has a complex tetraploidized genome, which may historically experience additional whole genome duplication than most other Cyprinids. Fine maps for female and male carp were constructed using a mapping panel containing one F1 family with 190 progeny. A total of 1,025 polymorphic markers were used to construct genetic maps. For the female map, 559 microsatellite markers in 50 linkage groups cover 3,468 cM of the genome. For the male map, 383 markers in 49 linkage groups cover 1,811 cM of the genome. The consensus map was constructed by integrating the new map with two published linkage maps, containing 732 markers and spanning 3,278 cM in 50 linkage groups. The number of consensus linkage groups corresponds to the number of common carp chromosomes. A significant difference on sex recombinant rate was observed that the ratio of female and male recombination rates was 4.2:1. Comparative analysis was performed between linkage map of common carp and genome of zebrafish (Danio rerio), which revealed clear 2:1 relationship of common carp linkage groups and zebrafish chromosomes. The results provided evidence that common carp did experienced a specific whole genome duplication event comparing with most other Cyprinids. The consensus linkage map provides an important tool for genetic and genome study of common carp and facilitates genetic selection and breeding for common carp industry.     

Recombination hotspots flank the Cryptococcus mating-type locus: implications for the evolution of a fungal sex chromosome

Recombination increases dramatically during meiosis to promote genetic exchange and generate recombinant progeny. Interestingly, meiotic recombination is unevenly distributed throughout genomes, and, as a consequence, genetic and physical map distances do not have a simple linear relationship. Recombination hotspots and coldspots have been described in many organisms and often reflect global features of chromosome structure. In particular, recombination frequencies are often distorted within or outside sex-determining regions of the genome. Here, we report that recombination is elevated adjacent to the mating-type locus (MAT) in the pathogenic basidiomycete Cryptococcus neoformans. Among fungi, C. neoformans has an unusually large MAT locus, and recombination is suppressed between the two >100-kilobase mating-type specific alleles. When genetic markers were introduced at defined physical distances from MAT, we found the meiotic recombination frequency to be ~20% between MAT and a flanking marker at 5, 10, 50, or 100 kilobases from the right border. As a result, the physical/genetic map ratio in the regions adjacent to MAT is distorted ~10- to 50-fold compared to the genome-wide average. Moreover, recombination frequently occurred on both sides of MAT and negative interference between crossovers was observed. MAT heterozygosity was not required for enhanced recombination, implying that this process is not due to a physical distortion from the two non-paired alleles and could also occur during same-sex mating. Sequence analysis revealed a correlation between high G + C content and these hotspot regions. We hypothesize that the presence of recombinational activators may have driven several key events during the assembly and reshaping of the MAT locus and may have played similar roles in the origins of both metabolic and biosynthetic gene clusters. Our findings suggest that during meiosis the MAT locus may be exchanged onto different genetic backgrounds and therefore have broad evolutionary implications with respect to mating-type switching in both model and pathogenic yeasts.     

Absolute Quantification of the Host-To-Parasite DNA Ratio in Theileria parva-Infected Lymphocyte Cell Lines

Theileria parva is a tick-transmitted intracellular apicomplexan pathogen of cattle in sub-Saharan Africa that causes East Coast fever (ECF). ECF is an acute fatal disease that kills over one million cattle annually, imposing a tremendous burden on African small-holder cattle farmers. The pathology and level of T. parva infections in its wildlife host, African buffalo (Syncerus caffer), and in cattle are distinct. We have developed an absolute quantification method based on quantitative PCR (qPCR) in which recombinant plasmids containing single copy genes specific to the parasite (apical membrane antigen 1 gene, ama1) or the host (hypoxanthine phosphoribosyltransferase 1, hprt1) are used as the quantification reference standards. Our study shows that T. parva and bovine cells are present in similar numbers in T. parva-infected lymphocyte cell lines and that consequently, due to its much smaller genome size, T. parva DNA comprises between 0.9% and 3% of the total DNA samples extracted from these lines. This absolute quantification assay of parasite and host genome copy number in a sample provides a simple and reliable method of assessing T. parva load in infected bovine lymphocytes, and is accurate over a wide range of host-to-parasite DNA ratios. Knowledge of the proportion of target DNA in a sample, as enabled by this method, is essential for efficient high-throughput genome sequencing applications for a variety of intracellular pathogens. This assay will also be very useful in future studies of interactions of distinct host-T. parva stocks and to fully characterize the dynamics of ECF infection in the field.     

A locus conferring tolerance to Theileria infection in African cattle

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Bos indicus (Boran) cattle demonstrating heritable tolerance to infection with T. parva (h² = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on bovine chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained variant in a paralogue of the FAF1 gene in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa’s most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.     

SNP Discovery and Chromosome Anchoring Provide the First Physically-Anchored Hexaploid Oat Map and Reveal Synteny with Model Species

A physically anchored consensus map is foundational to modern genomics research; however, construction of such a map in oat (Avena sativa L., 2n = 6x = 42) has been hindered by the size and complexity of the genome, the scarcity of robust molecular markers, and the lack of aneuploid stocks. Resources developed in this study include a modified SNP discovery method for complex genomes, a diverse set of oat SNP markers, and a novel chromosome-deficient SNP anchoring strategy. These resources were applied to build the first complete, physically-anchored consensus map of hexaploid oat. Approximately 11,000 high-confidence in silico SNPs were discovered based on nine million inter-varietal sequence reads of genomic and cDNA origin. GoldenGate genotyping of 3,072 SNP assays yielded 1,311 robust markers, of which 985 were mapped in 390 recombinant-inbred lines from six bi-parental mapping populations ranging in size from 49 to 97 progeny. The consensus map included 985 SNPs and 68 previously-published markers, resolving 21 linkage groups with a total map distance of 1,838.8 cM. Consensus linkage groups were assigned to 21 chromosomes using SNP deletion analysis of chromosome-deficient monosomic hybrid stocks. Alignments with sequenced genomes of rice and Brachypodium provide evidence for extensive conservation of genomic regions, and renewed encouragement for orthology-based genomic discovery in this important hexaploid species. These results also provide a framework for high-resolution genetic analysis in oat, and a model for marker development and map construction in other species with complex genomes and limited resources.     

A microsatellite-based genetic linkage map of the cichlid fish, Astatotilapia burtoni (Teleostei): a comparison of genomic architectures among rapidly speciating cichlids

Cichlid fishes compose an astonishingly large number of species and formed species flocks in record-breaking time. To facilitate efficient genome scans and comparisons of cichlid genomes, we constructed a medium-density genetic linkage map of microsatellite markers of Astatotilapia burtoni. The mapping cross was derived from two inbred laboratory lines to obtain F₂ progeny by intercrossing. The map revealed 25 linkage groups spanning 1249.3 cM of the genome (size ~950 Mb) with an average marker spacing of 6.12 cM. The seven Hox clusters, ParaHox C1, and two paralogs of Pdgfrβ were mapped to different linkage groups, thus supporting the hypothesis of a teleost-specific genome duplication. The A. burtoni linkage map was compared to the other two available maps for cichlids using shared markers that showed conservation and synteny among East African cichlid genomes. Interesting candidate genes for cichlid speciation were mapped using SNP markers.     

An AFLP-based genetic linkage map of Plasmodium chabaudi chabaudi

Background

Plasmodium chabaudi chabaudi can be considered as a rodent model of human malaria parasites in the genetic analysis of important characters such as drug resistance and immunity. Despite the availability of some genome sequence data, an extensive genetic linkage map is needed for mapping the genes involved in certain traits.

Methods

The inheritance of 672 Amplified Fragment Length Polymorphism (AFLP) markers from two parental clones (AS and AJ) of P. c. chabaudi was determined in 28 independent recombinant progeny clones. These, AFLP markers and 42 previously mapped Restriction Fragment Length Polymorphism (RFLP) markers (used as chromosomal anchors) were organized into linkage groups using Map Manager software.

Results

614 AFLP markers formed linkage groups assigned to 10 of 14 chromosomes, and 12 other linkage groups not assigned to known chromosomes. The genetic length of the genome was estimated to be about 1676 centiMorgans (cM). The mean map unit size was estimated to be 13.7 kb/cM. This was slightly less then previous estimates for the human malaria parasite, Plasmodium falciparum

Conclusion

The P. c. chabaudi genetic linkage map presented here is the most extensive and highly resolved so far available for this species. It can be used in conjunction with the genome databases of P. c chabaudi, P. falciparum and Plasmodium yoelii to identify genes underlying important phenotypes such as drug resistance and strain-specific immunity.     

Loop-mediated isothermal amplification (LAMP) assays for detection of Theileria parva infections targeting the PIM and p150 genes

We have developed two loop-mediated isothermal amplification (LAMP) assays for the detection of Theileria parva, the causative agent of East Coast fever (ECF), an economically important cattle disease in eastern, central and southern Africa. These assays target the polymorphic immunodominant molecule (PIM) and p150 LAMP genes. The primer set for each gene target consists of six primers, and each set recognises eight distinct regions on the target gene to give highly specific detection of T. parva. The detection limit of each primer set is 1 fg, which is equivalent to one copy of the PIM and p150 T. parva genes. These PIM and p150 LAMP primer sets amplify DNA of T. parva isolates from cattle and buffalo from different countries including Kenya, South Africa, Tanzania, Rwanda, Uganda and Burundi, indicating their ability to detect T. parva from different countries. With the advantages of simplicity, rapidity and cost effectiveness, these LAMP assays are good candidates for molecular epidemiology studies and for monitoring control programs in ECF-endemic, resource poor countries.     

Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.     

Composite genome map and recombination parameters derived from three archetypal lineages of Toxoplasma gondii

Toxoplasma gondii is a highly successful protozoan parasite in the phylum Apicomplexa, which contains numerous animal and human pathogens. T.gondii is amenable to cellular, biochemical, molecular and genetic studies, making it a model for the biology of this important group of parasites. To facilitate forward genetic analysis, we have developed a high-resolution genetic linkage map for T.gondii. The genetic map was used to assemble the scaffolds from a 10X shotgun whole genome sequence, thus defining 14 chromosomes with markers spaced at ~300 kb intervals across the genome. Fourteen chromosomes were identified comprising a total genetic size of ~592 cM and an average map unit of ~104 kb/cM. Analysis of the genetic parameters in T.gondii revealed a high frequency of closely adjacent, apparent double crossover events that may represent gene conversions. In addition, we detected large regions of genetic homogeneity among the archetypal clonal lineages, reflecting the relatively few genetic outbreeding events that have occurred since their recent origin. Despite these unusual features, linkage analysis proved to be effective in mapping the loci determining several drug resistances. The resulting genome map provides a framework for analysis of complex traits such as virulence and transmission, and for comparative population genetic studies.     

Transmission ratio distortion in an interspecific cross between Fusarium circinatum and Fusarium subglutinans

Previously, an interspecific cross between Fusarium circinatum and Fusarium subglutinans was used to generate a genetic linkage map. A ca. 55 % of genotyped markers displayed transmission ratio distortion (TRD) that demonstrated a genome-wide distribution. The working hypothesis for this study was that TRD would be non-randomly distributed throughout the genetic linkage map. This would indicate the presence of distorting loci. Using a genome-wide threshold of α = 0.01, 79 markers displaying TRD were distributed on all 12 linkage groups (LGs). Eleven putative transmission ratio distortion loci (TRDLs), spanning eight LGs, were identified in regions containing three or more adjacent markers displaying distortion. No epistatic interactions were observed between these TRDLs. Thus, it is uncertain whether the genome-wide TRD was due to linkage between markers and genomic regions causing distortion. The parental origins of markers followed a non-random distribution throughout the linkage map, with LGs containing stretches of markers originating from only one parent. Thus, due to the nature of the interspecific cross, the current hypothesis to explain these observations is that the observed genome-wide segregation was caused by the high level of genomic divergence between the parental isolates. Therefore, homologous chromosomes do not align properly during meiosis, resulting in aberrant transmission of markers. This also explains previous observations of the preferential transmission of F. subglutinans alleles to the F₁ progeny.