Novel putative causal mutations associated with fat traits in Nellore cattle uncovered by eQTLs located in open chromatin regions

Intramuscular fat (IMF) and backfat thickness (BFT) are critical economic traits impacting meat quality. However, the genetic variants controlling these traits need to be better understood. To advance knowledge in this area, we integrated RNA-seq and single nucleotide polymorphisms (SNPs) identified in genomic and transcriptomic data to generate a linkage disequilibrium filtered panel of 553,581 variants. Expression quantitative trait loci (eQTL) analysis revealed 36,916 cis-eQTLs and 14,408 trans-eQTLs. Association analysis resulted in three eQTLs associated with BFT and 24 with IMF. Functional enrichment analysis of genes regulated by these 27 eQTLs revealed noteworthy pathways that can play a fundamental role in lipid metabolism and fat deposition, such as immune response, cytoskeleton remodeling, iron transport, and phospholipid metabolism. We next used ATAC-Seq assay to identify and overlap eQTL and open chromatin regions. Six eQTLs were in regulatory regions, four in predicted insulators and possible CCCTC-binding factor DNA binding sites, one in an active enhancer region, and the last in a low signal region. Our results provided novel insights into the transcriptional regulation of IMF and BFT, unraveling putative regulatory variants.

Ruminal and feces metabolites associated with feed efficiency, water intake and methane emission in Nelore bulls

The objectives of this study were twofold: (1) to identify potential differences in the ruminal and fecal metabolite profiles of Nelore bulls under different nutritional interventions; and (2) to identify metabolites associated with cattle sustainability related-traits. We used different nutritional interventions in the feedlot: conventional (Conv; n = 26), and by-product (ByPr, n = 26). Thirty-eight ruminal fluid and 27 fecal metabolites were significantly different (P < 0.05) between the ByPr and Conv groups. Individual dry matter intake (DMI), residual feed intake (RFI), observed water intake (OWI), predicted water intake (WI), and residual water intake (RWI) phenotypes were lower (P < 0.05) in the Conv group, while the ByPr group exhibited lower methane emission (ME) (P < 0.05). Ruminal fluid dimethylamine was significantly associated (P < 0.05) with DMI, RFI, FE (feed efficiency), OWI and WI. Aspartate was associated (P < 0.05) with DMI, RFI, FE and WI. Fecal C22:1n9 was significantly associated with OWI and RWI (P < 0.05). Fatty acid C14:0 and hypoxanthine were significantly associated with DMI and RFI (P < 0.05). The results demonstrated that different nutritional interventions alter ruminal and fecal metabolites and provided new insights into the relationship of these metabolites with feed efficiency and water intake traits in Nelore bulls.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR &lt; 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

Simulation of group testing scenarios can boost COVID-19 screening power

The COVID-19 has severely affected economies and health systems around the world. Mass testing could work as a powerful alternative to restrain disease dissemination, but the shortage of reagents is a limiting factor. A solution to optimize test usage relies on ‘grouping’ or ‘pooling’ strategies, which combine a set of individuals in a single reaction. To compare different group testing configurations, we developed the poolingr package, which performs an innovative hybrid in silico/in vitro approach to search for optimal testing configurations. We used 6759 viral load values, observed in 2389 positive individuals, to simulate a wide range of scenarios. We found that larger groups (>100) framed into multi-stage setups (up to six stages) could largely boost the power to detect spreaders. Although the boost was dependent on the disease prevalence, our method could point to cheaper grouping schemes to better mitigate COVID-19 dissemination through identification and quarantine recommendation for positive individuals.

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.

Stool and Ruminal Microbiome Components Associated With Methane Emission and Feed Efficiency in Nelore Beef Cattle

Background: The impact of extreme changes in weather patterns on the economy and human welfare is one of the biggest challenges our civilization faces. From anthropogenic contributions to climate change, reducing the impact of farming activities is a priority since it is responsible for up to 18% of global greenhouse gas emissions. To this end, we tested whether ruminal and stool microbiome components could be used as biomarkers for methane emission and feed efficiency in bovine by studying 52 Brazilian Nelore bulls belonging to two feed intervention treatment groups, that is, conventional and by-product-based diets.

