The present study was aimed at identifying causative hub genes within modules formed by co-expression and protein–protein interaction (PPI) networks, followed by Bayesian network (BN) construction in ...the liver transcriptome of starved zebrafish. To this end, the GSE11107 and GSE112272 datasets from the GEO databases were downloaded and meta-analyzed using the MetaDE package, an add-on R package. Differentially expressed genes (DEGs) were identified based upon expression intensity N(µ = 0.2, σ2 = 0.4). Reconstruction of BNs was performed by the bnlearn R package on genes within modules using STRINGdb and CEMiTool. ndufs5 (shared among PPI, BN and COEX), rps26, rpl10, sdhc (shared between PPI and BN), ndufa6, ndufa10, ndufb8 (shared between PPI and COEX), skp1, atp5h, ndufb10, rpl5b, zgc:193613, zgc:123327, zgc:123178, wu:fc58f10, zgc:111986, wu:fc37b12, taldo1, wu:fb62f08, zgc:64133 and acp5a (shared between COEX and BN) were identified as causative hub genes affecting gene expression in the liver of starving zebrafish. Future work will shed light on using integrative analyses of miRNA and DNA microarrays simultaneously, and performing in silico and experimental validation of these hub-causative (CST) genes affecting starvation in zebrafish.
During the last few decades, the effective population size of indigenous zebu cattle breeds has declined drastically, resulting in the classification of some of them into the vulnerable, endangered, ...or critically endangered category. Drastic reductions in the effective size of a population may result in genetic bottlenecks and can affect within-breed genetic variability and its viability. The present study was undertaken with the objective of evaluating South Indian zebu cattle populations for mutation drift equilibrium and to detect the occurrence of recent genetic bottleneck events. A total of 293 cattle from eight indigenous breeds were genotyped at 27 FAO/ISAG-recommended microsatellite marker loci. Three different statistical tests, viz., the sign test, standardized differences test, and Wilcoxon sign rank test were performed using allele frequency data to detect loci with heterozygosity excess under the infinite alleles, stepwise, and two-phase mutation models. Under the infinite alleles model, the observed number of loci with heterozygosity excess (He > Heq) ranged between 10 and 19 among the investigated cattle breeds. However, the observed heterozygosity excess was not statistically significant (p > 0.05) in any of the studied breeds. Similarly, the standardized differences test and Wilcoxon sign rank test revealed no concrete evidence for the occurrence of a recent genetic bottleneck in South Indian zebu cattle breeds. The qualitative test for mode-shift distortion revealed a normal L-shaped distribution of allele frequencies, suggesting a lack of evidence for the loss of low-frequency alleles in all the investigated South Indian zebu cattle breeds.
Bayesian gene networks are powerful for modelling causal relationships and incorporating prior knowledge for making inferences about relationships. We used three algorithms to construct Bayesian gene ...networks around genes expressed in the bovine uterus and compared the efficacies of the algorithms. Dataset GSE33030 from the Gene Expression Omnibus (GEO) repository was analyzed using different algorithms for hub gene expression due to the effect of progesterone on bovine endometrial tissue following conception. Six different algorithms (grow-shrink, max-min parent children, tabu search, hill-climbing, max-min hill-climbing and restricted maximum) were compared in three higher categories, including constraint-based, score-based and hybrid algorithms. Gene network parameters were estimated using the bnlearn bundle, which is a Bayesian network structure learning toolbox implemented in R. The results obtained indicated the tabu search algorithm identified the highest degree between genes (390), Markov blankets (25.64), neighborhood sizes (8.76) and branching factors (4.38). The results showed that the highest number of shared hub genes (e.g., proline dehydrogenase 1 (
), Sam-pointed domain containing Ets transcription factor (
), monocyte-to-macrophage differentiation associated 2 (
), semaphorin 3E (
), solute carrier family 27 member 6 (
) and actin gamma 2 (
)) was seen between the hybrid and the constraint-based algorithms, and these genes could be recommended as central to the GSE33030 data series. Functional annotation of the hub genes in uterine tissue during progesterone treatment in the pregnancy period showed that the predicted hub genes were involved in extracellular pathways, lipid and protein metabolism, protein structure and post-translational processes. The identified hub genes obtained by the score-based algorithms had a role in 2-arachidonoylglycerol and enzyme modulation. In conclusion, different algorithms and subsequent topological parameters were used to identify hub genes to better illuminate pathways acting in response to progesterone treatment in the bovine uterus, which should help with our understanding of gene regulatory networks in complex trait expression.
Microscopy and polymerase chain reaction (PCR) were used to survey pathogenic trypanosome infection in naturally infected Nigerian cattle. In 411 animals sampled, microscopy detected 15.1% positive ...infection of at least one of Trypanosoma brucei, Trypanosoma congolense or Trypanosoma vivax, while PCR detected 63.7% positive infections of at least one of those species and Trypanosoma evansi. PCR detected 4.4%, 48.7%, 26.0% and 0.5% respectively of T. brucei, T. congolense, T. vivax and T. evansi infections. All of the T. congolense detected were savannah-type, except for two forest-type infections. Prevalence of mixed infections was 13.9%, being primarily co-infection by T. congolense and T. vivax while prevalence of mixed infections by T. evansi, T. vivax and T. congolense was 1.5%. Microscopy showed poor sensitivity but specificity greater than 94%. Infection rates were much higher in Southern than in Northern Nigeria. Infections were lowest in N’dama compared to Muturu, Sokoto Gudali and White Fulani breeds. Animals with T. vivax monoinfection and mixed infections showed significantly lower packed cell volume (PCV) values. Those infected with any Trypanosoma species with <200 parasites/μl showed higher PCV values than those infected with >200 parasites/μl. The new finding of savannah- and forest- type T. congolense in Nigeria and the relatively high abundance of mixed infections are of significant clinical relevance. This study also suggests that T. congolense is the most prevalent species in Nigeria.
We applied an untargeted metabolomics technique to analyze the plasma carboxyl-metabolome of beef steers with divergent average daily gain (ADG). Forty-eight newly weaned Angus crossbred beef steers ...were fed the same total mixed ration ad libitum for 42 days. On day 42, the steers were divided into two groups of lowest (LF:
= 8) and highest ADG (HF:
= 8), and blood samples were obtained from the two groups for plasma preparation. Relative quantification of carboxylic-acid-containing metabolites in the plasma samples was determined using a metabolomics technique based on chemical isotope labeling liquid chromatography mass spectrometry. Metabolites that differed (fold change (FC) ≥ 1.2 or ≤ 0.83 and FDR ≤ 0.05) between LF and HF were identified using a volcano plot. Metabolite set enrichment analysis (MSEA) of the differential metabolites was done to determine the metabolic pathways or enzymes that were potentially altered. In total, 328 metabolites were identified. Volcano plot analysis revealed 43 differentially abundant metabolites; several short chain fatty acids and ketone bodies had greater abundance in HF steers. Conversely, several long chain fatty acids were greater in LF steers. Five enzymatic pathways, such as fatty acyl CoA elongation and fatty-acid CoA ligase were altered based on MSEA. This study demonstrated that beef steers with divergent ADG had altered plasma carboxyl-metabolome, which is possibly caused by altered abundances and/or activities of enzymes involved in fatty acid oxidation and biosynthesis in the liver.
Data of 475 goats (242 Zaraibi, 84 Baladi, and 149 Damascus) kept at Sakha Farm, belonging to Animal Production Research Institute (APRI), Ministry of Agriculture, Dokki, Cairo, Egypt, from the ...period of 2011 to 2017 were used to estimate direct (additive) and maternal genetic effects for litter size at birth (LSB). Three animal models were used. Model 1 includes the fixed effects of breed, season, year of kidding, age of doe, and direct genetic effect. Model 2 is similar to model 1 beside the permanent environmental effect. Model 3 is similar to model 2 and includes the maternal genetic, covariance between direct and maternal genetic effects. The effects of breed, season, year of kidding, and age of doe were significant on LSB. Heritability estimates for LSB are 0.19±0.050, 0.19±0.160, and 0.06±0.007, as estimated from models 1, 2, and 3, respectively. Maternal heritability estimate using model 3 was 0.05±0.006. The removal of additive genetic maternal effect and covariance between direct and maternal effects from the model increased heritability of direct genetic effect by 0.13. Annual phenotypic trends for LSB in Zaraibi, Baladi, and Damascus goat breeds are negative. Annual genetic trends for LSB for the three breeds of goats are positive, significant, and being 0.10±0.01, 0.25±0.02, and 0.01±0.001 litter/year for Zaraibi, Baladi, and Damascus, respectively. The present results indicated that selection for LSB will take a long time. The speed and efficiency of selection is expected to increase by use of molecular markers in selection.
Type II melanoma-associated antigens (MAGE) are a subgroup of about a dozen proteins found in various locations in the genome and expressed in normal tissues, thus are not related to cancer as the ...type I MAGE genes. This gene family exists as a single copy in non-mammals and monotremata, but found as two copies in metatherians and occur as a diverse group in all eutherians. Our studies suggest MAGED2 as the ancestor of this subfamily and the most likely evolutionary history of eutherian type II MAGE genes is hereby proposed based on synteny conservation, phylogenetic relations, genome location, homology conservation, and the protein and gene structures. Type II genes can be divided into two: those with 13 exons (MAGED1, MAGED2, TRO, and MAGED4) and those with only one exon (MAGEE1, MAGEE2, MAGEF1, NSMCE3, MAGEH1, MAGEL2, and NDN) with different evolutionary patterns. Our results suggest a need to change the gene nomenclature to MAGE1 (the ancestral gene), currently designated as LOC103095671 and LOC100935086, in opossum and Tasmanian devil, respectively, and MAGE2 (the duplicated one), currently designated as LOC100617402 and NDNL2, respectively, to avoid confusion. We reconstructed the phylogenetic relationships among 23 mammalian species using the combined sequences of MAGED1, MAGED2, MAGEL2, and NDN, because of their high divergence, and found high levels of support, being able to resolve the phylogenetic relationships among Euarchontoglires, Laurasiatheria, Afrotheria, and Xenarthra, as an example that small, but phylogenetically informative sequences, can be very useful for resolving basal mammalian clades.
Data of 651 lambs (68 Romanov, 49 Rahmani, 151 ♀1/2 Rahmani X ♂1/2 Romanov) and 383 (♀3/4 Rahmani and 1/4♂ Romanov) were collected from Mehalet Mousa Farm, belonging to Animal Production Research ...Institute from the period of 2009 to 2016 to estimate phenotypic and genetic parameters. The traits studied were birth weight (BW), body weight at four week (BW4), body weight at eight weeks (BW8) and body weight at twelve weeks (BW12) or weaning weight. Least squares analysis of variance shows significance of the effects of breed groups, gender of lambs, birth type; month of birth and year of birth on all traits studied. Rahmani lambs had heavier BW, BW4, BW8 and BW12 while Romanov lambs had the lowest ones. The first generation (♀1/2 Rhamani X ♂1/2 Romanov) had heavier body weights than Romanov and the second generation (♀3/4 Rahmani X ♂1/4 Roamnov). Gender of lambs had highly significant effect on body weights. Males were significantly (p < 0.01) heavier than females for all traits studied. Least square means of BW, BW4, BW8 and BW12 for single lambs were 2.69, 10.43, 13.53 and 16.10 kg, respectively. Least square means of BW, BW4, BW8 and BW12 for twin lambs were 2.50, 9.37, 12.5 and 15.16 kg, respectively, while least square means of BW, BW4, BW8 and BW12 for triple lambs were 2.09, 7.86, 10.83 and 13.67 kg, respectively. Estimates of direct heritability measured by single trait animal model were 0.14, 0.23, 0.25 and 0.26 for BW, BW4, BW8 and BW12, respectively, and the corresponding measured by multi trait animal model were 0.17, 0.24, 0.32 and 0.36 for the same traits, respectively. All genetic and phenotypic correlations among different traits studied are positive and significant.
Background: Sheep is important in the socio-economic lives of people around the world. It is estimated that more than half of our once common livestock breeds are now endangered. Since genetic ...characterization of Nigerian sheep is still lacking, we analyzed ten morphological traits on 402 animals and 15 microsatellite DNA markers in 384 animals of the 4 Nigerian sheep breeds to better understand genetic diversity for breeding management and germplasm conservation. Results: Morphological traits of Uda and Balami were significantly (P 〈 0.05) higher than Yankasa, which were both higher than West African Dwarf (WAD) sheep. Stepwise discriminant analysis showed tail length, rump height, chest girth, ear length and chest depth as the most discriminating variables for classification. Mahalanobis distances show the least differentiation between Uda and Balami and the largest between WAD and Balami sheep. While 93.3% of WAD sheep were correctly assigned to their source genetic group, 63.9% of Yankasa, 61.2% of Balami and 45.2% of Uda were classified correctly by nearest neighbour discriminant analysis. The overall high Polymorphism Information Content (PIC) of all microsatellite markers ranged from 0.751 to 0.927 supporting their use in genetic characterization. Expected heterozygosity was high for all loci (0.783 to 0.93). Mean heterozygote deficiency across all populations (0.17 to 0.534) possibly indicate significant inbreeding (P 〈 0.05). Mean values for FST, FIT and F~s statistics across all loci were 0.088, 0.394 and 0.336 respectively. Yankasa and Balami are the most closely related breeds (DA = 0.184) while WAD and Balami are the farthest apart breeds (DA-- 0.665), which is coincident with distance based on morphological analysis and population structure assessed by STRUCTURE. Conclusions: These results suggest that within-breed genetic variation in Nigerian sheep is higher than between-breeds and may be a valuable tool for genetic improvement and conservation. The higher genetic variability in Yankasa suggests the presence of unique ancestral alleles reflecting the presence of certain functional genes which may result in better adaptability in characteristics are potentially useful in planning indigenous sheep. more agro-ecological zones of Nigeria. These genetic mprovement and conservation strategies in Nigerian
Skin is a major thermoregulatory organ in the body controlling homeothermy, a critical function for climate adaptation. We compared genes expressed between tropical- and temperate-adapted cattle to ...better understand genes involved in climate adaptation and hence thermoregulation. We profiled the skin of representative tropical and temperate cattle using RNA-seq. A total of 214,754,759 reads were generated and assembled into 72,993,478 reads and were mapped to unique regions in the bovine genome. Gene coverage of unique regions of the reference genome showed that of 24,616 genes, only 13,130 genes (53.34%) displayed more than one count per million reads for at least two libraries and were considered suitable for downstream analyses. Our results revealed that of 255 genes expressed differentially, 98 genes were upregulated in tropically-adapted White Fulani (WF;
) and 157 genes were down regulated in WF compared to Angus, AG (
). Fifteen pathways were identified from the differential gene sets through gene ontology and pathway analyses. These include the significantly enriched melanin metabolic process, proteinaceous extracellular matrix, inflammatory response, defense response, calcium ion binding and response to wounding. Quantitative PCR was used to validate six representative genes which are associated with skin thermoregulation and epithelia dysfunction (mean correlation 0.92;
< 0.001). Our results contribute to identifying genes and understanding molecular mechanisms of skin thermoregulation that may influence strategic genomic selection in cattle to withstand climate adaptation, microbial invasion and mechanical damage.