Alternative Traits for Genetic Evaluation of Mastitis Based on Lifetime Merit

Genetic selection has achieved little progress in reducing mastitis incidence. Mastitis traits are problematic due to the lack of sensitivity of the data and reliance on clinical diagnosis, often missing subclinical cases, and/or on monthly somatic cell count (SCC) measurements. The current measure for mastitis is the lactation average of the somatic cells score (LSCS). We studied two datasets: (1) 148 heifers divided into non-intramammary infected, sub-clinically infected and clinical mastitis groups; (2) data from 89,601 heifers from Israeli Holsteins through the same period divided into "udder healthy" (UH) and "non-healthy" (UNH) by a threshold of SCC 120,000 cells/mL in all nine monthly milk recordings. In study 1, non-infected heifers had significantly (p < 0.05) more partum, production days and overall lifetime milk production compared to clinical and sub-clinically infected. In study 2, UH heifers (20.3%) had significantly higher (p < 0.01) lifetime milk, production days, and lactations. Subdividing datasets by sires, the same analyses detected differences in percentages of UH daughters between the sire groups. Lifetime milk production correlated (r = +0.83, p < 0.001) with udder health status. SCC threshold of less than 120,000 cells/mL during all first lactation measurements indicated healthy udder, providing a valuable insight that this dichotomous trait is advantageous for calculating lifetime net-merit index (NM$) over LSCS.

Genetic and Genomic Analysis of Cow Mortality in the Israeli Holstein Population

“Livability” was defined as the inverse of the probability of death. The objectives of this study were to estimate the heritability, genetic and phenotypic trends for the livability of Israeli Holstein cows; estimate the genetic and environmental correlations between livability and the nine traits included in the Israeli breeding index; estimate the effect of the inclusion of livability in the Israeli breeding index on expected genetic gains; and compute a genome-wide association study (GWAS) for livability. Seven data sets were analyzed. All data were derived from the database of the Israeli dairy cattle herd-book. The mean livability for the complete data set of 523,954 cows born from 2000 through 2016 was 89.6%. Pregnancy reduced livability by 15%. Livability generally increased with parity and days in milk within parity. Heritability of livability was 0.0082. Phenotypic and genetic trends over the 14-year period from 2000 through 2013 were −0.42% and −0.22% per year. If livability is included in the Israeli breeding index, accounting for 9% of the index, livability would increase by 1.3% and protein production would decrease by 11 kg over the next decade, as compared to the current index. A marker in proximity to the oxytocin–vasopressin locus had the greatest effect in the GWAS. Oxytocin activity in cattle affects calving-associated pathologies and maternal death. Inclusion of livability in the Israeli breeding index is not recommended.

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.

Towards the detection of copy number variation from single sperm sequencing in cattle

BACKGROUND

Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms.

RESULTS

We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull's family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats.

CONCLUSION

Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.

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⁻⁸), 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.

Comparing BeadChip and WGS Genotyping: Non-Technical Failed Calling Is Attributable to Additional Variation within the Probe Target Sequence

Microarray-based genomic selection is a central tool to increase the genetic gain of economically significant traits in dairy cattle. Yet, the effectivity of this tool is slightly limited, as estimates based on genotype data only partially explain the observed heritability. In the analysis of the genomes of 17 Israeli Holstein bulls, we compared genotyping accuracy between whole-genome sequencing (WGS) and microarray-based techniques. Using the standard GATK pipeline, the short-variant discovery within sequence reads mapped to the reference genome (ARS-UCD1.2) was compared to the genotypes from Illumina BovineSNP50 BeadChip and to an alternative method, which computationally mimics the hybridization procedure by mapping reads to 50 bp spanning the BeadChip source sequences. The number of mismatches between the BeadChip and WGS genotypes was low (0.2%). However, 17,197 (40% of the informative SNPs) had extra variation within 50 bp of the targeted SNP site, which might interfere with hybridization-based genotyping. Consequently, with respect to genotyping errors, BeadChip varied significantly and systematically from WGS genotyping, introducing null allele-like effects and Mendelian errors (<0.5%), whereas the GATK algorithm of local de novo assembly of haplotypes successfully resolved the genotypes in the extra-variable regions. These findings suggest that the microarray design should avoid polymorphic genomic regions that are prone to extra variation and that WGS data may be used to resolve erroneous genotyping, which may partially explain missing heritability.

Genome-wide recombination map construction from single sperm sequencing in cattle

Background

Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle.

Results

To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle.

Conclusions

To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08415-w.

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.

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.

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 ; R²  = 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.

