Identifying Hub Genes for Heat Tolerance in Water Buffalo (Bubalus bubalis) Using Transcriptome Data

Estimation of the threshold for heat stress and genetic features for milk yield in Mehsana buffaloes in India

Heat stress is one of the primary environmental factors that harm both the productivity and health of buffaloes. The current study was conducted to estimate the threshold of temperature humidity index (THI)1 and genetic features for milk yield of first-lactation Mehsana buffaloes using an univariate repeatability test-day model. The data included 130,475 first lactation test-day milk yield (FLTDMY) records of 13,887 Mehsana buffaloes and the daily temperature and humidity. The statistical model included herd test day as fixed effects, days-in-milk (DIM) classes, age of the animal, as well as random factors such as the additive genetic effect (AGE) of animal in general conditions (intercept), AGE of the buffaloes subjected to heat stress (slope), permanent environmental effect of animal in general conditions (intercept), permanent environmental effect of animal under heat stress conditions (slope) and random residual effect. It was expected that the general effects and the heat-tolerance effects would be correlated, represented by the present investigation's repeatability models. The variance components of FLTDMY in the present study were computed using the REML method. The threshold for THI was 78. At the THI below the threshold, the heritability estimated for the FLTDMY trait was 0.29, and the additive genetic variance (AGV) for heat stress conditions was 0. At THI of 83, AGV for heat stress conditions was highest for FLTDMY. The genetic correlation of general AGE to heat-tolerant AGE was -0.40. The results indicated that a consistent selection for milk production, avoiding the thermal tolerance, may diminish the thermal tolerance capacity of Mehsana buffaloes.

Skin transcriptome analysis in Brangus cattle under heat stress

Heat stress is a major factor that negatively affects animal welfare and production systems. Livestock should adapt to tropical and subtropical areas and to meet this, composite breeds have been developed. This work aimed to evaluate gene expression profiles in the skin of Brangus cattle under heat stress using a case-control design, and to correlate this with skin histological characteristics. Two groups of bulls were set using rectal temperature as a criterion to define stress conditions: stressed (N = 5) and non-stressed (N = 5) groups. Skin transcriptomics was performed and correlations between breed composition, phenotypic and skin histological traits were evaluated. Results showed 4309 differentially expressed genes (P < 0.01), 2113 downregulated and 2196 upregulated. Enrichment and ontology analyses revealed 132 GO terms and 67 pathways (P < 0.01), including thermogenesis, glycolysis, gluconeogenesis, mitochondrial activity, antioxidant and immune response, and apoptosis. The identity of the terms and pathways indicated the diversity of mechanisms directed to relieve the animals' suffering, acting from simple passive mechanisms (conduction, convection and radiation) to more complex active ones (behavioural changes, evaporation, vasodilation and wheezing). Furthermore, significant differences between phenotypic and skin histological traits and correlations between pairs of traits suggested a direction towards heat dissipation processes. In this sense, number of vessels was positively correlated with number of sweat glands (P < 0.001) and both were positively correlated with zebuine genetic content (P < 0.05 and P < 0.01, respectively), gland size was positively correlated with epidermal thickness and negatively with hair length (P < 0.05), and epidermal thickness was negatively correlated with gland-epidermis distance (P < 0.0005). These results support the notion that response to heat stress is physiologically complex, producing significant changes in the expression of genes involved in several biological pathways, while the animal's ability to face it depends greatly on their skin features.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

Heat stress adaptation in cows - Physiological responses and underlying molecular mechanisms

