All-male production by marker-assisted selection for sex determining loci of admixed Oreochromis niloticus and Oreochromis aureus stocks.

Crossing Oreochromis niloticus (On) females with Oreochromis aureus (Oa) males results in all-male progeny that are essential for effective tilapia farming. However, a reproductive barrier between these species limits mating and mass-fry production. One approach to overcoming this barrier is to select parental stocks of mixed genetic backgrounds, which allow interspecific reproductive recognition, while closely maintaining the genetic profiles for sex-determination (SD) of the respective purebred species. Here, we test this approach in a data set of 160 On × Oa spawns of 109 male and 100 female parents randomly collected from admixed stocks, and genotyped for microsatellite markers representing the known SD loci on linkage groups (LGs) 1, 3, and 23. Following crossbreeding, the most significant paternal effects on male proportions in progeny were found for LG1-BYL018 (P < 2 × 10-32 ) and for LG3-UNH168 × LG23-UNH898 interaction (P < 1 × 10-17 ; R2  = 0.98). Furthermore, a maternal effect for LG3-UNH168 (P < 9 × 10-7 ) was associated with low female proportions in progeny (<7%), indicating a non-Mendelian effect on SD. Eighty-four males (77%) and 30 females (30%) were selected as parents, based on their genetic profiles for the SD loci that were associated with male production. Of these, 51 of 53 crosses produced all-male progeny, while two crosses had low female proportions in their progeny (<4%). This suggests that selection could be improved using the causative sequence variation underlying SD on LG3, since the large non-recombining block of the SD region in purebred Oa readily breaks down in hybrids. Nevertheless, marker-assisted selection for sex determining loci of admixed parental stocks may be used for all-male production.

A duplication of the Anti-Müllerian hormone gene is associated with genetic sex determination of different Oreochromis niloticus strains.

Sex determination (SD) mechanisms are ancient and conserved, yet much diversity is exhibited in primary sex-determining signals that trigger male or female development. In O. niloticus, SD is associated with a male-specific locus on linkage group (LG) 23 which harbors the Y-linked Anti-Müllerian hormone (amh) gene, and a truncated duplication, denoted amhΔy. We have evaluated the possible role of identified indels and SNPs in the amh gene on SD, based on conservation in different O. niloticus strains. A fluorescent assay for the detection of a 5 bp insertion in amhΔy exon VI, efficiently discriminated between XX, XY, and YY genotypes. Concordance rate between amhΔy and sex varied in six Oreochromis strains, from 100% (Ghana) through 90% (Swansea) to 85% (Thai-Chitralada). The association of amhΔy with sex was found to be conserved in all tested O. niloticus strains, and thus supports its key role in SD. However, the previously identified missense SNP (C/T) in amh exon II was found only in the Swansea strain, thus excluding its candidacy for the causal variation of SD across all strains. Effects of markers on LGs 1, 3, and 23 (amhΔy) fully explained sex distribution in one Thai-Chitralada family (R2 = 1.0), whereas in another family only the major effect of LG23 (amhΔy) was significant (R2 = 0.37). Thus, amhΔy on LG23 is associated with genetic SD, either as a single causal gene in different O. niloticus strains, or in combination with segregating genes on LGs 1 and 3 in the Thai-Chitralada hybrid strain.

Identification of the sex-determining region in flathead grey mullet (Mugil cephalus).

Elucidation of the sex-determination mechanism in flathead grey mullet (Mugil cephalus) is required to exploit its economic potential by production of genetically determined monosex populations and application of hormonal treatment to parents rather than to the marketed progeny. Our objective was to construct a first-generation linkage map of the M. cephalus in order to identify the sex-determining region and sex-determination system. Deep-sequencing data of a single male was assembled and aligned to the genome of Nile tilapia (Oreochromis niloticus). A total 245 M. cephalus microsatellite markers were designed, spanning the syntenic tilapia genome assembly at intervals of 10 Mb. In the mapping family of full-sib progeny, 156 segregating markers were used to construct a first-generation linkage map of 24 linkage groups (LGs), corresponding to the number of chromosomes. The linkage map spanned approximately 1200 cM with an average inter-marker distance of 10.6 cM. Markers segregating on LG9 in two independent mapping families showed nearly complete concordance with gender (R2  = 0.95). The sex determining locus was fine mapped within an interval of 8.6 cM on LG9. The sex of offspring was determined only by the alleles transmitted from the father, thus indicating an XY sex-determination system.

