BMC Genomics最新文献

Background: Staphylococcus aureus (S. aureus) mastitis results in economic losses during dairy production. Understanding the biological progression of bovine S. aureus mastitis is vital for its prevention. Lipoteichoic acid is a key virulence factor of S. aureus (aLTA), but the main biological pathways involved in its effect on bovine mammary epithetionallial cells (Mac-T) apoptosis and necrosis have not been fully explored. Folic acid (FA) has anti-inflammatory and anti-apoptotic effects. However, the role of FA in mediating the effects of aLTA on apoptosis and necrosis remains unknown.

Results: We found that low concentration of aLTA inhibited apoptosis and necrosis and that high concentration promoted the apoptosis and necrosis of Mac-T. FA pretreatment alleviated high concentration of aLTA induced apoptosis. Through transcriptomic analysis, we found that nuclear receptor subfamily 4 group A (NR4A), which alters the expression of downstream genes involved in apoptosis, proliferation, and inflammation, decreased under stimulation with a low concentration of aLTA and increased under stimulation with a high concentration of aLTA. Under stimulation with a high concentration of aLTA, the expression of the NR4A subfamily could be inhibited by FA. The results showed that aLTA may affect apoptosis and necrosis through the NR4A subfamily by targeting genes involved in bacterial invasion of epithelial cells, the IL-17 signaling pathway, DNA replication, longevity regulation, the cell cycle, and tight junction pathways. We further found that the expression trends of NR4A1 and the target genes of the NR4A subfamily (PTGS2, ESPL1, MCM5, and BUB1B) in the blood of healthy cows (Healthy), subclinical mastitis cows (SCM), and SCM supplemented with FA (SCM_FA) were consistent with those observed at the cellular level in this study.

Conclusions: Our study revealed that low and high concentrations of aLTA have opposite effects on apoptosis and necrosis of Mac-T and that FA can alleviate the apoptosis induced by high concentration of aLTA. Transcriptome analysis revealed that the NR4A subfamily play a role in the ability of FA to alleviate the apoptosis and necrosis induced by high concentration of aLTA.

Background: Although metabarcoding and metagenomic approaches have generated large datasets on worldwide phytoplankton species diversity, the intraspecific genetic diversity underlying the genetic adaptation of marine phytoplankton to specific environmental niches remains largely unexplored. This is mainly due to the lack of biological resources and tools for monitoring the dynamics of this diversity in space and time.

Results: To gain insight into population diversity, a novel method based on INDEL markers was developed on Bathycoccus prasinos (Mamiellophyceae), an abundant and cosmopolitan species with strong seasonal patterns. Long read sequencing was first used to characterize structural variants among the genomes of six B. prasinos strains sampled from geographically distinct regions in the world ocean. Markers derived from identified insertions/deletions were validated by PCR then used to genotype 55 B. prasinos strains isolated during the winter bloom 2018-2019 in the bay of Banyuls-sur-Mer (Mediterranean Sea, France). This led to their classification into eight multi-loci genotypes and the sequencing of strains representative of local diversity, further improving the available genetic diversity of B. prasinos. Finally, selected markers were directly tracked on environmental DNA sampled during 3 successive blooms from 2018 to 2021, showcasing a fast and cost-effective approach to follow local population dynamics.

Conclusions: This method, which involves (i) pre-identifying the genetic diversity of B. prasinos in environmental samples by PCR, (ii) isolating cells from selected environmental samples and (iii) identifying genotypes representative of B. prasinos diversity for sequencing, can be used to comprehensively describe the diversity and population dynamics not only in B. prasinos but also potentially in other generalist phytoplankton species.

Background: The placenta is essential for nutrient exchange and hormone production between the mother and the developing fetus and serves as an invaluable model for epigenetic research. Most epigenetic studies of the human placenta have used whole placentas from term pregnancies and have identified the presence of partially methylated domains (PMDs). However, the origin of these domains, which are typically absent in most somatic cells, remains unclear in the placental context.