Estimating Copy-Number Proportions: The Comeback of Sanger Sequencing

Determination of the relative copy numbers of mixed molecular species in nucleic acid samples is often the objective of biological experiments, including Single-Nucleotide Polymorphism (SNP), indel and gene copy-number characterization, and quantification of CRISPR-Cas9 base editing, cytosine methylation, and RNA editing. Standard dye-terminator chromatograms are a widely accessible, cost-effective information source from which copy-number proportions can be inferred. However, the rate of incorporation of dye terminators is dependent on the dye type, the adjacent sequence string, and the secondary structure of the sequenced strand. These variable rates complicate inferences and have driven scientists to resort to complex and costly quantification methods. Because these complex methods introduce their own biases, researchers are rethinking whether rectifying distortions in sequencing trace files and using direct sequencing for quantification will enable comparable accurate assessment. Indeed, recent developments in software tools (e.g., TIDE, ICE, EditR, BEEP and BEAT) indicate that quantification based on direct Sanger sequencing is gaining in scientific acceptance. This commentary reviews the common obstacles in quantification and the latest insights and developments relevant to estimating copy-number proportions based on direct Sanger sequencing, concluding that bidirectional sequencing and sophisticated base calling are the keys to identifying and avoiding sequence distortions.

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.

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.

Basal Levels of CD18 Antigen Presenting Cells in Cow Milk Associate with Copy Number Variation of Fc Gamma Receptors

Differentiation of cells by flow cytometry provides informative somatic cell counts (SCCs) that allow analyzing leukocyte population patterns in udder infections of different etiologies. Postulating that this approach also enhances the statistical power to detect genetic variants linked to cell levels in milk of healthy mammary glands, we used monoclonal antibodies anti-CD18, anti-CD4, anti--CD14, and anti-PMN to count cells presenting these surface antigens, and performed a genome-wide association study of these counts in 125 Israeli Holsteins genotyped using SNP BeadChips. We identified an informative haplotype of 15 SNPs in the centromeric end of BTA3 that was strongly associated with CD18 cells (p < 2.3 × 10-9). Within this region, examination of the network of genes interacting with ITGB2 (CD18) indicated an Fc-γ-receptor gene cluster, including FCGR2A (CD32). Sanger-sequence analysis of FCGR2s-linked exon 3 variation to CD18 counts. Meta-analysis of RNA-Seq data revealed a significant negative correlation (R = -0.51) between expression of CD32 and CD18 in milk. Assembly of DNA-Seq reads uncovered FCGR copy-number variation and a variant, designated V7, was abundant in dairy cattle, probably reflecting adaptation to selection pressure for low SCC in Holstein milk.

Avian Leptin: Bird's-Eye View of the Evolution of Vertebrate Energy-Balance Control

Discovery of the satiety hormone leptin in 1994 and its characterization in mammals provided a key tool to deciphering the complex mechanism governing adipose tissue regulation of appetite and energy expenditure. Surprisingly, despite the perfectly logical notion of an energy-storing tissue announcing the amount of fat stores using leptin signaling, alternate mechanisms were chosen in bird evolution. This conclusion emerged based on the recent discovery and characterization of genuine avian leptin - after it had been assumed missing by some, and erroneously identified by others. Critical evaluation of the past and present indications of the role of leptin in Aves provides a new perspective on the evolution of energy-balance control in vertebrates; proposing a regulation strategy alternative to the adipostat mechanism.

Bos taurus-indicus hybridization correlates with intralocus sexual-conflict effects of PRDM9 on male and female fertility in Holstein cattle

BACKGROUND

Crossover localization during meiotic recombination is mediated by the fast-evolving zinc-finger (ZnF) domain of gene PRDM9. To study its impact on dairy cattle performance, we compared its genetic variation between the relatively small Israeli (IL) Holsteins and the North American (US) Holsteins that count millions.

RESULTS

Initially, we analyzed the major BTA1 haplotypes present in IL Holsteins based on the 10 most telomeric SNPs of the BovineSNP50 BeadChip. Sequencing of representative haplotype carriers indicated that for all frequent haplotypes (> 6%), the variable PRDM9 ZnF array consisted of seven tandem ZnF repeats. Two rare haplotypes (frequency < 4%) carried an indicine PRDM9, whereas all others were variants of the taurine type. These two haplotypes included the minor SNP allele, which was perfectly linked with a previously described PRDM9 allele known to induce unique localization of recombination hotspots. One of them had a significant (p = 0.03) negative effect on IL sire fertility. This haplotype combined the rare minor alleles of the only SNPs with significant (p < 0.05) negative substitution effects on US sire fertility (SCR). Analysis of telomeric SNPs indicated general agreement of allele frequencies (R = 0.95) and of the substitution effects on sire fertility (SCR, R = 0.6) between the US and IL samples. Surprisingly, the alleles that had a negative impact on male fertility had the most positive substitution effects on female fertility traits (DPR, CCR and HCR).

