Feet and legs malformation in Nellore cattle: genetic analysis and prioritization of GWAS results.

Beef cattle affected by feet and legs malformations (FLM) cannot perform their productive and reproductive functions satisfactorily, resulting in significant economic losses. Accelerated weight gain in young animals due to increased fat deposition can lead to ligaments, tendon and joint strain and promote gene expression patterns that lead to changes in the normal architecture of the feet and legs. The possible correlated response in the FLM due to yearling weight (YW) selection suggest that this second trait could be used as an indirect selection criterion. Therefore, FLM breeding values and the genetic correlation between FLM and yearling weight (YW) were estimated for 295,031 Nellore animals by fitting a linear-threshold model in a Bayesian approach. A genome-wide association study was performed to identify genomic windows and positional candidate genes associated with FLM. The effects of single nucleotide polymorphisms (SNPs) on FLM phenotypes (affected or unaffected) were estimated using the weighted single-step genomic BLUP method, based on genotypes of 12,537 animals for 461,057 SNPs. Twelve non-overlapping windows of 20 adjacent SNPs explaining more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of candidate genes identified six genes (ATG7, EXT1, ITGA1, PPARD, SCUBE3, and SHOX) that may play a role in FLM expression due to their known role in skeletal muscle development, aberrant bone growth, lipid metabolism, intramuscular fat deposition and skeletogenesis. Identifying genes linked to foot and leg malformations enables selective breeding for healthier herds by reducing the occurrence of these conditions. Genetic markers can be used to develop tests that identify carriers of these mutations, assisting breeders in making informed breeding decisions to minimize the incidence of malformations in future generations, resulting in greater productivity and animal welfare.

Phenotypic variation in milk fatty acid composition and its association with stearoyl-CoA desaturase 1 (SCD1) gene polymorphisms in Gir cows.

Individual variation in milk fatty acid (FA) composition has been partially attributed to stearoyl-CoA desaturase 1 (SCD1) gene polymorphisms in taurine breeds, but much less is known for Zebu breeds. This study investigated the phenotypic variation in milk FA composition, and the influence of SCD1 variants on this trait and on milk fat desaturase indices (DI) in Gir cows. The functional impact of SCD1 variants was predicted using bioinformatics tools. Milk and blood samples were collected from 312 cows distributed in 10 herds from five states of Brazil. SCD1 variants were identified through target sequencing, and milk FA composition was determined by gas chromatography. Phenotypic variation in milk FA composition fell within the range reported for taurine breeds, with SCD18 index showing the lowest variation among the DI. Fourteen SCD1 variants were identified, six of which not previously described. Regarding the A293V polymorphism, all cows were homozygous for the C allele (coding for alanine), whereas all genotypes were detected for the second SNP affecting the 293 codon (G > A), with compelling evidence for functional effects. Significant associations (based on raw p-values) were found between this SNP and C12:0, cis-9, trans-11 CLA and short-chain FA, and between another SNP (rs523411937) and C15:0 and odd-chain linear FA. A new SNP on Chr26:21277069 was associated with trans-11 C18:1, cis-9, trans-11 CLA, C18:3 n-3 and n-3 FA. These findings indicate that SCD1 polymorphisms also contributes to the phenotypic variation in milk FA composition of Gir cows, with potential use in their breeding programmes.

Testicular hypoplasia in Nellore Cattle: Genetic analysis and functional analysis of genome-wide association study results.

Characterized by the incomplete development of the germinal epithelium of the seminiferous tubules, Testicular hypoplasia (TH) leads to decreased sperm concentration, increased morphological changes in sperm and azoospermia. Economic losses resulting from the disposal of affected bulls reduce the efficiency of meat production systems. A genome-wide association study and functional analysis were performed to identify genomic windows and the underlying positional candidate genes associated with TH in Nellore cattle. Phenotypic and pedigree data from 207,195 animals and genotypes (461,057 single nucleotide polymorphism, SNP) from 17,326 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. A possible correlated response on TH resulting from the selection for scrotal circumference was evaluated by using a two-trait analysis. Thus, estimated breeding values were calculated by fitting a linear-threshold animal model in a Bayesian approach. The SNP effects were estimated using the weighted single-step genomic BLUP method. Twelve non-overlapping windows of 20 adjacent SNP that explained more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of the candidate genes identified genes (KHDRBS3, GPX5, STAR, ERLIN2), which might play an important role in the expression of TH due to their known roles in the spermatogenesis process, synthesis of steroids and lipid metabolism.

Working memory and arithmetic impairments in children with FMR1 premutation and gray zone alleles.

