Genome Scan for Parent-of-Origin QTL Effects on Bovine Growth and Carcass Traits

Parent-of-origin effects (POE) such as genomic imprinting influence growth and body composition in livestock, rodents, and humans. Here, we report the results of a genome scan to detect quantitative trait loci (QTL) with POE on growth and carcass traits in Angus × Brahman cattle crossbreds. We identified 24 POE–QTL on 15 Bos taurus autosomes (BTAs) of which six were significant at 5% genome-wide (GW) level and 18 at the 5% chromosome-wide (CW) significance level. Six QTL were paternally expressed while 15 were maternally expressed. Three QTL influencing post-weaning growth map to the proximal end of BTA2 (linkage region of 0–9 cM; genomic region of 5.0–10.8 Mb), for which only one imprinted ortholog is known so far in the human and mouse genomes, and therefore may potentially represent a novel imprinted region. The detected QTL individually explained 1.4 ∼ 5.1% of each trait’s phenotypic variance. Comparative in silico analysis of bovine genomic locations show that 32 out of 1,442 known mammalian imprinted genes from human and mouse homologs map to the identified QTL regions. Although several of the 32 genes have been associated with quantitative traits in cattle, only two (GNAS and PEG3) have experimental proof of being imprinted in cattle. These results lend additional support to recent reports that POE on quantitative traits in mammals may be more common than previously thought, and strengthen the need to identify and experimentally validate cattle orthologs of imprinted genes so as to investigate their effects on quantitative traits

Constraint-Based, Score-Based and Hybrid Algorithms to Construct Bayesian Gene Networks in the Bovine Transcriptome

SIMPLE SUMMARY: In this study, we investigated and compared six different Bayesian network algorithms from three different categories to identify hub genes critical to gene expression networks activated in response to progesterone in the bovine uterus. We observed many common hub genes identified between constraint-based algorithms (CBAs) and hybrid algorithms (HAs), while it appeared that score-based algorithm (SBA) methods led to more accurate and relevant predictions of core genes. The results revealed that the identification of hub genes was affected by the type of network reconstruction and by the subsequently used topological parameters. Two identified genes known to have roles during pregnancy are ISG15 and DGAT2. The identified hub genes are associated with biological processes such as amino acid metabolism, hormonal signaling pathways and the immune system. Our analysis revealed a role for miRNAs in the regulation of this system. The biological and physiological roles (enzymatic and hormonal effects) of unannotated identified hub genes should be functionally validated by further studies. ABSTRACT: 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 (PRODH), Sam-pointed domain containing Ets transcription factor (SPDEF), monocyte-to-macrophage differentiation associated 2 (MMD2), semaphorin 3E (SEMA3E), solute carrier family 27 member 6 (SLC27A6) and actin gamma 2 (ACTG2)) 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

An Integrated Bioinformatics Approach to Identify Network Derived Hub Genes in Starving Zebrafish

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

Conservation of Repeats at the Mammalian KCNQ1OT1-CDKN1C Region Suggests a Role in Genomic Imprinting

KCNQ1OT1 is located in the region with the highest number of genes showing genomic imprinting, but the mechanisms controlling the genes under its influence have not been fully elucidated. Therefore, we conducted a comparative analysis of the KCNQ1/KCNQ1OT1-CDKN1C region to study its conservation across the best assembled eutherian mammalian genomes sequenced to date and analyzed potential elements that may be implicated in the control of genomic imprinting in this region. The genomic features in these regions from human, mouse, cattle, and dog show a higher number of genes and CpG islands (detected using cpgplot from EMBOSS), but lower number of repetitive elements (including short interspersed nuclear elements and long interspersed nuclear elements), compared with their whole chromosomes (detected by RepeatMasker). The KCNQ1OT1-CDKN1C region contains the highest number of conserved noncoding sequences (CNS) among mammals, where we found 16 regions containing about 38 different highly conserved repetitive elements (using mVista), such as LINE1 elements: L1M4, L1MB7, HAL1, L1M4a, L1Med, and an LTR element: MLT1H. From these elements, we found 74 CNS showing high sequence identity (>70%) between human, cattle, and mouse, from which we identified 13 motifs (using Multiple Em for Motif Elicitation/Motif Alignment and Search Tool) with a significant probability of occurrence, 3 of which were the most frequent and were used to find transcription factor–binding sites. We detected several transcription factors (using JASPAR suite) from the families SOX, FOX, and GATA. A phylogenetic analysis of these CNS from human, marmoset, mouse, rat, cattle, dog, horse, and elephant shows branches with high levels of support and very similar phylogenetic relationships among these groups, confirming previous reports. Our results suggest that functional DNA elements identified by comparative genomics in a region densely populated with imprinted mammalian genes may be related to the regulation of imprinted gene expression