Results: We identified a total of 5,693 amplicon sequence variants (ASVs) in the Nelore bulls’ microbiomes. A Differential abundance analysis with the ANCOM approach identified 30 bacterial and 15 archaeal ASVs as differentially abundant (DA) among treatment groups. An association analysis using Maaslin2 software and a linear mixed model indicated that bacterial ASVs are linked to the host’s residual methane emission (RCH₄) and residual feed intake (RFI) phenotype variation, suggesting their potential as targets for interventions or biomarkers.

Conclusion: The feed composition induced significant differences in both abundance and richness of ruminal and stool microbial populations in ruminants of the Nelore breed. The industrial by-product-based dietary treatment applied to our experimental groups influenced the microbiome diversity of bacteria and archaea but not of protozoa. ASVs were associated with RCH₄ emission and RFI in ruminal and stool microbiomes. While ruminal ASVs were expected to influence CH₄ emission and RFI, the relationship of stool taxa, such as Alistipes and Rikenellaceae (gut group RC9), with these traits was not reported before and might be associated with host health due to their link to anti-inflammatory compounds. Overall, the ASVs associated here have the potential to be used as biomarkers for these complex phenotypes.

Genetic effects of heat stress on milk fatty acids in a 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.

A Missense Mutation in the MYBPH Gene Is Associated With Abdominal Fat Traits in Meat-Type Chickens

Chicken is an important source of protein for human nutrition and a model system for growth and developmental biology. Although the genetic architecture of quantitative traits in meat-type chickens has been the subject of ongoing investigation, the identification of mutations associated with carcass traits of economic interest remains challenging. Therefore, our aim was to identify predicted deleterious mutation, which potentially affects protein function, and test if they were associated with carcass traits in chickens. For that, we performed a genome-wide association analysis (GWAS) for breast, thigh and drumstick traits in meat-type chickens and detected 19 unique quantitative trait loci (QTL). We then used: (1) the identified windows; (2) QTL for abdominal fat detected in a previous study with the same population and (3) previously obtained whole genome sequence data, to identify 18 predicted deleterious single nucleotide polymorphisms (SNPs) in those QTL for further association with breast, thigh, drumstick and abdominal fat traits. Using the additive model, a predicted deleterious SNP c.482C > T (SIFT score of 0.4) was associated (p-value < 0.05) with abdominal fat weight and percentage. This SNP is in the second exon of the MYBPH gene, and its allele frequency deviates from Hardy–Weinberg equilibrium. In conclusion, our study provides evidence that the c.482C > T SNP in the MYBPH gene is a putative causal mutation for fat deposition in meat-type chickens.

Genetic parameters for milk yield and quality traits of Brazilian Holstein cows as a function of temperature and humidity index

Measurements of milk yield (MY), somatic cell score (SCS), percentage of fat (FP), protein (PP), lactose (LP), casein (CP) and percentage of palmitic (C16:0), stearic (C18:0), oleic (C18:1), total saturated (SFA), unsaturated (UFA), monounsaturated (MUFA) and polyunsaturated (PUFA) fatty acids in milk from 5,224 Holstein cows were evaluated as a function of a temperature and humidity index (THI). Legendre orthogonal polynomials from second to seventh order were tested. The best fit order for MY, PP and C18:0 was the third, whereas the second for all other traits. The heritability estimates decreased for MY (0.31 to 0.14), FP (0.28 to 0.16), LP (0.43 to 0.30), SCS (0.14 to 0.09), SFA (0.33 to 0.22) and C16:0 (0.31 to 0.26), whereas increased for CP (0.32 to 0.42), MUFA (0.08 to 0.13), UFA (0.07 to 0.11) and C18:1 (0.07 to 0.11) as the THI level increased. For PP, heritabilities (0.26 to 0.39) presented larger values in intermediate THI. For PUFA and C18:0, heritabilities were approximately constant (0.13 to 0.14 and 0.15, respectively). However, the greatest variations may have been the result of the limitations of Legendre polynomials at the extreme points of the curve, and the pattern of heritabilities curves was approximately constant for the evaluated traits. Spearman's rank correlations between breeding values in extreme THI levels were greater than 0.80 for all traits considering all animals, only cows and only bulls. When considering the top 1% and the top 50% animals (only cows, only bulls and all), Spearman correlations smaller than 0.70 were found, suggesting reranking of the animals. Although there was little variation in the variance components over THI, it is possible that there is no heat stress in the animals studied, because, on average, there was no great impact of the thermal load on the traits. One possible explanation is the use of herds with little climatic difference among herds, as well as the use of fans and sprinklers into the barns. However, the THI levels may be important factors in the selection process, as reranking of animals was verified.

