Genome-wide association studies for the concentrations of insulin, triiodothyronine, and thyroxine in Chinese Holstein cattle

Genome wide association study of thyroid hormone levels following challenge with porcine reproductive and respiratory syndrome virus

Introduction: Porcine reproductive and respiratory syndrome virus (PRRSV) causes respiratory disease in piglets and reproductive disease in sows. Piglet and fetal serum thyroid hormone (i.e., T3 and T4) levels decrease rapidly in response to Porcine reproductive and respiratory syndrome virus infection. However, the genetic control of T3 and T4 levels during infection is not completely understood. Our objective was to estimate genetic parameters and identify quantitative trait loci (QTL) for absolute T3 and/or T4 levels of piglets and fetuses challenged with Porcine reproductive and respiratory syndrome virus. Methods: Sera from 5-week-old pigs (N = 1792) at 11 days post inoculation (DPI) with Porcine reproductive and respiratory syndrome virus were assayed for T3 levels (piglet_T3). Sera from fetuses (N = 1,267) at 12 or 21 days post maternal inoculation (DPMI) with Porcine reproductive and respiratory syndrome virus of sows (N = 145) in late gestation were assayed for T3 (fetal_T3) and T4 (fetal_T4) levels. Animals were genotyped using 60 K Illumina or 650 K Affymetrix single nucleotide polymorphism (SNP) panels. Heritabilities, phenotypic correlations, and genetic correlations were estimated using ASREML; genome wide association studies were performed for each trait separately using Julia for Whole-genome Analysis Software (JWAS). Results: All three traits were low to moderately heritable (10%-16%). Phenotypic and genetic correlations of piglet_T3 levels with weight gain (0-42 DPI) were 0.26 ± 0.03 and 0.67 ± 0.14, respectively. Nine significant quantitative trait loci were identified for piglet_T3, on Sus scrofa chromosomes (SSC) 3, 4, 5, 6, 7, 14, 15, and 17, and collectively explaining 30% of the genetic variation (GV), with the largest quantitative trait loci identified on SSC5, explaining 15% of the genetic variation. Three significant quantitative trait loci were identified for fetal_T3 on SSC1 and SSC4, which collectively explained 10% of the genetic variation. Five significant quantitative trait loci were identified for fetal_T4 on SSC1, 6, 10, 13, and 15, which collectively explained 14% of the genetic variation. Several putative immune-related candidate genes were identified, including CD247, IRF8, and MAPK8. Discussion: Thyroid hormone levels following Porcine reproductive and respiratory syndrome virus infection were heritable and had positive genetic correlations with growth rate. Multiple quantitative trait loci with moderate effects were identified for T3 and T4 levels during challenge with Porcine reproductive and respiratory syndrome virus and candidate genes were identified, including several immune-related genes. These results advance our understanding of growth effects of both piglet and fetal response to Porcine reproductive and respiratory syndrome virus infection, revealing factors associated with genomic control of host resilience.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Genome-Wide Association Analysis Reveals Novel Loci Related with Visual Score Traits in Nellore Cattle Raised in Pasture-Based Systems

Body conformation traits assessed based on visual scores are widely used in Zebu cattle breeding programs. The aim of this study was to identify genomic regions and biological pathways associated with body conformation (CONF), finishing precocity (PREC), and muscling (MUSC) in Nellore cattle. The measurements based on visual scores were collected in 20,807 animals raised in pasture-based systems in Brazil. In addition, 2775 animals were genotyped using a 35 K SNP chip, which contained 31,737 single nucleotide polymorphisms after quality control. Single-step GWAS was performed using the BLUPF90 software while candidate genes were identified based on the Ensembl Genes 69. PANTHER and REVIGO platforms were used to identify key biological pathways and STRING to create gene networks. Novel candidate genes were revealed associated with CONF, including ALDH9A1, RXRG, RAB2A, and CYP7A1, involved in lipid metabolism. The genes associated with PREC were ELOVL5, PID1, DNER, TRIP12, and PLCB4, which are related to the synthesis of long-chain fatty acids, lipid metabolism, and muscle differentiation. For MUSC, the most important genes associated with muscle development were SEMA6A, TIAM2, UNC5A, and UIMC1. The polymorphisms identified in this study can be incorporated in commercial genotyping panels to improve the accuracy of genomic evaluations for visual scores in beef cattle.

Diversity of Mitochondrial DNA Haplogroups and Their Association with Bovine Antral Follicle Count

