Genome-wide association study of copy number variations with shank traits in a F2 crossbred chicken population.

Copy number variations (CNVs) are large-scale changes in the DNA sequence that can affect the genetic structure and phenotype of an organism. The purpose of this study was to investigate the existing CNVs and their associations with the shank diameter (ShD) and shank length (ShL) traits using data from an F2 crossbred chicken population. To carry out the study, 312 chickens were genotyped using the Illumina 60k SNP Beadchip. The shank traits of the birds were measured from day 1 to 12 weeks of age. penncnv and cnvruler tools were used to find copy numbers and regions with copy number changes (CNVR), respectively. The CNVRanger package was used to perform a genome-wide association study between shank traits and CNVs. Gene ontology research in CNVRs was carried out using the david database. In this investigation, 966 CNVs and 606 regions with copy number changes were discovered. The copy number states and variations were randomly distributed along the length of the autosomal chromosomes. Weeks 1-4, 9 and 12 of growth revealed a significant association of copy number variations with shank traits, false discovery rate (FDR-corrected p-value < 0.01), and the majority of CNVs that were statistically significant were found on chromosomes 1-3. These CNV segments are nearby genes such as KCNJ12, FGF6 and MYF5, which are fundamental to growth and development. In addition, gene set analyses revealed terms related to muscle physiology, regulation of cellular processes and potassium channels.

Gene expression pattern and molecular mechanisms involved in Shal and Sangsari sheep fertility using RNA-Seq.

Ovulation rate and litter size are the main reproductive traits with high economic value in the sheep breeding industry. In this study, three Shal ewes (multiparous) and three Sangsari ewes (uniparous) at the age of 5 were used. The live weight was between 45 and 50 kg at an extremely body condition score of 3. These breeds are marked seasonal reproduction activity and are often bred in semi-closed breeding systems. Total RNAs were extracted from the ovarian tissues, and RNA sequencing was carried out. The DAVID (Database for Annotation, Visualization and Integrated Discovery) database was then used to annotate genes, and the string database and the Cytoscape software were used to investigate their interactions. Then path-act network analysis and gene-act network analysis were investigated. The results indicated that 19,932 genes were differentially expressed. The 5968 differentially expressed genes were identified in Shal ewe's ovarian tissue compared to Sangsari ewes (FDR < 0.05), of which 2921 genes were up-regulated and 3047 genes were down-regulated. Bioinformatics analysis exhibited that most of the biological processes and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways associated with significant DEGs (Differentially Expressed Genes) in the two studied breeds are associated with oocyte maturation and metabolism. MAPK signalling pathways and Ubiquitin-mediated proteolysis are the most important biological pathways associated with reproductive and fertility traits in the Shal breed. AKT3, MAPK8, MAPK9 and RELA genes are also important genes related to the fertility of multiparous sheep. Analysis of ovarian RNA-seq data identified that most of the differentially expressed genes were involved in various reproductive processes including folliculogenesis, ovulation, ovarian and embryonic development. The MAPK signalling pathway had the most interaction with other pathways, and the AKT3 gene could be a powerful candidate gene in the reproduction and fertility of Shal sheep. These results could pave the way for future efforts to address sheep prolificacy barriers.

Correction: Salimiyekta et al. Random Regression Model for Genetic Evaluation and Early Selection in the Iranian Holstein Population. Animals 2021, 11, 3492.

There was an error in the original publication [...].

An association of CEP78, MEF2C, VPS13A and ARRDC3 genes with survivability to heat stress in an F2 chicken population.

Heat stress is a serious problem in the poultry industry. An effective tool for improving heat tolerance can be genomic selection based on single nucleotide polymorphisms. This study was performed to identify genomic regions controlling survivability to heat stress in a population of F2 chickens that accidentally experienced acute heat stress, using Illumina 60K Chicken SNP Bead Chip. After quality control in markers, 47,730 SNPs remained for genome-wide association study (GWAS). The GWAS results indicated that markers Gga_rs16111480 (p = 8.503e-08), GGaluGA354375 (p = 5.99e-07) and Gga_rs14748694 (p = 7.085e-07) located on Z chromosome showed significant association with heat stress tolerance trait. The Gga_rs16111480 marker was located inside the CEP78 gene. The marker GGaluGA354375 was located inside the LOC101752071 gene and next to the MEF2C gene. The Gga_rs14748694 marker was adjacent to LOC101752071 and MEF2C genes. Moreover, the SNP maker of Gga_rs16111480 was located on 243 kb downstream of the VPS13A gene, and the GGaluGA354375 and Gga_rs14748694 SNPs were located on 947 kb and 888 kb downstream of the ARRDC3 gene, respectively. The results of this study suggest that apart from the gene LOC101752071, which its function was unknown, each of the two MEF2C and CEP78 genes were found to be closely related to heat stress resistance in bird.