Heat stress is a key abiotic stressor for dairy production in the tropics which is further compounded by the ongoing climate change. Heat stress not only adversely impacts the production and welfare of dairy cows but severely impacts the economics of dairying due to production losses and increased cost of rearing. Over the years, selection has ensured development of high producing breeds, however, the thermotolerance ability of animals has been largely overlooked. In the past decade, the ill effects of climate change have made it pertinent to rethink the selection strategies to opt for climate resilient breeds, to ensure optimum production and reproduction. This has led to renewed interest in evaluation of the impacts of heat stress on cows and the underlying mechanisms that results in their acclimatization and adaptation to varied thermal ambience. The understanding of heat stress and associated responses at various level of animal is crucial to device amelioration strategies to secure optimum production and welfare of cows. With this review, an effort has been made to provide an overview on temperature humidity index as an important indicator of heat stress, general effect of heat stress in dairy cows, and impact of heat stress and subsequent response at physiological, haematological, molecular and genetic level of dairy cows.

Transcriptome and coexpression network analysis reveals properties and candidate genes associated with grape (Vitis vinifera L.) heat tolerance

Temperature is one of the most important environmental factors affecting grape season growth and geographical distribution. With global warming and the increasing occurrence of extreme high-temperature weather, the impact of high temperatures on grape production has intensified. Therefore, identifying the molecular regulatory networks and key genes involved in grape heat tolerance is crucial for improving the resistance of grapes and promoting sustainable development in grape production. In this study, we observed the phenotypes and cellular structures of four grape varieties, namely, Thompson Seedless (TS), Brilliant Seedless (BS), Jumeigui (JMG), and Shine Muscat (SM), in the naturally high-temperature environment of Turpan. Heat tolerance evaluations were conducted. RNA-seq was performed on 36 samples of the four varieties under three temperature conditions (28°C, 35°C, and 42°C). Through differential expression analysis revealed the fewest differentially expressed genes (DEGs) between the heat-tolerant materials BS and JMG, and the DEGs common to 1890 were identified among the four varieties. The number of differentially expressed genes within the materials was similar, with a total of 3767 common DEGs identified among the four varieties. KEGG enrichment analysis revealed that fatty acid metabolism, starch and sucrose metabolism, plant hormone signal transduction, the MAPK signaling pathway, and plant-pathogen interactions were enriched in both between different temperatures of the same material, and between different materials of the same temperature. We also conducted statistical and expression pattern analyses of differentially expressed transcription factors. Based on Weighted correlation network analysis (WGCNA), four specific modules highly correlated with grape heat tolerance were identified by constructing coexpression networks. By calculating the connectivity of genes within the modules and expression analysis, six candidate genes (VIT_04s0044g01430, VIT_17s0000g09190, VIT_01s0011g01350, VIT_01s0011g03330, VIT_04s0008g05610, and VIT_16s0022g00540) related to heat tolerance were discovered. These findings provide a theoretical foundation for further understanding the molecular mechanisms of grape heat tolerance and offer new gene resources for studying heat tolerance in grapes.

The Effects of Temperature and Humidity Index on Growth Performance, Colon Microbiota, and Serum Metabolome of Ira Rabbits

This study investigates the effects of different THI values on growth performance, intestinal microbes, and serum metabolism in meat rabbits. The results showed that there were significant differences in THI in different location regions of the rabbit house. The high-THI group (HG) could significantly reduce average daily gain and average daily feed intake in Ira rabbits (p < 0.05). The low-THI group (LG) significantly increased the relative abundance of Blautia (p < 0.05). The HG significantly increased the relative abundance of Lachnospiraceae NK4A136 group and reduced bacterial community interaction (p < 0.05). The cytokine-cytokine receptor interactions, nuclear factor kappa B signaling pathway, and toll-like receptor signaling pathway in each rabbit's gut were activated when the THI was 26.14 (p < 0.05). Metabolic pathways such as the phenylalanine, tyrosine, and tryptophan biosynthesis and phenylalanine metabolisms were activated when the THI was 27.25 (p < 0.05). Meanwhile, the TRPV3 and NGF genes that were associated with heat sensitivity were significantly upregulated (p < 0.05). In addition, five metabolites were found to be able to predict THI levels in the environment with an accuracy of 91.7%. In summary, a THI of 26.14 is more suitable for the growth of meat rabbits than a THI of 27.25, providing a reference for the efficient feeding of meat rabbits.

