Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review.

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.

Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System.

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene's size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix-loop-helix domain that is typical of the ovary development regulator-factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.

Cross-species conservation of a transposase-linked element enables genetic sexing of commercial populations of Russian sturgeon (Acipenser gueldenstaedtii).

All-female culture of sturgeon is essential for efficient caviar production. However, Russian sturgeon (Acipenser gueldenstaedtii) does not exhibit external sexual dimorphism, and therefore, commercial farms apply gonadal endoscopy or ultrasound at the earliest age of 4-5 years to separate the sexes, with ~90% accuracy. Recently, a dominant genomic marker (AllWSEX2) has been found with association to femaleness in sturgeons. We developed a duplex PCR (dAllWSEX2) with the adjacent bmp7 gene as an internal control, to validate an effective PCR. Robust amplification of control fragments was observed for all samples of our commercial A. gueldenstaedtii stock (n = 337). The dAllWSEX2 assay was significantly associated with sex (n = 43, p < 1.6 × 10-8 ), yet four (18%) of the endoscopy-determined females were genetic males. To examine whether some females display a male genetic profile, we tested 96 egg-producing females, which were all verified as genetic females, indicating that the observed mismatches may be attributed to wrong sexing by endoscopy. Application of dAllWSEX2 on 100 7-month-old fish showed no sex-dependent differences in body weight, indicating that weighing is not an applicable tool for sorting females at a young age. Sanger sequencing of the bmp7 fragment revealed octaploidy and sex-independent variation, suggesting that the critical sex-determining region harboring AllWSEX2 is small. In keeping with a model of a single-ploidy encoding female determination, AllWSEX2 showed no variation despite being a transposase-linked repetitive element. Cross-species conservation of AllWSEX2, and absence of annotated sex-determination genes in this region suggests that, in sturgeons, the sex-determining mechanism is different from mechanisms identified in other fish.

Gene Variant of Barrier to Autointegration Factor 2 (Banf2w) Is Concordant with Female Determination in Cichlids.

Oreochromis fishes exhibit variability of sex-determination (SD) genes whose characterization contributes to understanding of the sex differentiation network, and to effective tilapia farming, which requires all-male culture. However, O. niloticus (On) amh is the only master-key regulator (MKR) of SD that has been mapped (XY/XX SD-system on LG23). In O. aureus (Oa), LG3 controls a WZ/ZZ SD-system that has recently been delimited to 9.2 Mbp, with an embedded interval rich with female-specific variation, harboring two paics genes and banf2. Developing genetic markers within this interval and using a hybrid Oa stock that demonstrates no recombination repression in LG3, we mapped the critical SD region to 235 Kbp on the orthologous On physical map (p < 1.5 × 10-26). DNA-seq assembly and peak-proportion analysis of variation based on Sanger chromatograms allowed the characterization of copy-number variation (CNV) of banf2. Oa males had three exons capable of encoding 90-amino-acid polypeptides, yet in Oa females, we found an extra copy with an 89-amino-acid polypeptide and three non-conservative amino acid substitutions, designated as banf2w. CNV analysis suggested the existence of two to five copies of banf2 in diploidic Cichlidae. Disrupting the Hardy-Weinberg equilibrium (p < 4.2 × 10-3), banf2w was concordant with female determination in Oa and in three cichlids with LG3 WZ/ZZ SD-systems (O. tanganicae, O. hornorum and Pelmatolapia mariae). Furthermore, exclusive RNA-seq expression in Oa females strengthened the candidacy of banf2w as the long-sought LG3 SD MKR. As banf genes mediate nuclear assembly, chromatin organization, gene expression and gonad development, banf2w may play a fundamental role inducing female nucleus formation that is essential for WZ/ZZ SD.

A novel c.1759T>G variant in follicle-stimulating hormone-receptor gene is concordant with male determination in the flathead grey mullet (Mugil cephalus).

