Phenotypic variation seems not to be associated with the genetic profile in Zygopetalum (Orchidaceae): a case study of a high-elevation rocky complex.

BACKGROUND: Hybridization associated with polyploidy studies is rare in the tropics. The genus Zygopetalum (Orchidaceae) was investigated here as a case study of Neotropical plants. In the rocky highlands of the Ibitipoca State Park (ISP), southeast Brazil, individuals with intermediate colors and forms between the species Z. maculatum and Z. triste were commonly identified. METHODS AND RESULTS: Chromosomal analysis and DNA quantity showed a uniform population. Regardless of the aspects related to the color and shape of floral structures, all individuals showed 2n = 96 chromosomes and an average of 14.05 pg of DNA. Irregularities in meiosis associated with chromosome number and C value suggest the occurrence of polyploidy. The genetic distance estimated using ISSR molecular markers revealed the existence of genetic variability not related to morphological clusters. Morphometric measurements of the flower pieces revealed that Z. maculatum shows higher variation than Z. triste although lacking a defined circumscription. CONCLUSION: The observed variation can be explained by the polyploid and phenotypic plasticity resulting from the interaction of the genotypes with the heterogeneous environments observed in this habitat.

Resequencing of Brazilian locally adapted cattle breeds revealed variants in candidate genes and transcription factors for meat fatty acid profile.

The beef cattle industry has experienced a shift driven by a market demand for healthier meat, cost efficiency and environmental sustainability in recent years. Consequently, there has been a growing focus on the fatty acids content and functions of meat in cattle breeding programmes. Besides, a deeper understanding of the biological mechanisms influencing the expression of different phenotypes related to fatty acid profiles is crucial. In this study, we aimed to identify Single-Nucleotide Variants (SNV) and Insertion/Deletion (InDels) DNA variants in candidate genes related to fatty acid profiles described in genomic, transcriptomic and proteomic studies conducted in beef cattle breeds. Utilizing whole-genome re-sequencing data from Brazilian locally adapted bovine breeds, namely Caracu and Pantaneiro, we identified SNVs and InDels associated with 23,947 genes. From these, we identified 318 candidate genes related to fatty acid profiles that contain variants. Subsequently, we select only genes with SNVs and InDels in their promoter, 5' UTR and coding region. Through the gene-biological process network, approximately 19 genes were highlighted. Furthermore, considering the studied trait and a literature review, we selected the main transcription factors (TF). Functional analysis via gene-TF network allowed us to identify the 30 most likely candidate genes for meat fatty acid profile in cattle. LIPE, MFSD2A and SREBF1 genes were highlighted in networks due to their biological importance. Further dissection of these genes revealed 15 new variants found in promoter regions of Caracu and Pantaneiro sequences. The gene networks facilitated a better functional understanding of genes and TF, enabling the identification of variants potentially related to the expression of candidate genes for meat fatty acid profiles in cattle.

Single-step genome-wide association studies and post-GWAS analyses for the number of oocytes and embryos in Gir cattle.

Genome-Wide Association Studies (GWAS) are used for identification of quantitate trait loci (QTL) and genes associated with several traits. We aimed to identify genomic regions, genes, and biological processes associated with number of total and viable oocytes, and number of embryos in Gir dairy cattle. A dataset with 17,526 follicular aspirations, including the following traits: number of viable oocytes (VO), number of total oocytes (TO), and number of embryos (EMBR) from 1641 Gir donors was provided by five different stock farms. A genotype file with 2093 animals and 395,524 SNP markers was used to perform a single-step GWAS analysis for each trait. The top 10 windows with the highest percentage of additive genetic variance explained by 100 adjacent SNPs were selected. The genomic regions identified in our work were overlapped with QTLs from QTL database on chromosomes 1, 2, 5, 6, 7, 8, 9, 13, 17, 18, 20, 21, 22, 24, and 29. These QTLs were classified as External, Health, Meat and carcass, Production or Reproduction traits, and about 38% were related to Reproduction. In total, 117 genes were identified, of which 111 were protein-coding genes. Exclusively associations were observed for 42 genes with EMBR, and 1 with TO. Also, 42 genes were in common between VO and TO, 28 between VO and EMBR, and four genes were in common among all traits. In conclusion, great part of the identified genes plays a functional role in initial embryo development or general cell functions. The protein-coding genes ARNT, EGR1, HIF1A, AHR, and PAX2 are good markers for the production of oocytes and embryos in Gir cattle.