Maternal origins based on the bovine mitochondrial D-loop region are proven to have two main origins: Bos taurus and Bos indicus. To examine the association between the maternal origins of bovine and reproductive traits, the complete mitochondrial D-loop region sequences from 501 Chinese Holstein cows and 94 individuals of other breeds were analyzed. Based on the results obtained from the haplotype analysis, 260 SNPs (single nucleotide polymorphism), 32 indels (insertion/deletion), and 219 haplotypes were identified. Moreover, the nucleotide diversity (π) and haplotype diversity (Hd) were 0.024 ± 0.001 and 0.9794 ± 0.003, respectively, indicating the abundance of genetic resources in Chinese Holstein cows. The results of the median-joining network analysis showed two haplogroups (HG, including HG1 and HG2) that diverged in genetic distance. Furthermore, the two haplogroups were significantly (p < 0.05) correlated with the antral follicle (diameter ≥ 8 mm) count, and HG1 individuals had more antral follicles than HG2 individuals, suggesting that these different genetic variants between HG1 and HG2 correlate with reproductive traits. The construction of a neighbor-joining phylogenetic tree and principal component analysis also revealed two main clades (HG1 and HG2) with different maternal origins: Bos indicus and Bos taurus, respectively. Therefore, HG1 originating from the maternal ancestors of Bos indicus may have a greater reproductive performance, and potential genetic variants discovered may promote the breeding process in the cattle industry.

Genetic parameters and genome-wide association for milk production traits and somatic cell score in different lactation stages of Shanghai Holstein population

The aim of this study was to investigate the genetic parameters and genetic architectures of six milk production traits in the Shanghai Holstein population. The data used to estimate the genetic parameters consisted of 1,968,589 test-day records for 305,031 primiparous cows. Among the cows with phenotypes, 3,016 cows were genotyped with Illumina Bovine SNP50K BeadChip, GeneSeek Bovine 50K BeadChip, GeneSeek Bovine LD BeadChip v4, GeneSeek Bovine 150K BeadChip, or low-depth whole-genome sequencing. A genome-wide association study was performed to identify quantitative trait loci and genes associated with milk production traits in the Shanghai Holstein population using genotypes imputed to whole-genome sequences and both fixed and random model circulating probability unification and a mixed linear model with rMVP software. Estimated heritabilities (h2) varied from 0.04 to 0.14 for somatic cell score (SCS), 0.07 to 0.22 for fat percentage (FP), 0.09 to 0.27 for milk yield (MY), 0.06 to 0.23 for fat yield (FY), 0.09 to 0.26 for protein yield (PY), and 0.07 to 0.35 for protein percentage (PP), respectively. Within lactation, genetic correlations for SCS, FP, MY, FY, PY, and PP at different stages of lactation estimated in random regression model were ranged from -0.02 to 0.99, 0.18 to 0.99, 0.04 to 0.99, 0.04 to 0.99, 0.01 to 0.99, and 0.33 to 0.99, respectively. The genetic correlations were highest between adjacent DIM but decreased as DIM got further apart. Candidate genes included those related to production traits (DGAT1, MGST1, PTK2, and SCRIB), disease-related (LY6K, COL22A1, TECPR2, and PLCB1), heat stress–related (ITGA9, NDST4, TECPR2, and HSF1), and reproduction-related (7SK and DOCK2) genes. This study has shown that there are differences in the genetic mechanisms of milk production traits at different stages of lactation. Therefore, it is necessary to conduct research on milk production traits at different stages of lactation as different traits. Our results can also provide a theoretical basis for subsequent molecular breeding, especially for the novel genetic loci.

Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo

Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (F_ROH) ranged from 0.0201 to 0.0371. Pearson correlations between F_ROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.

The Association of an SNP in the EXOC4 Gene and Reproductive Traits Suggests Its Use as a Breeding Marker in Pigs

In mammals, the exocyst complex component 4 (EXOC4) gene has often been reported to be involved in vesicle transport. The SNP rs81471943 (C/T) is located in the intron of porcine EXOC4, while six quantitative trait loci (QTL) within 5-10 Mb around EXOC4 are associated with ovary weight, teat number, total offspring born alive, and corpus luteum number. However, the molecular mechanisms between EXOC4 and the reproductive performance of pigs remains to be elucidated. In this study, rs81471943 was genotyped from a total of 994 Duroc sows, and the genotype and allele frequency of SNP rs81471943 (C/T) were statistically analyzed. Then, the associations between SNP rs81471943 and four reproductive traits, including number of piglets born alive (NBA), litter weight at birth (LWB), number of piglets weaned (NW), and litter weight at weaning (LWW), were determined. Sanger sequencing and PCR restriction fragment length polymorphism (PCR-RFLP) were utilized to identify the rs81471943 genotype. We found that the genotype frequency of CC was significantly higher than that of CT and TT, and CC was the most frequent genotype for NBA, LWB, NW, and LWW. Moreover, 5'-deletion and luciferase assays identified a positive transcription regulatory element in the EXOC4 promoter. After exploring the EXOC4 promoter, SNP -1781G/A linked with SNP rs81471943 (C/T) were identified by analysis of the transcription activity of the haplotypes, and SNP -1781 G/A may influence the potential binding of P53, E26 transformation specific sequence -like 1 transcription factor (ELK1), and myeloid zinc finger 1 (MZF1). These findings provide useful information for identifying a molecular marker of EXOC4-assisted selection in pig breeding.