Various master key regulators (MKRs) that control a binary switch of sex determination (SD) have been found in fish; these provide an excellent model for the study of vertebrate genetic SD. The SD region in flathead grey mullet has been previously mapped to a 1 Mbp region harboring 27 genes, of which one is follicle-stimulating hormone receptor (fshr). Although this gene is involved in gonad differentiation and function, it has not been considered as an MKR of SD. We systematically investigated polymorphism in mullet fshr using DNA shotgun sequences, and compared them between males and females. Capable of encoding nonconservative amino acid substitutions, c.1732G>A and c.1759T>G exhibited association with sex on a population level (N = 83; P ≤ 6.7 × 10-19). Hence, 1732 A and 1759 G represent a male-specific haplotype of the gene, designated as "fshry." Additional flanking SNPs showed a weaker degree of association with sex, delimiting the SD critical region to 143 nucleotides on exon 14. Lack of homozygotes for fshry, and the resulting divergence from Hardy-Weinberg equilibrium (N = 170; P ≤ 3.9 × 10-5), were compatible with a male heterogametic model (XY/XX). Capable of replacing a phenylalanine with valine, c.1759T>G alters a conserved position across the sixth transmembrane domain of vertebrate FSHRs. Amino acid substitutions in this position in vertebrates are frequently associated with constant receptor activation and consequently with FSH/FSHR signaling alteration; thus, indicating a potential role of fshr as an MKR of SD.

Preferential Mapping of Sex-Biased Differentially-Expressed Genes of Larvae to the Sex-Determining Region of Flathead Grey Mullet (Mugil cephalus).

Flathead gray mullet (Mugil cephalus) is a cosmopolitan mugilid species popular in fishery and aquaculture with an economic preference for all-female population. However, it displays neither sexual dimorphisms nor heteromorphic sex chromosomes. We have previously presented a microsatellite-based linkage map for this species locating a single sex determination region (SDR) on linkage group 9 (LG9) with evidence for XX/XY sex determination (SD) mechanism. In this work, we refine the critical SDR on LG9, and propose positional- and functional- candidate genes for SD. To elucidate the genetic mechanism of SD, we assembled and compared male and female genomic sequences of 19 syntenic genes within the putative SDR on mullet's LG9, based on orthology to tilapia's LG8 (tLG8) physical map. A total of 25 sequence-based markers in 12 genes were developed. For all markers, we observed association with sex in at least one of the two analyzed M. cephalus full-sib families, but not in the wild-type population. Recombination events were inferred within families thus setting the SDR boundaries to a region orthologous to ∼0.9 Mbp with 27 genes on tLG8. As the sexual phenotype is evident only in adults, larvae were assigned into two putative sex-groups according to their paternal haplotypes, following a model of XY/XX SD-system. A total of 107 sex-biased differentially expressed genes in larvae were observed, of which 51 were mapped to tLG8 (48% enrichment), as compared to 5% in random control. Furthermore, 23 of the 107 genes displayed sex-specific expression; and 22 of these genes were positioned to tLG8, indicating 96% enrichment. Of the 27 SDR genes, BCCIP, DHX32A, DOCK1, and FSHR (GTH-RI) are suggested as positional and functional gene candidates for SD.

Mapping of the Sex Determining Region on Linkage Group 12 of Guppy (Poecilia reticulata).

Poecilia reticulata is one of the most popular ornamental fish species with a higher demand for males due to their large colorful fins. The objectives of this study were mapping of the sex-determining (SD) region on linkage group 12 of guppy, and identification of a sex specific marker. We generated eight polymorphic microsatellite markers distributed along the distal 5.4 Mbp sequence of the previously identified SD region on linkage group (LG) 12. The markers were tested for association with sex in 156 individuals of the Red Blonde and Flame strains, and 126 progeny of four full-sibs Red Blonde families. A male-specific allele was found for microsatellite gu1066 at position of 25.3 Mbp on LG12 for both strains, and gu832 at position 24.4 Mbp for the Flame strain. Thus, a region of 0.9 Mbp between these markers, harboring 27 annotated genes, was selected for analysis. Based on association of copy number variation and sex determination we mapped a duplicated region between LGs 9 and 12, of 1.26 Mbp, containing 59 genes on LG12. The common interval between the segment bounded by gu1066 and gu832, and the duplicated region of 0.43 Mbp on LG12 harbors 11 genes of which 6 have multiple copies (54%). Growth arrest and DNA damage inducible gamma-like (GADD45G-like) is a plausible positional and functional candidate gene for its role in male fertility. We characterized the genomic structure of the gene in guppy, and found two isoforms; but no sex-biased differences were evident in genomic sequence and gene expression.

