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.

Estimating the Effect of the Kappa Casein Genotype on Milk Coagulation Properties in Israeli Holstein Cows.

In Israel, about 26% of produced milk is used to produce hard cheeses and 29% for soft cheeses. Milk with preferred coagulation properties requires a shorter coagulation time and yields a higher curd firmness than milk with inferior coagulation properties. Studies have shown that milk from cows with the B allele of kappa casein (κ-CN) produces more cheese than milk from those with A and E alleles. There is evidence that milk from AE or EE genotype cows is unsuitable for cheese production. In the early 1990s, the proportion of the B allele in Israeli Holstein cattle was about 17%, similar to its prevalence in the Holstein population worldwide. In recent years, however, its proportion has increased to about 40%. We analyzed milk coagulation properties as a function of the cow's κ-CN genotype, including time in minutes until the beginning of coagulation and curd firmness after 60 min-measured in volts via an optigraph device and scored on a scale of 0-4 by a laboratory technician. Cow selection was based on their sire's genotype, so that there would be sufficient genotypes that include the rare E allele. A total of 359 cows were sampled from 15 farms: 64 with genotype AA, 142 with AB, 41 with AE, 65 with BB, and 47 with BE. Data were analyzed via the general linear model procedure of SAS. We found the following: (a) There were significant differences between genotypes for optigraph-measured curd firmness. In a multi-comparison test, the BB genotype gave the highest curd firmness, and AB and BE showed a significant advantage compared to AA and AE (9.4, 8.6, 8.4, 6.9, 6.8 V, respectively). Assuming a frequency of about 55% for the A allele, about 30% of the milk delivered to dairy plants comes from AA cows. (b) There was a significant difference between the genotypes in technician-observed curd firmness, with BB scoring significantly higher than AA and AE. (c) The optigraph-measured curd firmness was significantly higher for milk from primiparous cows as compared to milk from second, third, or fourth lactation cows (8.9, 7.8, 7.9, 7.7 V, respectively). The technician-observed curd firmness was significantly higher for primiparous vs. multiparous cows. There was a clear advantage in curd firmness for genotypes that included the B allele compared to those with AA and AE genotypes. We can increase the proportion of the B allele in the population by insemination of cows using bulls with the genotypes AB and BB. This factor should therefore be included in the selection index.

Analysis of copy loss and gain variations in Holstein cattle autosomes using BeadChip SNPs.

BACKGROUND: Copy number variation (CNV) has been recently identified in human and other mammalian genomes, and there is a growing awareness of CNV's potential as a major source for heritable variation in complex traits. Genomic selection is a newly developed tool based on the estimation of breeding values for quantitative traits through the use of genome-wide genotyping of SNPs. Over 30,000 Holstein bulls have been genotyped with the Illumina BovineSNP50 BeadChip, which includes 54,001 SNPs (~SNP/50,000 bp), some of which fall within CNV regions. RESULTS: We used the BeadChip data obtained for 912 Israeli bulls to investigate the effects of CNV on SNP calls. For each of the SNPs, we estimated the frequencies of occurrence of loss of heterozygosity (LOH) and of gain, based either on deviation from the expected Hardy-Weinberg equilibrium (HWE) or on signal intensity (SI) using the PennCNV "detect" option. Correlations between LOH/CNV frequencies predicted by the two methods were low (up to r = 0.08). Nevertheless, 418 locations displayed significantly high frequencies by both methods. Efficiency of designating large genomic clusters of olfactory receptors as CNVs was 29%. Frequency values for copy loss were distinguishable in non-autosomal regions, indicating misplacement of a region in the current BTA7 map. Analysis of BTA18 placed major quantitative trait loci affecting net merit in the US Holstein population in regions rich in segmental duplications and CNVs. Enrichment of transporters in CNV loci suggested their potential effect on milk-production traits. CONCLUSIONS: Expansion of HWE and PennCNV analyses allowed estimating LOH/CNV frequencies, and combining the two methods yielded more sensitive detection of inherited CNVs and better estimation of their possible effects on cattle genetics. Although this approach was more effective than methodologies previously applied in cattle, it has severe limitations. Thus the number of CNVs reported here for the Holstein breed may represent as little as one-tenth of inherited common structural variation.

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.

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.

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.