Genetic regulators of mineral amount in Nelore cattle muscle predicted by a new co-expression and regulatory impact factor approach

Mineral contents in bovine muscle can affect meat quality, growth, health, and reproductive traits. To better understand the genetic basis of this phenotype in Nelore (Bos indicus) cattle, we analysed genome-wide mRNA and miRNA expression data from 114 muscle samples. The analysis implemented a new application for two complementary algorithms: the partial correlation and information theory (PCIT) and the regulatory impact factor (RIF), in which we included the estimated genomic breeding values (GEBVs) for the phenotypes additionally to the expression levels, originally proposed for these methods. We used PCIT to determine putative regulatory relationships based on significant associations between gene expression and GEBVs for each mineral amount. Then, RIF was adopted to determine the regulatory impact of genes and miRNAs expression over the GEBVs for the mineral amounts. We also investigated over-represented pathways, as well as pieces of evidences from previous studies carried in the same population and in the literature, to determine regulatory genes for the mineral amounts. For example, NOX1 expression level was positively correlated to Zinc and has been described as Zinc-regulated in humans. Based on our approach, we were able to identify genes, miRNAs and pathways not yet described as underlying mineral amount. The results support the hypothesis that extracellular matrix interactions are the core regulator of mineral amount in muscle cells. Putative regulators described here add information to this hypothesis, expanding the knowledge on molecular relationships between gene expression and minerals.

Potential Biomarkers for Feed Efficiency-Related Traits in Nelore Cattle Identified by Co-expression Network and Integrative Genomics Analyses

Feed efficiency helps to reduce environmental impacts from livestock production, improving beef cattle profitability. We identified potential biomarkers (hub genes) for feed efficiency, by applying co-expression analysis in Longissimus thoracis RNA-Seq data from 180 Nelore steers. Six co-expression modules were associated with six feed efficiency-related traits (p-value ≤ 0.05). Within these modules, 391 hub genes were enriched for pathways as protein synthesis, muscle growth, and immune response. Trait-associated transcription factors (TFs) ELF1, ELK3, ETS1, FLI1, and TCF4, were identified with binding sites in at least one hub gene. Gene expression of CCDC80, FBLN5, SERPINF1, and OGN was associated with multiple feed efficiency-related traits (FDR ≤ 0.05) and were previously related to glucose homeostasis, oxidative stress, fat mass, and osteoblastogenesis, respectively. Potential regulatory elements were identified, integrating the hub genes with previous studies from our research group, such as the putative cis-regulatory elements (eQTLs) inferred as affecting the PCDH18 and SPARCL1 hub genes related to immune system and adipogenesis, respectively. Therefore, our analyses contribute to a better understanding of the biological mechanisms underlying feed efficiency in bovine and the hub genes disclosed can be used as biomarkers for feed efficiency-related traits in Nelore cattle.

Impact of inbreeding on milk fatty acids of a Brazilian Holstein cattle