Quantitative trait loci on LGs 9 and 14 affect the reproductive interaction between two Oreochromis species, O. niloticus and O. aureus.

Effective farming of tilapia requires all-male culture, characterized by uniformity and high growth rate. Males of O. aureus (Oa) and females of O. niloticus (On) produce all-male offspring, but there is a behavioral reproductive barrier between the two species that prevents mass production. In crosses between Oa and On broodstocks, few hybrid females are attracted to the Oa male nests (denoted responders), and if they harbor the On alleles for the sex determination (SD) sites on linkage groups (LGs) 1, 3, and 23, all-male progeny are produced. Yet, without controlling for the alleles underlying SD, the parental stocks gradually lose their capability for all-male production. Hypothesizing that marker-assisted selection for female responders would allow production of sustainable broodstocks, we applied genotyping-by-sequencing to generate 4983 informative SNPs from 13 responding and 28 non-responding females from two full-sib families. Accounting for multiple comparisons in a genome-wide association study, seven SNPs met a false discovery rate of 0.061. Lowest nominal probabilities were on LGs 9 and 14, for which microsatellite DNA markers were designed within the candidate genes PTGDSL and CASRL, respectively. By increasing the sample size to 22 responders and 47 non-responders and by genotyping additional established microsatellites, we confirmed the association of these LGs with female responsiveness. The combined effects of microsatellites GM171 and CARSL-LOC100690618 on LGs 9 and 14 explained 37% of the phenotypic variance of reproductive interaction (p < 0.0001). Based on these findings, we propose a strategy for mass production of all-male tilapia hybrids through selection for genomic loci affecting SD and female responsiveness.

Identification of male-specific amh duplication, sexually differentially expressed genes and microRNAs at early embryonic development of Nile tilapia (Oreochromis niloticus).

BACKGROUND: The probable influence of genes and the environment on sex determination in Nile tilapia suggests that it should be regarded as a complex trait. Detection of sex determination genes in tilapia has both scientific and commercial importance. The main objective was to detect genes and microRNAs that were differentially expressed by gender in early embryonic development. RESULTS: Artificial fertilization of Oreochromis niloticus XX females with either sex-reversed ΔXX males or genetically-modified YY 'supermales' resulted in all-female and all-male embryos, respectively. RNA of pools of all-female and all-male embryos at 2, 5 and 9 dpf were used as template for a custom Agilent eArray hybridization and next generation sequencing. Fifty-nine genes differentially expressed between genders were identified by a false discovery rate of p < 0.05. The most overexpressed genes were amh and tspan8 in males, and cr/20ß-hsd, gpa33, rtn4ipl and zp3 in females (p < 1 × 10-9). Validation of gene expression using qPCR in embryos and gonads indicated copy number variation in tspan8, gpa33, cr/20ß-hsd and amh. Sequencing of amh identified a male-specific duplication of this gene, denoted amhy, differing from the sequence of amh by a 233 bp deletion on exonVII, hence lacking the capability to encode the protein motif that binds to the transforming growth factor beta receptor (TGF-ß domain). amh and amhy segregated in the mapping family in full concordance with SD-linked marker on LG23 signifying the QTL for SD. We discovered 831 microRNAs in tilapia embryos of which nine had sexually dimorphic expression patterns by a false discovery rate of p < 0.05. An up-regulated microRNA in males, pma-mir-4585, was characterized with all six predicted target genes including cr/20ß-hsd, down-regulated in males. CONCLUSIONS: This study reports the first discovery of sexually differentially expressed genes and microRNAs at a very early stage of tilapia embryonic development, i.e. from 2 dpf. Genes with sexually differential expression patterns are enriched for copy number variation. A novel male-specific duplication of amh, denoted amhy, lacking the TGF-ß domain was identified and mapped to the QTL region on LG23 for SD, thus indicating its potential role in SD.

