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

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 recombinant fungal defensin-like peptide-P2 combats Streptococcus dysgalactiae and biofilms

Streptococcus dysgalactiae, considered one of the main pathogens that causes bovine mastitis, is a serious threat to humans and animals. However, the excessive use of antibiotics and the characteristic of S. dysgalactiae forming biofilms in mastitic teat canal have serious clinical implications. In this study, in vivo and in vitro multiple mechanisms of action of P2, a mutant of fungal defensin plectasin, against S. dysgalactiae were systematically and comprehensively investigated for the first time. P2 showed potent antibacterial activity against S. dysgalactiae (minimum inhibitory concentration, MIC = 0.23-0.46 μM) and rapid bactericidal action by 3.0 lg units reduction in 2-4 h. No resistant mutants appeared after 30-d serial passage of S. dysgalactiae in the presence of P2. The results of electron microscopy and flow cytometer showed that P2 induced membrane damage of S. dysgalactiae, causing the leakage of cellular content and eventually cell death. Besides, P2 effectively inhibited early biofilm formation, eradicated mature biofilms, and killed 99.9% persisters which were resistant to 100 × MIC vancomycin; and confocal laser scanning microscopy (CLSM) also revealed the potent antibacterial and antibiofilm activity of P2 (the thickness of biofilm reduced from 18.82 to 7.94 μm). The in vivo therapeutic effect of P2 in mouse mastitis model showed that it decreased the number of mammary bacteria and alleviated breast inflammation by regulating cytokines and inhibiting bacterial proliferation, which were superior to vancomycin. These data indicated that P2 maybe a potential candidate peptide for mastitis treatment of S. dysgalactiae infections. KEY POINTS: •P2 showed potential in vitro antibacterial characteristics towards S. dysgalactiae. •P2 eradicated biofilms, killed persisters, and induced cell death of S. dysgalactiae. •P2 could effectively protect mice from S. dysgalactiae infection in gland.

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.

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.

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.

Exogenous phospholipase A2 affects inflammatory gene expression in primary bovine mammary epithelial cells

This Research Communication addresses the hypothesis that exogenously administered phospholipase A2 (PLA2) affects the inflammatory responses of bovine mammary epithelial cells (bMEC) in vitro with the aim of providing preliminary justification of investigation into the uses of exogenously administered PLA2 to manage or treat bovine mastitis. Primary bMEC lines from 11 lactating Holstein dairy cows were established and the expression of 14 pro-inflammatory genes compared under unchallenged and lipopolysaccharide (LPS)-challenged conditions, with and without concurrent treatment with bovine pancreatic PLA2G1B, a secreted form of PLA2. No differences in the expression of these genes were noted between PLA2-treated and untreated bMEC under unchallenged conditions. Following LPS challenge, untreated bMEC exhibited significant downregulation of CXCL8, IL1B, CCL20, and CXCL1. In contrast, PLA2-treated bMEC exhibited significant downregulation of IL1B and CCL20 only. These findings indicate that exogenous PLA2 affects the expression of some pro-inflammatory factors in immune-stimulated bMEC, but does not influence the constitutive expression of these factors. Further investigation of the influence of exogenous PLA2 in the bovine mammary gland is justified.

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.

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.

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.

Genomic and Phenomic Study of Mammary Pathogenic Escherichia coli

