Genomic prediction based on a joint reference population for the Xinjiang Brown cattle.

Introduction: Xinjiang Brown cattle constitute the largest breed of cattle in Xinjiang. Therefore, it is crucial to establish a genomic evaluation system, especially for those with low levels of breed improvement. Methods: This study aimed to establish a cross breed joint reference population by analyzing the genetic structure of 485 Xinjiang Brown cattle and 2,633 Chinese Holstein cattle (Illumina GeneSeek GGP bovine 150 K chip). The Bayes method single-step genome-wide best linear unbiased prediction was used to conduct a genomic evaluation of the joint reference population for the milk traits of Xinjiang Brown cattle. The reference population of Chinese Holstein cattle was randomly divided into groups to construct the joint reference population. By comparing the prediction accuracy, estimation bias, and inflation coefficient of the validation population, the optimal number of joint reference populations was determined. Results and Discussion: The results indicated a distinct genetic structure difference between the two breeds of adult cows, and both breeds should be considered when constructing multi-breed joint reference and validation populations. The reliability range of genome prediction of milk traits in the joint reference population was 0.142-0.465. Initially, it was determined that the inclusion of 600 and 900 Chinese Holstein cattle in the joint reference population positively impacted the genomic prediction of Xinjiang Brown cattle to certain extent. It was feasible to incorporate the Chinese Holstein into Xinjiang Brown cattle population to form a joint reference population for multi-breed genomic evaluation. However, for different Xinjiang Brown cattle populations, a fixed number of Chinese Holstein cattle cannot be directly added during multi-breed genomic selection. Pre-evaluation analysis based on the genetic structure, kinship, and other factors of the current population is required to ensure the authenticity and reliability of genomic predictions and improve estimation accuracy.

Genome-wide association analysis of resistance to frogeye leaf spot China race 7 in soybean based on high-throughput sequencing.

KEY MESSAGE: FLS is a disease that causes severe yield reduction in soybean. In this study, four genes (Glyma.16G176800, Glyma.16G177300, Glyma.16G177400 and Glyma.16G182300) were tentatively confirmed to play an important role in the resistance of soybean to FLS race 7. Frogeye leaf spot (FLS) causes severe yield loss in soybean and has been found in several countries worldwide. Therefore, it is necessary to select and utilize FLS-resistant varieties for the management of FLS. In the present study, 335 representative soybean materials were assessed for partial resistance to FLS race 7. Quantitative trait nucleotide (QTN) and FLS race 7 candidate genes were identified using genome-wide association analysis (GWAS) based on a site-specific amplified fragment sequencing (SLAF-seq) approach. A total of 23,156 single-nucleotide polymorphisms (SNPs) were used to evaluate the level of linkage disequilibrium with a minor allele frequency ≥ 5 and deletion data < 3%. These SNPs covered about 947.01 MBP, nearly 86.09% of the entire soybean genome. In addition, a compressed mixed linear model was utilized to identify association signals for partial resistance to FLS race 7. A total of 15 QTNs associated with resistance were found to be novel for FLS race 7 resistance. A total of 217 candidate genes located in the 200-kb genomic region of these peak SNPs were identified. Based on gene association analysis, qRT-PCR, haplotype analysis and virus-induced gene silencing (VIGS) systems were used to further verify candidate genes Glyma.16G176800, Glyma.16G177300, Glyma.16G177400 and Glyma.16G182300. This indicates that these four candidate genes may participate in FLS race 7 resistance responses.

Genome-wide association studies reveal novel QTLs, QTL-by-environment interactions and their candidate genes for tocopherol content in soybean seed.

Genome-wide association studies (GWAS) is an efficient method to detect quantitative trait locus (QTL), and has dissected many complex traits in soybean [Glycine max (L.) Merr.]. Although these results have undoubtedly played a far-reaching role in the study of soybean biology, environmental interactions for complex traits in traditional GWAS models are frequently overlooked. Recently, a new GWAS model, 3VmrMLM, was established to identify QTLs and QTL-by-environment interactions (QEIs) for complex traits. In this study, the GLM, MLM, CMLM, FarmCPU, BLINK, and 3VmrMLM models were used to identify QTLs and QEIs for tocopherol (Toc) content in soybean seed, including δ-Tocotrienol (δ-Toc) content, γ-Tocotrienol (γ-Toc) content, α-Tocopherol (α-Toc) content, and total Tocopherol (T-Toc) content. As a result, 101 QTLs were detected by the above methods in single-environment analysis, and 57 QTLs and 13 QEIs were detected by 3VmrMLM in multi-environment analysis. Among these QTLs, some QTLs (Group I) were repeatedly detected three times or by at least two models, and some QTLs (Group II) were repeatedly detected only by 3VmrMLM. In the two Groups, 3VmrMLM was able to correctly detect all known QTLs in group I, while good results were achieved in Group II, for example, 8 novel QTLs were detected in Group II. In addition, comparative genomic analysis revealed that the proportion of Glyma_max specific genes near QEIs was higher, in other words, these QEIs nearby genes are more susceptible to environmental influences. Finally, around the 8 novel QTLs, 11 important candidate genes were identified using haplotype, and validated by RNA-Seq data and qRT-PCR analysis. In summary, we used phenotypic data of Toc content in soybean, and tested the accuracy and reliability of 3VmrMLM, and then revealed novel QTLs, QEIs and candidate genes for these traits. Hence, the 3VmrMLM model has broad prospects and potential for analyzing the genetic structure of complex quantitative traits in soybean.