Construction of a microsatellites-based linkage map for the white grouper (Epinephelus aeneus).

The white grouper (Epinephelus aeneus) is a promising candidate for domestication and aquaculture due to its fast growth, excellent taste, and high market price. A linkage map is an essential framework for mapping quantitative trait loci for economic traits and the study of genome evolution. DNA of a single individual was deep-sequenced, and microsatellite markers were identified in 177 of the largest scaffolds of the sequence assembly. The success rate of developing polymorphic homologous markers was 94.9% compared with 63.1% of heterologous markers from other grouper species. Of the 12 adult mature fish present in the broodstock tank, two males and two females were identified as parents of the assigned offspring by parenthood analysis using 34 heterologous markers. A single full-sib family of 48 individuals was established for the construction of first-generation linkage maps based on genotyping data of 222 microsatellites. The markers were assigned to 24 linkage groups in accordance to the 24 chromosomal pairs. The female and male maps consisting of 203 and 202 markers spanned 1053 and 886 cM, with an average intermarker distance of 5.8 and 5.0 cM, respectively. Mapping of markers to linkage groups ends was enriched by using markers originating from scaffolds harboring telomeric repeat-containing RNA. Comparative mapping showed high synteny relationships among the white grouper, kelp grouper (E. bruneus), orange-spotted grouper (E. coioides), and Nile tilapia (Oreochromis niloticus). Thus, it would be useful to integrate the markers that were developed for different groupers, depending on sharing of sequence data, into a comprehensive consensus map.

Gender and sexual behavior modulate the composition of serum lipocalins in Nile tilapia (Oreochromis niloticus).

In tilapia species, plasma lipoproteins with high electrophoretic mobility function in intra- and intergender communication. Blood samples taken at onset and peak of daily sexual activity from dominant and subordinate Oreochromis niloticus males and females were fractionated by native gel electrophoresis and the fast-migrating proteins were subjected to mass spectrometry. Mining the sequence data of the Cichlid Genome Consortium, we identified 11 proteins from the lipocalin super-family and characterized their genes' structures. Phylogenetic and structural analyses subdivided these genes into two classes: (I) 3-coding-exon apolipoproteins and (II) more complex 6-coding-exon sulfide-bond-containing lipocalins. Five apolipoproteins and PTGDSL1, TBTBP, and MSP proteins were modulated by gender and sexual behavior. PTGDSL1 protein was only observed in the plasma serum of dominant males. However, the cysteine residue in the position that is crucial for synthetase activity in mammalian prostaglandin D synthetases was not conserved in PTGDSL1 or PTGDSL2 proteins. In line with previous reports suggesting their involvement in male functions as pheromone transporters, TBTBP and MSP proteins were not detected in females at the onset of daily activity. Their increasing amount in males was concordant with the increase in apolipoproteins AFP4L, APOA4a, APOA4b, APO14kD and APOC2, which were detected exclusively in dominant males, indicating a possible role in mobilization of the energy required to maintain their social hierarchy.

Comparative time-course study on pituitary and branchial response to salinity challenge in Mozambique tilapia (Oreochromis mossambicus) and Nile tilapia (O. niloticus).