Context Reproductive technologies may either contribute to the increase of genetic gains or inbreeding levels in dairy cattle breeding programs and herds. In general, studies have evaluated the impacts of inbreeding on traditional traits, such as production and reproduction. However, studies on novel traits, such as milk fatty acids profile, are scarce in the literature. Aim We aimed to evaluate the impact of inbreeding levels on milk fatty acids profile measured in a Brazilian Holstein cattle population. Methods Monthly records of the percentages (%) of palmitic, stearic, oleic, total saturated, unsaturated, monounsaturated and polyunsaturated fatty acids in milk were collected from 2047 Holstein cows between May and December 2012. The pedigree file contained 7963 animals over seven generations. The inbreeding trend was evaluated by linear regression of inbreeding coefficients on generations. The inbreeding impacts on breeding and phenotypic values were evaluated by linear regression of breeding values and phenotypic values on generations respectively. Breeding values were obtained by a Bayesian approach based on an animal model that included the systematic effect of contemporary group, cubic effect of days in milk as a covariate, and the random direct additive genetic and animal permanent environmental effects. Key results The average inbreeding coefficient was 0.28% for all animals, and 1.39% for inbred animals. Inbreeding tended to increase 0.24% per generation (P-value &lt;0.0001). The regression coefficients of phenotypic values on inbreeding were not significant for saturated, unsaturated, monounsaturated, palmitic, stearic and oleic fatty acids (P-values &gt;0.18, 0.10, 0.16, 0.26, 0.51 and 0.09 respectively); whereas for polyunsaturated fatty acids, they were significant (P-value &lt;0.04), with an increase of –0.000978% on phenotypic values for each 1% increase in inbreeding. The regression coefficients of breeding values on inbreeding were not significant for unsaturated, monounsaturated, polyunsaturated, stearic and oleic fatty acids (P-values &gt;0.33, 0.78, 0.48, 0.32 and 0.07 respectively). However, they were significant for saturated and palmitic fatty acids (P-values &lt;0.03 and 0.0001 respectively), presenting an increase of 0.000689% and 0.000283% respectively on breeding values for each 1% increase in inbreeding. Conclusions There was a significant increase in inbreeding coefficients over the studied generations; however, the impacts on both breeding and phenotypic values were small, when significant. Implications Our study may contribute to the understanding of how inbreeding levels may affect the fatty acids profile in a Brazilian Holstein cattle population, as well as encourage future research on this subject.

Statistical methods for identification of golden ratio

Several biological systems such as the biomechanics of human heart, locomotion, and phyllotaxis of plants present a harmonic behavior because their fractal structure are associated to the golden ratio. The golden ratio (Φ = 1.618033988749…), also known as Phi, golden mean, golden section or divine proportion, is an irrational constant found in various forms in nature and recently has been used in many health areas. However, there is no literature on a specific statistical test to identify the golden ratio structures. To validate the results from each survey, it is necessary that statistical techniques be correctly selected and implemented, and the absence of a test to identify the golden ratio may undermines the scientific papers which have this goal. Since the golden number is a ratio, some tests have been wrongly applied in its identification. The objective of this paper is to present and to evaluate methods for identification of golden ratio. Four tests were evaluated: t-Student with ratio statistic (TR), with delta statistic (TΔ), with difference statistic (TED), and Wilcoxon test with statistic difference (WD). Data simulating different samples sizes (n = 2-200) and variability scenarios were used. The tests were assessed regarding type I error rate and power. For TΔ, type I error rate increased along with sample size and variability, achieving 50% in the scenario of relative standard deviation of 12.5% and 20.0% for line segments of lengths a and b, and sample size equal 200. This test also showed lower power when compared to the others in all scenarios. Similarly, for TR, the type I error rate was sensitive to the increasing in sample size, varying from 5 to 60%. On the other hand, WD and TED were associated to low type I error rates (around 5%) and high power (6.1% for sample size equal 2-100% for sample size equal 200). The TΔ and TR were inadequate to identify the golden ratio, since they did not controlled the type I error rate and/or presented low power, leading to possible erroneous conclusions. Therefore WD and TED, both with statistical of difference, appeared as the most appropriate methods to test golden ratio structures.

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.