Genomic Selection for Milk Production Traits in Xinjiang Brown Cattle.

One-step genomic selection is a method for improving the reliability of the breeding value estimation. This study aimed to compare the reliability of pedigree-based best linear unbiased prediction (PBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP), single-trait and multitrait models, and the restricted maximum likelihood (REML) and Bayesian methods. Data were collected from the production performance records of 2207 Xinjiang Brown cattle in Xinjiang from 1983 to 2018. A cross test was designed to calculate the genetic parameters and reliability of the breeding value of 305 daily milk yield (305 dMY), milk fat yield (MFY), milk protein yield (MPY), and somatic cell score (SCS) of Xinjiang Brown cattle. The heritability of 305 dMY, MFY, MPY, and SCS estimated using the REML and Bayesian multitrait models was approximately 0.39 (0.02), 0.40 (0.03), 0.49 (0.02), and 0.07 (0.02), respectively. The heritability and estimated breeding value (EBV) and the reliability of milk production traits of these cattle calculated based on PBLUP and ssGBLUP using the multitrait model REML and Bayesian methods were higher than those of the single-trait model REML method; the ssGBLUP method was significantly better than the PBLUP method. The reliability of the estimated breeding value can be improved from 0.9% to 3.6%, and the reliability of the genomic estimated breeding value (GEBV) for the genotyped population can reach 83%. Therefore, the genetic evaluation of the multitrait model is better than that of the single-trait model. Thus, genomic selection can be applied to small population varieties such as Xinjiang Brown cattle, in improving the reliability of the genomic estimated breeding value.

HandyCNV: Standardized Summary, Annotation, Comparison, and Visualization of Copy Number Variant, Copy Number Variation Region, and Runs of Homozygosity.

Detection of CNVs (copy number variants) and ROH (runs of homozygosity) from SNP (single nucleotide polymorphism) genotyping data is often required in genomic studies. The post-analysis of CNV and ROH generally involves many steps, potentially across multiple computing platforms, which requires the researchers to be familiar with many different tools. In order to get around this problem and improve research efficiency, we present an R package that integrates the summarization, annotation, map conversion, comparison and visualization functions involved in studies of CNV and ROH. This one-stop post-analysis system is standardized, comprehensive, reproducible, timesaving, and user-friendly for researchers in humans and most diploid livestock species.

Discovering Copy Number Variation in Dual-Purpose XinJiang Brown Cattle.

Copy number variants (CNVs), which are a class of structural variant, can be important in relating genomic variation to phenotype. The primary aims of this study were to discover the common CNV regions (CNVRs) in the dual-purpose XinJiang-Brown cattle population and to detect differences between CNVs inferred using the ARS-UCD 1.2 (ARS) or the UMD 3.1 (UMD) genome assemblies based on the 150K SNP (Single Nucleotide Polymorphisms) Chip. PennCNV and CNVPartition methods were applied to calculate the deviation of the standardized signal intensity of SNPs markers to detect CNV status. Following the discovery of CNVs, we used the R package HandyCNV to generate and visualize CNVRs, compare CNVs and CNVRs between genome assemblies, and identify consensus genes using annotation resources. We identified 38 consensus CNVRs using the ARS assembly with 1.95% whole genome coverage, and 33 consensus CNVRs using the UMD assembly with 1.46% whole genome coverage using PennCNV and CNVPartition. We identified 37 genes that intersected 13 common CNVs (>5% frequency), these included functionally interesting genes such as GBP4 for which an increased copy number has been negatively associated with cattle stature, and the BoLA gene family which has been linked to the immune response and adaption of cattle. The ARS map file of the GGP Bovine 150K Bead Chip maps the genomic position of more SNPs with increased accuracy compared to the UMD map file. Comparison of the CNVRs identified between the two reference assemblies suggests the newly released ARS reference assembly is better for CNV detection. In spite of this, different CNV detection methods can complement each other to generate a larger number of CNVRs than using a single approach and can highlight more genes of interest.