The physiological response of Mozambique and Nile tilapia transferred from fresh to brackish (15 ppt) water was compared during a one-week time course. Response in the pituitary was measured by the gene expression pattern of prolactin (PRL I), growth hormone (GH), and calcium-sensing receptor (CaSR), while the response in the gills was measured by the gene expression pattern of the prolactin receptor (PRL-R), Na(+)/K(+)/2Cl(-) cotransporter (NKCC) and Na(+)/Cl(-) cotransporter (NCC), and by activity and expression of Na(+)/K(+)-ATPase (NKA). The time-course curves of plasma osmolality levels indicate a rapid elevation 24 h after transfer, which later decreased and maintained at stable level. PRL I expression decreased in both species, but with stronger response in the Nile tilapia, while no differences were found in the slightly elevated levels of GH mRNA. The branchial response demonstrated a faster up-regulation of NKA and NKCC in the Mozambique tilapia, but similar levels after a week, while Nile tilapia had stronger and constant down-regulation of NCC. The time-course response of the measured osmoregulatory parameters indicate that 24 h after transfer is a critical time point for brackish-water adaptation. The differences in responses to saltwater challenge between Mozambique and Nile tilapia shown in this study may be associated with the differences in saltwater tolerance between these two tilapiine species.

Fine Mapping of a QTL for Fertility on BTA7 and Its Association With a CNV in the Israeli Holsteins.

A quantitative trait locus (QTL) affecting female fertility, scored as the inverse of the number of inseminations to conception, on Bos taurus chromosome 7 was detected by a daughter design analysis of the Israeli Holstein population (P < 0.0003). Sires of five of the 10 families analyzed were heterozygous for the QTL. The 95% confidence interval of the QTL spans 27 cM from the centromere. Seven hundred and four SNP markers on the Illumina BovineSNP50 BeadChip within the QTL confidence interval were tested for concordance. A single SNP, NGS-58779, was heterozygous for all the five QTL heterozygous patriarchs, and homozygous for the remaining five QTL homozygous sires. A significant effect on fertility was associated with this marker in the sample of 900 sires genotyped (P < 10(-6)). Haplotype phase was the same for four of the five segregating sires. Thus concordance was obtained in nine of the ten families. We identified a common haplotype region associated with the rare and economically favorable allele of the SNP, spanning 270 kbp on BTA7 upstream to 4.72 Mbp. Eleven genes found in the common haplotype region should be considered as positional candidates for the identification of the causative quantitative trait nucleotide. Copy number variation was found in one of these genes, KIAA1683. Four gene variants were identified, but only the number of copies of a specific variant (V(1)) was significantly associated with breeding values of sires for fertility.

Identification of repetitive elements in the genome of Oreochromis niloticus: tilapia repeat masker.

The large-scale bacterial artificial chromosome-end sequencing project of Nile tilapia (Oreochromis niloticus) has generated extensive sequence data that allowed the examination of the repeat content in this fish genome and building of a repeat library specific for this species. This library was established based on Tilapiini repeat sequences from GenBank, sequences orthologous to the repeat library of zebrafish in Repbase, and novel repeats detected by genome analysis using MIRA assembler. We estimate that repeats constitute about 14% of the tilapia genome and also give estimates for the occurrence of the different repeats based on the Basic Local Alignment Search Tool searches within the database of known tilapia sequences. The frequent occurrence of novel repeats in the tilapia genome indicates the importance of using the species-specific repeat masker prior to sequence analyses. A web tool based on the RepeatMasker software was designed to assist tilapia genomics.

Combining mouse mammary gland gene expression and comparative mapping for the identification of candidate genes for QTL of milk production traits in cattle.