Linkage disequilibrium in Brazilian Santa Inês breed, Ovis aries

For genomic selection to be successful, there must be sufficient linkage disequilibrium between the markers and the causal mutations. The objectives of this study were to evaluate the extent of LD in ovine using the Santa Inês breed and to infer the minimum number of markers required to reach reasonable prediction accuracy. In total, 38,168 SNPs and 395 samples were used. The mean LD between adjacent marker pairs measured by r² and |D′| were 0.166 and 0.617, respectively. LD values between adjacent marker pairs ranged from 0.135 to 0.194 and from 0.568 to 0.650 for r² for |D′| across all chromosomes. The average r² between all pairwise SNPs on each chromosome was 0.018. SNPs separated by between 0.10 to 0.20 Mb had an estimated average r² equal to 0.1033. The identified haplotype blocks consisted of 2 to 21 markers. Moreover, estimates of average coefficients of inbreeding and effective population size were 0.04 and 96, respectively. LD estimated in this study was lower than that reported in other species and was characterized by short haplotype blocks. Our results suggest that the use of a higher density SNP panel is recommended for the implementation of genomic selection in the Santa Inês breed.

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 growth characteristics of free-range chickens under univariate random regression models

Repeated measures from the same individual have been analyzed by using repeatability and finite dimension models under univariate or multivariate analyses. However, in the last decade, the use of random regression models for genetic studies with longitudinal data have become more common. Thus, the aim of this research was to estimate genetic parameters for body weight of four experimental chicken lines by using univariate random regression models. Body weight data from hatching to 84 days of age (n = 34,730) from four experimental free-range chicken lines (7P, Caipirão da ESALQ, Caipirinha da ESALQ and Carijó Barbado) were used. The analysis model included the fixed effects of contemporary group (gender and rearing system), fixed regression coefficients for age at measurement, and random regression coefficients for permanent environmental effects and additive genetic effects. Heterogeneous variances for residual effects were considered, and one residual variance was assigned for each of six subclasses of age at measurement. Random regression curves were modeled by using Legendre polynomials of the second and third orders, with the best model chosen based on the Akaike Information Criterion, Bayesian Information Criterion, and restricted maximum likelihood. Multivariate analyses under the same animal mixed model were also performed for the validation of the random regression models. The Legendre polynomials of second order were better for describing the growth curves of the lines studied. Moderate to high heritabilities (h(2) = 0.15 to 0.98) were estimated for body weight between one and 84 days of age, suggesting that selection for body weight at all ages can be used as a selection criteria. Genetic correlations among body weight records obtained through multivariate analyses ranged from 0.18 to 0.96, 0.12 to 0.89, 0.06 to 0.96, and 0.28 to 0.96 in 7P, Caipirão da ESALQ, Caipirinha da ESALQ, and Carijó Barbado chicken lines, respectively. Results indicate that genetic gain for body weight can be achieved by selection. Also, selection for body weight at 42 days of age can be maintained as a selection criterion.

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.

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.

Correlated responses and genetic parameters for performance and carcass traits in a broiler line

The growth rate of broilers has triplicated in the last decades. The body weight is used as one of the selection criteria whereas the carcass traits are valuable market requirements. Thus, the meat industry like animals with high weights at slaughter and better carcass traits. However, the genetic relation of carcass traits with several body weights is unknown. Therefore, we established genetic associations among performance and carcass traits in a broiler chicken line and estimated genetic gain and trends. We also evaluated what age of selection would lead to a more efficient indirect selection of carcass traits. The data set with information of weights in different ages and carcass traits of 128,459 chickens was used. The pedigree data used contained 132,442 chickens. Genetic analysis were realized using ASREML® software applied a restricted maximum likelihood method. Heritability estimates ranged from moderate to high, which indicates that these traits can have high selection response. Genetic correlations between performance and carcass traits varied from moderate to high, which indicates the presence of a genetic association whereas genetic trends indicated that direct selection is occurring for body weight at different ages. Theselection at 30 and 38 days should be considered instead of the slaughter weight, as it anticipates selection in around 12 days.