Sequence variation of necdin gene in Bovidae

Abstract Background Necdin (NDN), a member of the melanoma antigen family showing imprinted pattern of expression, has been implicated as causing Prader-Willi symptoms, and known to participate in cellular growth, cellular migration and differentiation. The region where NDN is located has been associated to QTLs affecting reproduction and early growth in cattle, but location and functional analysis of the molecular mechanisms have not been established. Methods Here we report the sequence variation of the entire coding sequence from 72 samples of cattle, yak, buffalo, goat and sheep, and discuss its variation in Bovidae. Median-joining network analysis was used to analyze the variation found in the species. Synonymous and non-synonymous substitution rates were determined for the analysis of all the polymorphic sites. Phylogenetic analysis were carried out among the species of Bovidae to reconstruct their relationships. Results From the phylogenetic analysis with the consensus sequences of the studied Bovidae species, we found that only 11 of the 26 nucleotide changes that differentiate them produced amino acid changes. All the SNPs found in the cattle breeds were novel and showed similar percentages of nucleotides with non-synonymous substitutions at the N-terminal, MHD and C-terminal (12.3, 12.8 and 12.5%, respectively), and were much higher than the percentage of synonymous substitutions (2.5, 2.6 and 4.9%, respectively). Three mutations in cattle and one in sheep, detected in heterozygous individuals were predicted to be deleterious. Additionally, the analysis of the biochemical characteristics in the most common form of the proteins in each species show very little difference in molecular weight, pI, net charge, instability index, aliphatic index and GRAVY (Table 4) in the Bovidae species, except for sheep, which had a higher molecular weight, instability index and GRAVY. Conclusions There is sufficient variation in this gene within and among the studied species, and because NDN carry key functions in the organism, it can have effects in economically important traits in the production of these species. NDN sequence is phylogenetically informative in this group, thus we propose this gene as a phylogenetic marker to study the evolution and conservation in Bovidae

Interethnic diversity of the CD209 (rs4804803) gene promoter polymorphism in African but not American sickle cell disease

Elucidating the genomic diversity of CD209 gene promoter polymorphism could assist in clarifying disease pathophysiology as well as contribution to co-morbidities. CD209 gene promoter polymorphism has been shown to be associated with susceptibility to infection. We hypothesize that CD209 mutant variants occur at a higher frequency among Africans and in sickle cell disease. We analyzed the frequency of the CD209 gene (rs4804803) in healthy control and sickle cell disease (SCD) populations and determined association with disease. Genomic DNA was extracted from blood samples collected from 145 SCD and 231 control Africans (from Mali), 331 SCD and 379 control African Americans and 159 Caucasians. Comparative analysis among and between groups was carried out by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Per ethnic diversification, we found significant disparity in genotypic (23.4% versus 16.9% versus 3.2%) and allelic frequencies (48.7% versus 42.1% versus 19.8%) of the homozygote mutant variant of the CD209 (snp 309A/G) gene promoter between Africans, African Americans and Caucasians respectively. Comparative evaluation between disease and control groups reveal a significant difference in genotypic (10.4% versus 23.4%; p = 0.002) and allelic frequencies (39.7% versus 48.7%; p = 0.02) of the homozygote mutant variant in African SCD and healthy controls respectively, an observation that is completely absent among Americans. Comparing disease groups, we found no difference in the genotypic (p = 0.19) or allelic (p = 0.72) frequencies of CD209 homozygote mutant variant between Africans and Americans with sickle cell disease. The higher frequency of CD209 homozygote mutant variants in the African control group reveals a potential impairment of the capacity to mount an immune response to infectious diseases, and possibly delineate susceptibility to or severity of infectious co-morbidities within and between groups

REVIEW: Sperm: Comparative morphology and function related to altered reproductive strategies and fertility in mammals

Reproduction and development are fascinating processes that exhibit significant differences among mammals. Successful propagation of species depends on male fertility and the ability of the sperm to fertilize and activate the egg and support early development. Remarkably, the specifics of how sperm structure and function have evolved and their crucial roles in fertility remain largely unknown. This synthesis provides a concise review of comparative anatomy and physiology of mammalian sperm structure and how inevitable changes regulate fertility

Analysis of culling records and estimation of genetic parameters for longevity and some production traits in Holstein dairy cattle