Methods for genomic evaluation of a relatively small genotyped dairy population and effect of genotyped cow information in multiparity analyses.

Methods for genomic prediction were evaluated for an Israeli Holstein dairy population of 713,686 cows and 1,305 progeny-tested bulls with genotypes. Inclusion of genotypes of 343 elite cows in an evaluation method that considers pedigree, phenotypes, and genotypes simultaneously was also evaluated. Two data sets were available: a complete data set with production records from 1985 through 2011, and a reduced data set with records after 2006 deleted. For each production trait, a multitrait animal model was used to compute traditional genetic evaluations for parities 1 through 3 as separate traits. Evaluations were calculated for the reduced and complete data sets. The evaluations from the reduced data set were used to calculate parent average for validation bulls, which was the benchmark for comparing gain in predictive ability from genomics. Genomic predictions for bulls in 2006 were calculated using a Bayesian regression method (BayesC), genomic BLUP (GBLUP), single-step GBLUP (ssGBLUP), and weighted ssGBLUP (WssGBLUP). Predictions using BayesC and GBLUP were calculated either with or without an index that included parent average. Genomic predictions that included elite cow genotypes were calculated using ssGBLUP and WssGBLUP. Predictive ability was assessed by coefficients of determination (R(2)) and regressions of predictions of 135 validation bulls with no daughters in 2006 on deregressed evaluations of those bulls in 2011. A reduction in R(2) and regression coefficients was observed from parities 1 through 3. Fat and protein yields had the lowest R(2) for all the methods. On average, R(2) was lowest for parent averages, followed by GBLUP, BayesC, ssGBLUP, and WssGBLUP. For some traits, R(2) for direct genomic values from BayesC and GBLUP were lower than those for parent averages. Genomic estimated breeding values using ssGBLUP were the least biased, and this method appears to be a suitable tool for genomic evaluation of a small genotyped population, as it automatically accounts for parental index, allows for inclusion of female genomic information without preadjustments in evaluations, and uses the same model as in traditional evaluations. Weighted ssGBLUP has the potential for higher evaluation accuracy.

Predictive ability of selected subsets of single nucleotide polymorphisms (SNPs) in a moderately sized dairy cattle population.

Several studies have shown that computation of genomic estimated breeding values (GEBV) with accuracies significantly greater than parent average (PA) estimated breeding values (EBVs) requires genotyping of at least several thousand progeny-tested bulls. For all published analyses, GEBV computed from the selected samples of markers have lower or equal accuracy than GEBV derived on the basis of all valid single nucleotide polymorphisms (SNPs). In the current study, we report on four new methods for selection of markers. Milk, fat, protein, somatic cell score, fertility, persistency, herd life and the Israeli selection index were analyzed. The 972 Israeli Holstein bulls genotyped with EBV for milk production traits computed from daughter records in 2012 were assigned into a training set of 844 bulls with progeny test EBV in 2008, and a validation set of 128 young bulls. Numbers of bulls in the two sets varied slightly among the nonproduction traits. In EFF12, SNPs were first selected for each trait based on the effects of each marker on the bulls' 2012 EBV corrected for effective relationships, as determined by the SNP matrix. EFF08 was the same as EFF12, except that the SNPs were selected on the basis of the 2008 EBV. In DIFmax, the SNPs with the greatest differences in allelic frequency between the bulls in the training and validation sets were selected, whereas in DIFmin the SNPs with the smallest differences were selected. For all methods, the numbers of SNPs retained varied over the range of 300 to 6000. For each trait, except fertility, an optimum number of markers between 800 and 5000 was obtained for EFF12, based on the correlation between the GEBV and current EBV of the validation bulls. For all traits, the difference between the correlation of GEBV and current EBV and the correlation of the PA and current EBV was >0.25. EFF08 was inferior to EFF12, and was generally no better than PA EBV. DIFmax always outperformed DIFmin and generally outperformed EFF08 and PA. Furthermore, GEBV based on DIFmax were generally less biased than PA. It is likely that other methods of SNP selection could improve upon these results.

