Analysis of protein-protein interaction network based on transcriptome profiling of ovine granulosa cells identifies candidate genes in cyclic recruitment of ovarian follicles

Deciphering sperm functions using biological networks

The global human population is exponentially increasing, which requires the production of quality food through efficient reproduction as well as sustainable production of livestock. Lack of knowledge and technology for assessing semen quality and predicting bull fertility is hindering advances in animal science and food animal production and causing millions of dollars of economic losses annually. The intent of this systemic review is to summarize methods from computational biology for analysis of gene, metabolite, and protein networks to identify potential markers that can be applied to improve livestock reproduction, with a focus on bull fertility. We provide examples of available gene, metabolic, and protein networks and computational biology methods to show how the interactions between genes, proteins, and metabolites together drive the complex process of spermatogenesis and regulate fertility in animals. We demonstrate the use of the National Center for Biotechnology Information (NCBI) and Ensembl for finding gene sequences, and then use them to create and understand gene, protein and metabolite networks for sperm associated factors to elucidate global cellular processes in sperm. This study highlights the value of mapping complex biological pathways among livestock and potential for conducting studies on promoting livestock improvement for global food security.

Unveiling the population genetic structure of Iranian horses breeds by whole-genome resequencing analysis

Preserving genetic diversity is pivotal for enhancing genetic improvement and facilitating adaptive responses to selection. This study focuses on identifying key genetic variants, including single nucleotide polymorphisms (SNPs), insertion/deletion polymorphisms (INDELs), and copy number variants (CNVs), while exploring the genomic evolutionary connectedness among seven Iranian horses representing five indigenous breeds: Caspian, Turkemen, DareShuri, Kurdish, and Asil. Using whole-genome resequencing, we generated 2.7 Gb of sequence data, with raw reads ranging from 1.2 Gb for Caspian horses to 0.38 Gb for Turkoman horses. Post-filtering, approximately 1.9 Gb of reads remained, with ~ 1.5 Gb successfully mapped to the horse reference genome (EquCab3.0), achieving mapping rates between 76.4% (Caspian) and 98.35% (Turkoman). We identified 2,909,816 SNPs in Caspian horses, constituting around 0.1% of the genome. Notably, 71% of these SNPs were situated in intergenic regions, while 8.5 and 6.8% were located upstream and downstream, respectively. A comparative analysis of SNPs between Iranian and non-Iranian horse breeds showed that Caspian horses had the lowest number of shared SNPs with Turkoman horses. Instead, they showed a closer genetic relationship with DareShuri, Quarter, Arabian, Standardbred, and Asil breeds. Hierarchical clustering highlighted Caspian horses as a distinct cluster, underscoring their distinctive genomic signature. Caspian horses exhibit a unique genetic profile marked by an enrichment of private mutations in neurological genes, influencing sensory perception and awareness. This distinct genetic makeup shapes mating preferences and signifies a separate evolutionary trajectory. Additionally, significant non-synonymous single nucleotide polymorphisms (nsSNPs) in reproductive genes offer intervention opportunities for managing Caspian horses. These findings reveal the population genetic structure of Iranian horse breeds, contributing to the advancement of knowledge in areas such as conservation, performance traits, climate adaptation, reproduction, and resistance to diseases in equine science.

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds

