The design and application of a 50 K SNP chip for a threatened Aotearoa New Zealand passerine,the hihi

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or “SNP chips”, enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50 K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10 WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536 hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species.     

Genetic linkage maps of Populus nigra L. including AFLPs, SSRs, SNPs, and sex trait

Black poplar (Populus nigra L.) is a tree of ecological and economic interest. A better knowledge of P. nigra genome is needed for an effective protection and use of its genetic resources. The main objective of this study is the construction of a highly informative genetic map of P. nigra species including genes of adaptive and economic interest. Two genotypes originated from contrasted natural Italian populations were crossed to generate a F₁ mapping pedigree of 165 individuals. Amplification fragment length polymorphism (AFLP), simple sequence repeat (SSR), and single nucleotide polymorphism (SNP) markers were used to genotype 92 F₁ individuals, and the pseudo-test-cross strategy was applied for linkage analysis. The female parent map included 368 markers (274 AFLPs, 91 SSRs, and 3 SNPs) and spanned 2,104 cM with 20 linkage groups, and the male parent map, including 317 markers (205 AFLPs, 106 SSRs, 5 SNPs, and sex trait), spanned 2,453 cM with 23 main linkage groups. The sex, as morphological trait, was mapped on the linkage group XIX of the male parent map. The generated maps are among the most informative in SSRs when compared to the Populus maps published so far and allow a complete alignment with the 19 haploid chromosomes of Populus sequence genome. These genetic maps provide informative tools for a better understanding of P. nigra genome structure and genetic improvement of this ecologically and economically important European tree species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Mating system and environmental variation drive patterns of adaptation in Boechera spatifolia (Brassicaceae)

Determining the relative contribution of population genetic processes to the distribution of natural variation is a major goal of evolutionary biology. Here, we take advantage of variation in mating system to test the hypothesis that local adaptation is constrained by asexual reproduction. We explored patterns of variation in ecological traits and genome‐wide molecular markers in Boechera spatifolia (Brassicaceae), a species that contains both apomictic (asexual) and sexual individuals. Using a combination of quantitative genetics, neutral genetic (SSR) and genome‐wide single nucleotide polymorphism, we assessed the hypothesis that asexual lineages should have reduced signatures of adaptation relative to sexual conspecifics. All three measures (traits, SSRs, SNPs) demonstrated that apomicts are genetically distinct from sexuals, regardless of population location. Additionally, phylogenetic clustering revealed that the apomictic group shared a single common ancestor. Across the landscape, sexual genome‐wide SNP variation was strongly associated with latitude (r² > 0.9), indicating that sexual populations have differentiated across an environmental gradient. Furthermore, flowering time and growth rate, as assessed in a common garden, strongly covary with the elevation and latitude of the source population. Despite a wide geographic distribution that largely overlaps with sexual populations, there was little evidence for differentiation in molecular markers or quantitative characters among apomictic populations. Combined, these data indicated that, in contrast to asexual populations, sexual populations show evidence of local adaptation.     

Identification, characterisation and application of single nucleotide polymorphisms for diversity assessment in cassava (Manihot esculenta Crantz)

To monitor genetic diversity in the field it is important that it is measured accurately. Here, we elucidate the potential of single nucleotide polymorphisms (SNPs) for measuring genetic diversity in cassava. The nature and frequency of SNPs was characterised and their utility in genetic diversity assessment compared to that of simple sequence repeats (SSRs). This was achieved by direct sequencing of amplicons in diverse cassava varieties. A total of 26 SNPs were identified from quality sequences of nine genes, giving an estimated frequency of one SNP every 121 nucleotides. Nucleotide diversity ranged from 7.8 × 10⁻⁴ to 5.6 × 10⁻³. Average haplotype-based polymorphic information content (PIC = 0.414) was higher than for individual SNPs (PIC = 0.228). The Mantel test indicated interdependence (r = 0.219; P < 0.001) between SNP and SSR genotypic data. Individual SNPs had lower PIC values than SSRs. For this reason larger numbers of SNPs may be necessary to achieve the same level of discrimination among genotypes provided by SSRs.     