Comparative expression analysis of water buffalo (Bubalus bubalis) to identify genes associated with economically important traits

The milk, meat, skins, and draft power of domestic water buffalo (Bubalus bubalis) provide substantial contributions to the global agricultural economy. The world's water buffalo population is primarily found in Asia, and the buffalo supports more people per capita than any other livestock species. For evaluating the workflow, output rate, and completeness of transcriptome assemblies within and between reference-free (RF) de novo transcriptome and reference-based (RB) datasets, abundant bioinformatics studies have been carried out to date. However, comprehensive documentation of the degree of consistency and variability of the data produced by comparing gene expression levels using these two separate techniques is lacking. In the present study, we assessed the variations in the number of differentially expressed genes (DEGs) attained with RF and RB approaches. In light of this, we conducted a study to identify, annotate, and analyze the genes associated with four economically important traits of buffalo, viz., milk volume, age at first calving, post-partum cyclicity, and feed conversion efficiency. A total of 14,201 and 279 DEGs were identified in RF and RB assemblies. Gene ontology (GO) terms associated with the identified genes were allocated to traits under study. Identified genes improve the knowledge of the underlying mechanism of trait expression in water buffalo which may support improved breeding plans for higher productivity. The empirical findings of this study using RNA-seq data-based assembly may improve the understanding of genetic diversity in relation to buffalo productivity and provide important contributions to answer biological issues regarding the transcriptome of non-model organisms.

Lysolecithin-derived feed additive improves feedlot performance, carcass characteristics, and muscle fatty acid profile of Bos indicus-influenced cattle fed in a tropical environment

Lysolecithin might increase ruminal and intestinal emulsification, leading to increased digestibility, but there is minimum information about which is the most appropriate phase to start supplementation and its impacts on feedlot performance and muscle fatty acid profile. Two experiments were conducted to evaluate the effects of phase-feeding of Lysoforte™ eXtend (LYSO). In the first experiment, 1,760 predominantly Bos indicus bullocks (initial body weight of 400 ± 0.561 kg) were allocated in a complete randomized block design. LYSO was supplemented at 1 g/1% of ether extract from the diet. Treatments were no LYSO supplementation (NON); LYSO starting during the growing period and continuing during the finishing period; LYSO starting during the finishing period (FIN); and LYSO during adaptation, growing, and finishing periods (ALL). In the second experiment, the same treatments were evaluated with 96 bullocks (64 Nellore and 32½ Nellore × ½ Angus) in a 4 × 2 factorial arrangement (treatments × genotype). For both studies, intake and average daily gain were accessed; carcass characteristics were evaluated in the first experiment, while digestibility of nutrients and profile of muscle fatty acids were measured in the second experiment. In the first experiment, LYSO increased final body weight (P &lt; 0.022) and average daily gain (GRO and FIN; P &lt; 0.05). In the second study, a treatment × breed × feeding phase interaction was observed with Nellore having a greater average daily gain (P &lt; 0.05) than crossbreds in every feeding phase that LYSO was introduced to the diet. A treatment × feeding phase interaction was observed for digestibility, such that LYSO increased total dry matter (P = 0.004), crude protein (P = 0.043), and NDF (P = 0.001) digestibility during the finishing period. A treatment × breed × day classification was observed (P &lt; 0.05). During the finishing phase, crossbreds treated with LYSO had greater DMI (P &lt; 0.05) on very hot days than NON. Also, animals treated with LYSO presented a greater C18:3 n3 concentration (P = 0.047) in the longissimus. Overall, feeding LYSO during GRO and FIN enhanced feedlot performance and should lead to higher intakes during very hot days of the finishing feeding period.

Weighted Gene Co-Expression Network Analysis to Explore Hub Genes of Resveratrol Biosynthesis in Exocarp and Mesocarp of 'Summer Black' Grape

