EEF1A1 transcription cofactor gene polymorphism is associated with muscle gene expression and residual feed intake in Nelore cattle.

Cis-acting effects of noncoding variants on gene expression and regulatory molecules constitute a significant factor for phenotypic variation in complex traits. To provide new insights into the impacts of single-nucleotide polymorphisms (SNPs) on transcription factors (TFs) and transcription cofactors (TcoF) coding genes, we carried out a multi-omic analysis to identify cis-regulatory effects of SNPs on these genes' expression in muscle and describe their association with feed efficiency-related traits in Nelore cattle. As a result, we identified one SNP, the rs137256008C > T, predicted to impact the EEF1A1 gene expression (ß = 3.02; P-value = 3.51E-03) and the residual feed intake trait (ß = - 3.47; P-value = 0.02). This SNP was predicted to modify transcription factor sites and overlaps with several QTL for feed efficiency traits. In addition, co-expression network analyses showed that animals containing the T allele of the rs137256008 SNP may be triggering changes in the gene network. Therefore, our analyses reinforce and contribute to a better understanding of the biological mechanisms underlying gene expression control of feed efficiency traits in bovines. The cis-regulatory SNP can be used as biomarker for feed efficiency in Nelore cattle.

Genetic effects of heat stress on milk fatty acids in Brazilian Holstein cattle.

The present study aimed to estimate covariance components of milk fatty acids (FA) and to compare the genomic estimated breeding values under general and heat-stress effects. Data consisted of 38,762 test-day records from 6,344 Holstein cows obtained from May 2012 through January 2018 on 4 dairy herds from Brazil. Single-trait repeatability test-day models with random regressions as a function of temperature-humidity index values were used for genetic analyses. The models included contemporary groups, parity order (1-6), and days in milk classes as fixed effects, and general and thermotolerance additive genetic and permanent environmental as random effects. Notably, differences in heritability estimates between environments (general and heat stress) increased (0.03 to 0.06) for unsaturated FA traits, such as unsaturated, monounsaturated, and polyunsaturated, at higher heat-stress levels. In contrast, heritability estimated between environments for saturated FA traits, including saturated FA, palmitic acid (C16:0), and stearic acid (C18:0) did not observe significant differences between environments. In addition, our study revealed negative genetic correlations between general and heat-stress additive genetic effects (antagonistic effect) for the saturated FA, C16:0, C18:0, and C18:1, which ranged from -0.007 to -0.32. Spearman's ranking correlation between genomic estimated breeding values ranged from -0.27 to 0.99. Results indicated a moderate to strong interaction of genotype by the environment for most FA traits comparing a heat-stress environment with thermoneutral conditions. Our findings point out novel opportunities to explore the use of FA milk profile and heat-stress models.

Genome-wide association study for milk production traits in a Brazilian Holstein population.

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.

Estimates of genomic heritability and genome-wide association study for fatty acids profile in Santa Inês sheep.

BACKGROUND: Despite the health concerns and nutritional importance of fatty acids, there is a relative paucity of studies in the literature that report genetic or genomic parameters, especially in the case of sheep populations. To investigate the genetic architecture of fatty acid composition of sheep, we conducted genome-wide association studies (GWAS) and estimated genomic heritabilities for fatty acid profile in Longissimus dorsi muscle of 216 male sheep. RESULTS: Genomic heritability estimates for fatty acid content ranged from 0.25 to 0.46, indicating that substantial genetic variation exists for the evaluated traits. Therefore, it is possible to alter fatty acid profiles through selection. Twenty-seven genomic regions of 10 adjacent SNPs associated with fatty acids composition were identified on chromosomes 1, 2, 3, 5, 8, 12, 14, 15, 16, 17, and 18, each explaining ≥0.30% of the additive genetic variance. Twenty-three genes supporting the understanding of genetic mechanisms of fat composition in sheep were identified in these regions, such as DGAT2, TRHDE, TPH2, ME1, C6, C7, UBE3D, PARP14, and MRPS30. CONCLUSIONS: Estimates of genomic heritabilities and elucidating important genomic regions can contribute to a better understanding of the genetic control of fatty acid deposition and improve the selection strategies to enhance meat quality and health attributes.

Genetic parameters for milk fatty acids, milk yield and quality traits of a Holstein cattle population reared under tropical conditions.