Maternity validation using sire-only BovineSNP50 BeadChip data.

Based on pairwise identity-by-state (IBS) distances and whole-genome SNP data, kinship was investigated in the Israeli Holstein population. A total of 789 bulls, including most of the artificial insemination sires in service since 1987, were genotyped by the BovineSNP50 BeadChip. This sample included up to five generations. For each bull-by-bull combination, three states are possible for each marker: no match, a single match and both alleles match. Summing over all markers, the 932 598 IBS scores (three match frequencies*310 866 bull-by-bull combinations) were visualized using three-dimensional coordinates that corresponded to the frequencies of the three possible states. Results were reduced to two dimensions using the transformations x' = 0.7071(1 + freq1-freq2) and y' = 1.2247freq0. Bull-by-bull pairs were grouped according to their level of kinship, and canonical scores were calculated using discriminant analysis and the x' and y' features. Of the 474 pairs of recorded maternal grandsire-grandson with both individuals genotyped, the probability for 28 pairs to belong to this level of kinship was low (P < 0.05), suggesting an error rate of around 3% per generation in pedigree determination.

A deleterious effect associated with UNH159 is attenuated in twin embryos of an inbred line of blue tilapia Oreochromis aureus.

Offspring of a highly inbred gynogenetic line of Oreochromis aureus displayed 12-fold increase in twinning rate compared to the outbred population. Asymmetric conjoined twins, which consist of a normal embryo attached to a malformed-atrophic twin, were frequently encountered in both gynogenetic (90·7%) and outbred (38·2%) embryos. The monozygotic origin of these twins was determined using five microsatellite markers. Progeny of heterozygous parents for the microsatellite UNH159 were separated into sub-sets of twins and normal full-sibs. Consistent with previous reports, the normal embryo sub-set exhibited elimination of both types of homozygotes for the UNH159 genetic marker at 2-8 days after fertilization. Unexpectedly, this elimination was less frequent in twins. The UNH159 marker as well as RNA-binding motif protein, X-linked (rbmx), SRY-box containing gene 3 (sox3) and alpha-thalassemia/mental retardation syndrome X-linked (atrx) genes were mapped to linkage group 2. These gene orthologues are all located on the mammalian X chromosome and atrx is necessary for the X-chromosome inactivation.

Signatures of contemporary selection in the Israeli Holstein dairy cattle.

Strong selection in the Israeli Holstein dairy cattle population over the last three decades should have left clear signatures of selection. Two experimental approaches were applied to detect evidence of contemporary selection based on the 54K BeadChip genotypes of ~1000 Israeli Holstein bulls: (i) the long-range haplotype test, which searches for structural evidence resulting from selective sweep, and (ii) direct analysis of the changes in haplotypes frequencies over time combined with linkage disequilibrium blocks haplotype-based association analysis. Ten traits were analyzed: the PD07 Israeli selection index, milk, milk fat, % fat, milk protein, % protein, somatic cell score, female fertility, milk production persistency and herd life. The long-range haplotype test detected ~15% of the 3288 haplotypes that showed significant positive frequency trends (P < 0.05) and was significantly correlated with the substitution effects of the haplotypes and the selection intensities for the different traits. Thirty signatures of recent selection, which correspond to both approaches and affect the Israeli PD07 selection index, were identified on 17 of the 29 autosomes. The second experimental approach also was used to estimate the selection intensity of the different traits. The correlation between the selection intensities for the traits analyzed, derived from changes in haplotype frequencies in the population of bulls, and those derived from trait-based analysis of the cow population was 0.93 over all traits. Thus, the changes in haplotypes frequencies in the bulls' population accurately estimate genetic trends in the general cow population and can be used to detect signatures of recent selection.