Resveratrol is a polyphenol compound beneficial to human health, and its main source is grapes. In the present study, the molecular regulation of resveratrol biosynthesis in developing grape berries was investigated using weighted gene co-expression network analysis (WGCNA). At the same time, the reason for the resveratrol content difference between grape exocarp (skin) and mesocarp (flesh) was explored. Hub genes (CHS, STS, F3'5'H, PAL, HCT) related to resveratrol biosynthesis were screened with Cytoscape software. The expression level of hub genes in the exocarp was significantly higher than that in the mesocarp, and the expressions of the hub genes and the content of resveratrol in exocarp peaked at the maturity stage. While the expression levels of PAL, CHS and STS in the mesocarp, reached the maximum at the maturity stage, and F3'5'H and HCT decreased. These hub genes likely play a key role in resveratrol biosynthesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis further indicated that resveratrol biosynthesis was related to flavonoid biosynthesis, phenylalanine metabolism, phenylpropanoid biosynthesis, and stilbene biosynthesis pathways. This study has theoretical significance for exploring genes related to resveratrol biosynthesis in the exocarp and mesocarp of grapes, and provides a theoretical basis for the subsequent function and regulatory mechanism of hub genes.

Economically sustainable shade design for feedlot cattle

Provision of shade reduces radiant heat load on feedlot cattle, thus reducing demand of water and energy for thermoregulation. While the positive effects of shade on animal welfare are widely known, the literature lacks data on the magnitude of its economic impacts. In this study, we propose the concept of novel shade design to prove that a correctly oriented and dimensioned roof structure, which optimizes shade to be displaced within the pens, motivates cattle to seek shade, protect them from short-wave solar radiation, and is resilient to counteract weather adverse conditions. The beneficial outcome is improvement in animal welfare and productive performance, as well as increments on financial return and sustainability. To attest these benefits, eight hundred B. indicus × Bos taurus bulls were randomly assigned in pens with or without shade from a galvanized steel-roof structure. Performance data (e.g., dry matter intake, body weight gain, feed efficiency and hot carcass weight) and heat stress indicators (e.g., subcutaneous temperature, body-surface temperature, respiratory rate and water intake) were assessed along the study period. The economic outcomes derived from shade implementation were determined using the net present value. Meteorological variables were also monitored every 1 min, and grouped in a thermal comfort index for feedlot cattle, the InComfort Index (InCI). The shade structure efficiently reduced radiant heat load on cattle in pens with shade. According to the classification of the InCI, during very hot days (InCI > 0.6; around noon with mean solar radiation above 800 W m^-2 and mean air temperature above 33°C), greater proportion (80%) of animals in shaded pens were using shade. Under such circumstances, cattle in shade had water intake reduced by 3.4 L per animal, body temperature was lower by 5°C, subcutaneous temperature was lower by 1°C and respiration rate was lower by 10 breaths min^-1 compared to animals in pens without shade (P = 0.0001). Although dry matter intake was similar (P = 0.6805), cattle in pens with shade had higher average daily gain reflected in a heavier hot carcass weight (8 kg animal^-1; P = 0.0002). Considering an initial investment of $90 per animal to build a structure that lasts 15 years, the expected payback time is four finishing cycles (~110 days per cycle). In conclusion, this study confirms that the proposed novel shade design is economically profitable, improves performance, and enhances animal welfare.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

Transcriptome and Metabolome Analyses Revealed the Response Mechanism of Sugar Beet to Salt Stress of Different Durations