Expansive mutations in familial mental retardation 1 (FMR1) gene have been associated with different phenotypes. Full mutations are associated with intellectual disability and autism spectrum disorder; premutations are associated with math learning difficulties and working memory impairments. In gray zone, neuropsychological development has not yet been described. Objectives: This study aimed to describe the frequency of FMR1 premutation and gray zone alleles in a school population sample representing a broad spectrum of variation in math achievement and detail school achievement and cognitive performance in the children identified with FMR1 premutation or gray zone alleles. Methods: We described a two-phase study. In the first phase, 2,195 school-age children were screened for math achievement. In the second phase, 378 children with normal intelligence were neuropsychologically assessed and genotyped for FMR1. Of these, 121 children (61 girls) performed below percentile 25 in mathematics (MD group) and 257 children (146 girls) performed above percentile 25 (control group). Results: Four pupils presented expanded alleles, one premutation and three gray zone alleles. The girl with the premutation and one boy with a gray zone allele presented impairments in working memory and arithmetic performance below percentile 6, compatible with the diagnosis of developmental dyscalculia. These children's difficulties were not associated with inaccuracy of nonsymbolic number representations or literacy impairments. Dyscalculia in these children seems to be associated mainly with working memory impairments. Conclusions: FMR1 expansions in the gray zone may contribute to dyscalculia in otherwise healthy and normally intelligent children.

Mutações expansivas no gene FMR1 têm sido associadas a diferentes fenótipos. Mutações completas estão associadas a deficiência intelectual e transtorno do espectro do autismo; pré-mutações, com dificuldades de aprendizagem de matemática e comprometimentos de memória de trabalho. Na zona cinzenta o desenvolvimento neuropsicológico ainda não foi descrito. Objetivos: Descrever a frequência de alelos pré-mutados e zona cinzenta em uma amostra escolar que representa amplo espectro de variação do desempenho em Matemática e detalhar o desempenho escolar e cognitivo em crianças identificadas com alelos pré-mutados ou zona cinzenta. Métodos: Aqui, descrevemos um estudo de duas fases. Na primeira fase, 2.195 crianças em idade escolar foram selecionadas para desempenho em Matemática. Na segunda fase, 378 crianças com inteligência normal foram avaliadas neuropsicologicamente e, em seguida, por genotipagem FMR1. Resultados: Tiveram desempenho abaixo do percentil 25 em Matemática (grupo DM) 121 crianças (61 meninas), e tiveram desempenho acima do percentil 25 (grupo controle) 257 crianças (146 meninas). Quatro alunos apresentaram alelos expandidos, sendo uma pré-mutação e três alelos da zona cinza. A menina com a pré-mutação e um menino com o alelo da zona cinza apresentaram prejuízos na memória de trabalho e desempenho aritmético abaixo do percentil 6, compatíveis com o diagnóstico de discalculia do desenvolvimento. As dificuldades dessas crianças não foram associadas à imprecisão de representações não simbólicas de números ou deficiências de alfabetização. A discalculia nessas crianças parece estar associada principalmente a deficiências da memória de trabalho. Conclusões: Em conclusão, expansões na zona cinzenta do FMR1 podem contribuir para a discalculia em crianças saudáveis com inteligência normal.

Genome-wide association study for meat tenderness in beef cattle identifies patterns of the genetic contribution in different post-mortem stages.

The beef tenderization process during the post-mortem period is one of the most important sensorial attributes and it is well-established. The aim of this study was to identify the genetic contribution pattern to meat tenderness at 7-(LMD7), 14-(LMD14), and 21-(LMD21) days post-mortem. The heritabilities for LMD7 (0.194), LMD14 (0.142) and LMD21 (0.048) are well established in the population evaluated here. However, its genetic contribution in terms of genomic candidate regions is still poorly understood. Tenderness was measured in the Longissiums thoracis using Warner-Bratzler shear force in the three post-mortem periods. A total of 4323 crossbred beef cattle were phenotyped and genotyped using the Illumina BovineSNP50K. The percentage of the total genetic variance was estimated using the weighted single-step genomic best linear unbiased prediction method. The main candidate windows for LMD7 were associated with proteolysis of myofibrillar structures and the weakening endomysium and perimysium. Candidate windows for LMD14 and LMD21 were mapped in bovine QTLs for body composition, height and growth. Results presented herein highlight, the largest contribution of proteolysis related processes before 14-days post-mortem and body composition characteristics in later stages for meat tenderness.

Putative Causal Variant on Vlgr1 for the Epileptic Phenotype in the Model Wistar Audiogenic Rat.

Wistar Audiogenic Rat is an epilepsy model whose animals are predisposed to develop seizures induced by acoustic stimulation. This model was developed by selective reproduction and presents a consistent genetic profile due to the several generations of inbreeding. In this study, we performed an analysis of WAR RNA-Seq data, aiming identified at genetic variants that may be involved in the epileptic phenotype. Seventeen thousand eighty-five predicted variants were identified as unique to the WAR model, of which 15,915 variants are SNPs and 1,170 INDELs. We filter the predicted variants by pre-established criteria and selected five for validation by Sanger sequencing. The genetic variant c.14198T>C in the Vlgr1 gene was confirmed in the WAR model. Vlgr1 encodes an adhesion receptor that is involved in the myelination process, in the development of stereocilia of the inner ear, and was already associated with the audiogenic seizures presented by the mice Frings. The transcriptional quantification of Vlgr1 revealed the downregulation this gene in the corpus quadrigeminum of WAR, and the protein modeling predicted that the mutated residue alters the structure of a domain of the VLGR1 receptor. We believe that Vlgr1 gene may be related to the predisposition of WAR to seizures and suggest the mutation Vlgr1/Q4695R as putative causal variant, and the first molecular marker of the WAR strain.