A simple method for genomic selection of moderately sized dairy cattle populations.

An efficient algorithm for genomic selection of moderately sized populations based on single nucleotide polymorphism chip technology is described. A total of 995 Israeli Holstein bulls with genetic evaluations based on daughter records were genotyped for either the BovineSNP50 BeadChip or the BovineSNP50 v2 BeadChip. Milk, fat, protein, somatic cell score, female fertility, milk production persistency and herd-life were analyzed. The 400 markers with the greatest effects on each trait were first selected based on individual analysis of each marker with the genetic evaluations of the bulls as the dependent variable. The effects of all 400 markers were estimated jointly using a 'cow model,' estimated from the data truncated to exclude lactations with freshening dates after September 2006. Genotype probabilities for each locus were computed for all animals with missing genotypes. In Method I, genetic evaluations were computed by analysis of the truncated data set with the sum of the marker effects subtracted from each record. Genomic estimated breeding values for the young bulls with genotypes, but without daughter records, were then computed as their parent averages combined with the sum of each animal's marker effects. Method II genomic breeding values were computed based on regressions of estimated breeding values of bulls with daughter record on their parent averages, sum of marker effects and birth year. Method II correlations of the current breeding values of young bulls without daughter records in the truncated data set were higher than the correlations of the current breeding values with the parent averages for fat and protein production, persistency and herd-life. Bias of evaluations, estimated as a difference between the mean of current breeding values of the young bulls and their genomic evaluations, was reduced for milk production traits, persistency and herd-life. Bias for milk production traits was slightly negative, as opposed to the positive bias of parent averages. Correlations of Method II with the means of daughter records adjusted for fixed effects were higher than parent averages for fat, protein, fertility, persistency and herd-life. Reducing the number of markers included in the analysis from 400 to 300 did not reduce correlations of genomic breeding values for protein with current breeding values, but did slightly reduce correlations with means of daughter records. Method II has the advantages as compared with the method of VanRaden in that genotypes of cows can be readily incorporated into the Method II analysis, and it is more effective for moderately sized populations.

Linkage and Physical Mapping of Sex Region on LG23 of Nile Tilapia (Oreochromis niloticus).

Evidence supports that sex determination (SD) in tilapia is controlled by major genetic factors that may interact with minor genetic as well as environmental factors, thus implying that SD should be analyzed as a quantitative trait. Quantitative trait loci (QTL) for SD in Oreochromis niloticus were previously detected on linkage groups (LG) 1 and 23. Twenty-one short single repeats (SSR) of >12 TGs and one single nucleotide polymorphism were identified using the unpublished tilapia genome sequence on LG23. All markers showed two segregating alleles in a mapping family that was obtained by a cross between O. niloticus male (XY) and sex-reversed female (ΔXY) yielding 29 females (XX) and 61 males (XY and YY). Interval mapping analysis mapped the QTL peak between SSR markers ARO172 and ARO177 with a maximum F value of 78.7 (P < 7.6 × 10(-14)). Twelve adjacent markers found in this region were homozygous in females and either homozygous for the alternative allele or heterozygous in males. This segment was defined as the sex region (SR). The SR encompasses 1.5 Mbp on a single tilapia scaffold (no. 101) harboring 51 annotated genes. Among 10 candidate genes for SD that were tested for gene expression, anti-Müllerian hormone (Amh), which is located in the center of the SR, showed the highest overexpression in male vs. female embryos at 3 to 7 days postfertilization.

Differential expression of ruminant ZNF496 variants: Association with quantitative trait locus affecting bovine milk concentration and fertility.