Salinity is one of the most serious threats to agriculture worldwide. Sugar beet is an important sugar-yielding crop and has a certain tolerance to salt; however, the genome-wide dynamic response to salt stress remains largely unknown in sugar beet. In the present study, physiological and transcriptome analyses of sugar beet leaves and roots were compared under salt stress at five time points. The results showed that different salt stresses influenced phenotypic characteristics, leaf relative water content and root activity in sugar beet. The contents of chlorophyll, malondialdehyde (MDA), the activities of peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) were also affected by different salt stresses. Compared with control plants, there were 7391 and 8729 differentially expressed genes (DEGs) in leaves and roots under salt stress, respectively. A total of 41 hub genes related to salt stress were identified by weighted gene co-expression network analysis (WGCNA) from DEGs, and a transcriptional regulatory network based on these genes was constructed. The expression pattern of hub genes under salt stress was confirmed by qRT-PCR. In addition, the metabolite of sugar beet was compared under salt stress for 24 h. A total of 157 and 157 differentially accumulated metabolites (DAMs) were identified in leaves and roots, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis further indicated that DEGs and DAMs act on the starch and sucrose metabolism, alpha-linolenic acid metabolism, phenylpropanoid biosynthesis and plant hormone signal transduction pathway. In this study, RNA-seq, WGCNA analysis and untargeted metabolomics were combined to investigate the transcriptional and metabolic changes of sugar beet during salt stress. The results provided new insights into the molecular mechanism of sugar beet response to salt stress, and also provided candidate genes for sugar beet improvement.

Genome-Wide Association Study for Lactation Performance in the Early and Peak Stages of Lactation in Holstein Dairy Cows

Simple Summary

Although genome-wide association studies (GWAS) have been carried out within a variety of cattle breeds, they are mainly based on the accumulated 305-day lactation yield traits estimated by summing the test-day recorded every day during the lactation period, or combining the weekly or monthly test-day records by linear interpolation. Since the additive genetic variance for milk yield and composition changes during lactation, the genetic effects of QTL related to these traits are not constant during the lactation period. Therefore, a better understanding of the genetic architecture of milk production traits in different lactation stages (e.g., beginning, peak, and end stages of lactation) is needed. The aim of this study was to detect genomic loci associated with lactation performance during 9 to 50 days in milk (DIM) in Holstein dairy cows. Candidate genes identified for milk production traits showed contrasting results between the EARLY and PEAK stages of lactation. Based on the results of this study, it can be concluded that in any genomic study it should be taken into account that the genetic effects of genes related to the lactation performance are not constant during the lactation period.

Abstract

This study aimed to detect genomic loci associated with the lactation performance during 9 to 50 days in milk (DIM) in Holstein dairy cows. Daily milk yield (MY), fat yield (FY), and protein yield (PY) during 9 to 50 DIM were recorded on 134 multiparous Holstein dairy cows distributed in four research herds. Fat- and protein-corrected milk (FPCM), fat-corrected milk (FCM), and energy-corrected milk (ECM) were predicted. The records collected during 9 to 25 DIM were put into the early stage of lactation (EARLY) and those collected during 26 to 50 DIM were put into the peak stage of lactation (PEAK). Then, the mean of traits in each cow included in each lactation stage (EARLY and PEAK) were estimated and used as phenotypic observations for the genome-wide association study. The included animals were genotyped with the Illumina BovineHD Genotyping BeadChip (Illumina Inc., San Diego, CA, USA) for a total of 777,962 single nucleotide polymorphisms (SNPs). After quality control, 585,109 variants were analyzed using GEMMA software in a mixed linear model. Although there was no SNP associated with traits included at the 5% genome-wide significance threshold, 18 SNPs were identified to be associated with milk yield and composition at the suggestive genome-wide significance threshold. Candidate genes identified for milk production traits showed contrasting results between the EARLY and PEAK stages of lactation. This suggests that differential sets of candidate genes underlie the phenotypic expression of the considered traits in the EARLY and PEAK stages of lactation. Although further functional studies are needed to validate our findings in independent populations, it can be concluded that in any genomic study it should be taken into account that the genetic effects of genes related to the lactation performance are not constant during the lactation period.

Whole-Genome-Based Web Genomic Resource for Water Buffalo (Bubalus bubalis)