Hu sheep (HS), a breed of sheep carrying the FecB mutation gene, is known for its "year-round estrus and multiple births" and is an ideal model for studying the high fecundity mechanisms of livestock. Through analyzing and comparing the genomic selection features of Hu sheep and other sheep breeds, we identified a series of candidate genes that may play a role in Hu sheep's high fecundity mechanisms. In this study, we conducted whole-genome resequencing on six breeds and screened key mutations significantly correlated with high reproductive traits in sheep. Notably, the CC2D1B gene was selected by the fixation index (FST) and the cross-population composite likelihood ratio (XP-CLR) methods in HS and other five breeds. It was worth noting that the CC2D1B gene in HS was different from that in other sheep breeds, and seven missense mutations have been identified. Furthermore, the linkage disequilibrium (LD) analysis revealed a strong linkage disequilibrium in this specific gene region. Subsequently, by performing different grouping based on FecB genotypes in Hu sheep, genome-wide selective signal analysis screened several genes related to reproduction, such as BMPR1B and PPM1K. Besides, FST analysis identified functional genes related to reproductive traits, including RHEB, HSPA2, PPP1CC, HVCN1, and CCDC63. Additionally, a missense mutation was found in the CCDC63 gene and the haplotype was different between the high reproduction (HR) group and low reproduction (LR) group in HS. In summary, we discovered genetic differentiation among six distinct breeding sheep breeds at the whole genome level. Additionally, we identified a set of genes which were associated with reproductive performance in Hu sheep and visualized how these genes differed in different breeds. These findings laid a theoretical foundation for understanding genetic mechanisms behind high prolific traits in sheep.

RNA-Seq-Based Transcriptome Analysis of Nitric Oxide Scavenging Response in Neurospora crassa

While the biological role of naturally occurring nitric oxide (NO) in filamentous fungi has been uncovered, the underlying molecular regulatory networks remain unclear. In this study, we conducted an analysis of transcriptome profiles to investigate the initial stages of understanding these NO regulatory networks in Neurospora crassa, a well-established model filamentous fungus. Utilizing RNA sequencing, differential gene expression screening, and various functional analyses, our findings revealed that the removal of intracellular NO resulted in the differential transcription of 424 genes. Notably, the majority of these differentially expressed genes were functionally linked to processes associated with carbohydrate and amino acid metabolism. Furthermore, our analysis highlighted the prevalence of four specific protein domains (zinc finger C2H2, PLCYc, PLCXc, and SH3) in the encoded proteins of these differentially expressed genes. Through protein-protein interaction network analysis, we identified eight hub genes with substantial interaction connectivity, with mss-4 and gel-3 emerging as possibly major responsive genes during NO scavenging, particularly influencing vegetative growth. Additionally, our study unveiled that NO scavenging led to the inhibition of gene transcription related to a protein complex associated with ribosome biogenesis. Overall, our investigation suggests that endogenously produced NO in N. crassa likely governs the transcription of genes responsible for protein complexes involved in carbohydrate and amino acid metabolism, as well as ribosomal biogenesis, ultimately impacting the growth and development of hyphae.

Transcriptome-wide N6-methyladenine methylation in granulosa cells of women with decreased ovarian reserve

Background

The emerging epitranscriptome plays an essential role in female fertility. As the most prevalent internal mRNA modification, N6-methyladenine (m6A) methylation regulate mRNA fate and translational efficiency. However, whether m6A methylation was involved in the aging-related ovarian reserve decline has not been investigated. Herein, we performed m6A transcriptome-wide profiling in the ovarian granulosa cells of younger women (younger group) and older women (older group).

Results

m6A methylation distribution was highly conserved and enriched in the CDS and 3’UTR region. Besides, an increased number of m6A methylated genes were identified in the older group. Bioinformatics analysis indicated that m6A methylated genes were enriched in the FoxO signaling pathway, adherens junction, and regulation of actin cytoskeleton. A total of 435 genes were differently expressed in the older group, moreover, 58 of them were modified by m6A. Several specific genes, including BUB1B, PHC2, TOP2A, DDR2, KLF13, and RYR2 which were differently expressed and modified by m6A, were validated using qRT-PCR and might be involved in the decreased ovarian functions in the aging ovary.

Conclusions

Hence, our finding revealed the transcriptional significance of m6A modifications and provide potential therapeutic targets to promote fertility reservation for aging women.

Characterization of cadmium-responsive transcription factors in wolf spider Pardosa pseudoannulata