Random Regression Model for Genetic Evaluation and Early Selection in the Iranian Holstein Population.

The objective of this study was to use a model to predict breeding values for sires and cows at an early stage of the first lactation of cows and progeny groups in the Iranian Holstein population to enable the early selection of sires. An additional objective was to estimate genetic and phenotypic parameters associated with this model. The accuracy of predicted breeding values was investigated using cross-validation based on sequential genetic evaluations emulating yearly evaluation runs. The data consisted of 2,166,925 test-day records from 456,712 cows calving between 1990 and 2015. (Co)-variance components and breeding values were estimated using a random regression test-day model and the average information (AI) restricted maximum likelihood method (REML). Legendre polynomial functions of order three were chosen to fit the additive genetic and permanent environmental effects, and a homogeneous residual variance was assumed throughout lactation. The lowest heritability of daily milk yield was estimated to be just under 0.14 in early lactation, and the highest heritability of daily milk yield was estimated to be 0.18 in mid-lactation. Cross-validation showed a highly positive correlation of predicted breeding values between consecutive yearly evaluations for both cows and sires. Correlation between predicted breeding values based only on records of early lactation (5-90 days) and records including late lactation (181-305 days) were 0.77-0.87 for cows and 0.81-0.94 for sires. These results show that we can select sires according to their daughters' early lactation information before they finish the first lactation. This can be used to decrease generation interval and to increase genetic gain in the Iranian Holstein population.

Genetic variance of metabolomic features and their relationship with body weight and body weight gain in Holstein cattle1.

In recent years, metabolomics has been used to clarify the biology underlying biological samples. In the field of animal breeding, investigating the magnitude of genetic control on the metabolomic profiles of animals and their relationships with quantitative traits adds valuable information to animal improvement schemes. In this study, we analyzed metabolomic features (MFs) extracted from the metabolomic profiles of 843 male Holstein calves. The metabolomic profiles were obtained using nuclear magnetic resonance (NMR) spectroscopy. We investigated 2 alternative methods to control for peak shifts in the NMR spectra, binning and aligning, to determine which approach was the most efficient for assessing genetic variance. Series of univariate analyses were implemented to elucidate the heritability of each MF. Furthermore, records on BW and ADG from 154 to 294 d of age (ADG154-294), 294 to 336 d of age (ADG294-336), and 154 to 336 d of age (ADG154-336) were used in a series of bivariate analyses to establish the genetic and phenotypic correlations with MFs. Bivariate analyses were only performed for MFs that had a heritability significantly different from zero. The heritabilities obtained in the univariate analyses for the MFs in the binned data set were low (<0.2). In contrast, in the aligned data set, we obtained moderate heritability (0.2 to 0.5) for 3.5% of MFs and high heritability (more than 0.5) for 1% of MFs. The bivariate analyses showed that ~12%, ~3%, ~9%, and ~9% of MFs had significant additive genetic correlations with BW, ADG154-294, ADG294-336, and ADG154-336, respectively. In all of the bivariate analyses, the percentage of significant additive genetic correlations was higher than the percentage of significant phenotypic correlations of the corresponding trait. Our results provided insights into the influence of the underlying genetic mechanisms on MFs. Further investigations in this field are needed for better understanding of the genetic relationship among the MFs and quantitative traits.

Differential Expression of the Alpha S1 Casein and Beta-Lactoglobulin Genes in Different Physiological Stages of the Adani Goats Mammary Glands.

