Animal & Avian Sciences Research Works
Permanent URI for this collectionhttp://hdl.handle.net/1903/1600
Browse
Item Allele-Specific Expression of CD4+ T Cells in Response to Marek’s Disease Virus Infection(MDPI, 2019-09-17) Bai, Hao; He, Yanghua; Ding, Yi; Carrillo, José A.; Selvaraj, Ramesh K.; Zhang, Huanmin; Chen, Jilan; Song, JiuzhouMarek’s disease (MD) is a T cell lymphoma disease induced by Marek’s disease virus (MDV), a highly oncogenic α herpesvirus primarily affecting chickens. MD is a chronic infectious disease that threatens the poultry industry. However, the mechanisms of genetic resistance for MD are complex and not completely understood. In this study, to identify high-confidence candidate genes of MD genetic resistance, high throughput sequencing (RNA-seq) was used to obtain transcriptomic data of CD4+ T cells isolated from MDV-infected and non-infected groups of two reciprocal crosses of individuals mating by two highly inbred chicken lines (63 MD-resistant and 72 MD-susceptible). After RNA-seq analysis with two biological replicates in each group, we identified 61 and 123 single nucleotide polymorphisms (SNPs) (false discovery rate (FDR) < 0.05) annotated in 39 and 132 genes in intercrosses 63 × 72 and 72 × 63, respectively, which exhibited allele-specific expression (ASE) in response to MDV infection. Similarly, we identified 62 and 79 SNPs annotated in 66 and 96 genes in infected and non-infected groups, respectively. We identified 534 and 1543 differentially expressed genes (DEGs) (FDR < 0.05) related to MDV infection in intercrosses 63 × 72 and 72 × 63, respectively. We also identified 328 and 20 DEGs in infected and non-infected groups, respectively. The qRT-PCR using seven DEGs further verified our results of RNA-seq analysis. The qRT-PCR of 11 important ASE genes was performed for gene functional validation in CD4+ T cells and tumors. Combining the analyses, six genes (MCL1, SLC43A2, PDE3B, ADAM33, BLB1, and DMB2), especially MCL1, were highlighted as the candidate genes with the potential to be involved in MDV infection. Gene-set enrichment analysis revealed that many ASE genes are linked to T cell activation, T cell receptor (TCR), B cell receptor (BCR), ERK/MAPK, and PI3K/AKT-mTOR signaling pathways, which play potentially important roles in MDV infection. Our approach underlines the importance of comprehensive functional studies for gaining valuable biological insight into the genetic factors behind MD and other complex traits, and our findings provide additional insights into the mechanisms of MD and disease resistance breeding in poultry.Item Analyses of inter-individual variations of sperm DNA methylation and their potential implications in cattle(Springer Nature, 2019-11-21) Liu, Shuli; Fang, Lingzhao; Zhou, Yang; Santos, Daniel J.A.; Xiang, Ruidong; Daetwyler, Hans D.; Chamberlain, Amanda J.; Cole, John B.; Li, Cong-jun; Yu, Ying; Ma, Li; Zhang, Shengli; Liu, George E.DNA methylation has been shown to be involved in many biological processes, including X chromosome inactivation in females, paternal genomic imprinting, and others.Item Analysis of recent segmental duplications in the bovine genome(Springer Nature, 2009-12-01) Liu, George E; Ventura, Mario; Cellamare, Angelo; Chen, Lin; Cheng, Ze; Zhu, Bin; Li, Congjun; Song, Jiuzhou; Eichler, Evan EDuplicated sequences are an important source of gene innovation and structural variation within mammalian genomes. We performed the first systematic and genome-wide analysis of segmental duplications in the modern domesticated cattle (Bos taurus). Using two distinct computational analyses, we estimated that 3.1% (94.4 Mb) of the bovine genome consists of recently duplicated sequences (≥ 1 kb in length, ≥ 90% sequence identity). Similar to other mammalian draft assemblies, almost half (47% of 94.4 Mb) of these sequences have not been assigned to cattle chromosomes. In this study, we provide the first experimental validation large duplications and briefly compared their distribution on two independent bovine genome assemblies using fluorescent in situ hybridization (FISH). Our analyses suggest that the (75-90%) of segmental duplications are organized into local tandem duplication clusters. Along with rodents and carnivores, these results now confidently establish tandem duplications as the most likely mammalian archetypical organization, in contrast to humans and great ape species which show a preponderance of interspersed duplications. A cross-species survey of duplicated genes and gene families indicated that duplication, positive selection and gene conversion have shaped primates, rodents, carnivores and ruminants to different degrees for their speciation and adaptation. We identified that bovine segmental duplications corresponding to genes are significantly enriched for specific biological functions such as immunity, digestion, lactation and reproduction. Our results suggest that in most mammalian lineages segmental duplications are organized in a tandem configuration. Segmental duplications remain problematic for genome and assembly and we highlight genic regions that require higher quality sequence characterization. This study provides insights into mammalian genome evolution and generates a valuable resource for cattle genomics research.Item Analysis of the fecal microbiota of fast- and slow-growing rainbow trout (Oncorhynchus mykiss)(Springer Nature, 2019-10-29) Chapagain, Pratima; Arivett, Brock; Cleveland, Beth M.; Walker, Donald M.; Salem, MohamedDiverse microbial communities colonizing the intestine of fish contribute to their growth, digestion, nutrition, and immune function. We hypothesized that fecal samples representing the gut microbiota of rainbow trout could be associated with differential growth rates observed in fish breeding programs. If true, harnessing the functionality of this microbiota can improve the profitability of aquaculture. The first objective of this study was to test this hypothesis if gut microbiota is associated with fish growth rate (body weight). Four full-sibling families were stocked in the same tank and fed an identical diet. Two fast-growing and two slow-growing fish were selected from each family for 16S rRNA microbiota profiling.Item Antimicrobial Effect and Probiotic Potential of Phage Resistant Lactobacillus plantarum and its Interactions with Zoonotic Bacterial Pathogens(MDPI, 2019-06-05) Nagarajan, Vinod; Peng, Mengfei; Tabashsum, Zajeba; Salaheen, Serajus; Padilla, Joselyn; Biswas, DebabrataDevelopment of phage-resistant probiotic particularly Lactobacillus is an alternative approach to enhance their beneficial effects as in animal feed supplements. In this study, we developed phage-resistant Lactobacillus plantarum (LP+PR) mutant and compared their antimicrobial effects and probiotic potential against zoonotic bacterial pathogens including Salmonella enterica serovar Typhimurium, enterohemorrhagic Escherichia coli (EHEC), Staphylococcus aureus, and Listeria monocytogenes with phage-sensitive L. plantarum (LP) strain. LP+PR strain showed markedly higher growth rate than wild-type LP strain. In co-culture with LP+PR and in the presence of cell-free cultural supernatants (CFCSs) of LP+PR, the growth of S. Typhimurium, EHEC, S. aureus, and L. monocytogenes were reduced significantly (P < 0.05). The adhesion ability of LP+PR was slightly higher than the LP on human epithelial INT-407 cells. Most importantly, LP+PR strain significantly inhibited the adhesive and invasive abilities of all four zoonotic pathogens to INT-407 cells (P < 0.05). Moreover, real-time qPCR revealed that in the presence of LP+PR strain or its CFCSs, expression of virulence genes of these zoonotic bacterial pathogens were suppressed significantly (P < 0.05). These findings suggest that the LP+PR strain is capable of inhibiting major zoonotic bacterial pathogens efficiently and would be a potential candidate for industrial usage in animal production or fermentation.Item Application of orange essential oil as an antistaphylococcal agent in a dressing model(Springer Nature, 2012-08-16) Muthaiyan, Arunachalam; Biswas, Debabrata; Crandall, Philip G; Wilkinson, Brian J; Ricke, Steven CStaphylococcus aureus is the pathogen most often and prevalently involved in skin and soft tissue infections. In recent decades outbreaks of methicillin-resistant S. aureus (MRSA) have created major problems for skin therapy, and burn and wound care units. Topical antimicrobials are most important component of wound infection therapy. Alternative therapies are being sought for treatment of MRSA and one area of interest is the use of essential oils. With the increasing interest in the use and application of natural products, we screened the potential application of terpeneless cold pressed Valencia orange oil (CPV) for topical therapy against MRSA using an in vitro dressing model and skin keratinocyte cell culture model. The inhibitory effect of CPV was determined by disc diffusion vapor assay for MRSA and vancomycin intermediate-resistant S. aureus (VISA) strains. Antistaphylococcal effect of CPV in an in vitro dressing model was tested on S. aureus inoculated tryptic soya agar plate. Bactericidal effect of CPV on MRSA and VISA infected keratinocyte cells was examined by enumeration of extra- and intra-cellular bacterial cells at different treatment time points. Cytotoxic effects on human skin cells was tested by adding CPV to the keratinocyte (HEK001) cells grown in serum free KSFM media, and observed by phase-contrast microscope. CPV vapour effectively inhibited the MRSA and VISA strains in both disc diffusion vapour assay and in vitro dressing model. Compared to untreated control addition of 0.1% CPV to MRSA infected keratinocyte decreased the viable MRSA cells by 2 log CFU/mL in 1 h and in VISA strain 3 log CFU/mL reduction was observed in 1 h. After 3 h viable S. aureus cells were not detected in the 0.2% CPV treatment. Bactericidal concentration of CPV did not show any cytotoxic effect on the human skin keratinocyte cells in vitro. At lower concentration addition of CPV to keratinocytes infected with MRSA and VISA rapidly killed the bacterial cells without causing any toxic effect to the keratinocytes. Therefore, the results of this study warrant further in vivo study to evaluate the potential of CPV as a topical antistaphylococcal agent.Item Are CD45RO+ and CD45RA- genuine markers for bovine memory T cells?(Springer Nature, 2022-10-11) Anmol, Kandel; Akanksha, Hada; Zhengguo, XiaoEffective vaccination induces memory T cells, which protect the host against pathogen re-infections. Therefore, detection of memory T cells is essential for evaluating vaccine efficacy, which was originally dependent on cytokine induction assays. Currently, two isoforms of CD45 tyrosine phosphatase, CD45RO expression and CD45RA exclusion (CD45RO+/ CD45RA-) are used extensively for detecting memory T cells in cattle. The CD45RO+/CD45RA- markers were first established in humans around three decades ago, and were adopted in cattle soon after. However, in the last two decades, some published data in humans have challenged the initial paradigm, and required multiple markers for identifying memory T cells. On the contrary, memory T cell detection in cattle still mostly relies on CD45RO+/CD45RA- despite some controversial evidence. In this review, we summarized the current literature to examine if CD45RO+/CD45RA- are valid markers for detecting memory T cells in cattle. It seems CD45RA and CD45RO (CD45RA/RO) as markers for identifying bovine memory T cells are questionable.Item Array CGH-based detection of CNV regions and their potential association with reproduction and other economic traits in Holsteins(Springer Nature, 2019-03-07) Liu, Mei; Fang, Lingzhao; Liu, Shuli; Pan, Michael G.; Seroussi, Eyal; Cole, John B.; Ma, Li; Chen, Hong; Liu, George E.Copy number variations (CNVs) are structural variants consisting of large-scale insertions and deletions of genomic fragments. Exploring CNVs and estimating their effects on phenotypes are useful for genome selection but remain challenging in the livestock.Item Avian reovirus: a furious and fast evolving pathogen(Microbiology Society, 2023-10) Egana-Labrin, Sofia; Broadbent, Andrew J.Avian reoviruses (ARVs) have a significant economic impact on the poultry industry, affecting commercial and backyard flocks. Spread feco-orally, or vertically, many do not cause morbidity, but pathogenic strains can contribute to several diseases, including tenosynovitis/arthritis, which is clinically the most significant. The last decade has seen a surge in cases in the US, and due to ongoing evolution, seven genotypic clusters have now been identified. Control efforts include strict biosecurity and vaccination with commercial and autogenous vaccines. Research priorities include improving understanding of pathogenesis and developing new vaccines guided by ongoing molecular and serologic surveillance.Item Bos taurus–indicus hybridization correlates with intralocus sexual-conflict effects of PRDM9 on male and female fertility in Holstein cattle(Springer Nature, 2019-08-28) Seroussi, Eyal; Shirak, Andrey; Gershoni, Moran; Ezra, Ephraim; Santos, Daniel Jordan de Abreu; Ma, Li; Liu, George E.Crossover localization during meiotic recombination is mediated by the fast-evolving zinc-finger (ZnF) domain of gene PRDM9. To study its impact on dairy cattle performance, we compared its genetic variation between the relatively small Israeli (IL) Holsteins and the North American (US) Holsteins that count millions.Item CeHR Axenic Liquid Medium for C. elegans growth(2005-02-24T14:24:16Z) Hamza, IqbalCeHR medium was originally formulated by Dr. Eric Clegg at the United States Army Center for Environmental Health Research to grow worms in an axenic liquid media. (http://usacehr.detrick.army.mil/CleggLab.html) We are thankful to Eric for his invaluable help and advice through all the modifications we have since made to optimize CeHR medium to suit our research goals. We modified CeHR medium (mCeHR) because of our interests in defining the individual micronutrient components in order to identify the cellular and molecular components of metal homeostasis using C. elegans as a model multicellular eukaryote.Item Changes in various metabolic parameters in blood and milk during experimental Escherichia coli mastitis for primiparous Holstein dairy cows during early lactation(Springer Nature, 2014-10-17) Moyes, Kasey M; Larsen, Torben; Sørensen, Peter; Ingvartsen, Klaus LThe objective of this study was to characterize the changes in various metabolic parameters in blood and milk during IMI challenge with Escherichia coli (E. coli) for dairy cows during early lactation. Thirty, healthy primiparous Holstein cows were infused (h = 0) with ~20-40 cfu of live E. coli into one front mammary quarter at ~4-6 wk in lactation. Daily feed intake and milk yield were recorded. At –12, 0, 3, 6, 12, 18, 24, 36, 48, 60, 72, 96, 108, 120, 132, 144, 156, 168, 180 and 192 h relative to challenge rectal temperatures were recorded and quarter foremilk was collected for analysis of shedding of E. coli. Composite milk samples were collected at -180, -132, -84, -36, -12, 12, 24, 36, 48, 60, 72, 84, 96, 132 and 180 h relative to challenge (h = 0) and analyzed for lactate dehydrogenase (LDH), somatic cell count, fat, protein, lactose, citrate, beta-hydroxybutyrate (BHBA), free glucose (fglu), and glucose-6-phosphate (G6P). Blood was collected at -12, 0, 3, 6, 12, 18, 24, 36, 60, 72, 84, 132 and 180 h relative to challenge and analyzed for plasma non-esterified fatty acids (NEFA), BHBA and glucose concentration. A generalized linear mixed model was used to determine the effect of IMI challenge on metabolic responses of cows during early lactation. By 12 h, E. coli was recovered from challenged quarters and shedding continued through 72 h. Rectal temperature peaked by 12 h post-challenge and returned to pre-challenge values by 36 h post-IMI challenge. Daily feed intake and milk yield decreased (P <0.05) by 1 and 2 d, respectively, after mastitis challenge. Plasma BHBA decreased (12 h; P <0.05) from 0.96 ± 1.1 at 0 h to 0.57 ± 0.64 mmol/L by 18 h whereas concentration of plasma NEFA (18 h) and glucose (24 h) were significantly greater, 11 and 27%, respectively, after challenge. In milk, fglu, lactose, citrate, fat and protein yield were lower whereas yield of BHBA and G6P were higher after challenge when compared to pre-challenge values. Changes in metabolites in blood and milk were most likely associated with drops in feed intake and milk yield. However, the early rise in plasma NEFA may also signify enhanced adipose tissue lipolysis. Lower concentrations of plasma BHBA may be attributed to an increase transfer into milk after IMI. Decreases in both milk lactose yield and % after challenge may be partly attributed to reduced conversion of fglu to lactose. Rises in G6P yield and concentration in milk after challenge (24 h) may signify increased conversion of fglu to G6P. Results identify changes in various metabolic parameters in blood and milk after IMI challenge with E. coli in dairy cows that may partly explain the partitioning of nutrients and changes in milk components after IMI for cows during early lactation.Item Characterisation of sRNAs enriched in outer membrane vesicles of pathogenic Flavobacterium psychrophilum causing Bacterial Cold Water Disease in rainbow trout(Wiley, 2024-06-19) Chapagain, Pratima; Ali, Ali; Kidane, Destaalem T.; Farone, Mary; Salem, MohamedFlavobacterium psychrophilum (Fp) causes Bacterial Cold Water Disease in salmonids. During host-pathogen interactions, gram-negative bacteria, such as Fp, release external membrane vesicles (OMVs) harbouring cargos, such as DNA, RNA and virulence factors. This study aimed to characterise the potential role of the OMVs’ small RNAs (sRNAs) in the Fp-rainbow trout host-pathogen interactions. sRNAs carried within OMVs were isolated from Fp. RNA-Seq datasets from whole-cell Fp and their isolated OMVs indicated substantial enrichment of specific sRNAs in the OMVs compared to the parent cell. Many of the OMV-packaged sRNAs were located in the pathogenicity islands of Fp. Conservation of sRNAs in 65 strains with variable degrees of virulence was reported. Dual RNA-Seq of host and pathogen transcriptomes on day 5 post-infection of Fp -resistant and -susceptible rainbow trout genetic lines revealed correlated expression of OMV-packaged sRNAs and their predicted host's immune gene targets. In vitro, treatment of the rainbow trout epithelial cell line RTgill-W1 with OMVs showed signs of cytotoxicity accompanied by dynamic changes in the expression of host genes when profiled 24 h following treatment. The OMV-treated cells, similar to the Fp -resistant fish, showed downregulated expression of the suppressor of cytokine signalling 1 (SOCS1) gene, suggesting induction of phagosomal maturation. Other signs of modulating the host gene expression following OMV-treatment include favouring elements from the phagocytic, endocytic and antigen presentation pathways in addition to HSP70, HSP90 and cochaperone proteins, which provide evidence for a potential role of OMVs in boosting the host immune response. In conclusion, the study identified novel microbial targets and inherent characteristics of OMVs that could open up new avenues of treatment and prevention of fish infections.Item Characterization of Copy Number Variation’s Potential Role in Marek’s Disease(MDPI, 2017-05-09) Xu, Lingyang; He, Yanghua; Ding, Yi; Sun, Guirong; Carrillo, Jose Adrian; Li, Yaokun; Ghaly, Mona M.; Ma, Li; Zhang, Huanmin; Liu, George E.; Song, JiuzhouMarek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.Item Characterization of recombination features and the genetic basis in multiple cattle breeds(Springer Nature, 2018-04-27) Shen, Botong; Jiang, Jicai; Seroussi, Eyal; Liu, George E.; Ma, LiCrossover generated by meiotic recombination is a fundamental event that facilitates meiosis and sexual reproduction. Comparative studies have shown wide variation in recombination rate among species, but the characterization of recombination features between cattle breeds has not yet been performed. Cattle populations in North America count millions, and the dairy industry has genotyped millions of individuals with pedigree information that provide a unique opportunity to study breed-level variations in recombination. Based on large pedigrees of Jersey, Ayrshire and Brown Swiss cattle with genotype data, we identified over 3.4 million maternal and paternal crossover events from 161,309 three-generation families. We constructed six breed- and sex-specific genome-wide recombination maps using 58,982 autosomal SNPs for two sexes in the three dairy cattle breeds. A comparative analysis of the six recombination maps revealed similar global recombination patterns between cattle breeds but with significant differences between sexes. We confirmed that male recombination map is 10% longer than the female map in all three cattle breeds, consistent with previously reported results in Holstein cattle. When comparing recombination hotspot regions between cattle breeds, we found that 30% and 10% of the hotspots were shared between breeds in males and females, respectively, with each breed exhibiting some breed-specific hotspots. Finally, our multiple-breed GWAS found that SNPs in eight loci affected recombination rate and that the PRDM9 gene associated with hotspot usage in multiple cattle breeds, indicating a shared genetic basis for recombination across dairy cattle breeds. Collectively, our results generated breed- and sex-specific recombination maps for multiple cattle breeds, provided a comprehensive characterization and comparison of recombination patterns between breeds, and expanded our understanding of the breed-level variations in recombination features within an important livestock species.Item Coding and noncoding genes involved in atrophy and compensatory muscle growth in Nile Tilapia(MDPI, 2022-08-12) Ali, Ali; Shaalan, Walaa M.; Al-Tobasei, Rafet; Salem, MohamedImprovements in growth-related traits reduce fish time and production costs to reach market size. Feed deprivation and refeeding cycles have been introduced to maximize aquaculture profits through compensatory growth. However, the molecular compensatory growth signature is still uncertain in Nile tilapia. In this study, fish were subjected to two weeks of fasting followed by two weeks of refeeding. The growth curve in refed tilapia was suggestive of a partial compensatory response. Transcriptome profiling of starved and refed fish was conducted to identify genes regulating muscle atrophy and compensatory growth. Pairwise comparisons revealed 5009 and 478 differentially expressed (differential) transcripts during muscle atrophy and recovery, respectively. Muscle atrophy appears to be mediated by the ubiquitin-proteasome and autophagy/lysosome systems. Autophagy-related 2A, F-box and WD repeat domain containing 7, F-box only protein 32, miR-137, and miR-153 showed exceptional high expression suggesting them as master regulators of muscle atrophy. On the other hand, the muscle compensatory growth response appears to be mediated by the continuous stimulation of muscle hypertrophy which exceeded normal levels found in control fish. For instance, genes promoting ribosome biogenesis or enhancing the efficiency of translational machinery were upregulated in compensatory muscle growth. Additionally, myogenic microRNAs (e.g., miR-1 and miR-206), and hypertrophy-associated microRNAs (e.g., miR-27a-3p, miR-29c, and miR-29c) were reciprocally expressed to favor hypertrophy during muscle recovery. Overall, the present study provided insights into the molecular mechanisms regulating muscle mass in fish. The study pinpoints extensive growth-related gene networks that could be used to inform breeding programs and also serve as valuable genomic resources for future mechanistic studies.Item Coding and Noncoding Genes Involved in Atrophy and Compensatory Muscle Growth in Nile Tilapia(MDPI, 2022-08-12) Ali, Ali; Shaalan, Walaa M.; Al-Tobasei, Rafet; Salem, MohamedImprovements in growth-related traits reduce fish time and production costs to reach market size. Feed deprivation and refeeding cycles have been introduced to maximize aquaculture profits through compensatory growth. However, the molecular compensatory growth signature is still uncertain in Nile tilapia. In this study, fish were subjected to two weeks of fasting followed by two weeks of refeeding. The growth curve in refed tilapia was suggestive of a partial compensatory response. Transcriptome profiling of starved and refed fish was conducted to identify genes regulating muscle atrophy and compensatory growth. Pairwise comparisons revealed 5009 and 478 differentially expressed (differential) transcripts during muscle atrophy and recovery, respectively. Muscle atrophy appears to be mediated by the ubiquitin-proteasome and autophagy/lysosome systems. Autophagy-related 2A, F-box and WD repeat domain containing 7, F-box only protein 32, miR-137, and miR-153 showed exceptional high expression suggesting them as master regulators of muscle atrophy. On the other hand, the muscle compensatory growth response appears to be mediated by the continuous stimulation of muscle hypertrophy which exceeded normal levels found in control fish. For instance, genes promoting ribosome biogenesis or enhancing the efficiency of translational machinery were upregulated in compensatory muscle growth. Additionally, myogenic microRNAs (e.g., miR-1 and miR-206), and hypertrophy-associated microRNAs (e.g., miR-27a-3p, miR-29c, and miR-29c) were reciprocally expressed to favor hypertrophy during muscle recovery. Overall, the present study provided insights into the molecular mechanisms regulating muscle mass in fish. The study pinpoints extensive growth-related gene networks that could be used to inform breeding programs and also serve as valuable genomic resources for future mechanistic studies.Item Comparative transcriptome in large-scale human and cattle populations(Springer Nature, 2022-08-22) Yao, Yuelin; Liu, Shuli; Xia, Charley; Gao, Yahui; Pan, Zhangyuan; Canela-Xandri, Oriol; Khamseh, Ava; Rawlik, Konrad; Wang, Sheng; Li, Bingjie; Zhang, Yi; Pairo-Castineira, Erola; D’Mellow, Kenton; Li, Xiujin; Yan, Ze; Li, Cong-jun; Yu, Ying; Zhang, Shengli; Ma, Li; Cole, John B.; Ross, Pablo J.; Zhou, Huaijun; Haley, Chris; Liu, George E.; Fang, Lingzhao; Tenesa, AlbertCross-species comparison of transcriptomes is important for elucidating evolutionary molecular mechanisms underpinning phenotypic variation between and within species, yet to date it has been essentially limited to model organisms with relatively small sample sizes. Here, we systematically analyze and compare 10,830 and 4866 publicly available RNA-seq samples in humans and cattle, respectively, representing 20 common tissues. Focusing on 17,315 orthologous genes, we demonstrate that mean/median gene expression, inter-individual variation of expression, expression quantitative trait loci, and gene co-expression networks are generally conserved between humans and cattle. By examining large-scale genome-wide association studies for 46 human traits (average n = 327,973) and 45 cattle traits (average n = 24,635), we reveal that the heritability of complex traits in both species is significantly more enriched in transcriptionally conserved than diverged genes across tissues. In summary, our study provides a comprehensive comparison of transcriptomes between humans and cattle, which might help decipher the genetic and evolutionary basis of complex traits in both species.Item Comparative whole genome DNA methylation profiling across cattle tissues reveals global and tissue-specific methylation patterns(Springer Nature, 2020-07-06) Zhou, Yang; Liu, Shuli; Hu, Yan; Fang, Lingzhao; Gao, Yahui; Xia, Han; Schroeder, Steven G.; Rosen, Benjamin D.; Connor, Erin E.; Li, Cong-jun; Baldwin, Ransom L.; Cole, John B.; Van Tassell, Curtis P.; Yang, Liguo; Ma, Li; Liu, George E.Efforts to improve animal health, and understand genetic bases for production, may benefit from a comprehensive analysis of animal genomes and epigenomes. Although DNA methylation has been well studied in humans and other model species, its distribution patterns and regulatory impacts in cattle are still largely unknown. Here, we present the largest collection of cattle DNA methylation epigenomic data to date. Using Holstein cattle, we generated 29 whole genome bisulfite sequencing (WGBS) datasets for 16 tissues, 47 corresponding RNA-seq datasets, and 2 whole genome sequencing datasets. We did read mapping and DNA methylation calling based on two different cattle assemblies, demonstrating the high quality of the long-read-based assembly markedly improved DNA methylation results. We observed large differences across cattle tissues in the methylation patterns of global CpG sites, partially methylated domains (PMDs), hypomethylated regions (HMRs), CG islands (CGIs), and common repeats. We detected that each tissue had a distinct set of PMDs, which showed tissue-specific patterns. Similar to human PMD, cattle PMDs were often linked to a general decrease of gene expression and a decrease in active histone marks and related to long-range chromatin organizations, like topologically associated domains (TADs). We tested a classification of the HMRs based on their distributions relative to transcription start sites (TSSs) and detected tissue-specific TSS-HMRs and genes that showed strong tissue effects. When performing cross-species comparisons of paired genes (two opposite strand genes with their TSS located in the same HMR), we found out they were more consistently co-expressed among human, mouse, sheep, goat, yak, pig, and chicken, but showed lower consistent ratios in more divergent species. We further used these WGBS data to detect 50,023 experimentally supported CGIs across bovine tissues and found that they might function as a guard against C-to-T mutations for TSS-HMRs. Although common repeats were often heavily methylated, some young Bov-A2 repeats were hypomethylated in sperm and could affect the promoter structures by exposing potential transcription factor binding sites. This study provides a comprehensive resource for bovine epigenomic research and enables new discoveries about DNA methylation and its role in complex traits.Item A computational model of the hypothalamic - pituitary - gonadal axis in female fathead minnows (Pimephales promelas) exposed to 17α-ethynylestradiol and 17β-trenbolone(Springer Nature, 2011-05-05) Li, Zhenhong; Kroll, Kevin J; Jensen, Kathleen M; Villeneuve, Daniel L; Ankley, Gerald T; Brian, Jayne V; Sepúlveda, María S; Orlando, Edward F; Lazorchak, James M; Kostich, Mitchell; Armstrong, Brandon; Denslow, Nancy D; Watanabe, Karen HEndocrine disrupting chemicals (e.g., estrogens, androgens and their mimics) are known to affect reproduction in fish. 17α-ethynylestradiol is a synthetic estrogen used in birth control pills. 17β-trenbolone is a relatively stable metabolite of trenbolone acetate, a synthetic androgen used as a growth promoter in livestock. Both 17α-ethynylestradiol and 17β-trenbolone have been found in the aquatic environment and affect fish reproduction. In this study, we developed a physiologically-based computational model for female fathead minnows (FHM, Pimephales promelas), a small fish species used in ecotoxicology, to simulate how estrogens (i.e., 17α-ethynylestradiol) or androgens (i.e., 17β-trenbolone) affect reproductive endpoints such as plasma concentrations of steroid hormones (e.g., 17β-estradiol and testosterone) and vitellogenin (a precursor to egg yolk proteins). Using Markov Chain Monte Carlo simulations, the model was calibrated with data from unexposed, 17α-ethynylestradiol-exposed, and 17β-trenbolone-exposed FHMs. Four Markov chains were simulated, and the chains for each calibrated model parameter (26 in total) converged within 20,000 iterations. With the converged parameter values, we evaluated the model's predictive ability by simulating a variety of independent experimental data. The model predictions agreed with the experimental data well. The physiologically-based computational model represents the hypothalamic-pituitary-gonadal axis in adult female FHM robustly. The model is useful to estimate how estrogens (e.g., 17α-ethynylestradiol) or androgens (e.g., 17β-trenbolone) affect plasma concentrations of 17β-estradiol, testosterone and vitellogenin, which are important determinants of fecundity in fish.