Transcription factors (TFs) act on the regulation of gene expression, which is prevalent in all organisms, and their characterization may provide important clues for understanding the regulatory mechanism of gene expression. In this research, inhibited growth (delayed developmental time and decreased body weight) and increased activities of antioxidant enzymes (peroxidase, superoxide dismutase, and catalase) were recorded in Pardosa pseudoannulata in response to cadmium burden. Expression profiles of TFs were analyzed based on the transcriptome profiling of P. pseudoannulata, and 1711 TFs genes were differentially expressed with 995 up-regulated and 716 down-regulated. Most of the differentially expressed TFs belonged to zf-C2H2, ZBTB, Homeobox, and bHLH families. Interestingly, hub genes smads, TCF7L2, EGR1, and GATA5 were identified to be the candidate Cd-responsive TFs related to growth of spider. The expression level of Sod2 (superoxide dismutase) was regulated by the up-regulated TF foxo3, implying its important role in the antioxidant defense of spider. Moreover, sequence analysis demonstrated that smads and foxo3 were conserved among spiders and insects. This study revealed for the first time the role of TFs in molecular response of P. pseudoannulata to Cd stress, providing the basis for the protection of tarantula under Cd stress.

Transcriptome Analysis of the Chicken Follicular Theca Cells with miR-135a-5p Suppressed

As a class of transcription regulators, numerous miRNAs have been verified to participate in regulating ovary follicular development in chickens (Gallus gallus). Previously we showed that gga-miR-135a-5p has significant differential expression between high and low-yield chicken ovaries, and the abundance of gga-miR-135a-5p is significantly higher in follicular theca cells than in granulosa cells. However, the exact role of gga-miR-135a-5p in chicken follicular theca cells is unclear. In this study, primary chicken follicular theca cells were isolated and then transfected with gga-miR-135a-5p inhibitor. Transcriptome sequencing was performed in chicken follicular theca cells with or without transfection. Differentially expressed genes (DEGs) were analyzed using bioinformatics. A dual-luciferase reporter assay was used to verify the target relationship between gga-miR-135a-5p and predicted targets within the DEGs. Compared with the normal chicken follicle theca cells, 953 up-regulated and 1060 down-regulated genes were detected in cells with gga-miR-135a-5p inhibited. The up-regulated genes were significantly enriched in Gene Ontology terms and pathways involved in cell proliferation and differentiation. In chicken follicular theca cells, Krüppel-like factor 4 (KLF4), ATPase phospholipid transporting 8A1 (ATP8A1), and Complexin-1 (CPLX1) were significantly up-regulated when the expression of gga-miR-135a-5p was inhibited. In addition, KLF4, ATP8A1, and CPLX1 confirmed as targets of gga-miR-135a-5p by using a dual-luciferase assay in vitro. The results suggest that gga-mir-135a-5p may involve in proliferation and differentiation in chicken ovarian follicular theca cells by targeting KLF4, ATP8A1, and CPLX1.

Genetic Analyses of Tanzanian Local Chicken Ecotypes Challenged with Newcastle Disease Virus

Newcastle Disease (ND) is a continuing global threat to domestic poultry, especially in developing countries, where severe outbreaks of velogenic ND virus (NDV) often cause major economic losses to households. Local chickens are of great importance to rural family livelihoods through provision of high-quality protein. To investigate the genetic basis of host response to NDV, three popular Tanzanian chicken ecotypes (regional populations) were challenged with a lentogenic (vaccine) strain of NDV at 28 days of age. Various host response phenotypes, including anti-NDV antibody levels (pre-infection and 10 days post-infection, dpi), and viral load (2 and 6 dpi) were measured, in addition to growth rate. We estimated genetic parameters and conducted genome-wide association study analyses by genotyping 1399 chickens using the Affymetrix 600K chicken SNP chip. Estimates of heritability of the evaluated traits were moderate (0.18–0.35). Five quantitative trait loci (QTL) associated with growth and/or response to NDV were identified by single-SNP analyses, with some regions explaining ≥1% of genetic variance based on the Bayes-B method. Immune related genes, such as ETS1, TIRAP, and KIRREL3, were located in regions associated with viral load at 6 dpi. The moderate estimates of heritability and identified QTL indicate that NDV response traits may be improved through selective breeding of chickens to enhance increased NDV resistance and vaccine efficacy in Tanzanian local ecotypes.