Escherichia coli is a major etiological agent of intra-mammary infections (IMI) in cows, leading to acute mastitis and causing great economic losses in dairy production worldwide. Particular strains cause persistent IMI, leading to recurrent mastitis. Virulence factors of mammary pathogenic E. coli (MPEC) involved pathogenesis of mastitis as well as those differentiating strains causing acute or persistent mastitis are largely unknown. This study aimed to identify virulence markers in MPEC through whole genome and phenome comparative analysis. MPEC strains causing acute (VL2874 and P4) or persistent (VL2732) mastitis were compared to an environmental strain (K71) and to the genomes of strains representing different E. coli pathotypes. Intra-mammary challenge in mice confirmed experimentally that the strains studied here have different pathogenic potential, and that the environmental strain K71 is non-pathogenic in the mammary gland. Analysis of whole genome sequences and predicted proteomes revealed high similarity among MPEC, whereas MPEC significantly differed from the non-mammary pathogenic strain K71, and from E. coli genomes from other pathotypes. Functional features identified in MPEC genomes and lacking in the non-mammary pathogenic strain were associated with synthesis of lipopolysaccharide and other membrane antigens, ferric-dicitrate iron acquisition and sugars metabolism. Features associated with cytotoxicity or intra-cellular survival were found specifically in the genomes of strains from severe and acute (VL2874) or persistent (VL2732) mastitis, respectively. MPEC genomes were relatively similar to strain K-12, which was subsequently shown here to be possibly pathogenic in the mammary gland. Phenome analysis showed that the persistent MPEC was the most versatile in terms of nutrients metabolized and acute MPEC the least. Among phenotypes unique to MPEC compared to the non-mammary pathogenic strain were uric acid and D-serine metabolism. This study reveals virulence factors and phenotypic characteristics of MPEC that may play a role in pathogenesis of E. coli mastitis.

Detection of copy number variations and their effects in Chinese bulls

Background

Copy number variations (CNVs) are a main source of genomic structural variations underlying animal evolution and production traits. Here, with one pure-blooded Angus bull as reference, we describe a genome-wide analysis of CNVs based on comparative genomic hybridization arrays in 29 Chinese domesticated bulls and examined their effects on gene expression and cattle growth traits.

Results

We identified 486 copy number variable regions (CNVRs), covering 2.45% of the bovine genome, in 24 taurine (Bos taurus), together with 161 ones in 2 yaks (Bos grunniens) and 163 ones in 3 buffaloes (Bubalus bubalis). Totally, we discovered 605 integrated CNVRs, with more “loss” events than both “gain” and “both” ones, and clearly clustered them into three cattle groups. Interestingly, we confirmed their uneven distributions across chromosomes, and the differences of mitochondrion DNA copy number (gain: taurine, loss: yak & buffalo). Furthermore, we confirmed approximately 41.8% (253/605) and 70.6% (427/605) CNVRs span cattle genes and quantitative trait loci (QTLs), respectively. Finally, we confirmed 6 CNVRs in 9 chosen ones by using quantitative PCR, and further demonstrated that CNVR22 had significantly negative effects on expression of PLA2G2D gene, and both CNVR22 and CNVR310 were associated with body measurements in Chinese cattle, suggesting their key effects on gene expression and cattle traits.

Conclusions

The results advanced our understanding of CNV as an important genomic structural variation in taurine, yak and buffalo. This study provides a highly valuable resource for Chinese cattle’s evolution and breeding researches.

Electronic supplementary material

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

UVB-induced gene expression in the skin of Xiphophorus maculatus Jp 163 B

Xiphophorus fish and interspecies hybrids represent long-standing models to study the genetics underlying spontaneous and induced tumorigenesis. The recent release of the Xiphophorus maculatus genome sequence will allow global genetic regulation studies of genes involved in the inherited susceptibility to UVB-induced melanoma within select backcross hybrids. As a first step toward this goal, we report results of an RNA-Seq approach to identify genes and pathways showing modulated transcription within the skin of X. maculatus Jp 163 B upon UVB exposure. X. maculatus Jp 163 B were exposed to various doses of UVB followed by RNA-Seq analysis at each dose to investigate overall gene expression in each sample. A total of 357 genes with a minimum expression change of 4-fold (p-adj<0.05) were identified as responsive to UVB. The molecular genetic response of Xiphophorus skin to UVB exposure permitted assessment of; (1) the basal expression level of each transcript for each skin sample, (2) the changes in expression levels for each gene in the transcriptome upon exposure to increasing doses of UVB, and (3) clusters of genes that exhibit similar patterns of change in expression upon UVB exposure. These data provide a foundation for understanding the molecular genetic response of fish skin to UVB exposure.