CONCLUSIONS

A negative genetic correlation between male and female fertility is encoded within the BTA1 telomere. Cloning the taurine PRDM9 gene, which is the common form carried by Holsteins, we encountered the infiltration of an indicine PRDM9 variant into this population. During meiosis, in heterozygous males, the indicine PRDM9 variant may induce incompatibility of recombination hotspots and male infertility. However, this variant is associated with favorable female fertility, which would explain its survival and the general negative correlation (R = - 0.3) observed between male and female fertility in US Holsteins. Further research is needed to explain the mechanism underlying this positive effect and to devise a methodology to unlink it from the negative effect on male fertility during breeding.

Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins

BACKGROUND

Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock.

RESULTS

We identified 1043 CNV regions (CNVRs) from array comparative genomic hybridization (CGH) data of 47 Holstein bulls. Using a probe-based CNV association approach, we detected 87 CNVRs significantly (Bonferroni-corrected P value < 0.05) associated with at least one out of 41 complex traits. Within them, 39 CNVRs were simultaneously associated with at least 2 complex traits. Notably, 24 CNVRs were markedly related to daughter pregnancy rate (DPR). For example, CNVR661 containing CYP4A11 and CNVR213 containing CTR9, respectively, were associated with DPR and other traits related to reproduction, production, and body conformation. CNVR758 was also significantly related to DPR, with a nearby gene CAPZA3, encoding one of F-actin-capping proteins which play a role in determining sperm architecture and male fertility. We corroborated these CNVRs by examining their overlapped quantitative trait loci and comparing with previously published CNV results.

CONCLUSION

To our knowledge, this is one of the first genome-wide association studies based on CNVs called by array CGH in Holstein cattle. Our results contribute substantial information about the potential CNV impacts on reproduction, health, production, and body conformation traits, which lay the foundation for incorporating CNV into the future dairy cattle breeding program.

Application of pancreatic phospholipase A2 for treatment of bovine mastitis

Recent findings have indicated that secreted phospholipases A2 (sPLA2s) have anti-inflammatory functions, including relief of symptoms in a mouse model of mastitis. This prompted us to investigate the therapeutic application of sPLA2, PLA2G1B, for bovine mastitis. Initial testing of PLA2G1B's effect on bovine mammary epithelial cell (bMEC) line PS revealed no changes in cell viability or cytokine-secretion pattern. However, when cells were first treated with lipopolysaccharide endotoxin (LPS) or live bacteria (Escherichia coli or Staphylococcus aureus), incubation with PLA2G1B significantly improved cell viability, suggesting involvement of sPLA2s in protecting membranes from lipid-peroxidation damage, rather than a bactericidal action. When PLA2G1B was applied simultaneously with LPS, a significant short-term reduction in interleukin-8 secretion was observed compared with bMECs treated only with LPS, supporting previous reports that PLA2G1B affects interleukin-8 signaling in similar cells. Following the favorable outcome of the in vitro experiments, we tested PLA2G1B in vivo by mammary infusion into infected glands. In one of a small sample (n = 4) of lactating cows chronically infected with Streptococcus dysgalactiae, a single PLA2G1B treatment completely cleared inflammation and bacteria, demonstrating its potential to cure subclinical mastitis. PLA2G1B treatment did not affect coagulase-negative staphylococci infection. These types of mastitis may involve formation of a resistant biofilm, and its elimination may relate to sPLA2s' characteristic ability to aggregate with cellular debris, facilitating their internalization by macrophages. In a bovine model of clinical mastitis based on introduction of E. coli via the streak canal, a single mammary infusion of PLA2G1B led to faster recovery to pre-infection milk-yield levels and decrease of somatic cell counts. In this case, all of sPLA2s' modes of resolving inflammation may apply, including competitive binding of the sPLA2s’ receptor, the inactivation of which confers resistance to endotoxic shock. Hence, this study strongly supports further research into PLA2G1B as a cure for bovine mastitis.

Characterization of recombination features and the genetic basis in multiple cattle breeds

BACKGROUND