Runs of homozygosity and signatures of selection for number of oocytes and embryos in the Gir Indicine cattle.

Runs of homozygosity (ROH) and signatures of selection are the results of selection processes in livestock species that have been shown to affect several traits in cattle. The aim of the current work was to verify the profile of ROH and inbreeding depression in the number of total (TO) and viable oocytes (VO) and the number of embryos (EMBR) in Gir Indicine cattle. In addition, we aim to identify signatures of selection, genes, and enriched regions between Gir subpopulations sorted by breeding value for these traits. The genotype file contained 2093 animals and 420,718 SNP markers. Breeding values used to sort Gir animals were previously obtained. ROH and signature of selection analyses were performed using PLINK software, followed by ROH-based (FROH) and pedigree-based inbreeding (Fped) and a search for genes and their functions. An average of 50 ± 8.59 ROHs were found per animal. ROHs were separated into classes according to size, ranging from 1 to 2 Mb (ROH1-2Mb: 58.17%), representing ancient inbreeding, ROH2-4Mb (22.74%), ROH4-8Mb (11.34%), ROH8-16Mb (5.51%), and ROH>16Mb (2.24%). Combining our results, we conclude that the increase in general FROH and Fped significantly decreases TO and VO; however, in different chromosomes traits can increase or decrease with FROH. In the analysis for signatures of selection, we identified 15 genes from 47 significant genomic regions, indicating differences in populations with high and low breeding value for the three traits.

Assessment of copy number variants in three Brazilian locally adapted cattle breeds using whole-genome re-sequencing data.

Further characterization of genetic structural variations should strongly focus on small and endangered local breeds given their role in unraveling genes and structural variants underlying selective pressures and phenotype variation. A comprehensive genome-wide assessment of copy number variations (CNVs) based on whole-genome re-sequencing data was performed on three Brazilian locally adapted cattle breeds (Caracu Caldeano, Crioulo Lageano, and Pantaneiro) using the ARS-UCD1.2 genome assembly. Data from 36 individuals with an average coverage depth of 14.07× per individual was used. A total of 24 945 CNVs were identified distributed among the breeds (Caracu Caldeano = 7285, Crioulo Lageano = 7297, and Pantaneiro = 10 363). Deletion events were 1.75-2.07-fold higher than duplications, and the total length of CNVs is composed mostly of a high number of segments between 10 and 30 kb. CNV regions (CNVRs) are not uniformly scattered throughout the genomes (n = 463), and 105 CNVRs were found overlapping among the studied breeds. Functional annotation of the CNVRs revealed variants with high consequence on protein sequence harboring relevant genes, in which we highlighted the BOLA-DQB, BOLA-DQA5, CD1A, ß-defensins, PRG3, and ULBP21 genes. Enrichment analysis based on the gene list retrieved from the CNVRs disclosed over-represented terms (p < 0.01) strongly associated with immunity and cattle resilience to harsh environments. Additionally, QTL associated with body conformation and dairy-related traits were also unveiled within the CNVRs. These results provide better understanding of the selective forces shaping the genome of such cattle breeds and identify traces of natural selection pressures by which these populations have been exposed to challenging environmental conditions.

Association between dietary creatine and visuospatial short-term memory in older adults.

Aims: The purpose was to examine the relationship between habitual dietary creatine intake obtained in food and visuospatial short-term memory (VSSM). Methods: Forty-two participants (32 females, 10 males; > 60 yrs of age) completed a 5-day dietary recall to estimate creatine intake and performed a cognitive assessment which included a visuospatial short-term memory test (forward and reverse corsi block test) and a mini-mental state examination (MMSE). Pearson correlation coefficients were determined. Further, cohorts were derived based on the median creatine intake. Results: There was a significant correlation between the forward Corsi (r = 0.703, P < 0.001), reverse Corsi (r = 0.715, P < 0.001), and the memory sub-component of the MMSE (r = 0.406, P = 0.004). A median creatine intake of 0.382 g/day was found. Participants consuming greater than the median had a significantly higher Corsi (P = 0.005) and reverse Corsi (P < 0.001) scores compared to participants ingesting less than the median. Conclusions: Dietary creatine intake is positively associated with measures of memory in older adults. Clinical Implications: Older adults should consider food sources containing creatine (i.e. red meat, seafood) due to the positive association with visuospatial short-term memory.

