SNPmplexViewer--toward a cost-effective traceability system

<p>Abstract</p> <p>Background</p> <p>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 <it>SNaPshot</it>, 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.</p> <p>Findings</p> <p>To further decrease <it>SNaPshot</it>'s cost, we introduced the Perl script <it>SNPmplexViewer</it>, which facilitates the analysis of trace files for reactions performed without the use of fluorescent size standards. <it>SNPmplexViewer </it>automatically aligns reference and target trace electropherograms, run with and without fluorescent size standards, respectively. <it>SNPmplexViewer </it>produces a modified target trace file containing a normalised trace in which the reference size standards are embedded. <it>SNPmplexViewer </it>also outputs aligned images of the two electropherograms together with a difference profile.</p> <p>Conclusions</p> <p>Modified trace files generated by <it>SNPmplexViewer </it>enable genotyping of <it>SnaPshot </it>reactions performed without fluorescent size standards, using common fragment-sizing software packages. <it>SNPmplexViewer</it>'s normalised output may also improve the genotyping software's performance. Thus, <it>SNPmplexViewer </it>is a general free tool enabling the reduction of <it>SNaPshot</it>'s cost as well as the fast viewing and comparing of trace electropherograms for fragment analysis. <it>SNPmplexViewer </it>is available at <url>http://cowry.agri.huji.ac.il/cgi-bin/SNPmplexViewer.cgi</url>.</p

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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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 <it>PennCNV </it>"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.</p> <p>Conclusions</p> <p>Expansion of HWE and <it>PennCNV </it>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.</p

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

<p>Abstract</p> <p>Background</p> <p>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.</p> <p>Results</p> <p>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.</p> <p>Conclusion</p> <p>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 (<it>ABCG2</it>, <it>DGAT1</it>, <it>GDF8, IGF2</it>) 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.</p

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.The study was supported by the Israel Academy of Sciences grants no. 876/ 14 and 1294/17, and Chief Scientist of the Israeli Ministry of Agriculture 0469/14 (to MFE and ES)

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

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.https://doi.org/10.1186/s12863-019-0773-

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