Crossover generated by meiotic recombination is a fundamental event that facilitates meiosis and sexual reproduction. Comparative studies have shown wide variation in recombination rate among species, but the characterization of recombination features between cattle breeds has not yet been performed. Cattle populations in North America count millions, and the dairy industry has genotyped millions of individuals with pedigree information that provide a unique opportunity to study breed-level variations in recombination.

RESULTS

Based on large pedigrees of Jersey, Ayrshire and Brown Swiss cattle with genotype data, we identified over 3.4 million maternal and paternal crossover events from 161,309 three-generation families. We constructed six breed- and sex-specific genome-wide recombination maps using 58,982 autosomal SNPs for two sexes in the three dairy cattle breeds. A comparative analysis of the six recombination maps revealed similar global recombination patterns between cattle breeds but with significant differences between sexes. We confirmed that male recombination map is 10% longer than the female map in all three cattle breeds, consistent with previously reported results in Holstein cattle. When comparing recombination hotspot regions between cattle breeds, we found that 30% and 10% of the hotspots were shared between breeds in males and females, respectively, with each breed exhibiting some breed-specific hotspots. Finally, our multiple-breed GWAS found that SNPs in eight loci affected recombination rate and that the PRDM9 gene associated with hotspot usage in multiple cattle breeds, indicating a shared genetic basis for recombination across dairy cattle breeds.

CONCLUSIONS

Collectively, our results generated breed- and sex-specific recombination maps for multiple cattle breeds, provided a comprehensive characterization and comparison of recombination patterns between breeds, and expanded our understanding of the breed-level variations in recombination features within an important livestock species.

Comparative omics and feeding manipulations in chicken indicate a shift of the endocrine role of visceral fat towards reproduction

Background

The mammalian adipose tissue plays a central role in energy-balance control, whereas the avian visceral fat hardly expresses leptin, the key adipokine in mammals. Therefore, to assess the endocrine role of adipose tissue in birds, we compared the transcriptome and proteome between two metabolically different types of chickens, broilers and layers, bred towards efficient meat and egg production, respectively.

Results

Broilers and layer hens, grown up to sexual maturation under free-feeding conditions, differed 4.0-fold in weight and 1.6-fold in ovarian-follicle counts, yet the relative accumulation of visceral fat was comparable. RNA-seq and mass-spectrometry (MS) analyses of visceral fat revealed differentially expressed genes between broilers and layers, 1106 at the mRNA level (FDR ≤ 0.05), and 203 at the protein level (P ≤ 0.05). In broilers, Ingenuity Pathway Analysis revealed activation of the PTEN-pathway, and in layers increased response to external signals. The expression pattern of genes encoding fat-secreted proteins in broilers and layers was characterized in the RNA-seq and MS data, as well as by qPCR on visceral fat under free feeding and 24 h-feed deprivation. This characterization was expanded using available RNA-seq data of tissues from red junglefowl, and of visceral fat from broilers of different types. These comparisons revealed expression of new adipokines and secreted proteins (LCAT, LECT2, SERPINE2, SFTP1, ZP1, ZP3, APOV1, VTG1 and VTG2) at the mRNA and/or protein levels, with dynamic gene expression patterns in the selected chicken lines (except for ZP1; FDR/P ≤ 0.05) and feed deprivation (NAMPT, SFTPA1 and ZP3) (P ≤ 0.05). In contrast, some of the most prominent adipokines in mammals, leptin, TNF, IFNG, and IL6 were expressed at a low level (FPKM/RPKM< 1) and did not show differential mRNA expression neither between broiler and layer lines nor between fed vs. feed-deprived chickens.

Conclusions

Our study revealed that RNA and protein expression in visceral fat changes with selective breeding, suggesting endocrine roles of visceral fat in the selected phenotypes. In comparison to gene expression in visceral fat of mammals, our findings points to a more direct cross talk of the chicken visceral fat with the reproductive system and lower involvement in the regulation of appetite, inflammation and insulin resistance.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4675-0) contains supplementary material, which is available to authorized users.

Mapping of leptin and its syntenic genes to chicken chromosome 1p

Background

Misidentification of the chicken leptin gene has hampered research of leptin signaling in this species for almost two decades. Recently, the genuine leptin gene with a GC-rich (~70%) repetitive-sequence content was identified in the chicken genome but without indicating its genomic position. This suggests that such GC-rich sequences are difficult to sequence and therefore substantial regions are missing from the current chicken genome assembly.

Results