A single nucleotide polymorphism in the intergenic region upstream of the ZNF496 gene on Bos taurus chromosome 7 displayed significant population-wide linkage disequilibrium with milk protein percentage in the Israeli Holstein population. The frequency of the allele associated with increased protein concentration was 10%. This single nucleotide polymorphism was located in the promoter region from which a 10-exon transcript of the bovine and the ovine ZNF496 genes are transcribed. The gene architecture was similar to the mouse ortholog Zkscan17. A 5-exon murine antisense transcript was complementary to the 5' untranslated Zkscan17 region that included a sequence domain conserved between mouse and ruminants, suggesting a regulatory function. In the bovine ZNF496 chromosomal region, segregation of a quantitative trait locus (QTL) for milk protein percentage was confirmed in a daughter design sire family. Concordance was not obtained between QTL status of bulls and any of the polymorphisms in the functional elements of ZNF496. This excludes these variations as the causative polymorphism under the assumption of no epigenetic effect for this locus. However, ZNF496 variants were differentially expressed in bovine ovaries, and only the paternal variant was expressed in liver and kidney in a sheep family with polymorphic ZNF496 sequence. Thus, the search for the mutation underlying the minor QTL allele, which is a top economically favorable allele in Israeli Holstein cattle, may be complicated by the presence of an imprinting center in this QTL confidence interval.

Fine-mapping of a locus on linkage group 23 for sex determination in Nile tilapia (Oreochromis niloticus).

Genetic markers in tilapia species associated with loci affecting sex determination (SD), sex-specific mortality or both were mapped to linkage groups (LG) 1, 2, 3, 6 and 23. The objective of this study was to use these markers to fine-map the locus with the greatest effect on SD in Oreochromis niloticus. Our parental stock, full-sibs of Nile tilapia (Swansea origin), were divided into three groups: (i) untreated, (ii) feminized by diethylstilbestrol and (iii) masculinized by 17α-methyltestosterone. We analysed the first group for association of microsatellite markers representing these five LGs. The strongest association with gender was found on LG23 for marker UNH898 (χ(2) ; P=8.6×10(-5) ). Allele 276 was found almost exclusively in males, and we hypothesized that this allele is a male-associated allele (MAA). Sex-reversed individuals were used for mating experiments with and without the segregating MAA. Mating of individuals lacking the MAA resulted in all-female progeny. Mating of two heterozygotes for MAA gave rise to 81 males and 30 females. Analysis of association between gender and genotypes identified the MAA in 98.6% of males as opposed to 8.0% of females (χ(2) ; P=2.5×10(-18) ). Eight markers that flank UNH898 were genotyped to map the locus on LG23 within a confidence interval of 16-21 cM. Mating of homozygous individuals for MAA is underway for production of all-male populations.

Development of a 25-plex SNP assay for traceability in cattle.

Single nucleotide polymorphisms (SNPs) are amenable to automation and therefore have become the marker of choice for DNA profiling. SNaPshot, a primer extension-based method, was used to multiplex 25 SNPs that have been previously validated as useful for identity control. Detection of extended products was based on four different fluorochromes and extension primers with oligonucleotide tails of differing lengths, thus controlling the concise length of the entire chromatogram to 81 bases. Allele frequencies for Holstein, Simmental, Limousin, Angus, Charolais and Tux Cattle were estimated and significant positive Pearson-correlation coefficients were obtained among the analysed breeds. The probability that two randomly unrelated individuals would share identical genotypes for all 25 loci varied from 10(-8) to 10(-10) for these breeds. For parentage control, the exclusion power was found to be 99.9% when the genotypes of both putative parents are known. A traceability test of duplicated samples indicated a high genotyping precision of >0.998. This was further corroborated by analysis of 60 cases of parent-sib pairs and trio families. The 25-plex SNaPshot assay is adapted for low- and high-throughput capacity and thus presents an alternative for DNA-based traceability in the major commercial cattle breeds.

Copy number variation of lipocalin family genes for male-specific proteins in tilapia and its association with gender.