Single‐nucleotide polymorphism discovery in Leptographium longiclavatum,a mountain pine beetle‐associated symbiotic fungus,using whole‐genome resequencing

Single‐nucleotide polymorphisms (SNPs) are rapidly becoming the standard markers in population genomics studies; however, their use in nonmodel organisms is limited due to the lack of cost‐effective approaches to uncover genome‐wide variation, and the large number of individuals needed in the screening process to reduce ascertainment bias. To discover SNPs for population genomics studies in the fungal symbionts of the mountain pine beetle (MPB), we developed a road map to discover SNPs and to produce a genotyping platform. We undertook a whole‐genome sequencing approach of Leptographium longiclavatum in combination with available genomics resources of another MPB symbiont, Grosmannia clavigera. We sequenced 71 individuals pooled into four groups using the Illumina sequencing technology. We generated between 27 and 30 million reads of 75 bp that resulted in a total of 1, 181 contigs longer than 2 kb and an assembled genome size of 28.9 Mb (N₅₀ = 48 kb, average depth = 125x). A total of 9052 proteins were annotated, and between 9531 and 17 266 SNPs were identified in the four pools. A subset of 206 genes (containing 574 SNPs, 11% false positives) was used to develop a genotyping platform for this species. Using this roadmap, we developed a genotyping assay with a total of 147 SNPs located in 121 genes using the Illumina^® Sequenom iPLEX Gold. Our preliminary genotyping (success rate = 85%) of 304 individuals from 36 populations supports the utility of this approach for population genomics studies in other MPB fungal symbionts and other fungal nonmodel species.     

SNP discovery using Paired‐End RAD‐tag sequencing on pooled genomic DNA of Sisymbrium austriacum (Brassicaceae)

Single nucleotide polymorphisms SNPs are rapidly replacing anonymous markers in population genomic studies, but their use in non model organisms is hampered by the scarcity of cost‐effective approaches to uncover genome‐wide variation in a comprehensive subset of individuals. The screening of one or only a few individuals induces ascertainment bias. To discover SNPs for a population genomic study of the Pyrenean rocket (Sisymbrium austriacum subsp. chrysanthum), we undertook a pooled RAD‐PE (Restriction site Associated DNA Paired‐End sequencing) approach. RAD tags were generated from the PstI‐digested pooled genomic DNA of 12 individuals sampled across the species distribution range and paired‐end sequenced using Illumina technology to produce ~24.5 Mb of sequences, covering ~7% of the specie's genome. Sequences were assembled into ~76 000 contigs with a mean length of 323 bp (N₅₀ = 357 bp, sequencing depth = 24x). In all, >15 000 SNPs were called, of which 47% were annotated in putative genic regions based on homology with the Arabidopsis thaliana genome. Gene ontology (GO) slim categorization demonstrated that the identified SNPs covered extant genic variation well. The validation of 300 SNPs on a larger set of individuals using a KASPar assay underpinned the utility of pooled RAD‐PE as an inexpensive genome‐wide SNP discovery technique (success rate: 87%). In addition to SNPs, we discovered >600 putative SSR markers.     

Comparison of SNPs and microsatellites for assessing the genetic structure of chicken populations