Aneuploids and its increment on diversity of Lippia alba polyploid complex: genetic aspects and origin.

BACKGROUND: Aneuploidy is associated with add or lack of individual chromosomes. The knowledge regarding aneuploidy is still rare in wild and tropical populations. Lippia alba is a tropical polyploid complex naturally formed, with 2x, 3x, 4x, 6x, and aneuploid individuals. The species presents pharmacological and medicinal importance, due to its essential oil compounds, which are related to the ploidal level. Considering the singularity of aneuploids emergence and stability, we proposed to investigate putative cytotypes involved in the aneuploids formation. METHODS AND RESULTS: Molecular, cytogenetic, reproductive, and chemical approaches were adopted. The results showed that the aneuploids possibly have independent origin considering the genetic, chemical and karyotypical profiles. The chemical composition of aneuploids is related to genetic similarity. The aneuploid origin may involve 2x and 3x cytotypes being possible to rise four scenarios of crosses to explain that. CONCLUSIONS: The results, in general, contribute to the comprehension of the origin of aneuploids and highlight the genetic profile of these accessions as a key element on the understanding of the chemical profile of L. alba accessions.

Could dietary creatine intake modulate overweight elderly's selective attention and inhibitory function?

Aims: There is evidence that both aging and increased adiposity may impact creatine levels in the brain, and brain creatine levels are important for cognition. The aim of this study was to assess correlation between dietary creatine intake and cognition in in elderly women with overweight. Methods: Twenty seven overweight women over 60 years of age who were part of a larger study participated in an Eriksen Flanker Task (EFT) to asssess cognitive performance. Additionally, diet was assessed over 5 days via daily diary nutritional recalls and the estimate of the daily amount of creatine was calculated. Results: In the EFT when incongruente stimulus were presented there was a significant diferences between those with low and high intake of creatine (-35.3 ± 5.84; p < 0.001). Similarly, reaction time to answer incongruent stimulus (r = -0.383; p = 0.004) and the percent of correct answers (r = 0.743; p < 0.001) showed weak to strong correlations with self-reported daily creatine intake. Conclusions: In conclusion, our results suggest that in elderly women with overweight that dietary intake of creatine may influence cognitive ability. Clinical Implications: Our findings support the idea that intake of dietary creatine may be an important factor for cognition in older adults.

Psychological effect of acute creatine pre-workout supplementation induces performance improvement in resistance exercise.

The purpose of this study was to test whether believed versus actual acute creatine ingestion impacted resistance exercise performance. Fifteen men (21.9 ± 2.7 years old) completed four bouts of three sets each of squat and bench press to volitional fatigue at a 10RM load with 1-min between-sets rest interval. Thirty minutes prior to each exercise bout, they received the following treatments in a randomized order: 1) nothing (CON); 2) 0.3 g·kg-1 dextrose placebo (PLC); 3) 0.3 g·kg-1 dextrose, identified as creatine (Cr-False); 4) 0.3 g·kg 20 -1 creatine, identified as creatine (CrTrue). Between-treatments comparisons included the total repetitions completed and the rate of perceived exertion. Results revealed (p < 0.05) higher repetitions performed for all treatments versus CON for both squat and bench press. In the squat, more repetitions were performed with Cr-True (p < 0.001) and CrFalse (p < 0.001) than with either CON or PLC. Bayes Factor analyses revealed strong (PLC to Cr-True BF = 19.1) and very strong (PLC to CrFalse BF = 45.3) posterior probability favouring positive effects for both "creatine" conditions over PLC for the squat. In conclusion, in acute measures, belief versus ingestion of creatine yields similar exercise performance.