BACKGROUND: Many studies have found segregating quantitative trait loci (QTL) for milk production traits in different dairy cattle populations. However, even for relatively large effects with a saturated marker map the confidence interval for QTL location by linkage analysis spans tens of map units, or hundreds of genes. Combining mapping and arraying has been suggested as an approach to identify candidate genes. Thus, gene expression analysis in the mammary gland of genes positioned in the confidence interval of the QTL can bridge the gap between fine mapping and quantitative trait nucleotide (QTN) determination. RESULTS: We hybridized Affymetrix microarray (MG-U74v2), containing 12,488 murine probes, with RNA derived from mammary gland of virgin, pregnant, lactating and involuting C57BL/6J mice in a total of nine biological replicates. We combined microarray data from two additional studies that used the same design in mice with a total of 75 biological replicates. The same filtering and normalization was applied to each microarray data using GeneSpring software. Analysis of variance identified 249 differentially expressed probe sets common to the three experiments along the four developmental stages of puberty, pregnancy, lactation and involution. 212 genes were assigned to their bovine map positions through comparative mapping, and thus form a list of candidate genes for previously identified QTLs for milk production traits. A total of 82 of the genes showed mammary gland-specific expression with at least 3-fold expression over the median representing all tissues tested in GeneAtlas. CONCLUSION: This work presents a web tool for candidate genes for QTL (cgQTL) that allows navigation between the map of bovine milk production QTL, potential candidate genes and their level of expression in mammary gland arrays and in GeneAtlas. Three out of four confirmed genes that affect QTL in livestock (ABCG2, DGAT1, GDF8, IGF2) were over expressed in the target organ. Thus, cgQTL can be used to determine priority of candidate genes for QTN analysis based on differential expression in the target organ.

Amh and Dmrta2 genes map to tilapia (Oreochromis spp.) linkage group 23 within quantitative trait locus regions for sex determination.

Recent studies have revealed that the major genes of the mammalian sex determination pathway are also involved in sex determination of fish. Several studies have reported QTL in various species and strains of tilapia, regions contributing to sex determination have been identified on linkage groups 1, 3, and 23. Genes contributing to sex-specific mortality have been detected on linkage groups 2, 6, and 23. To test whether the same genes might control sex determination in mammals and fishes, we mapped 11 genes that are considered putative master key regulators of sex determination: Amh, Cyp19, Dax1, Dmrt2, Dmrta2, Fhl3l, Foxl2, Ixl, Lhx9, Sf1, and Sox8. We identified polymorphisms in noncoding regions of these genes and genotyped these sites for 90 individuals of an F2 mapping family. Mapping of Dax1 joined LG16 and LG21 into a single linkage group. The Amh and Dmrta2 genes were mapped to two distinct regions of LG23. The Amh gene was mapped 5 cM from UNH879 within a QTL region for sex determination and 2 cM from UNH216 within a QTL region for sex-specific mortality. Dmrta2 was mapped 4 cM from UNH848 within another QTL region for sex determination. Cyp19 was mapped to LG1 far from a previously reported QTL region for sex determination on this chromosome. Seven other candidate genes mapped to LG4, -11, -12, -14, and -17.

Accelerated expansion of group IID-like phospholipase A2 genes in Bos taurus.

Low-molecular-weight, calcium-dependent phospholipase A2 genes (PLA2s) that belong to the secreted type of PLA2s are clustered within a syntenic group on human 1p35-p36 and mouse 4qD3. We reassembled trace files available from the Whole Genome Sequencing (WGS) Project, obtaining an 86-kb contig with three tandem PLA2G2D duplications in the Hereford strain. We used mate-pair data to monitor the assembly and to exclude chimeric clones, demonstrating that the current WGS data may be assembled even in a highly repetitive region with a coverage exceeding fivefold. The genomic structure indicated that most of the PLA2G2D transcripts are formed by four exons. Two alternative first exons were present in all duplications. In two duplications insertions of satellite DNA in the third intron created a novel exon that gave rise to a two-exon product. Linkage and comparative mapping placed the bovine PLA2G2 locus on BTA2, indicating that it evolved from an ancestral PLA2G2D locus common to human, cattle, and rodents. Bovine PLA2G2D variants were capable of encoding 147-amino-acid polypeptides that consisted of putative signal peptide and metal-binding domains. Cysteine residues were conserved in positions analogous to those forming the seven disulfide bonds characteristic of PLA2G2 genes. Quantitative PCR analysis of bovine PLA2G2D transcripts indicated that their expression levels varied between the dry period and lactation in the mammary gland samples and that their expression was polymorphic in liver tissue. The recent burst of duplication and divergence of the bovine PLA2G2D genes and their polymorphic nature are typical of innate immune response genes.