Many studies in human genetics compare informativeness of single‐nucleotide polymorphisms (SNPs) and microsatellites (single sequence repeats; SSR) in genome scans, but it is difficult to transfer the results directly to livestock because of different population structures. The aim of this study was to determine the number of SNPs needed to obtain the same differentiation power as with a given standard set of microsatellites. Eight chicken breeds were genotyped for 29 SSRs and 9216 SNPs. After filtering, only 2931 SNPs remained. The differentiation power was evaluated using two methods: partitioning of the Euclidean distance matrix based on a principal component analysis (PCA) and a Bayesian model‐based clustering approach. Generally, with PCA‐based partitioning, 70 SNPs provide a comparable resolution to 29 SSRs. In model‐based clustering, the similarity coefficient showed significantly higher values between repeated runs for SNPs compared to SSRs. For the membership coefficients, reflecting the proportion to which a fraction segment of the genome belongs to the ith cluster, the highest values were obtained for 29 SSRs and 100 SNPs respectively. With a low number of loci (29 SSRs or ≤100 SNPs), neither marker types could detect the admixture in the Gödöllö Nhx population. Using more than 250 SNPs allowed a more detailed insight into the genetic architecture. Thus, the admixed population could be detected. It is concluded that breed differentiation studies will substantially gain power even with moderate numbers of SNPs.     

Oceanographic variation influences spatial genomic structure in the sea scallop,Placopecten magellanicus

Environmental factors can influence diversity and population structure in marine species and accurate understanding of this influence can both improve fisheries management and help predict responses to environmental change. We used 7163 SNPs derived from restriction site‐associated DNA sequencing genotyped in 245 individuals of the economically important sea scallop, Placopecten magellanicus, to evaluate the correlations between oceanographic variation and a previously identified latitudinal genomic cline. Sea scallops span a broad latitudinal area (>10 degrees), and we hypothesized that climatic variation significantly drives clinal trends in allele frequency. Using a large environmental dataset, including temperature, salinity, chlorophyll a, and nutrient concentrations, we identified a suite of SNPs (285–621, depending on analysis and environmental dataset) potentially under selection through correlations with environmental variation. Principal components analysis of different outlier SNPs and environmental datasets revealed similar northern and southern clusters, with significant associations between the first axes of each (R²_adj = .66–.79). Multivariate redundancy analysis of outlier SNPs and the environmental principal components indicated that environmental factors explained more than 32% of the variance. Similarly, multiple linear regressions and random‐forest analysis identified winter average and minimum ocean temperatures as significant parameters in the link between genetic and environmental variation. This work indicates that oceanographic variation is associated with the observed genomic cline in this species and that seasonal periods of extreme cold may restrict gene flow along a latitudinal gradient in this marine benthic bivalve. Incorporating this finding into management may improve accuracy of management strategies and future predictions.     

Development of genomic resources for the tick Ixodes ricinus: isolation and characterization of single nucleotide polymorphisms

Assessing the genetic variability of the tick Ixodes ricinus—an important vector of pathogens in Europe—is an essential step for setting up antitick control methods. Here, we report the first identification of a set of SNPs isolated from the genome of I. ricinus, by applying a reduction in genomic complexity, pyrosequencing and new bioinformatics tools. Almost 1.4 million of reads (average length: 528 nt) were generated with a full Roche 454 GS FLX run on two reduced representation libraries of I. ricinus. A newly developed bioinformatics tool (DiscoSnp), which isolates SNPs without requiring any reference genome, was used to obtain 321 088 putative SNPs. Stringent selection criteria were applied in a bioinformatics pipeline to select 1768 SNPs for the development of specific primers. Among 384 randomly SNPs tested by Fluidigm genotyping technology on 464 individuals ticks, 368 SNPs loci (96%) exhibited the presence of the two expected alleles. Hardy–Weinberg equilibrium tests conducted on six natural populations of ticks have shown that from 26 to 46 of the 384 loci exhibited significant heterozygote deficiency.     

Population structure and genetic diversity in a commercial maize breeding program assessed with SSR and SNP markers