<p>We used data collected on 48,739 cows from 1982 to 2013 from Foka farm in the Isfahan province of Iran to assess the reasons for culling in Iranian Holstein dairy cattle in the context of breeding values of economic traits, and also to study the relationship between longevity and production traits. Estimation of (co)variance components and genetic parameters for studied traits including MY, FY, herd life (HL) and length of production life (LPL) was by AI-REML algorithm in WOMBAT software. The results showed that 27.11% of cows were culled voluntarily under the farmers’ command. The number and reasons for culling animals varied considerably over the years. Also, the results showed that involuntarily culled animals were valuable with significantly higher breeding values for MY and FY (<i>P</i> < .0001) than those of voluntarily culled individuals. The estimates of heritability for longevity traits (HL and LPL) were low (0.11 and 0.09), but their genetic correlation with production traits were high and positive (0.63–0.70). Our results suggest that some opportunity may exist for genetic improvements in fat and milk yield in Iranian Holsteins.</p

Nucleotide sequence variability analysis of Major Histocompatibility Complex Class II DQA1 gene in Nigerian goats

Major Histocompatibility Complex (MHC) molecules loaded with peptides derived from invading pathogens are recognised by the immune system to produce a highly effective and specific response against foreign pathogens. A 310-bp fragment of exon 2 of the MHC Class II DQA1 gene was amplified in 27 animals made up of three major Nigerian goat breeds [West African Dwarf (WAD), Red Sokoto (RS) and Sahel (SH)]. Twenty amino acid polymorphic sites were found in Nigerian goats. Comparison of predicted amino acid residues of DQA1 exon 2 alleles of Nigerian goats with similar alleles from other caprine species revealed considerable congruence in amino acid substitution pattern. A significant positive selection signature was detected at the DQA1 locus of Nigerian goats in that non-synonymous substitutions occurred at a faster rate compared to synonymous substitutions (dN:dS ratio = 1.28 ; Z-Statistics= 1.634; P<0.05). The evolutionary tree constructed using UPGMA, revealed that the southern WAD goat appeared to be more related to the northern RS than SH goat at the DQA1 locus. It will be interesting therefore, for future studies to investigate the association of the genetic variants in DQA1 gene of Nigerian goats with resistance/susceptiblity to diseases in order to conserve these precious animal genetic resources

Extensive Genomic Variability of Knops Blood Group Polymorphisms is Associated with Sickle Cell Disease in Africa

Sickle cell disease (SCD) is a multisystem disorder characterized by chronic hemolytic anemia, vaso-occlusive crises, and marked variability in disease severity. Patients require transfusions to manage disease complications, with complements, directed by complement regulatory genes (CR1) and its polymorphisms, implicated in the development of alloantibodies. We hypothesize that CR1 polymorphisms affect complement regulation and function, leading to adverse outcome in SCD. To this end, we determined the genomic diversity of complement regulatory genes by examining single nucleotide polymorphisms associated with Knops blood group antigens. Genomic DNA samples from 130 SCD cases and 356 control Africans, 331 SCD cases and 497 control African Americans, and 254 Caucasians were obtained and analyzed, utilizing a PCR–-RFLP (polymerase chain reaction-restriction fragment length polymorphism) assay. Analyzing for ethnic diversity, we found significant differences in the genotypic and allelic frequencies of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between Africans, African Americans, and Caucasians ( P < 0.05). The homozygote mutant variants had significantly higher frequencies in Africans and African Americans but were insignificant in Caucasians (80.2% and 59.6% vs 5.9% for Sl1/2; and 36% and 24% vs 1.8% for McCa/b). With SCD, we did not detect any difference among cases and controls either in Africa or in the United States. However, we found significant difference in genotypic ( P < 0.0001) and allelic frequencies ( P < 0.0001) of Sl1/Sl2 (rs17047661) and McCa/b (rs17047660) polymorphisms between SCD groups from Africa and the United States. There was no difference in haplotype frequencies of these polymorphisms among or between groups. The higher frequency of CR1 homozygote mutant variants in Africa but not United States indicates a potential pathogenic role, possibly associated with complicated disease pathophysiology in the former and potentially protective in the latter. The difference between sickle cell groups suggests potential genetic drift or founder effect imposed on the disease in the United States, but not in Africa, and a possible confirmation of the ancestral susceptibility hypothesis. The lower haplotype frequencies among sickle cell and control populations in the United States may be due to the admixture and the dilution of African genetic ancestry in the African American population