Concordance analysis between estimation methods of milk fatty acid content

Considering the milk fatty acid influence on human health, the aim of this study was to compare gas chromatography (GC) and Fourier transform infrared (FTIR) spectroscopy for the determination of these compounds. Fatty acid content (g/100g of fat) were obtained by both methods and compared through Pearson's correlation, linear Bayesian regression, and the Bland-Altman method. Despite the high correlations between the measurements (r=0.60-0.92), the regression coefficient values indicated higher measures for palmitic acid, oleic acid, unsaturated and monounsaturated fatty acids and lower values for stearic acid, saturated and polyunsaturated fatty acids estimated by GC in comparison to FTIR results. This inequality was confirmed in the Bland-Altman test, with an average bias varying from -8.65 to 6.91g/100g of fat. However, the inclusion of 94% of the samples into the concordance limits suggested that the variability of the differences between the methods was constant throughout the range of measurement. Therefore, despite the inequality between the estimates, the methods displayed the same pattern of milk fat composition, allowing similar conclusions about the milk samples under evaluation.

Definição de grupos genéticos aditivos visando melhor predição de valores genéticos em bovinos de corte

O objetivo do estudo foi avaliar a estrutura de grupos genéticos aditivos como uma alternativa para melhor ajuste de animais com paternidade desconhecida nas avaliações genéticas. As características estudadas foram: peso a desmama; peso ao sobreano; ganho de peso pós-desmama em 345 dias; perímetro escrotal ao sobreano e o escore visual de musculosidade ao sobreano. Um banco em que todos os animais tinham paternidade conhecida foi utilizado como grupo controle. A partir deste, três cenários foram criados em que 30; 50 ou 70% dos indivíduos foram aleatoriamente assumidos como animais com paternidade desconhecida. As estratégias consideradas de grupos genéticos aditivos foram: ano de nascimento do animal com paternidade desconhecida; fazenda de nascimento do animal com paternidade desconhecida; ano de nascimento e fazenda de nascimento do animal e um grupo controle negativo que manteve os animais com pais desconhecidos. A estratégia adequada de grupo genético aditivo foi escolhida como sendo aquela que resultou em um maior coeficiente de regressão e com valor genético aditivo mais próximos ao predito para os animais do grupo controle. Os resultados que incluíram a estratégia de grupo genético aditivo nos modelos de predição dos valores genéticos mostraram maior coerência com o grupo controle, frente a não inclusão destes nas análises. Dentre as estratégias propostas, ano de nascimento e fazenda de nascimento do animal e apenas ano de nascimento do animal com paternidade desconhecida obtiveram os melhores resultados.

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.

The effect of using "race of child" instead of "race of mother" on the black-white gap in infant mortality due to birth defects

OBJECTIVE

For at least 20 years, birth defects have been the leading cause of infant mortality in the United States. Some studies have reported higher rates for black infants than white infants of mortality due to birth defects, while other studies have reported no black-white differences. The authors analyzed the effect on these rates of a change in the way the National Center for Health Statistics (NCHS) tabulates "race" for newborns.

METHODS

The authors calculated infant mortality rates due to birth defects for 1980-1993 using two standard methods of assigning newborns to "racial" categories: a "race of child" algorithm and the "race of mother" approach currently used by NCHS.

RESULTS

From 1980 through 1993, birth defect-specific infant mortality rates (BD-IMRs) were significantly higher for black infants than white infants 12 of the 14 years by "race of mother" and only 5 of 14 years by "race of child." Calculation of BD-IMRs by "race of mother" reduced the rate for white infants and increased the rate for black infants in each of the 14 years. The choice of method for assigning newborns to "racial" categories had a progressively greater effect over time on the black-white gap in BD-IMRs.

CONCLUSIONS