Water buffalo (Bubalus bubalis), belonging to the Bovidae family, is an economically important animal as it is the major source of milk, meat, and drought in numerous countries. It is mainly distributed in tropical and subtropical regions with a global population of approximately 202 million. The advent of low cost and rapid sequencing technologies has opened a new vista for global buffalo researchers. In this study, we utilized the genomic data of five commercially important buffalo breeds, distributed globally, namely, Mediterranean, Egyptian, Bangladesh, Jaffrarabadi, and Murrah. Since there is no whole-genome sequence analysis of these five distinct buffalo breeds, which represent a highly diverse ecosystem, we made an attempt for the same. We report the first comprehensive, holistic, and user-friendly web genomic resource of buffalo (BuffGR) accessible at http://backlin.cabgrid.res.in/buffgr/, that catalogues 6028881 SNPs and 613403 InDels extracted from a set of 31 buffalo tissues. We found a total of 7727122 SNPs and 634124 InDels distributed in four breeds of buffalo (Murrah, Bangladesh, Jaffarabadi, and Egyptian) with reference to the Mediterranean breed. It also houses 4504691 SSR markers from all the breeds along with 1458 unique circRNAs, 37712 lncRNAs, and 938 miRNAs. This comprehensive web resource can be widely used by buffalo researchers across the globe for use of markers in marker trait association, genetic diversity among the different breeds of buffalo, use of ncRNAs as regulatory molecules, post-transcriptional regulations, and role in various diseases/stresses. These SNPs and InDelscan also be used as biomarkers to address adulteration and traceability. This resource can also be useful in buffalo improvement programs and disease/breed management.

Genome-wide association analyses identify genotype-by-environment interactions of growth traits in Simmental cattle

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.

Water buffalo production in the Brazilian Amazon Basin: a review

The Brazilian Amazon has witnessed, in the last decades, an increase in the water buffalo (Bubalus bubalis) inventory, with interesting productivity results. As the Brazilian Amazon contains the main water buffalo population in the Americas, the aim of this work is to review its most relevant production systems and some peculiarities about meat and milk production in this territory. The opening section describes the Amazon Basin, the most common water buffalo breeds, a brief history of the local livestock farming beginning in 1644. Also, it presents how water buffaloes gradually replaced bovine herds, especially where the latter had a lower productive performance. The use of extensive or more intensified models is pointed out and the ecosystems in which buffaloes are raised are detailed since native or cultivated pastures can be used in floodplains or drylands. Buffalo raising is favored in the Amazon due to the climate, soil, genetic variability of forages, animal adaptability, and physical space. Thus, it is clear that buffaloes have a high potential for meat and milk production and are an alternative in the use of altered areas of the Amazon; and, in the recent past, the low profitability of buffalo farming in traditional production systems in the Amazon was the reason which made this activity economically unattractive. Most recent technologies as outdoor confinements and silvopastoral systems are pointed out as more suitable regarding land-use policies, and buffalo farming for meat and milk production fits perfectly in this context, with productivity and beneficial socioeconomic.

Differential expression of miRNAs and related mRNAs during heat stress in buffalo heifers