Information about the genetic diversity and population structure in elite breeding material is of fundamental importance for the improvement of crops. The objectives of our study were to (a) examine the population structure and the genetic diversity in elite maize germplasm based on simple sequence repeat (SSR) markers, (b) compare these results with those obtained from single nucleotide polymorphism (SNP) markers, and (c) compare the coancestry coefficient calculated from pedigree records with genetic distance estimates calculated from SSR and SNP markers. Our study was based on 1,537 elite maize inbred lines genotyped with 359 SSR and 8,244 SNP markers. The average number of alleles per locus, of group specific alleles, and the gene diversity (D) were higher for SSRs than for SNPs. Modified Roger’s distance (MRD) estimates and membership probabilities of the STRUCTURE matrices were higher for SSR than for SNP markers but the germplasm organization in four heterotic pools was consistent with STRUCTURE results based on SSRs and SNPs. MRD estimates calculated for the two marker systems were highly correlated (0.87). Our results suggested that the same conclusions regarding the structure and the diversity of heterotic pools could be drawn from both markers types. Furthermore, although our results suggested that the ratio of the number of SSRs and SNPs required to obtain MRD or D estimates with similar precision is not constant across the various precision levels, we propose that between 7 and 11 times more SNPs than SSRs should be used for analyzing population structure and genetic diversity.     

Using SNP markers to dissect linkage disequilibrium at a major quantitative trait locus for resistance to the potato cyst nematode Globodera pallida on potato chromosome V

The damage caused by the parasitic root cyst nematode Globodera pallida is a major yield-limiting factor in potato cultivation . Breeding for resistance is facilitated by the PCR-based marker ‘HC’, which is diagnostic for an allele conferring high resistance against G. pallida pathotype Pa2/3 that has been introgressed from the wild potato species Solanum vernei into the Solanum tuberosum tetraploid breeding pool. The major quantitative trait locus (QTL) controlling this nematode resistance maps on potato chromosome V in a hot spot for resistance to various pathogens including nematodes and the oomycete Phytophthora infestans. An unstructured sample of 79 tetraploid, highly heterozygous varieties and breeding clones was selected based on presence (41 genotypes) or absence (38 genotypes) of the HC marker. Testing the clones for resistance to G. pallida confirmed the diagnostic power of the HC marker. The 79 individuals were genotyped for 100 single nucleotide polymorphisms (SNPs) at 10 loci distributed over 38 cM on chromosome V. Forty-five SNPs at six loci spanning 2 cM in the interval between markers GP21-GP179 were associated with resistance to G. pallida. Based on linkage disequilibrium (LD) between SNP markers, six LD groups comprising between 2 and 18 SNPs were identified. The LD groups indicated the existence of multiple alleles at a single resistance locus or at several, physically linked resistance loci. LD group C comprising 18 SNPs corresponded to the ‘HC’ marker. LD group E included 16 SNPs and showed an association peak, which positioned one nematode resistance locus physically close to the R1 gene family. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Tuber Blight Development in Potato Cultivars in Response to Different Genotypes of Phytophthora infestans

Migrations or introduction of new genotypes of Phytophthora infestans to a specific region imposes a different perspective for potato production. During 2009–2010, a late blight epidemic affected the Northeastern United States, which quickly spread through several states. The epidemic was characterized by the appearance of a new genotype of P. infestans designated US‐22, which was isolated from tomato and potato. Potato tubers are an essential component of late blight epidemics where the pathogen cannot overwinter on Solanaceous plants. Six potato cultivars were inoculated with 12 isolates of P. infestans (five different genotypes), including isolates of the genotype US‐22. Tuber blight development was characterized in terms of tissue darkening expressed as area under the disease progress curve values and lenticel infection. The responses indicated that US‐8 was more aggressive than US‐22, but US‐22 isolates obtained from potato were more aggressive on potato than those acquired from tomato. Tuber periderm responses to infection were limited, yet US‐8 isolates infected the periderm more often than US‐22 isolates. There were significant differences among the cultivars tested but cv. Jacqueline Lee was the most resistant overall. Although isolates of P. infestans genotype US‐22 were less aggressive in comparison with US‐8 isolates, US‐22 isolates still infected potato tubers and were as aggressive us US‐8 isolates on some cultivars. Management of late blight caused by isolates of US‐22 through host resistance may be feasible but imposes a different set of criteria for consideration from those that US‐8 imposed.     

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome‐wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly‐differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs’ allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16–0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species – pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.     