Information about genetic parameters is essential for selection decisions and genetic evaluation. These estimates are population specific; however, there are few studies with dairy cattle populations reared under tropical and sub-tropical conditions. Thus, the aim was to obtain estimates of heritability and genetic correlations for milk yield and quality traits using pedigree and genomic information from a Holstein population maintained in a tropical environment. Phenotypic records (n = 36 457) of 4203 cows as well as the genotypes for 57 368 single nucleotide polymorphisms from 755 of these cows were used. Covariance components were estimated using the restricted maximum likelihood method under a mixed animal model, considering a pedigree-based relationship matrix or a combined pedigree-genomic matrix. High heritabilities (around 0.30) were estimated for lactose and protein content in milk whereas moderate values (between 0.19 and 0.26) were obtained for percentages of fat, saturated fatty acids and palmitic acid in milk. Genetic correlations ranging from -0.38 to -0.13 were determined between milk yield and composition traits. The smaller estimates compared to other similar studies can be due to poor environmental conditions, which may reduce genetic variability. These results highlight the importance in using genetic parameters estimated in the population under evaluation for selection decisions.

An alternative mode of translation permits production of a variant NBS1 protein from the common Nijmegen breakage syndrome allele.

Nijmegen breakage syndrome (NBS) is a rare chromosomal-instability syndrome associated with cancer predisposition, radiosensitivity and radioresistant DNA synthesis-S phase checkpoint deficiency, which results in the failure to suppress DNA replication origins following DNA damage. Approximately 90% of NBS patients are homozygous for the 657del5 allele, a truncating mutation of NBS1 that causes premature termination at codon 219. Because null mutations in MRE11 and RAD50, which encode binding partners of NBS1, are lethal in vertebrates, and mouse Nbs1-null mutants are inviable, we tested the hypothesis that the NBS1 657del5 mutation was a hypomorphic defect. We showed that NBS cells contain the predicted 26-kD amino-terminal protein fragment, NBS1p26, and a 70-kD NBS1 protein (NBS1p70) lacking the native N terminus. The NBSp26 protein is not physically associated with the MRE11 complex, whereas the p70 species is physically associated with it. NBS1p70 is produced by internal translation initiation within the NBS1 mRNA using an open reading frame generated by the 657del5 frameshift. We propose that the common NBS1 allele encodes a partially functional protein that diminishes the severity of the NBS phenotype.

Cell-cycle-regulated association of RAD50/MRE11/NBS1 with TRF2 and human telomeres.

Telomeres allow cells to distinguish natural chromosome ends from damaged DNA and protect the ends from degradation and fusion. In human cells, telomere protection depends on the TTAGGG repeat binding factor, TRF2 (refs 1-4), which has been proposed to remodel telomeres into large duplex loops (t-loops). Here we show by nanoelectrospray tandem mass spectrometry that RAD50 protein is present in TRF2 immunocomplexes. Protein blotting showed that a small fraction of RAD50, MRE11 and the third component of the MRE11 double-strand break (DSB) repair complex, the Nijmegen breakage syndrome protein (NBS1), is associated with TRF2. Indirect immunofluorescence demonstrated the presence of RAD50 and MRE11 at interphase telomeres. NBS1 was associated with TRF2 and telomeres in S phase, but not in G1 or G2. Although the MRE11 complex accumulated in irradiation-induced foci (IRIFs) in response to gamma-irradiation, TRF2 did not relocate to IRIFs and irradiation did not affect the association of TRF2 with the MRE11 complex, arguing against a role for TRF2 in DSB repair. Instead, we propose that the MRE11 complex functions at telomeres, possibly by modulating t-loop formation.

The Mre11 complex and ATM: collaborating to navigate S phase.

Recently, findings regarding a group of cancer predisposition and chromosome instability syndromes, Nijmegen breakage syndrome (NBS), the ataxia-telangiectasia-like disorder (A-TLD) and ataxia telangiectasia have shed light on the unexpected role of recombinational DNA repair proteins in DNA-damage-dependent cell-cycle regulation. Mutations in the Mre11 complex cause A-TLD and NBS. In addition, functions of the Mre11 complex have been biochemically linked to ATM, the large protein kinase that is defective in ataxia-telangiectasia cells by the observation that Nbs1 is a bona fide substrate of the ATM kinase.

In situ visualization of DNA double-strand break repair in human fibroblasts.