A radiation hybrid panel of chicken-hamster Wg3hCl2 cells was used to map the genome location of the chicken leptin gene. Contrary to our expectations, based on comparative genome mapping and sequence characteristics, the chicken leptin was not located on a microchromosome, which are known to contain GC-rich and repetitive regions, but at the distal tip of the largest chromosome (1p). Following conserved synteny with other vertebrates, we also mapped five additional genes to this genomic region (ARF5, SND1, LRRC4, RBM28, and FLNC), bridging the genomic gap in the current Galgal5 build for this chromosome region. All of the short scaffolds containing these genes were found to consist of GC-rich (54 to 65%) sequences comparing to the average GC-content of 40% on chromosome 1. In this syntenic group, the RNA-binding protein 28 (RBM28) was in closest proximity to leptin. We deduced the full-length of the RBM28 cDNA sequence and profiled its expression patterns detecting a negative correlation (R = − 0.7) between the expression of leptin and of RBM28 across tissues that expressed at least one of the genes above the average level. This observation suggested a local regulatory interaction between these genes. In adipose tissues, we observed a significant increase in RBM28 mRNA expression in breeds with lean phenotypes.

Conclusion

Mapping chicken leptin together with a cluster of five syntenic genes provided the final proof for its identification as the true chicken ortholog. The high GC-content observed for the chicken leptin syntenic group suggests that other similar clusters of genes in GC-rich genomic regions are missing from the current genome assembly (Galgal5), which should be resolved in future assemblies of the chicken genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0543-1) contains supplementary material, which is available to authorized users.

Gene Augmentation Therapy for a Missense Substitution in the cGMP-Binding Domain of Ovine CNGA3 Gene Restores Vision in Day-Blind Sheep

Purpose

Applying CNGA3 gene augmentation therapy to cure a novel causative mutation underlying achromatopsia (ACHM) in sheep.

Methods

Impaired vision that spontaneously appeared in newborn lambs was characterized by behavioral, electroretinographic (ERG), and histologic techniques. Deep-sequencing reads of an affected lamb and an unaffected lamb were compared within conserved genomic regions orthologous to human genes involved in similar visual impairment. Observed nonsynonymous amino acid substitutions were classified by their deleteriousness score. The putative causative mutation was assessed by producing compound CNGA3 heterozygotes and applying gene augmentation therapy using the orthologous human cDNA.

Results

Behavioral assessment revealed day blindness, and subsequent ERG examination showed attenuated photopic responses. Histologic and immunohistochemical examination of affected sheep eyes did not reveal degeneration, and cone photoreceptors expressing CNGA3 were present. Bioinformatics and sequencing analyses suggested a c.1618G>A, p.Gly540Ser substitution in the GMP-binding domain of CNGA3 as the causative mutation. This was confirmed by genetic concordance test and by genetic complementation experiment: All five compound CNGA3 heterozygotes, carrying both p.Arg236* and p.Gly540Ser mutations in CNGA3, were day-blind. Furthermore, subretinal delivery of the intact human CNGA3 gene using an adeno-associated viral vector (AAV) restored photopic vision in two affected p.Gly540Ser homozygous rams.

Conclusions

The c.1618G>A, p.Gly540Ser substitution in CNGA3 was identified as the causative mutation for a novel form of ACHM in Awassi sheep. Gene augmentation therapy restored vision in the affected sheep. This novel mutation provides a large-animal model that is valid for most human CNGA3 ACHM patients; the majority of them carry missense rather than premature-termination mutations.

Correspondence on Lovell et al.: identification of chicken genes previously assumed to be evolutionarily lost

Through RNA-Seq analyses, we identified 137 genes that are missing in chicken, including the long-sought-after nephrin and tumor necrosis factor genes. These genes tended to cluster in GC-rich regions that have poor coverage in genome sequence databases. Hence, the occurrence of syntenic groups of vertebrate genes that have not been observed in Aves does not prove the evolutionary loss of such genes.Please see related Research article: http://dx.doi.org/10.1186/s13059-014-0565-1 and Please see response from Lovell et al: https://www.dx.doi.org/10.1186/s13059-017-1234-y.

Sequence motifs capable of forming DNA stem-loop structures act as a replication diode