Calculations of trends in "race"-specific BD-IMRs by may vary substantially by whether "race of mother" or "race of child" is used. Identifying the method of tabulation is imperative for appropriate comparisons and interpretations.

Signaling function of reconstituted CD16: zeta: gamma receptor complex isoforms

Natural killer cells express an Fc receptor for IgG (CD16) in association with disulfide-linked dimers composed of two homologous subunits: the zeta chain of the T cell antigen receptor complex and the gamma chain of the mast cell/basophil Fc receptor for IgE. The ability of zeta and gamma to transduce CD16-mediated activation signals was compared by reconstituting distinct CD16 receptor isoforms composed of various combinations of zeta- and gamma-containing dimers. Stably transformed non-hematopoietic and hematopoietic cell lines were established that expressed chimeric molecules comprising the extracellular domain of CD16 joined to the transmembrane and intracellular domains of zeta or gamma. Reconstituted CD16 receptor complexes triggered Ca2+ influx, tyrosine phosphorylation, and IL-2 production in stable transformants of the Jurkat T cell line. However, cross-linking of the CD16/gamma chimera induced a specific pattern of tyrosine phosphorylation and was more efficient at signal transduction than a CD16, zeta-zeta complex, suggesting that zeta and gamma cytoplasmic domains may be coupled to distinct tyrosine kinase pathways that differentially regulate CD16-mediated activation signals. By contrast, both CD16/zeta and CD16/gamma chimeric molecules were not functional in stable transformants of the fibroblast Chinese Hamster Ovary cell line, indicating a requirement for downstream signaling components present in hematopoietic cells. Finally, the zeta transmembrane domain appears to preferentially associate with CD16 rather than the CD3:TCR complex, suggesting that a hierarchy of molecular interactions governs NK and T cell differentiation.

Patterns and extent of isotype-specificity in the murine H chain switch DNA rearrangement

We have analyzed the configuration of the H chain locus of 41 hybridomas by Southern blot analysis. Each H chain switch region was determined to be germ line, rearranged, or deleted. Including 13 previously analyzed hybridomas, 60% of those with rearrangements on both alleles showed a correlation of the two alleles, i.e., both the expressed and the nonexpressed alleles have rearranged to the same H chain constant region gene segment. When the two H chain alleles did not rearrange to the same gene, they often rearranged to neighboring H chain genes. These results support a role for isotype-specific factors in H chain switch recombination. The action of these isotype-specific factors may be propagated to some extent along the chromosome, which would lead to rearrangements to neighboring genes.

Patterns and extent of isotype-specificity in the murine H chain switch DNA rearrangement

We have analyzed the configuration of the H chain locus of 41 hybridomas by Southern blot analysis. Each H chain switch region was determined to be germ line, rearranged, or deleted. Including 13 previously analyzed hybridomas, 60% of those with rearrangements on both alleles showed a correlation of the two alleles, i.e., both the expressed and the nonexpressed alleles have rearranged to the same H chain constant region gene segment. When the two H chain alleles did not rearrange to the same gene, they often rearranged to neighboring H chain genes. These results support a role for isotype-specific factors in H chain switch recombination. The action of these isotype-specific factors may be propagated to some extent along the chromosome, which would lead to rearrangements to neighboring genes.

Products and implied mechanism of H chain switch recombination

The Ig H chain switch is a DNA recombination event. The recombination occurs between two or more switch regions, areas of tandem sequence duplication that lie upstream of the corresponding H chain C region genes. We have determined the DNA sequence at four recombination sites in three molecularly cloned, rearranged switch regions. All eight donor and recipient recombination sites are at the common pentamers GGGGT, GAGCT, and GGTGG. One of the switch recombination events is an inversion of S gamma 3 sequences. Another of the recombinational events is an internal S gamma 1 deletion, which may be switch enzyme mediated. These results, together with other switch recombination site sequences, suggest that switch recombination is mediated by cutting enzymes with modest specificity and religation enzymes with no specificity.