A method was developed to examine DNA repair within the intact cell. Ultrasoft x-rays were used to induce DNA double-strand breaks (DSBs) in defined subnuclear volumes of human fibroblasts and DNA repair was visualized at those sites. The DSBs remained in a fixed position during the initial stages of DNA repair, and the DSB repair protein hMre11 migrated to the sites of damage within 30 minutes. In contrast, hRad51, a human RecA homolog, did not localize at sites of DNA damage, a finding consistent with the distinct roles of these proteins in DNA repair.

DNA damage-dependent nuclear dynamics of the Mre11 complex.

The Mre11 complex has been implicated in diverse aspects of the cellular response to DNA damage. We used in situ fractionation of human fibroblasts to carry out cytologic analysis of Mre11 complex proteins in the double-strand break (DSB) response. In situ fractionation removes most nucleoplasmic protein, permitting immunofluorescent localization of proteins that become more avidly bound to nuclear structures after induction of DNA damage. We found that a fraction of the Mre11 complex was bound to promyelocyte leukemia protein bodies in undamaged cells. Within 10 min after gamma irradiation, nuclear retention of the Mre11 complex in small granular foci was observed and persisted until 2 h postirradiation. In light of the previous demonstration that the Mre11 complex associated with ionizing radiation (IR)-induced DSBs, we infer that the protein retained under these conditions was associated with DNA damage. We also observed increased retention of Rad51 following IR treatment, although IR induced Rad51 foci were distinct from Mre11 foci. The ATM kinase, which phosphorylates Nbs1 during activation of the S-phase checkpoint, was not required for the Mre11 complex to associate with DNA damage. These data suggest that the functions of the Mre11 complex in the DSB response are implicitly dependent upon its ability to detect DNA damage.

DNA ligase I mediates essential functions in mammalian cells.

DNA replication, repair, and recombination are essential processes in mammalian cells. Hence, the application of gene targeting to the study of these DNA metabolic pathways requires the creation of nonnull mutations. We have developed a method for introducing partially defective mutants in murine embryonic stem cells that circumvents the problem of cellular lethality of targeted mutations at essential loci. Using this approach, we have determined that mammalian DNA ligase I is essential for cell viability. Thus, DNA ligases II and III are not redundant with DNA ligase I for the function(s) associated with cell proliferation. Partial complementation of the lethal DNA ligase I null mutation allowed the creation of deficient embryonic stem cell lines. We found that a wild-type DNA ligase I cDNA, as well as a variant DNA ligase I cDNA, was able to rescue the lethality of the homozygous null mutation, whereas an N-terminal deletion mutant consisting of the minimal DNA ligase I catalytic domain was not. This observation demonstrates that sequences outside the DNA ligase I catalytic domain are essential for DNA ligase I function in vivo.

The Mre11-Rad50-Xrs2 protein complex facilitates homologous recombination-based double-strand break repair in Saccharomyces cerevisiae.

Saccharomyces cerevisiae mre11Delta mutants are profoundly deficient in double-strand break (DSB) repair, indicating that the Mre11-Rad50-Xrs2 protein complex plays a central role in the cellular response to DNA DSBs. In this study, we examined the role of the complex in homologous recombination, the primary mode of DSB repair in yeast. We measured survival in synchronous cultures following irradiation and scored sister chromatid and interhomologue recombination genetically. mre11Delta strains were profoundly sensitive to ionizing radiation (IR) throughout the cell cycle. Mutant strains exhibited decreased frequencies of IR-induced sister chromatid and interhomologue recombination, indicating a general deficiency in homologous recombination-based DSB repair. Since a nuclease-deficient mre11 mutant was not impaired in these assays, it appears that the role of the S. cerevisiae Mre11-Rad50-Xrs2 protein complex in facilitating homologous recombination is independent of its nuclease activities.

hMre11 and hRad50 nuclear foci are induced during the normal cellular response to DNA double-strand breaks.