Caffeine Supplementation or Carbohydrate Mouth Rinse Improves Performance.

Training volume is one of the critical variables required to promote resistance training benefits (e. g., hypertrophy, muscular strength). Thus, strategies to improve training volume are required. We tested the hypothesis that there is an increase in performance and reduction in the rate of perceived exertion in strength training with caffeine supplementation, carbohydrate mouth rinse, and a synergistic effect of caffeine supplementation plus carbohydrate mouth rinse. We recruited 29 physically active women: 24±4 years, 60.0±7.9 kg, 161.0±6.0 cm. This study was a double-blind, randomized, placebo-controlled, and crossover one. The subjects performed all sessions of strength training under different conditions: caffeine (6.5 mg·kg-1 body mass); carbohydrate (6 g of maltodextrin rinsed for 10 s); caffeine plus carbohydrate; or placebo. We applied the ANOVA for repeated measures through the null-hypothesis statistical test, and the Bayes factors analyses approach. The subjects showed significant improvement in the total repetitions (P<0.01; BF10 ~99%) for both lower and upper limb exercises in all conditions compared to placebo without difference among interventions. The rate of perceived exertion (P>0.05) did not differ among interventions. However, BF10 analyses showed a higher probability of lower RPE for CAF intervention. We conclude that either caffeine supplementation or carbohydrate mouth rinse can improve performance in resistance training. There is no additional effect of caffeine plus carbohydrate.

Going further post-RNA-seq: In silico functional analyses revealing candidate genes and regulatory elements related to mastitis in dairy cattle.

This study aimed to obtain a better understanding of the regulatory genes and molecules involved in the development of mastitis. For this purpose, the transcription factors (TF) and MicroRNAs (miRNA) related to differentially expressed genes previously found in extracorporeal udders infected with Streptococcus agalactiae were investigated. The Gene-TF network highlighted LOC515333, SAA3, CD14, NFKBIA, APOC2 and LOC100335608 and genes that encode the most representative transcription factors STAT3, PPARG, EGR1 and NFKB1 for infected udders. In addition, it was possible to highlight, through the analysis of the gene-miRNA network, genes that could be post-transcriptionally regulated by miRNAs, such as the relationship between the CCL5 gene and the miRNA bta-miR-363. Overall, our data demonstrated genes and regulatory elements (TF and miRNA) that can play an important role in mastitis resistance. The results provide new insights into the first functional pathways and the network of genes that orchestrate the innate immune responses to infection by Streptococcus agalactiae. Our results will increase the general knowledge about the gene networks, transcription factors and miRNAs involved in fighting intramammary infection and maintaining tissue during infection and thus enable a better understanding of the pathophysiology of mastitis.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr.

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data.

BACKGROUND: The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics. RESULTS: A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds. CONCLUSION: The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.

High-fat diet-induced obesity worsens TH2 immune response and immunopathologic characteristics in murine model of eosinophilic oesophagitis.

BACKGROUND: Eosinophilic oesophagitis (EoE) is an emergent chronic immune-mediated disease of the oesophagus, which affects both children and adults. It is clinically characterized by dysphagia, food impaction and oesophageal eosinophilia. Epidemiological studies indicate that obesity can worsen allergic symptoms; however, its effect on EoE immunopathological response has not been evaluated yet. This study aimed to assess the effect of obesity on allergic inflammation and T helper-2 profile in an EoE experimental model. METHODS: Obesity was induced by high-fat feeding. After 7 weeks of diet, male BALB/c mice were subcutaneously sensitized and orally challenged with OVA. RESULTS: Obesity itself induced a significant mast cell and eosinophil accumulation in the oesophagus, trachea, gut and lung. After allergy induction, this number was higher, when compared to lean-allergic mice. Moreover, obese-allergic mice showed higher remodelling area, in the oesophagus, associated with higher IL-5 and TSLP mRNA expression. In contrast, FoxP3 and IL-10 were less expressed in comparison with lean-allergic mice. In addition, the amount of CD11c+ MHCII+ PDL1+ dendritic cells was reduced, while the number of CD11c+ MHCII+ CD80+ DCs and CD3+ CD4+ GATA3+ IL-4+ cells was increased in obese-allergic mice in the spleen and lymph nodes when compared to lean-allergic mice. CONCLUSION: Obesity aggravated the immune histopathological characteristics in the EoE experimental model, which was associated with the reduction in the regulatory profile, and the increased inflammatory cells influx, related to the TH 2 profile. Altogether, the data provide new knowledge about obesity as a risk factor, worsening EoE symptoms, and contribute for future treatment strategies for this specific profile.