Lipocalins are involved in the binding of small molecules like sex steroids. We show here that the previously reported tilapia male-specific protein (MSP) is a lipocalin encoded by a variety of paralogous and homologous genes in different tilapia species. Exon-intron boundaries of MSP genes were typical of the six-exon genomic structure of lipocalins, and the transcripts were capable of encoding 200 amino-acid polypeptides that consisted of a putative signal peptide and a lipocalin domain. Cysteine residues are conserved in positions analogous to those forming the three disulfide bonds characteristic of the ligand pocket. The calculated molecular mass of the secreted MSP (20.4 kDa) was less than half of that observed, suggesting that it is highly glycosylated like its homologue tributyltin-binding protein. Analysis of sequence variations revealed three types of paralogs MSPA, MSPB and MSPC. Expression of both MSPA and MSPB was detected in testis. In haploid Oreochromis niloticus embryos, each of these types consisted of two closely related paralogs, and asymmetry between MSP copy numbers on the maternal (six copies) and the paternal (three copies) chromosomes was observed. Using this polymorphism we mapped MSPA and MSPC to linkage group 12 of an F(2) mapping family derived from a cross between O. niloticus and Oreochromis aureus. Females with high MSP copy number were more frequent by more than twofold than males. Gender-MSPC combinations showed significant deviation from expected Mendelian segregation (P=0.009) suggesting elimination of males with MSPC copies. We discuss different hypotheses to explain this elimination, including possibility for allelic conflict resulted by the hybridization.

Genetics of sex determination in tilapiine species.

We identified DNA markers linked to sex determining genes in six closely related species of tilapiine fishes. The mode of sex determination differed among species. In Oreochromis karongae and Tilapia mariae the sex-determining locus is on linkage group (LG) 3 and the female is heterogametic (WZ-ZZ system). In O. niloticus and T. zillii the sex-determining locus is on LG1 and the male is heterogametic (XX-XY system). A more complex pattern was observed in O. aureus and O. mossambicus, in which markers on both LG1 and LG3 were associated with sex. We found evidence for sex-linked lethal effects on LG1, as well as interactions between loci in the two linkage groups. Comparison of genetic and physical maps demonstrated a broad region of recombination suppression harboring the sex-determining locus on LG3. Sex-specific recombination suppression was found in the female heterogametic sex. Sequence analysis showed the accumulation of repetitive elements in this region. Phylogenetic analysis suggests that at least two transitions in the mode of sex determination have occurred in this clade. This variation in sex determination mechanisms among closely related species makes tilapias an excellent model system for studying the evolution of sex chromosomes in vertebrates.

Association between loci with deleterious alleles and distorted sex ratios in an inbred line of tilapia (Oreochromis aureus).

Three microsatellite markers (UNH159, UNH231, and UNH216) were examined for association with both deleterious genes and sex-ratio distortions in a full-sib family of 222 progeny from the fourth generation of a meiogynogenetic tilapia line (Oreochromis aureus). The three markers were mapped previously to different linkage groups and were shown to be associated with genes with deleterious alleles in this line. A restricted maximum likelihood model was used for analysis of major effects and their interactions on sex ratio and viability. This model was based on selective mortality of genders, ignoring effects of possible sex-determining genes. The results showed that deleterious genes linked to UNH216 and UNH231 exert higher lethality in females than in males (P < .0005 and P < .05, respectively). UNH159 was not associated directly with sex ratio distortion, but acts strongly as a modifier of sex ratio in combination with UNH216 and UNH231. Each of the three loci was found to have a significant effect on viability (P < .05) in the maximum likelihood analysis. The deleterious single-locus effects act strongly against females, while most of the epistatic interactions exert higher lethality in males. This contradiction results in a close to 1:1 sex ratio at maturity. The genetic mechanism and significance of such a balance between genders are still unknown. A detailed analysis of sex-specific lethality may be applied by screening in appropriate series of matings and fine mapping with additional markers. Our data suggest that UNH216 and UNH231 are linked to sex ratio distortion genes and that UNH159 may be linked to a modifier of these genes.