The present experiment was aimed to study differential expression of miRNAs and related mRNAs during heat stress (HS) in buffalo heifers. Twelve Murrah buffalo heifers aged between 1.5 and 2.0 years, weighting between 250 and 300 Kg were randomly assigned into two equal groups. The animals were kept in the psychrometric chamber under Thermo-neutral (TN; THI = 72) and HS (THI = 87-90) conditions for 6 h every day between 1000 and 1600 h for 21 days. The blood sampling was done at 1500 h on 15th day of the experiment and physiological parameters viz. pulse rate (PR), respiratory rate (RR) and rectal temperature (RT) were recorded at 1500 h on day -5, -3, -1, 0, +1, +3, +5 with respect to blood sampling. PBMCs were used for extraction of miRNAs and total RNA; and first strand cDNA was synthesized. qPCR was performed for relative gene expression studies. Physiological, hematological (erythrocytic indices), biochemical (triglycerides, urea, ALT, AST, LDH), redox (SOD, ROS) and endocrine parameters (T₄) altered significantly (P < 0.05) during HS as compared to TN. Out of eight targeted miRNAs only four were expressed in buffalo heifers. The relative expression of bta-mir-142, bta-mir-1248 and bta-mir-2332 was significantly (P < 0.05) up-regulated whereas expression of bta-mir-2478 was significantly (P < 0.05) down-regulated during HS as compared to TN. The relative expression of the predicted target genes i.e. HSF1, HSP60, HSP70, HSPA8 and HSP90 were significantly (P < 0.05) up-regulated whereas HSF4 expression was significantly (P < 0.05) down-regulated during HS as compared to TN. It can be concluded that a THI of 87-90 could lead to a moderate HS in buffalo heifers. Differential expression studies of miRNAs and related mRNAs in present study deciphers the role of miRNAs in the heat tolerance in buffalo heifers.

Thermoregulatory responses in riverine buffaloes against heat stress: An updated review

High heat and humidity stress have been a perpetual perilous for the buffalo's production and productivity in tropics and subtropics including India. Productive potential of livestock's species including buffaloes is maximum with in thermo-neutral zone (TNZ) and if ambient temperature exceeds TNZ and upper critical temperature expose livestock's to heat stress conditions. For decades, heat stress has been the prime factor to plummet buffalo's growth, development, reproduction and production in tropics and subtropics including India. In general, buffaloes are homeotherms and known as temperature regulators as they resist the variations in ambient temperatures. Generally, buffaloes like other livestock's display amalgamation of thermoregulatory responses to withstand the changes occurred in their micro and macro environment. These thermoregulatory responses are behavioural, physiological, neuro-endocrine and molecular responses acting synergistically to counteract the deleterious effects of heat stress. Amidst all responses, molecular responses play major role to confer thermo-tolerance through expression of highly conserved family of proteins known as heat shock proteins (HSPs). Despite of these thermoregulatory responses, heat stress prodigiously muddles buffalo's production and productivity. The present review highlights the thermoregulatory responses manifested by riverine buffaloes against heat stress.

Seasonal effects on miRNA and transcriptomic profile of oocytes and follicular cells in buffalo (Bubalus bubalis)

Season clearly influences oocyte competence in buffalo (Bubalus bubalis); however, changes in the oocyte molecular status in relation to season are poorly understood. This study characterizes the microRNA (miRNA) and transcriptomic profiles of oocytes (OOs) and corresponding follicular cells (FCs) from buffalo ovaries collected in the breeding (BS) and non-breeding (NBS) seasons. In the BS, cleavage and blastocyst rates are significantly higher compared to NBS. Thirteen miRNAs and two mRNAs showed differential expression (DE) in FCs between BS and NBS. DE-miRNAs target gene analysis uncovered pathways associated with transforming growth factor β (TGFβ) and circadian clock photoperiod. Oocytes cluster in function of season for their miRNA content, showing 13 DE-miRNAs between BS and NBS. Between the two seasons, 22 differentially expressed genes were also observed. Gene Ontology (GO) analysis of miRNA target genes and differentially expressed genes (DEGs) in OOs highlights pathways related to triglyceride and sterol biosynthesis and storage. Co-expression analysis of miRNAs and mRNAs revealed a positive correlation between miR-296-3p and genes related to metabolism and hormone regulation. In conclusion, season significantly affects female fertility in buffalo and impacts on oocyte transcriptomic of genes related to folliculogenesis and acquisition of oocyte competence.

Reproduction in female wild cattle: Influence of seasonality on ARTs