Genetic relationships and polyploid origins in the Lippia alba complex.

PREMISE: Plant genomes vary in size and complexity due in part to polyploidization. Latitudinal analyses of polyploidy are biased toward floras of temperate regions, with much less research done in the tropics. Lippia alba has been described as a tropical polyploid complex with diploid, triploid, tetraploid, and hexaploid accessions. However, no data regarding relationships among the ploidal levels and their origins have been reported. Our goals are to clarify the relationships among accessions of Lippia alba and the origins of each ploidal level. METHODS: We investigated 98 samples representing all five geographical regions of Brazil and all ploidal levels using microsatellite (SSR) allelic variation and DNA sequences of ITS and trnL-F. Nine morphological structures were analyzed from 33 herbarium samples, and the chemical compounds of 78 accessions were analyzed by GC-MS. RESULTS: Genetic distance analysis, the alignment block pattern, as well as RAxML and Bayesian trees showed that accessions grouped by ploidal level. The triploids form a well-defined group that originated from a single group of diploids. The tetraploids and hexaploid grouped together in SSR and trnL-F analyses. The recovered groups agree with chemical data and morphology. CONCLUSIONS: The accessions grouped by ploidal level. Only one origin of triploids from a single group of diploids was observed. The tetraploid origin is uncertain; however, it appears to have contributed to the origin of the hexaploid. This framework reveals linkages among the ploidal levels, providing new insights into the evolution of a polyploid complex of tropical plants.

Development of microsatellite markers for Lippia alba and related Lippia species.

Microsatellite primers were developed in Lippia alba complex to better understanding the origins and evolution of the species. We sought to increase the numbers of available simple sequence repeat (SSR) markers. We performed low-coverage (~ twofold) genomic DNA sequencing of a diploid accession and generated a de novo assembly comprising 175,572 contigs. Sixteen SSR loci were selected and of these 13 SSR loci were successfully amplified in 20 L. alba tetraploid accessions and in 12 other Lippia species. Only one SSR locus was monomorphic, whereas 12 loci were polymorphic, yielding one to nine alleles. The heterozygosity was similar among markers, with values of 0.274-0.485; the polymorphism information content values varied from 0.237 to 0.367. These markers were successfully amplified in related species with 85% of transferability on average. Thus, we demonstrate the utility of including a de novo assembly step to obtain SSR markers from low-coverage genomic datasets.

Single-step genome-wide association studies (GWAS) and post-GWAS analyses to identify genomic regions and candidate genes for milk yield in Brazilian Girolando cattle.

Milk production is economically important to the Brazilian agribusiness, and the majority of the country's milk production derives from Girolando (Gir × Holstein) cows. This study aimed to identify quantitative trait loci (QTL) and candidate genes associated with 305-d milk yield (305MY) in Girolando cattle. In addition, we investigated the SNP-specific variances for Holstein and Gir breeds of origin within the sequence of candidate genes. A single-step genomic BLUP procedure was used to identify QTL associated with 305MY, and the most likely candidate genes were identified through follow-up analyses. Genomic breeding values specific for Holstein and Gir were estimated in the Girolando animals using a model that uses breed-specific partial relationship matrices, which were converted to breed of origin SNP effects. Differences between breed of origin were evaluated by comparing estimated SNP variances between breeds. From 10 genome regions explaining most additive genetic variance for 305MY in Girolando cattle, 7 candidate genes were identified on chromosomes 1, 4, 6, and 26. Within the sequence of these 7 candidate genes, Gir breed of origin SNP alleles showed the highest genetic variance. These results indicated QTL regions that could be further explored in genomic selection panels and which may also help in understanding the gene mechanisms involved in milk production in the Girolando breed.