Genome scans on experimentally evolved populations reveal candidate regions for adaptation to plant resistance in the potato cyst nematode Globodera pallida

Improving resistance durability involves to be able to predict the adaptation speed of pathogen populations. Identifying the genetic bases of pathogen adaptation to plant resistances is a useful step to better understand and anticipate this phenomenon. Globodera pallida is a major pest of potato crop for which a resistance QTL, GpaV_vrn, has been identified in Solanum vernei. However, its durability is threatened as G. pallida populations are able to adapt to the resistance in few generations. The aim of this study was to investigate the genomic regions involved in the resistance breakdown by coupling experimental evolution and high‐density genome scan. We performed a whole‐genome resequencing of pools of individuals (Pool‐Seq) belonging to G. pallida lineages derived from two independent populations having experimentally evolved on susceptible and resistant potato cultivars. About 1.6 million SNPs were used to perform the genome scan using a recent model testing for adaptive differentiation and association to population‐specific covariables. We identified 275 outliers and 31 of them, which also showed a significant reduction in diversity in adapted lineages, were investigated for their genic environment. Some candidate genomic regions contained genes putatively encoding effectors and were enriched in SPRYSECs, known in cyst nematodes to be involved in pathogenicity and in (a)virulence. Validated candidate SNPs will provide a useful molecular tool to follow frequencies of virulence alleles in natural G. pallida populations and define efficient strategies of use of potato resistances maximizing their durability.     

Demographic consequences of greater clonal than sexual reproduction in Dicentra canadensis

Clonality is a widespread life history trait in flowering plants that may be essential for population persistence, especially in environments where sexual reproduction is unpredictable. Frequent clonal reproduction, however, could hinder sexual reproduction by spatially aggregating ramets that compete with seedlings and reduce inter‐genet pollination. Nevertheless, the role of clonality in relation to variable sexual reproduction in population dynamics is often overlooked. We combined population matrix models and pollination experiments to compare the demographic contributions of clonal and sexual reproduction in three Dicentra canadensis populations, one in a well‐forested landscape and two in isolated forest remnants. We constructed stage‐based transition matrices from 3 years of census data to evaluate annual population growth rates, λ. We used loop analysis to evaluate the relative contribution of different reproductive pathways to λ. Despite strong temporal and spatial variation in seed set, populations generally showed stable growth rates. Although we detected some pollen limitation of seed set, manipulative pollination treatments did not affect population growth rates. Clonal reproduction contributed significantly more than sexual reproduction to population growth in the forest remnants. Only at the well‐forested site did sexual reproduction contribute as much as clonal reproduction to population growth. Flowering plants were more likely to transition to a smaller size class with reduced reproductive potential in the following year than similarly sized nonflowering plants, suggesting energy trade‐offs between sexual and clonal reproduction at the individual level. Seed production had negligible effects on growth and tuber production of individual plants. Our results demonstrate that clonal reproduction is vital for population persistence in a system where sexual reproduction is unpredictable. The bias toward clonality may be driven by low fitness returns for resource investment in sexual reproduction at the individual level. However, chronic failure in sexual reproduction may exacerbate the imbalance between sexual and clonal reproduction and eventually lead to irreversible loss of sex in the population.     

Microsatellite markers reveal two admixed genetic groups and an ongoing displacement within the French population of the invasive plant pathogen Phytophthora infestans