We previously identified a conserved multiprotein complex that includes hMre11 and hRad50. In this study, we used immunofluorescence to investigate the role of this complex in DNA double-strand break (DSB) repair. hMre11 and hRad50 form discrete nuclear foci in response to treatment with DSB-inducing agents but not in response to UV irradiation. hMre11 and hRad50 foci colocalize after treatment with ionizing radiation and are distinct from those of the DSB repair protein, hRad51. Our data indicate that an irradiated cell is competent to form either hMre11-hRad50 foci or hRad51 foci, but not both. The multiplicity of hMre11 and hRad50 foci is much higher in the DSB repair-deficient cell line 180BR than in repair-proficient cells. hMre11-hRad50 focus formation is markedly reduced in cells derived from ataxia-telangiectasia patients, whereas hRad51 focus formation is markedly increased. These experiments support genetic evidence from Saccharomyces cerevisiae indicating that Mre11-Rad50 have roles distinct from that of Rad51 in DSB repair. Further, these data indicate that hMre11-hRad50 foci form in response to DNA DSBs and are dependent upon a DNA damage-induced signaling pathway.

Human Rad50 is physically associated with human Mre11: identification of a conserved multiprotein complex implicated in recombinational DNA repair.

In this report, we describe the identification and molecular characterization of a human RAD50 homolog, hRAD50. hRAD50 was included in a collection of cDNAs which were isolated by a direct cDNA selection strategy focused on the chromosomal interval spanning 5q23 to 5q31. Alterations of the 5q23-q31 interval are frequently observed in myelodysplasia and myeloid leukemia. This strategy was thus undertaken to create a detailed genetic map of that region. Saccharomyces cerevisiae RAD50 (ScRAD50) is one of three yeast RAD52 epistasis group members (ScRAD50, ScMRE11, and ScXRS2) in which mutations eliminate meiotic recombination but confer a hyperrecombinational phenotype in mitotic cells. The yeast Rad50, Mre11, and Xrs2 proteins appear to act in a multiprotein complex, consistent with the observation that the corresponding mutants confer essentially identical phenotypes. In this report, we demonstrate that the human Rad50 and Mre11 proteins are stably associated in a protein complex which may include three other proteins. hRAD50 is expressed in all tissues examined, but mRNA levels are significantly higher in the testis. Other human RAD52 epistasis group homologs exhibit this expression pattern, suggesting the involvement of human RAD52 epistasis group proteins in meiotic recombination. Human RAD52 epistasis group proteins are highly conserved and act in protein complexes that are analogous to those of their yeast counterparts. These findings indicate that the function of the RAD52 epistasis group is conserved in human cells.

Mre11 complex and DNA replication: linkage to E2F and sites of DNA synthesis.

We show that the Mre11 complex associates with E2F family members via the Nbs1 N terminus. This association and Nbs1 phosphorylation are correlated with S-phase checkpoint proficiency, whereas neither is sufficient individually for checkpoint activation. The Nbs1 E2F interaction occurred near the Epstein-Barr virus origin of replication as well as near a chromosomal replication origin in the c-myc promoter region and was restricted to S-phase cells. The Mre11 complex colocalized with PCNA at replication forks throughout S phase, both prior to and coincident with the appearance of nascent DNA. These data suggest that the Mre11 complex suppresses genomic instability through its influence on both the regulation and progression of DNA replication.

Radiosensitivity in Nijmegen Breakage Syndrome cells is attributable to a repair defect and not cell cycle checkpoint defects.

Cells derived from Nijmegen Breakage Syndrome (NBS) patients display radiosensitivity and cell cycle checkpoint defects. Here, we examine whether the radiosensitivity of NBS cells is the result of a repair defect or whether it can be attributed to impaired checkpoint arrest. We report a small increased fraction of unrejoined double strand breaks and, more significantly, increased chromosome breaks in noncycling NBS cells at 24 h after irradiation. One of the NBS lines examined (347BR) was atypical in showing a nearly normal checkpoint response. In contrast to the mild checkpoint defect, 347BR displays marked y-ray sensitivity similar to that shown by other NBS lines. Thus, the gamma-ray sensitivity correlates with the repair defect rather than impaired checkpoint control. Taken together, the results provide direct evidence for a repair defect in NBS cells and are inconsistent with the suggestion that the radiosensitivity is attributable only to impaired checkpoint arrest. 347BR also displays elevated spontaneous damage that cannot be attributed to impaired G2-M arrest, suggesting a function of Nbsl in decreasing or limiting the impact of spontaneously arising double strand breaks.

Alteration of N-terminal phosphoesterase signature motifs inactivates Saccharomyces cerevisiae Mre11.