Spatial dynamics of the Girolando breed in Brazil: analysis of genetic integration and environmental factors.

Brazil is one of the world's largest milk producers. Several scientific studies have been developed related to landscape analyses that combine genetic with landscape structure data. In the present study, we aimed to analyze the relationship between genetic, environmental, and socioeconomic aspects of production in Girolando cattle in Brazil, as well as verify the spatial patterns of its genetic diversity. Genetic values and accuracy of 46,289 animals were used as well as information from DNA of 310 Girolando animals. Canonic, discriminant, and cluster analyses were conducted in SAS® and K-means method in ArcGIS 10.3 software. The relationship between genetic and geographic distance was analyzed using different methods in software Alleles in Space®. Clusters with animals with higher genetic values for milk production are located in municipalities with lower gross domestic product, fewer family-based establishments, and lower human development index. These clusters are associated with regions with higher area planted with crops, lower percentage of pastures that were less degraded, higher humidity, lower temperature range, and lower normalized difference vegetation index (NDVI) values. The greater the geographical distance between groups of animals, the greater the genetic distance between them with a significant distinction over 504 km. There is high genetic heterogeneity among animals. From these results, it will be possible to develop methodologies for better evaluation of the animals within the production systems.

Genome-wide association studies for heat stress response in Bos taurus × Bos indicus crossbred cattle.

Heat stress is an important issue in the global dairy industry. In tropical areas, an alternative to overcome heat stress is the use of crossbred animals or synthetic breeds, such as the Girolando. In this study, we performed a genome-wide association study (GWAS) and post-GWAS analyses for heat stress in an experimental Gir × Holstein F2 population. Rectal temperature (RT) was measured in heat-stressed F2 animals, and the variation between 2 consecutive RT measurements (ΔRT) was used as the dependent variable. Illumina BovineSNP50v1 BeadChip (Illumina Inc., San Diego, CA) and single-SNP approach were used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene-transcription factor (TF) networks, generated from enriched TF. The breed origin of marker alleles in the F2 population was assigned using the breed of origin of alleles (BOA) approach. Heritability and repeatability estimates (± standard error) for ΔRT were 0.13 ± 0.08 and 0.29 ± 0.06, respectively. Association analysis revealed 6 SNP significantly associated with ΔRT. Genes involved with biological processes in response to heat stress effects (LIF, OSM, TXNRD2, and DGCR8) were identified as putative candidate genes. After performing the BOA approach, the 10% of F2 animals with the lowest breeding values for ΔRT were classified as low-ΔRT, and the 10% with the highest breeding values for ΔRT were classified as high-ΔRT. On average, 49.4% of low-ΔRT animals had 2 alleles from the Holstein breed (HH), and 39% had both alleles from the Gir breed (GG). In high-ΔRT animals, the average proportion of animals for HH and GG were 1.4 and 50.2%, respectively. This study allowed the identification of candidate genes for ΔRT in Gir × Holstein crossbred animals. According to the BOA approach, Holstein breed alleles could be associated with better response to heat stress effects, which could be explained by the fact that Holstein animals are more affected by heat stress than Gir animals and thus require a genetic architecture to defend the body from the deleterious effects of heat stress. Future studies can provide further knowledge to uncover the genetic architecture underlying heat stress in crossbred cattle.

Assessment of runs of homozygosity islands and estimates of genomic inbreeding in Gyr (Bos indicus) dairy cattle.

BACKGROUND: Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (FROH), genomic relationship matrix approach (FGRM) and based on the observed versus expected number of homozygous genotypes (FHOM), and from pedigree-based coefficient (FPED). RESULTS: ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH1-2 Mb) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of FROH, FGRM, FHOM, and FPED approaches. FPED estimates ranged from 0.00 to 0.327 and FROH from 0.001 to 0.201. Low to moderate correlations were observed between FPED-FROH and FGRM-FROH, with values ranging from -0.11 to 0.51. Low to high correlations were observed between FROH-FHOM and moderate between FPED-FHOM and FGRM-FHOM. Correlations between FROH from different lengths and FPED gradually increased with ROH length. CONCLUSIONS: Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low FPED-FROH correlations for small segments indicate that FPED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that FROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.