Calculating peak‐height ratios between single‐nucleotide polymorphisms (SNP) alleles in sequencing chromatograms is a practical method for estimating their copy number proportions (CNPs). However, it is surprising that sequencing DNA from different directions might yield different results. We analyzed three adjacent SNPs within the ovine period circadian‐clock 2 (PER2) gene that displayed such behavior. We compared Sanger and DNA‐seq sequencing for this locus and applied high‐resolution melt and MFOLD analyses to point to the DNA secondary structure that underlined this phenomenon. A synthetic system of oligonucleotides cloned into plasmids was used to further test the effect of such structures on sequencing. Our analyses indicated that a stem–loop structure capable of G–T pairing mediated the orientation bias by stabilizing this structure for specific alleles in heterozygous situations. We propose that this wobble‐like pairing hinders DNA polymerase passage on one strand while, on the complementary strand, the nonpaired A–C nucleotide counterparts allow unobstructed replication. Experimentation with synthetic amplicons that form similar stem–loop structures supported our hypothesis. We coined the term ‘replication diode’ for this effect and demonstrated that we can minimize it by lowering DNA and salt concentration. We also demonstrated that common genomic palindromes might induce the replication diode effect by applying bidirectional sequencing to an amplicon containing the palindrome within the human miRNA 1‐1 gene. Hence, to obtain reliable peak‐height ratios, bidirectional sequencing should be practiced at the lowest possible ionic strength whenever estimating CNPs. Further research is needed to determine whether the observed variable stem–loop structures affect PER2 regulation in vivo.

Unveiling genomic regions that underlie differences between Afec-Assaf sheep and its parental Awassi breed

BACKGROUND

Sheep production in Israel has improved by crossing the fat-tailed local Awassi breed with the East Friesian and later, with the Booroola Merino breed, which led to the formation of the highly prolific Afec-Assaf strain. This strain differs from its parental Awassi breed in morphological traits such as tail and horn size, coat pigmentation and wool characteristics, as well as in production, reproductive and health traits. To identify major genes associated with the formation of the Afec-Assaf strain, we genotyped 41 Awassi and 141 Afec-Assaf sheep using the Illumina Ovine SNP50 BeadChip array, and analyzed the results with PLINK and EMMAX software. The detected variable genomic regions that differed between Awassi and Afec-Assaf sheep (variable genomic regions; VGR) were compared to selection signatures that were reported in 48 published genome-wide association studies in sheep. Because the Afec-Assaf strain, but not the Awassi breed, carries the Booroola mutation, association analysis of BMPR1B used as the test gene was performed to evaluate the ability of this study to identify a VGR that includes such a major gene.

RESULTS

Of the 20 detected VGR, 12 were novel to this study. A ~7-Mb VGR was identified on Ovies aries chromosome OAR6 where the Booroola mutation is located. Similar to other studies, the most significant VGR was detected on OAR10, in a region that contains candidate genes affecting horn type (RXFP2), climate adaptation (ALOX5AP), fiber diameter (KATNAl1), coat pigmentation (FRY) and genes associated with fat distribution. The VGR on OAR2 included BNC2, which is also involved in controlling coat pigmentation in sheep. Six other VGR contained genes that were shown to be involved in coat pigmentation by analyzing their mammalian orthologues. Genes associated with fat distribution in humans, including GRB14 and COBLL1, were located in additional VGR. Sequencing DNA from Awassi and Afec-Assaf individuals revealed non-synonymous mutations in some of these candidate genes.

CONCLUSIONS

Our results highlight VGR that differentiate the Awassi breed from the Afec-Assaf strain, some of which may include genes that confer an advantage to Afec-Assaf and Assaf over Awassi sheep with respect to intensive sheep production under Mediterranean conditions.

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 (R²  = 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.

Identification of the Long-Sought Leptin in Chicken and Duck: Expression Pattern of the Highly GC-Rich Avian leptin Fits an Autocrine/Paracrine Rather Than Endocrine Function

More than 20 years after characterization of the key regulator of mammalian energy balance, leptin, we identified the leptin (LEP) genes of chicken (Gallus gallus) and duck (Anas platyrhynchos). The extreme guanine-cytosine content (∼70%), the location in a genomic region with low-complexity repetitive and palindromic sequence elements, the relatively low sequence conservation, and low level of expression have hampered the identification of these genes until now. In vitro-expressed chicken and duck leptins specifically activated signaling through the chicken leptin receptor in cell culture. In situ hybridization demonstrated expression of LEP mRNA in granular and Purkinje cells of the cerebellum, anterior pituitary, and in embryonic limb buds, somites, and branchial arches, suggesting roles in adult brain control of energy balance and during embryonic development. The expression patterns of LEP and the leptin receptor (LEPR) were explored in chicken, duck, and quail (Coturnix japonica) using RNA-sequencing experiments available in the Short Read Archive and by quantitative RT-PCR. In adipose tissue, LEP and LEPR were scarcely transcribed, and the expression level was not correlated to adiposity. Our identification of the leptin genes in chicken and duck genomes resolves a long lasting controversy regarding the existence of leptin genes in these species. This identification was confirmed by sequence and structural similarity, conserved exon-intron boundaries, detection in numerous genomic, and transcriptomic datasets and characterization by PCR, quantitative RT-PCR, in situ hybridization, and bioassays. Our results point to an autocrine/paracrine mode of action for bird leptin instead of being a circulating hormone as in mammals.

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⁻⁹). 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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-774) contains supplementary material, which is available to authorized users.