Results: Using whole-genome bisulfite sequencing and analysis of histone H3 modifications, we generated epigenetic profiles of human cytotrophoblasts during the first trimester and at term, as well as human trophoblast stem cells. Our study focused specifically on PMDs. We found that genomic regions likely to form PMDs are resistant to global DNA demethylation during trophectoderm reprogramming, and that PMDs arise through a slow methylation process within condensed chromatin near the nuclear lamina. In addition, we found significant differences in histone H3 modifications between PMDs in cytotrophoblasts and trophoblast stem cells.

Conclusions: Our findings suggest that spatiotemporal genomic features shape megabase-scale DNA methylation patterns, including PMDs, in the human placenta and highlight distinct differences in PMDs between human cytotrophoblasts and trophoblast stem cells. These findings advance our understanding of placental biology and provide a basis for further research into human development and related diseases.

Background: Echolocation was a key development in toothed whale evolution, enabling their adaptation and diversification across various environments. Previous bioacoustic and morphological studies suggest that environmental pressures have influenced the evolution of echolocation in toothed whales. This hypothesis demands further investigation, especially regarding the molecular mechanisms involved in the adaptive radiation of toothed whales across multiple habitats. Here we show that the coding sequences of four hearing genes involved in echolocation (CDH23, prestin, TMC1, and CLDN14) have different signatures of molecular evolution among riverine, coastal, and oceanic dolphins, suggesting that the evolutionary constraints of these habitats shaped the underlying genetic diversity of the toothed whale sonar.

Results: Our comparative analysis across 37 odontocete species revealed patterns of accelerated evolution within coastal and riverine lineages, supporting the hypothesis that shallow habitats pose specific selective pressures to sonar propagation, which are not found in the deep ocean. All toothed whales with genes evolving under positive selection are shallow coastal species, including three species that have recently diverged from freshwater lineages (Cephalorhynchus commersonii, Sotalia guianensis, and Orcaella heinsohni - CDH23), and three species that operate specialized Narrow Band High Frequency (NBHF) Sonars (Phocoena sinus - prestin, Neophocaena phocaenoides - TMC1 and Cephalorhynchus commersonii - CDH23). For river dolphins and deep-diving toothed whales, we found signatures of positive selection and molecular convergence affecting specific sites on CDH23, TMC1, and prestin. Positively selected sites (PSS) were different in number, identity, and substitution rates (dN/dS) across riverine, coastal, and oceanic toothed whales.

Conclusion: Here we shed light on potential molecular mechanisms underlying the diversification of toothed whale echolocation. Our results suggest that toothed whale hearing genes changed under different selective pressures in coastal, riverine, and oceanic environments.

Brassinosteroids (BRs), one of the major classes of phytohormones are essential for various processes of plant growth, development, and adaptations to biotic and abiotic stresses. In Arabidopsis, AtCYP90D1 acts as a bifunctional cytochrome P450 monooxygenase, catalyzing C-23 hydroxylation in the brassinolide biosynthetic pathway. The present study reports the functional characterizations of PtoCYP90D1, one of the AtCYP90D1 homologous genes from Populus tomentosa. The qRT-PCR analysis showed that PtoCYP90D1 was highly expressed in roots and old leaves. Overexpression of PtoCYP90D1 (PtoCYP90D1-OE) in poplar promoted growth and biomass yield, as well as increased xylem area and cell layers. Transgenic plants exhibited a significant increase in plant height and stem diameter as compared to the wild type. In contrast, the CRISPR/Cas9-generated mutation of PtoCYP90D1 (PtoCYP90D1-KO) resulted in significantly decreased biomass production in transgenic plants. Further studies revealed that cell wall components increased significantly in PtoCYP90D1-OE lines but not in PtoCYP90D1-KO lines, as compared to wild-type plants. Overall, the findings indicate a positive role of PtoCYP90D1 in improving growth rate and elevating biomass production in poplar, which will have positive implications for its versatile industrial or agricultural applications.

Background: Odontotermes formosanus (Shiraki) is a highly damaging agroforestry pest. Serratia marcescens is a broad-spectrum insecticidal pathogen and is highly lethal to O. formosanus. However, little is known about the mechanism between them. To improve the biological control of pests, a more in-depth analysis of the interactions between the pests and the pathogens is essential.