Genetic mechanisms underlying spermatic and testicular traits within and among cattle breeds: systematic review and prioritization of GWAS results.

Reduced bull fertility imposes economic losses in bovine herds. Specifically, testicular and spermatic traits are important indicators of reproductive efficiency. Several genome-wide association studies (GWAS) have identified genomic regions associated with these fertility traits. The aims of this study were as follows: 1) to perform a systematic review of GWAS results for spermatic and testicular traits in cattle and 2) to identify key functional candidate genes for these traits. The identification of functional candidate genes was performed using a systems biology approach, where genes shared between traits and studies were evaluated by a guilt by association gene prioritization (GUILDify and ToppGene software) in order to identify the best functional candidates. These candidate genes were integrated and analyzed in order to identify overlapping patterns among traits and breeds. Results showed that GWAS for testicular-related traits have been developed for beef breeds only, whereas the majority of GWAS for spermatic-related traits were conducted using dairy breeds. When comparing traits measured within the same study, the highest number of genes shared between different traits was observed, indicating a high impact of the population genetic structure and environmental effects. Several chromosomal regions were enriched for functional candidate genes associated with fertility traits. Moreover, multiple functional candidate genes were enriched for markers in a species-specific basis, taurine (Bos taurus) or indicine (Bos indicus). For the different candidate regions identified in the GWAS in the literature, functional candidate genes were detected as follows: B. Taurus chromosome X (BTX) (TEX11, IRAK, CDK16, ATP7A, ATRX, HDAC6, FMR1, L1CAM, MECP2, etc.), BTA17 (TRPV4 and DYNLL1), and BTA14 (MOS, FABP5, ZFPM2). These genes are responsible for regulating important metabolic pathways or biological processes associated with fertility, such as progression of spermatogenesis, control of ciliary activity, development of Sertoli cells, DNA integrity in spermatozoa, and homeostasis of testicular cells. This study represents the first systematic review on male fertility traits in cattle using a system biology approach to identify key candidate genes for these traits.

Combining multi-OMICs information to identify key-regulator genes for pleiotropic effect on fertility and production traits in beef cattle.

The identification of biological processes related to the regulation of complex traits is a difficult task. Commonly, complex traits are regulated through a multitude of genes contributing each to a small part of the total genetic variance. Additionally, some loci can simultaneously regulate several complex traits, a phenomenon defined as pleiotropy. The lack of understanding on the biological processes responsible for the regulation of these traits results in the decrease of selection efficiency and the selection of undesirable hitchhiking effects. The identification of pleiotropic key-regulator genes can assist in developing important tools for investigating biological processes underlying complex traits. A multi-breed and multi-OMICs approach was applied to study the pleiotropic effects of key-regulator genes using three independent beef cattle populations evaluated for fertility traits. A pleiotropic map for 32 traits related to growth, feed efficiency, carcass and meat quality, and reproduction was used to identify genes shared among the different populations and breeds in pleiotropic regions. Furthermore, data-mining analyses were performed using the Cattle QTL database (CattleQTLdb) to identify the QTL category annotated in the regions around the genes shared among breeds. This approach allowed the identification of a main gene network (composed of 38 genes) shared among breeds. This gene network was significantly associated with thyroid activity, among other biological processes, and displayed a high regulatory potential. In addition, it was possible to identify genes with pleiotropic effects related to crucial biological processes that regulate economically relevant traits associated with fertility, production and health, such as MYC, PPARG, GSK3B, TG and IYD genes. These genes will be further investigated to better understand the biological processes involved in the expression of complex traits and assist in the identification of functional variants associated with undesirable phenotypes, such as decreased fertility, poor feed efficiency and negative energetic balance.

Identification of Candidate Genes for Reactivity in Guzerat (Bos indicus) Cattle: A Genome-Wide Association Study.

Temperament is fundamental to animal production due to its direct influence on the animal-herdsman relationship. When compared to calm animals, the aggressive, anxious or fearful ones exhibit less weight gain, lower reproductive efficiency, decreased milk production and higher herd maintenance costs, all of which contribute to reduced profits. However, temperament is a trait that is complex and difficult to assess. Recently, a new quantitative system, REATEST®, for assessing reactivity, a phenotype of temperament, was developed. Herein, we describe the results of a Genome-wide association study for reactivity, assessed using REATEST® with a sample of 754 females from five dual-purpose (milk and meat production) Guzerat (Bos indicus) herds. Genotyping was performed using a 50k SNP chip and a two-step mixed model approach (Grammar-Gamma) with a one-by-one marker regression was used to identify QTLs. QTLs for reactivity were identified on chromosomes BTA1, BTA5, BTA14, and BTA25. Five intronic and two intergenic markers were significantly associated with reactivity. POU1F1, DRD3, VWA3A, ZBTB20, EPHA6, SNRPF and NTN4 were identified as candidate genes. Previous QTL reports for temperament traits, covering areas surrounding the SNPs/genes identified here, further corroborate these associations. The seven genes identified in the present study explain 20.5% of reactivity variance and give a better understanding of temperament biology.