Products and implied mechanism of H chain switch recombination

The Ig H chain switch is a DNA recombination event. The recombination occurs between two or more switch regions, areas of tandem sequence duplication that lie upstream of the corresponding H chain C region genes. We have determined the DNA sequence at four recombination sites in three molecularly cloned, rearranged switch regions. All eight donor and recipient recombination sites are at the common pentamers GGGGT, GAGCT, and GGTGG. One of the switch recombination events is an inversion of S gamma 3 sequences. Another of the recombinational events is an internal S gamma 1 deletion, which may be switch enzyme mediated. These results, together with other switch recombination site sequences, suggest that switch recombination is mediated by cutting enzymes with modest specificity and religation enzymes with no specificity.

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.

Palliative bipolar electrocoagulation therapy of obstructing esophageal cancer

In a multicenter pilot study, a prototype bipolar electrocoagulation tumor probe was employed for palliation of obstructing circumferential esophageal cancer in 20 patients. Mean number of initial treatment sessions was 1.7. Dysphagia and tumor channel size improved significantly after treatment. Mean treatment interval before repeat treatment was 7.6 weeks. Major complications included delayed hemorrhage (two patients) and esophageal-pulmonary fistula (two patients). This new device may provide a less expensive alternative to laser or surgery for palliation of malignant esophageal obstruction.

The immunoglobulin heavy chain switch: structural features of gamma 1 recombinant switch regions

The immunoglobulin heavy chain isotype switch is mediated by a DNA rearrangement involving specific genomic segments referred to as switch regions. Switch regions are composed of tandemly repeated simple sequences. The role of the tandemly repeated structure of switch regions in the switch recombination process is not understood. We mapped eight recombination sites--six in the gamma 1 and two in the gamma 3 tandem arrays. In addition, we obtained molecular clones representing three of the six gamma 1 rearrangements, and determined the nucleotide sequences of the recombination sites in each. In general, the rearrangements are confined to the tandem repeat units, and are not clustered in a particular portion of either the gamma 3 or gamma 1 switch region. Nucleotide sequence analysis of one of the recombinant clones, gamma M35, reveals evidence for a successive switch event wherein a recombination between S mu and S gamma 3 was followed by recombination 57 bp downstream with S gamma 1. gamma 1 sequence data from the molecular clones we obtained, together with similar data from other investigators regarding the gamma 1, gamma 2b, and gamma 2a switch regions, reveals that recombinations tend to occur at homologous positions of the respective gamma-unit repeats, adjacent to the elements AGCT and GGGG found in each. This finding suggests that the cutting and religation step of the recombination process is mediated by a recombinase common to the four gamma-isotypes.

The immunoglobulin heavy chain switch: structural features of gamma 1 recombinant switch regions

The immunoglobulin heavy chain isotype switch is mediated by a DNA rearrangement involving specific genomic segments referred to as switch regions. Switch regions are composed of tandemly repeated simple sequences. The role of the tandemly repeated structure of switch regions in the switch recombination process is not understood. We mapped eight recombination sites--six in the gamma 1 and two in the gamma 3 tandem arrays. In addition, we obtained molecular clones representing three of the six gamma 1 rearrangements, and determined the nucleotide sequences of the recombination sites in each. In general, the rearrangements are confined to the tandem repeat units, and are not clustered in a particular portion of either the gamma 3 or gamma 1 switch region. Nucleotide sequence analysis of one of the recombinant clones, gamma M35, reveals evidence for a successive switch event wherein a recombination between S mu and S gamma 3 was followed by recombination 57 bp downstream with S gamma 1. gamma 1 sequence data from the molecular clones we obtained, together with similar data from other investigators regarding the gamma 1, gamma 2b, and gamma 2a switch regions, reveals that recombinations tend to occur at homologous positions of the respective gamma-unit repeats, adjacent to the elements AGCT and GGGG found in each. This finding suggests that the cutting and religation step of the recombination process is mediated by a recombinase common to the four gamma-isotypes.