Potato late blight is an example of a re‐emerging disease of plants. Phytophthora infestans was first introduced into Europe during the 19th century, where it caused the Irish potato famine. During the 20th century several additional introduction events have been suspected, especially in the mid‐70s due to the import of large quantities of potato needed after the shortage caused by drought in 1976. Here, we investigate the genetic population structure of Phytophthora infestans, at the first stages of a recent invasion process in France. A total of 220 isolates was collected from 20 commercial fields of the potato susceptible cultivar Bintje, during two consecutive years (2004 and 2005). Clustering analyses based on eight recently developed microsatellite markers reveal that French P. infestans populations are made of two differentiated genetic clusters of isolates (F_ST = 0.19). This result suggests multiple introductions of P. infestans into France, either through the introduction of a composite population of isolates or through the successive introduction of isolates having differentiated genetic backgrounds. Both clusters identified have a strong clonal structure and are similar regarding genetic diversity and mating type composition. The maintenance of differentiation between the two genetic clusters should result from the low or non‐existent contribution of sexual reproduction in French P. infestans populations.     

Causes and consequences of large clonal assemblies in a poplar hybrid zone

Asexual reproduction is a common and fundamental mode of reproduction in plants. Although persistence in adverse conditions underlies most known cases of clonal dominance, proximal genetic drivers remain unclear, in particular for populations dominated by a few large clones. In this study, we studied a clonal population of the riparian tree Populus alba in the Douro river basin (northwestern Iberian Peninsula) where it hybridizes with Populus tremula, a species that grows in highly contrasted ecological conditions. We used 73 nuclear microsatellites to test whether genomic background (species ancestry) is a relevant cause of clonal success, and to assess the evolutionary consequences of clonal dominance by a few genets. Additional genotyping‐by‐sequencing data were produced to estimate the age of the largest clones. We found that a few ancient (over a few thousand years old) and widespread genets dominate the population, both in terms of clone size and number of sexual offspring produced. Interestingly, large clones possessed two genomic regions introgressed from P. tremula, which may have favoured their spread under stressful environmental conditions. At the population level, the spread of large genets was accompanied by an overall ancient (>0.1 Myr) but soft decline of effective population size. Despite this decrease, and the high clonality and dominance of sexual reproduction by large clones, the Douro hybrid zone still displays considerable genetic diversity and low inbreeding. This suggests that even in extreme cases as in the Douro, asexual and sexual dominance of a few large, geographically extended individuals does not threaten population survival.     

Mapping and Validation of the Major Sex-Determining Region in Nile Tilapia (Oreochromis niloticus L.) Using RAD Sequencing

Sex in Oreochromis niloticus (Nile tilapia) is principally determined by an XX/XY locus but other genetic and environmental factors also influence sex ratio. Restriction Associated DNA (RAD) sequencing was used in two families derived from crossing XY males with females from an isogenic clonal line, in order to identify Single Nucleotide Polymorphisms (SNPs) and map the sex-determining region(s). We constructed a linkage map with 3,802 SNPs, which corresponded to 3,280 informative markers, and identified a major sex-determining region on linkage group 1, explaining nearly 96% of the phenotypic variance. This sex-determining region was mapped in a 2 cM interval, corresponding to approximately 1.2 Mb in the O. niloticus draft genome. In order to validate this, a diverse family (4 families; 96 individuals in total) and population (40 broodstock individuals) test panel were genotyped for five of the SNPs showing the highest association with phenotypic sex. From the expanded data set, SNPs Oni23063 and Oni28137 showed the highest association, which persisted both in the case of family and population data. Across the entire dataset all females were found to be homozygous for these two SNPs. Males were heterozygous, with the exception of five individuals in the population and two in the family dataset. These fish possessed the homozygous genotype expected of females. Progeny sex ratios (over 95% females) from two of the males with the “female” genotype indicated that they were neomales (XX males). Sex reversal induced by elevated temperature during sexual differentiation also resulted in phenotypic males with the “female” genotype. This study narrows down the region containing the main sex-determining locus, and provides genetic markers tightly linked to this locus, with an association that persisted across the population. These markers will be of use in refining the production of genetically male O. niloticus for aquaculture.