Discovery and characterization of the first genuine avian leptin gene in the rock dove (Columba livia)

Leptin, the key regulator of mammalian energy balance, has been at the center of a great controversy in avian biology for the last 15 years since initial reports of a putative leptin gene (LEP) in chickens. Here, we characterize a novel LEP in rock dove (Columba livia) with low similarity of the predicted protein sequence (30% identity, 47% similarity) to the human ortholog. Searching the Sequence-Read-Archive database revealed leptin transcripts, in the dove's liver, with 2 noncoding exons preceding 2 coding exons. This unusual 4-exon structure was validated by sequencing of a GC-rich product (76% GC, 721 bp) amplified from liver RNA by RT-PCR. Sequence alignment of the dove leptin with orthologous leptins indicated that it consists of a leader peptide (21 amino acids; aa) followed by the mature protein (160 aa), which has a putative structure typical of 4-helical-bundle cytokines except that it is 12 aa longer than human leptin. Extra residues (10 aa) were located within the loop between 2 5'-helices, interrupting the amino acid motif that is conserved in tetrapods and considered essential for activation of leptin receptor (LEPR) but not for receptor binding per se. Quantitative RT-PCR of 11 tissues showed highest (P < .05) expression of LEP in the dove's liver, whereas the dove LEPR peaked (P < .01) in the pituitary. Both genes were prominently expressed in the gonads and at lower levels in tissues involved in mammalian leptin signaling (adipose; hypothalamus). A bioassay based on activation of the chicken LEPR in vitro showed leptin activity in the dove's circulation, suggesting that dove LEP encodes an active protein, despite the interrupted loop motif. Providing tools to study energy-balance control at an evolutionary perspective, our original demonstration of leptin signaling in dove predicts a more ancient role of leptin in growth and reproduction in birds, rather than appetite control.

Quack leptin

BACKGROUND

A LEP transcript up-regulated in lungs of ducks (Anas platyrhynchos) infected by avian influenza A virus was recently described in the Nature Genetics manuscript that reported the duck genome. In vertebrates, LEP gene symbol is reserved for leptin, the key regulator of energy balance in mammals.

RESULTS

Launching an extensive search for this gene in the genome data that was submitted to the public databases along with duck genome manuscript and extending this search to all avian genomes in the whole-genome shotgun-sequencing database, we were able to report the first identification of coding sequences capable of encoding the full leptin protein precursor in wild birds. Gene structure, synteny and sequence-similarity (up to 54% identity and 68% similarity) to reptilian leptin evident in falcons (Falco peregrinus and cherrug), tits (Pseudopodoces humilis), finches (Taeniopygia guttata) and doves (Columba livia) confirmed that the bird leptin was a true ortholog of its mammalian form. Nevertheless, in duck, like other domestic fowls the LEP gene was not identifiable.

CONCLUSION

Lack of the LEP gene in poultry suggests that birds that have lost it are particularly suited to domestication. Identification of an intact avian gene for leptin in wild birds might explain in part the evolutionary conservation of its receptor in leptin-less fowls.

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.

Nonbactericidal secreted phospholipase A2s are potential anti-inflammatory factors in the mammary gland

The recent burst of duplication and divergence of the bovine PLA2G2D genes is considered typical of immune response genes, and it was recently shown that PLA2G2D is abundantly expressed in mouse leukocytes and acts as an immunosuppressive phospholipase. Analysis of 1,143 Holstein bulls indicated that the four common haplotypes spanning PLA2G2D display copy number variation ranging from 1 to 4 per haploid genome. Association of the fourth haplotype with negative total merit remained significant (P < 0.002) when corrected for population relatedness. We compared the lipase and bactericidal activities of bovine pancreatic PLA2G1B with human PLA2G2A and G2D and bovine PLA2G2D1 and G2D4 proteins, which had been subcloned, expressed, and refolded by us, and the impact of point mutations in the calcium binding site was investigated. All tested phospholipases were ineffective bactericides of Escherichia coli isolated from bovine mastitis. However, in lactating mice treated with E. coli or lipopolysaccharide (LPS), intramammary injection of bovine PLA2G1B relieved visual and histological inflammation and reduced blood levels of infiltrating lactose. Further studies are warranted to determine whether the observed anti-inflammatory effect involves competitive binding of the receptor Pla2r1 which may mimic the LPS resistance effect in Pla2r1-deficient mice.

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.