Wild cattle species, often considered less alluring than certain conservation-dependent species, have not attracted the same level of interest as the charismatic megafauna from the general public, private or corporate donors, and other funding agencies. Currently, most wild cattle populations are vulnerable or threatened with extinction. The implementation of reproductive technologies to maintain genetically healthy cattle populations in situ and ex situ has been considered for more than 30 years. Protocols developed for domestic cattle breeds have been used with some success in various wild cattle species. However, inherent differences in the natural life history of these species makes extrapolation of domestic cattle protocols difficult, and in some cases, minimally effective. Reproductive seasonality, driven by either photoperiod or nutritional resource availability, has significant influence on the success of assisted reproductive technologies (ARTs). This review focuses on the physiological processes that differ in breeding (ovulatory) and non-breeding (anovulatory) seasons in female cattle, and the potential methods used to overcome these challenges. Techniques to be discussed within the context of seasonality include: estrus synchronization and ovulation induction, ovarian superstimulation, artificial insemination (AI), multiple ovulation embryo transfer (MOET), and ovum pick-up (OPU) with in vitro fertilization (IVF) and embryo transfer (ET).

Asian water buffalo: domestication, history and genetics

The domestic Asian water buffalo (Bubalus bubalis) is found on all five continents, with a global population of some 202 million. The livelihoods of more people depend on this species than on any other domestic animal. The two distinct types (river and swamp) descended from different wild Asian water buffalo (Bubalus arnee) populations that diverged some 900 kyr BP and then evolved in separate geographical regions. After domestication in the western region of the Indian subcontinent (ca. 6300 years BP), the river buffalo spread west as far as Egypt, the Balkans and Italy. Conversely, after domestication in the China/Indochina border region ca. 3000-7000 years BP, swamp buffaloes dispersed through south-east Asia and China as far as the Yangtze River valley. Molecular and morphological evidence indicates that swamp buffalo populations have strong geographic genetic differentiation and a lack of gene flow, but strong phenotypic uniformity. In contrast, river buffalo populations show a weaker phylogeographic structure, but higher phenotypic diversity (i.e. many breeds). The recent availability of a high-quality reference genome and of a medium-density marker panel for genotyping has triggered a number of genome-wide investigations on diversity, evolutionary history, production traits and functional elements. The growing molecular knowledge combined with breeding programmes should pave the way to improvements in production, environmental adaptation and disease resistance in water buffalo populations worldwide.

Transcriptome Analysis Reveals Potential Regulatory Genes Related to Heat Tolerance in Holstein Dairy Cattle

Heat stress affects the physiology and production performance of Chinese Holstein dairy cows. As such, the selection of heat tolerance in cows and elucidating its underlying mechanisms are vital to the dairy industry. This study aimed to investigate the heat tolerance associated genes and molecular mechanisms in Chinese Holstein dairy cows using a high-throughput sequencing approach and bioinformatics analysis. Heat-induced physiological indicators and milk yield changes were assessed to determine heat tolerance levels in Chinese Holstein dairy cows by Principal Component Analysis method following Membership Function Value Analysis. Results indicated that rectal temperature (RT), respiratory rate (RR), and decline in milk production were significantly lower (p < 0.05) in heat tolerant (HT) cows while plasma levels of heat shock protein (HSP: HSP70, HSP90), and cortisol were significantly higher (p < 0.05) when compared to non-heat tolerant (NHT) Chinese Holstein dairy cows. By applying RNA-Seq analysis, we identified 200 (81 down-regulated and 119 up-regulated) significantly (|log2fold change| ≥ 1.4 and p ≤ 0.05) differentially expressed genes (DEGs) in HT versus NHT Chinese Holstein dairy cows. In addition, 14 of which were involved in protein–protein interaction (PPI) network. Importantly, several hub genes (OAS2, MX2, IFIT5 and TGFB2) were significantly enriched in immune effector process. These findings might be helpful to expedite the understanding for the mechanism of heat tolerance in Chinese Holstein dairy cows.