BACKGROUND: Milk proteins genes have been the focus of the researches as the candidate target genes that play a decisive role when animal breeding is desired. OBJECTIVES: In the present study, the transcriptional levels of Beta-lactoglobulin (BLG) and Alpha S1 casein (CSN1S1) genes were investigated during prenatal, milking and drying times in mammary glands of the Adani goats which showed high and low breeding values. MATERIALS AND METHODS: The breeding values of the animals were estimated first by applying multi-trait random regression model. Using the biopsy gun, the mammary gland samples were taken and real-time PCR was applied to search the expression of the genes. Fixed factors of the model were the breeding value groups, sampling times and their interactions. RESULTS: The interactions were significant for both genes. At milking time, the high breeding value group exhibited more transcriptional levels for BLG and less transcriptional levels for CSN1S1 gene compared with the low breeding value group. The expression patterns of these genes were also different between the two breeding value groups. The maximum level of BLG and CSN1S1 transcriptions were found to occur at drying time. CONCLUSIONS: A difference in the gene expression was observed between the two groups which indicate the change in the nucleotide sequence for transcription factor binding sites, or miRNA binding sites, otherwise in the coding regions. Therefore, the variations in the coding and promoter regions of this gene should be investigated in the further studies.

Identification of learning and memory genes in canine; promoter investigation and determining the selective pressure.

One of the important behaviors of dogs is trainability which is affected by learning and memory genes. These kinds of the genes have not yet been identified in dogs. In the current research, these genes were found in animal models by mining the biological data and scientific literatures. The proteins of these genes were obtained from the UniProt database in dogs and humans. Not all homologous proteins perform similar functions, thus comparison of these proteins was studied in terms of protein families, domains, biological processes, molecular functions, and cellular location of metabolic pathways in Interpro, KEGG, Quick Go and Psort databases. The results showed that some of these proteins have the same performance in the rat or mouse, dog, and human. It is anticipated that the protein of these genes may be effective in learning and memory in dogs. Then, the expression pattern of the recognized genes was investigated in the dog hippocampus using the existing information in the GEO profile. The results showed that BDNF, TAC1 and CCK genes are expressed in the dog hippocampus, therefore, these genes could be strong candidates associated with learning and memory in dogs. Subsequently, due to the importance of the promoter regions in gene function, this region was investigated in the above genes. Analysis of the promoter indicated that the HNF-4 site of BDNF gene and the transcription start site of CCK gene is exposed to methylation. Phylogenetic analysis of protein sequences of these genes showed high similarity in each of these three genes among the studied species. The dN/dS ratio for BDNF, TAC1 and CCK genes indicates a purifying selection during the evolution of the genes.

Association of single nucleotide polymorphism of GHSR and TGFB2 genes with growth and body composition traits in sire and dam lines of a broiler chicken.

Growth hormone secretagogue receptor (GHSR) modulates many physiological processes by binding to its ligand, as well as transforming growth factor-beta 2 (TGFB2) regulates cell growth and development in animals and, therefore, are candidate genes for performance in chickens. In the current study, single nucleotide polymorphisms of GHSR C3286 > T and TGFB2 T(-640) > C were genotyped in sire and dam lines of a broiler chicken to evaluate the association with the growth and body composition traits. Least squares means analysis showed that the GHSR C3286 > T SNP was significantly (P < 0.01) associated with growth (DFI and ADG) and body composition traits (AFW and %AFW). In addition, the TGFB2 T(-640) > C SNP was associated with ADG (P < 0.05) and DFI and body composition traits (DW, LBW, BAKWT, %BMW, %HNDWT and %CW) (P < 0.01). Significant associations of the single nucleotide polymorphisms (SNPs) on the traits reported in the present study might be the distinct usage of codons in avian, or relating to an enhancer element and modulating the expression of the gene in chicken. The data indicated that these SNPs could be valuable genetic elements for selection of chickens for better performance in the population.

Sequence and phylogenetic analysis of the non-structural 3A and 3B protein-coding regions of foot-and-mouth disease virus subtype A Iran 05.

The A Iran 05 foot-and-mouth disease virus (FMDV) subtype was detected in Iran during 2005 and has proven to be highly virulent. This study was undertaken to focus on molecular and phylogenetic analysis of 3A and 3B coding-regions in the A Iran 05 field isolate. To assess the genetic relatedness of A Iran 05 isolate the nucleotide and predicted amino acid sequences of the 3AB region of type A FMDV isolates were compared with twenty previously described type A FMDV isolates. The phylogenetic tree based on the 672 bp 3AB gene sequences of type A FMDV from thirteen different locations clustered them into five distinct lineages. The A Iran 05 isolate clustered in lineage A along with four type A variants and was closely matched with viruses isolated in Turkey and Pakistan during 2005~2006. The number of protein sequence differences exhibited by each of the isolates revealed that A Iran 05 isolate contains three amino acid substitutions at positions 47 and 119 of 3A and 27 of the 3B coding region. The nucleotide identity between A Iran 05 and the other four isolates of lineage A was estimated to be 98%.