Pegylated leptin antagonist with strong orexigenic activity in mice is not effective in chickens

A chicken gene orthologous to human leptin receptor (LEPR) has been characterized and found to be active in leptin signaling in vitro in response to a variety of recombinant leptins and leptin-containing blood samples. However, the endogenous ligand of chicken LEPR (cLEPR) - the putative chicken leptin - has been reported by us and others to be undetectable at the DNA, mRNA, protein and activity levels. These reports have raised questions as to cLEPR's role. Here we analyzed the effects of a pegylated superactive mouse leptin antagonist (PEG-SMLA) in chicken. We showed that the leptin antagonist efficiently and specifically blocks leptin signaling through the cLEPR in vitro. The effect of the leptin antagonist was then studied in vivo by daily administration of 10 mg/kg for 10 consecutive days to White Leghorn female chickens (G. gallus), at the age of two weeks. Despite the efficient attenuation of the cLEPR in vitro, no effect was observed on body weight, feed intake, feed efficiency or fat accumulation in the treated birds. Since similar treatment in rodents leads to a highly pronounced increase in appetite and body weight that are observed from the first day of treatment, it is concluded that the cLEPR is not implicated in the control of appetite or adipose homeostasis in chickens.

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.

SNPmplexViewer--toward a cost-effective traceability system

Background

Beef traceability has become mandatory in many regions of the world and is typically achieved through the use of unique numerical codes on ear tags and animal passports. DNA-based traceability uses the animal's own DNA code to identify it and the products derived from it. Using SNaPshot, a primer-extension-based method, a multiplex of 25 SNPs in a single reaction has been practiced for reducing the expense of genotyping a panel of SNPs useful for identity control.

Findings

To further decrease SNaPshot's cost, we introduced the Perl script SNPmplexViewer, which facilitates the analysis of trace files for reactions performed without the use of fluorescent size standards. SNPmplexViewer automatically aligns reference and target trace electropherograms, run with and without fluorescent size standards, respectively. SNPmplexViewer produces a modified target trace file containing a normalised trace in which the reference size standards are embedded. SNPmplexViewer also outputs aligned images of the two electropherograms together with a difference profile.

Conclusions

Modified trace files generated by SNPmplexViewer enable genotyping of SnaPshot reactions performed without fluorescent size standards, using common fragment-sizing software packages. SNPmplexViewer's normalised output may also improve the genotyping software's performance. Thus, SNPmplexViewer is a general free tool enabling the reduction of SNaPshot's cost as well as the fast viewing and comparing of trace electropherograms for fragment analysis. SNPmplexViewer is available at http://cowry.agri.huji.ac.il/cgi-bin/SNPmplexViewer.cgi.

Biochemical and in vitro biological significance of natural sequence variation in the ovine leptin gene

The hormone leptin is involved in diverse biological processes, including regulation of food intake, body-weight homeostasis and energy balance. Sequence variation in the bovine leptin gene has been found to be associated with variations in carcass fat content and average daily gain, as well as in milk yield, milk somatic cell count and several traits governing reproduction. We sequenced genomic DNA and cDNA samples of individuals from three divergent sheep breeds and revealed synonymous as well as novel non-synonymous allelic variation at the third exon of the ovine leptin gene (oLEP) as compared to the sequence published at Accession No. U84247 (reference sequence). In addition, two alternatively spliced oLEP transcripts were found in the abdominal fat tissue. The biochemical and the in vitro biological significance of the sequence variation in the oLEP was examined by generating recombinant oLEP-protein variants namely: p.Q28del, p.N78S, p.R84Q, p.P99Q, p.V123L and p.R138Q, carrying the corresponding sequence variation. Surface plasmon resonance experiments revealed, in most cases, reduced affinity of the oLEP protein variants examined, to human leptin-binding domain (hLBD), relative to the reference variant, being 0.75, 0.60, 0.60, 0.89, 0.92 and 1.03, respectively. In competitive binding assays between biotinylated oLEP and the recombinant leptin protein variants, p.N78S and p.R84Q variants exhibited the lowest affinity to hLBD (0.18 and 0.41, respectively) as compared to the reference hormone. We then tested the protein variants' ability to induce proliferation in Baf-3 cells stably expressing the long form of the human leptin receptor: significant differences in proliferative activity were only found for p.N78S (1.8-fold higher) and p.R138Q (4.2-fold lower) relative to the reference oLEP variant.

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.

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.