Results: We used RNA-seq, enzyme activity assays and real-time fluorescent quantitative PCR (qPCR) to explore the defense responses of O. formosanus against SM1. RNA-seq results showed that 1,160, 2,531 and 4,536 genes were differentially expressed at 3, 6 and 12 h after SM1 infection, and Kyoto Encyclopedia of Genes and Genomes (KEGG) results indicated that immune response and energy metabolism were involved in the defense of O. formosanus against SM1. Reactive oxygen species (ROS) levels and ROS synthesis genes were significantly elevated, and the antioxidant system were induced in O. formosanus after SM1 infection. In addition, the cellular immune genes were affected, and the Toll, immune deficiency (Imd), Janus kinase/signal transducer and activator of transcription (JAK/STAT), c-Jun N-terminal Kinase (JNK) and melanization pathways were activated. In vitro, Oftermicin, an antimicrobial peptide, had a significantly inhibitory effect on SM1. Furthermore, the expression levels and enzyme activities of phosphofructokinase (PFK), lactate dehydrogenase (LDH), succinate dehydrogenase (SDH) and isocitrate dehydrogenase (IDH) in glycolysis and tricarboxylic acid (TCA) cycles were increased.

Conclusions: Our results clearly demonstrated that O. formosanus defended against SM1 by activating the antioxidant system, innate immunity and energy metabolism. This study would provide useful information for the development of biological controls of O. formosanus.

Genetic profiling of the biodiversity in cultivated crop plants is necessary to preserve important genes and utilize them in a breeding program. Cucumber is used as a model plant to study various characteristics of Cucurbitaceae. Its adaptation to a wide range of climatic conditions suggested analyzing the landraces. The present study was conducted to evaluate the differences, at the genetic level, among landraces spanning five continents. DNA extracted from fifty-six landraces selected from USDA germplasm bank to cover a global representative sample of world cucumber landraces was used for polymerase chain reaction using twenty-eight polymorphic expressed sequence tags simple sequence repeat (EST-SSR) markers. Twenty-eight EST-SSR markers covering all seven chromosomes yielded 98 bands with an average of 3.42 bands per marker. Polymorphic information content ranged from 0.00 (EC35) to 0.74 (EC17) with an average of 0.34. Six clusters provided an appropriate summary of the variation among the landraces, with the two largest groups including 32 (Asiatic) and 17 (European and American) landraces, respectively. Four small groups, three with two members, and one with one member (PI 525155-Egypt) were dissimilar to the two main groups. Landraces from the same region were often clustered together. Genetic similarity of the landraces was revealed by marker banding patterns. The locations of genetic diversity for cucumber landraces can be identified from this study.

Background: Embryonic genome activation (EGA) is a critical step in early embryonic development, as it marks the transition from relying on maternal factors to the initiation of transcription from embryo's own genome. The factors associated with EGA are not well understood and need further investigation. PRD-like (PRDL) homeodomain transcription factors (TFs) are considered to play crucial roles in this early event during development but these TFs have evolved differently, even within mammalian lineages. Different numbers of PRDL TFs have been predicted in bovine (Bos taurus); however, their divergent evolution requires species-specific confirmation and functional investigations.

Results: In this study, we conducted molecular cloning of mRNAs for the PRDL TFs ARGFX, DUXA, LEUTX, NOBOX, TPRX1, TPRX2, and TPRX3 in bovine oocytes or in vitro fertilized (IVF) preimplantation embryos. Our results confirmed the expression of PRDL TF genes in early bovine development at the cDNA level and uncovered their structures. For each investigated PRDL TF gene, we isolated at least one homeodomain-encoding cDNA fragment, indicative of DNA binding and thus potential role in transcriptional regulation in developing bovine embryos. Additionally, our cDNA cloning approach allowed us to reveal breed-related differences in bovine, as evidenced by the identification of a high number of single nucleotide variants (SNVs) across the PRDL class homeobox genes. Subsequently, we observed the prediction of the 9aa transactivation domain (9aaTAD) motif in the putative protein sequence of TPRX3 leading us to conduct functional analysis of this gene. We demonstrated that the TPRX3 overexpression in bovine fibroblast induces not only protein-coding genes but also short noncoding RNAs involved in splicing and RNA editing. We supported this finding by identifying a shared set of genes between our and published bovine early embryo development datasets.