Efficient inference of haplotypes from genotypes on a large animal pedigree.

We present a simple algorithm for reconstruction of haplotypes from a sample of multilocus genotypes. The algorithm is aimed specifically for analysis of very large pedigrees for small chromosomal segments, where recombination frequency within the chromosomal segment can be assumed to be zero. The algorithm was tested both on simulated pedigrees of 155 individuals in a family structure of three generations and on real data of 1149 animals from the Israeli Holstein dairy cattle population, including 406 bulls with genotypes, but no females with genotypes. The rate of haplotype resolution for the simulated data was >91% with a standard deviation of 2%. With 20% missing data, the rate of haplotype resolution was 67.5% with a standard deviation of 1.3%. In both cases all recovered haplotypes were correct. In the real data, allele origin was resolved for 22% of the heterozygous genotypes, even though 70% of the genotypes were missing. Haplotypes were resolved for 36% of the males. Computing time was insignificant for both data sets. Despite the intricacy of large-scale real pedigree genotypes, the proposed algorithm provides a practical rule-based solution for resolving haplotypes for small chromosomal segments in commercial animal populations.

Altered self-erythrocyte recognition and destruction in an inbred line of tilapia (Oreochromis aureus).

Carboxyfluorescein diacetate (cFDA)-stained autologous and syngeneic tilapia (Oreochromis aureus) erythrocytes are recognized by effector peripheral blood leukocytes and lysed after a short culture period of 4 h. The hemolysis level was evaluated by measuring the fluorescence of the released cFDA. The degree of lysis of stained target erythrocytes of 60 individuals revealed a trimodal distribution statistically stratified into three groups of low (LR), intermediate (IR), and high (HR) responders. Depletion of the majority of phagocytes from leukocytes lowered the lysis level of HR to that of LR. A highly significant increase of LR cytotoxicity was obtained after the addition of conditioned medium from HR but only in the presence of phagocytes. Genetic analysis of offspring from four crosses (IR x HR, IR x LR, HR x LR, and LR x LR) revealed a quantitative trait locus (QTL) segregating for the level of response linked to markers UNH207 and UNH231 on linkage group 6 of tilapia. Based on segregation analysis of 58 gynogenetic BIU-1 offspring, the distances from the centromere were estimated as 21.5, 11.5, and 9.0 cM for UNH207, UNH231, and the QTL, respectively. It is suggested that 1) self-target recognition and destruction requires both cFDA-altered self-erythrocyte membrane and membrane structures normally present in autologous, syngeneic, and xenogeneic targets; 2) natural cytotoxic cells and/or macrophages are involved in erythrocyte lysis; and 3) the lysis level is codominantly inherited by a QTL segregating on tilapia linkage group 6.

Correcting for bias in estimation of quantitative trait loci effects.

Estimates of quantitative trait loci (QTL) effects derived from complete genome scans are biased, if no assumptions are made about the distribution of QTL effects. Bias should be reduced if estimates are derived by maximum likelihood, with the QTL effects sampled from a known distribution. The parameters of the distributions of QTL effects for nine economic traits in dairy cattle were estimated from a daughter design analysis of the Israeli Holstein population including 490 marker-by-sire contrasts. A separate gamma distribution was derived for each trait. Estimates for both the alpha and beta parameters and their SE decreased as a function of heritability. The maximum likelihood estimates derived for the individual QTL effects using the gamma distributions for each trait were regressed relative to the least squares estimates, but the regression factor decreased as a function of the least squares estimate. On simulated data, the mean of least squares estimates for effects with nominal 1% significance was more than twice the simulated values, while the mean of the maximum likelihood estimates was slightly lower than the mean of the simulated values. The coefficient of determination for the maximum likelihood estimates was five-fold the corresponding value for the least squares estimates.