Current and future developments in patents for quantitative trait loci in dairy cattle.

Many studies have proposed that rates of genetic gain in dairy cattle can be increased by direct selection on the individual quantitative loci responsible for the genetic variation in these traits, or selection on linked genetic markers. The development of DNA-level genetic markers has made detection of QTL nearly routine in all major livestock species. The studies that attempted to detect genes affecting quantitative traits can be divided into two categories: analysis of candidate genes, and genome scans based on within-family genetic linkage. To date, 12 patent cooperative treaty (PCT) and US patents have been registered for DNA sequences claimed to be associated with effects on economic traits in dairy cattle. All claim effects on milk production, but other traits are also included in some of the claims. Most of the sequences found by the candidate gene approach are of dubious validity, and have been repeated in only very few independent studies. The two missense mutations on chromosomes 6 and 14 affecting milk concentration derived from genome scans are more solidly based, but the claims are also disputed. A few PCT in dairy cattle are commercialized as genetic tests where commercial dairy farmers are the target market.

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.

Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle.

We previously localized a quantitative trait locus (QTL) on chromosome 6 affecting milk fat and protein concentration to a 4-cM confidence interval, centered on the microsatellite BM143. We characterized the genes and sequence variation in this region and identified common haplotypes spanning five polymorphic sites in the genes IBSP, SPP1, PKD2, and ABCG2 for two sires heterozygous for this QTL. Expression of SPP1 and ABCG2 in the bovine mammary gland increased from parturition through lactation. SPP1 and all the coding exons of ABCG2 and PKD2 were sequenced for these two sires. The single nucleotide change capable of encoding a substitution of tyrosine-581 to serine (Y581S) in the ABCG2 transporter was the only polymorphism corresponding to the segregation status of all 3 heterozygous and 15 homozygous sires for the QTL in the Israeli and U.S. Holstein populations. The allele substitution fixed effects on the genetic evaluations of 335 Israeli sires were -341 kg milk, +0.16% fat, and +0.13% protein (F-value = 200). No other polymorphism gave significant effect for fat and protein concentration in models that also included Y581S. The allele substitution effects on the genetic evaluations of 670 cows, daughters of two heterozygous sires, were -226 kg milk, 0.09% fat, and 0.08% protein (F-value = 394), with partial dominance towards the 581S homozygotes. We therefore propose that Y581S in ABCG2 is the causative site for this QTL.

Cloning and characterization of FAM13A1--a gene near a milk protein QTL on BTA6: evidence for population-wide linkage disequilibrium in Israeli Holsteins.

A cluster of genes coding for proteins of the extracellular matrix (ECM) containing sequence motifs essential for integrin-receptor interactions is located on HSA4q21 and on BTA6, within the critical region of a quantitative trait locus (QTL) affecting milk protein production. Genes within this cluster are involved in the formation of bone and lobuloalveolar structures in mammary gland and in kidney function. We cloned a bovine gene neighboring this ECM cluster, termed FAM13A1, the first member of a novel gene family (FAM13). A short predominant 5.1-kb mRNA variant capable of encoding 697 amino acids is transcribed from 18-exon orthologous genes in human and mouse. All putative protein orthologs contained a bipartite nuclear-localization signal and two coil-coiled domains. We detected two other FAM13 paralogs, C5ORF5 and FAM13C1, on HSA5q31 and HSA10q21, respectively. All FAM13 paralogs produce transcripts that are complementary to adjacent genes, suggesting that antisense transcription may regulate their functions. The structure of a longer 5.9-kb variant, FAM13A1_v2, that has an extra 7-exon putative RhoGAP domain at the N-terminus provides further clues to FAM13 function. Analysis of 400 bulls revealed a population-wide linkage disequilibrium between FAM13A1 polymorphisms and the QTL.