Conclusions: Providing full-length cDNA evidence for previously predicted homeobox genes that belong to PRDL class improves the annotation of the bovine genome. Updating the annotation with seven developmentally-important genes will enhance the accuracy of RNAseq analysis with datasets derived from bovine preimplantation embryos. In addition, the absence of TPRX3 in humans highlights the species-specific and TF-specific regulation of biological processes during early embryo development.

Background: The Japanese cedar (Cryptomeria japonica D. Don) is one of the most important Japanese forest trees, occupying approximately 44% of artificial forests and planted in East Asia, the Azores Archipelago, and certain islands in the Indian Ocean. Although the huge genome of the species (ca. 9 Gbp) with abundant repeat elements may have represented an obstacle for genetic analysis, this species is easily propagated by cutting, flowered by gibberellic acid, transformed by Agrobacterium, and edited by CRISPR/Cas9. These characteristics of C. japonica recommend it as a model conifer species for which reference genome sequences are necessary.

Results: Herein, we report the first chromosome-level assembly of C. japonica (2n = 22) using third-generation selfed progeny (estimated homozygosity rate = 0.96). Young leaf tissue was used to extract high molecular weight DNA (> 50 kb) for HiFi PacBio long-read sequencing and to construct an Hi-C/Omni-C library for Illumina short-read sequencing. The 29× and 26× genome coverage of HiFi and Illumina reads, respectively, for de novo assembly yielded 2,651 contigs (9.1 Gbp, N50 contig size 12.0 Mbp). Hi-C analysis mapped 97% of the nucleotides on 11 chromosomes. The assembly was verified through comparison with a consensus linkage map comprising 7,781 markers. BUSCO analysis identified ∼ 91% conserved genes.

Conclusions: Annotations of genes and comparisons of repeat elements with other Cupressaceae and Pinaceae species provide a fundamental resource for conifer research.

The current investigation was undertaken to elucidate the population-stratifying and ancestry-informative markers in Indian, Chinese, and wild yak populations using whole genome resequencing (WGS) analysis while employing various selection strategies (Delta, Pairwise Wright's Fixation Index-F_ST, and Informativeness of Assignment) and marker densities (5-25 thousand). The study used WGS data on 105 individuals from three separate yak cohorts i.e., Indian yak (n = 29), Chinese yak (n = 61), and wild yak (n = 15). Variant calling in the GATK program with strict quality control resulted in 1,002,970 high-quality and independent (LD-pruned) SNP markers across the yak autosomes. Analysis was undertaken in toolbox for ranking and evaluation of SNPs (TRES) program wherein three different criteria i.e., Delta, Pairwise Wright's Fixation Index-F_ST, and Informativeness of Assignment were employed to identify population-stratifying and ancestry-informative markers across various datasets. The top-ranked 5,000 (5K), 10,000 (10K), 15,000 (15K), 20,000 (20K), and 25,000 (25K) SNPs were identified from each dataset while their composition and performance was assessed using different criteria. The average genomic breed clustering of Indian, Chinese, and wild yak cohorts with full density dataset (105 individuals with 1,002,970 markers) was 81.74%, 80.02%, and 83.62%, respectively. Informativeness of Assignment criterion with 10K density emerged as the best combination for three yak cohorts with 86.94%, 96.46%, and 98.20% clustering for Indian, Chinese, and wild yak, respectively. There was an average increase of 7.56%, 22.72%, and 30.35% in genomic breed clustering scores of Indian, Chinese, and wild yak cohorts over the estimates of the original dataset. The selected markers showed overlap multiple protein-coding genes within a 10 kb window including ADGRB3, ANK1, CACNG7, CALN1, CHCHD2, CREBBP, GLI3, KHDRBS2, and OSBPL10. This is the first report ever on elucidating low-density SNP marker sets with population-stratifying and ancestry-informative properties in three yak groups using WGS data. The results gain significance for application of genomic selection using cost-effective low-density SNP panels in global yak species.