GWAS-Based Identification of New Loci for Milk Yield, Fat, and Protein in Holstein Cattle.

High-yield and high-quality of milk are the primary goals of dairy production. Understanding the genetic architecture underlying these milk-related traits is beneficial so that genetic variants can be targeted toward the genetic improvement. In this study, we measured five milk production and quality traits in Holstein cattle population from China. These traits included milk yield, fat, and protein. We used the estimated breeding values as dependent variables to conduct the genome-wide association studies (GWAS). Breeding values were estimated through pedigree relationships by using a linear mixed model. Genotyping was carried out on the individuals with phenotypes by using the Illumina BovineSNP150 BeadChip. The association analyses were conducted by using the fixed and random model Circulating Probability Unification (FarmCPU) method. A total of ten single-nucleotide polymorphisms (SNPs) were detected above the genome-wide significant threshold (p < 4.0 × 10-7), including six located in previously reported quantitative traits locus (QTL) regions. We found eight candidate genes within distances of 120 kb upstream or downstream to the associated SNPs. The study not only identified the effect of DGAT1 gene on milk fat and protein, but also discovered novel genetic loci and candidate genes related to milk traits. These novel genetic loci would be an important basis for molecular breeding in dairy cattle.

Genome-wide association study of milk and reproductive traits in dual-purpose Xinjiang Brown cattle.

BACKGROUND: Dual-purpose cattle are more adaptive to environmental challenges than single-purpose dairy or beef cattle. Balance among milk, reproductive, and mastitis resistance traits in breeding programs is therefore more critical for dual-purpose cattle to increase net income and maintain well-being. With dual-purpose Xinjiang Brown cattle adapted to the Xinjiang Region in northwestern China, we conducted genome-wide association studies (GWAS) to dissect the genetic architecture related to milk, reproductive, and mastitis resistance traits. Phenotypic data were collected for 2410 individuals measured during 1995-2017. By adding another 445 ancestors, a total of 2855 related individuals were used to derive estimated breeding values for all individuals, including the 2410 individuals with phenotypes. Among phenotyped individuals, we genotyped 403 cows with the Illumina 150 K Bovine BeadChip. RESULTS: GWAS were conducted with the FarmCPU (Fixed and random model circulating probability unification) method. We identified 12 markers significantly associated with six of the 10 traits under the threshold of 5% after a Bonferroni multiple test correction. Seven of these SNPs were in QTL regions previously identified to be associated with related traits. One identified SNP, BovineHD1600006691, was significantly associated with both age at first service and age at first calving. This SNP directly overlapped a QTL previously reported to be associated with calving ease. Within 160 Kb upstream and downstream of each significant SNP identified, we speculated candidate genes based on functionality. Four of the SNPs were located within four candidate genes, including CDH2, which is linked to milk fat percentage, and GABRG2, which is associated with milk protein yield. CONCLUSIONS: These findings are beneficial not only for breeding through marker-assisted selection, but also for genome editing underlying the related traits to enhance the overall performance of dual-purpose cattle.

The study of blood transcriptome profiles in Holstein cows with miscarriage during peri-implantation.

OBJECTIVE: In this study, the transcriptome profile of cow experiencing miscarriage during peri-implantation was investigated. METHODS: Total transcriptomes were checked by RNA sequencing, and the analyzed by bioinformatics methods, the differentially expressed genes (DEGs) were analysed with hierarchical clustering and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis. RESULTS: The results suggested that serum progesterone levels were significantly decreased in cows that miscarried as compared to the pregnant cows at 18, 21, 33, 39, and 51 days after artificial insemination. The RNA sequencing results suggested that 32, 176, 5, 10, and 2 DEGs were identified in the pregnant cows and miscarried cows at 18, 21, 33, 39, and 51 d after artificial insemination. And 15, 101, 1, 2, and 2 DEGs were upregulated, and 17, 74, 4, and 8 DEGs were downregulated in the cows in the pregnant and miscarriage groups, respectively at 18, 21, 33, and 39, but no gene was downregulated at 51 d after artificial insemination. These DEGs were distributed to 13, 20, 3, 6, and 20 pathways, and some pathway essential for pregnancy, such as cell adhesion molecules, tumor necrosis factor signaling pathway and PI3K-Akt signaling pathway. CONCLUSION: This analysis has identified several genes and related pathways crucial for pregnancy and miscarriage in cows, as well as these genes supply molecular markers to predict the miscarriage in cows.