Saccharomyces cerevisiae Mre11, Rad50, and Xrs2 function in a protein complex that is important for nonhomologous recombination. Null mutants of MRE11, RAD50, and XRS2 are characterized by ionizing radiation sensitivity and mitotic interhomologue hyperrecombination. We mutagenized the four highly conserved phosphoesterase signature motifs of Mre11 to create mre11-11, mre11-2, mre11-3, and mre11-4 and assessed the functional consequences of these mutant alleles with respect to mitotic interhomologue recombination, chromosome loss, ionizing radiation sensitivity, double-strand break repair, and protein interaction. We found that mre11 mutants that behaved as the null were sensitive to ionizing radiation and deficient in double-strand break repair. We also observed that these null mutants exhibited a hyperrecombination phenotype in mitotic cells, consistent with previous reports, but did not exhibit an increased frequency of chromosome loss. Differential ionizing radiation sensitivities among the hypomorphic mre11 alleles correlated with the trends observed in the other phenotypes examined. Two-hybrid interaction testing showed that all but one of the mre11 mutations disrupted the Mre11-Rad50 interaction. Mutagenesis of the phosphoesterase signatures in Mre11 thus demonstrated the importance of these conserved motifs for recombinational DNA repair.

T cell leakiness in scid mice.

FACS analysis showed that the incidence of leaky T cells increases with age, such that virtually all old scid mice (greater than 1 year) contain detectable CD3+ cells. The number of detectable T cells remained very low; individual old scid mice generally contained less than 10(5) CD3+ cells. When CD3+ populations in individual leaky mice were analyzed for expression of the T cell subset markers, CD4 and CD8, the ratios of CD4/CD8 were found to be markedly skewed relative to normal mice. This suggested the presence of very few T cell clones. Indeed, the analysis of TCR gene rearrangements in polyclonally stimulated T cell cultures revealed only 1-5 clones in the pooled spleen and lymph nodes of individual old scid mice. These studies also indicated that TCR gene rearrangements in the majority of the stimulated T cell cultures did not contain abnormal J-associated deletions that are characteristic of antigen receptor genes of scid lymphomas. Four of five alloreactive T cell clones from leaky scid mice also apparently lacked abnormal J-associated deletions in their rearranged TCR alleles. Therefore, most leaky lymphocytes appear to derive from progenitors with normal or near-normal scid recombinase activity. However, one of five leaky T cell clones (S1233) and one Con A stimulated monoclonal culture (8706) contained both normally and abnormally rearranged TCR genes. The configuration of TCR loci in such clones may reflect the ability of the defective scid recombinase to mediate normal rearrangements at a low frequency.

Genetic grouping strategies in selection efficiency of composite beef cattle ( × ).

The inclusion of genetic groups in sire evaluation has been widely used to represent genetic differences among animals not accounted for by the absence of parentage data. However, the definition of these groups is still arbitrary, and studies assessing the effects of genetic grouping strategies on the selection efficiency are rare. Therefore, the aim in this study was to compare genetic grouping strategies for animals with unknown parentage in prediction of breeding values (EBV). The total of 179,302 records of weaning weight (WW), 29,825 records of scrotal circumference (SC), and 70,302 records of muscling score (MUSC) from Montana Tropical animals, a Brazilian composite beef cattle population, were used. Genetic grouping strategies involving year of birth, sex of the unknown parent, birth farm, breed composition, and their combinations were evaluated. Estimated breeding values were predicted for each approach simulating a loss of genealogy data. Thereafter, these EBV were compared to those obtained in an analysis involving a real relationship matrix to estimate selection efficiency and correlations between EBV and animal rankings. The analysis model included the fixed effects of contemporary groups and class of the dam age at calving, the covariates of additive and nonadditive genetic effects, and age, and the additive genetic effect of animal as random effects. A second model also included the fixed effects of genetic group. The use of genetic groups resulted in means of selection efficiency and correlation of 70.4 to 97.1% and 0.51 to 0.94 for WW, 85.8 to 98.8% and 0.82 to 0.98 for SC, and 85.1 to 98.6% and 0.74 to 0.97 for MUSC, respectively. High selection efficiencies were observed for year of birth and breed composition strategies. The maximum absolute difference in annual genetic gain estimated through the use of complete genealogy and genetic groups were 0.38 kg for WW, 0.02 cm for SC, and 0.01 for MUSC, with lower differences obtained when year of birth was adopted as a genetic group criterion. Grouping strategy must consider selection decisions and the number of genetic groups formed, in the way that genetic groups represent the genetic differences in population and allow an adequate prediction of EBV.