Genome-wide association study of <Emphasis Type="Italic">Mycobacterium avium</Emphasis> subspecies <Emphasis Type="Italic">Paratuberculosis</Emphasis> infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50?K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50?K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P?<?5?×?10^??5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1?Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5?Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

     

CYP1A1 based on metabolism of xenobiotics by cytochrome P450 regulates chicken male germ cell differentiation

This study aimed to explore the regulatory mechanism of metabolism of xenobiotics by cytochrome P450 during the differentiation process of chicken embryonic stem cells (ESCs) into spermatogonial stem cells (SSCs) and consummate the induction differentiation system of chicken embryonic stem cells (cESCs) into SSCs in vitro. We performed RNA-Seq in highly purified male ESCs, male primordial germ cells (PGCs), and SSCs that are associated with the male germ cell differentiation. Thereinto, the metabolism of xenobiotics by cytochrome P450 was selected and analyzed with Venny among male ESC vs male PGC, male PGC vs SSC, and male ESC vs SSC groups and several candidates differentially expressed genes (DEGs) were excavated. Finally, quantitative real-time PCR (qRT-PCR) detected related DEGs under the condition of retinoic acid (RA) induction in vitro, and the expressions were compared with RNA-Seq. By knocking down CYP1A1, we detected the effect of CYP1A1-mediated metabolism of xenobiotics by cytochrome P450 on male germ cell differentiation by qRT-PCR and immunocytochemistry. Results showed that 17,742 DEGs were found during differentiation of ESCs into SSCs and enriched in 72 differently significant pathways. Thereinto, the metabolism of xenobiotics by cytochrome P450 was involved in the whole differentiation process of ESCs into SSCs and several candidate DEGs: CYP1A1, CYP3A4, CYP2D6, ALDH3B1, and ALDH1A3 were expressed with the same trend with RNA-Seq. Knockdown of CYP1A1 caused male germ cell differentiation under restrictions. Our findings showed that the metabolism of xenobiotics by cytochrome P450 was significantly different during the process of male germ cell differentiation and was persistently activated when we induced cESCs to differentiate into SSCs with RA in vitro, which illustrated that the metabolism of xenobiotics by cytochrome P450 played a crucial role in the differentiation process of ESCs into SSCs.     

Developmental transcriptome analysis of floral transition in <Emphasis Type="Italic">Rosa odorata</Emphasis> var. <Emphasis Type="Italic">gigantea</Emphasis>

Key message

Expression analyses revealed that floral transition of Rosa odorata var. gigantea is mainly regulated by VRN1, COLs, DELLA and KSN, with contributions by the effects of phytohormone and starch metabolism.

Abstract

Seasonal plants utilize changing environmental and developmental cues to control the transition from vegetative growth to flowering at the correct time of year. This study investigated global gene expression profiles at different developmental stages of Rosa odorata var. gigantea by RNA-sequencing, combined with phenotypic characterization and physiological changes. Gene ontology enrichment analysis of the differentially expressed genes (DEGs) between four different developmental stages (vegetative meristem, pre-floral meristem, floral meristem and secondary axillary buds) indicated that DNA methylation and the light reaction played a large role in inducing the rose floral transition. The expression of SUF and FLC, which are known to play a role in delaying flowering until vernalization, was down-regulated from the vegetative to the pre-floral meristem stage. In contrast, the expression of VRN1, which promotes flowering by repressing FLC expression, increased. The expression of DELLA proteins, which function as central nodes in hormone signaling pathways, and probably involve interactions between GA, auxin, and ABA to promote the floral transition, was well correlated with the expression of floral integrators, such as AGL24, COL4. We also identified DEGs associated with starch metabolism correlated with SOC1, AGL15, SPL3, AGL24, respectively. Taken together, our results suggest that vernalization and photoperiod are prominent cues to induce the rose floral transition, and that DELLA proteins also act as key regulators. The results summarized in the study on the floral transition of the seasonal rose lay a foundation for further functional demonstration, and have profound economic and ornamental values.

     

Bioinformatic prediction and analysis of glucolipid metabolic regulation by miR-34a in <Emphasis Type="Italic">Megalobrama amblycephala</Emphasis>

The objective of this study was to analyze the target genes and regulatory function of miR-34a in Megalobrama amblycephala using second-generation high-throughput sequencing and bioinformatic tools. Functional enrichment analysis was performed by gene ontology. MiR-34a and target gene expression levels were measured in M. amblycephala fed normal and high-carbohydrate diets. The results revealed that miR-34a was highly conserved in several species, and miR-34a of M. amblycephala has a close evolutionary relationship to that of zebrafish and common carp. miRanda, TargetScan, RNAhybrid predicted 5,185, 6,282 and 2,168 target genes, respectively, and 645 target genes were in common. According to annotation information, the target genes were enriched in phosphate metabolism, glycerophospholipid metabolism, Golgi vesicle transport, cell division, and other biological processes (P?<?0.05). Pathway enrichment analysis revealed that these target genes were mainly enriched in alpha-linolenic acid and linoleic acid metabolism, ether lipid metabolism, VEGF signaling pathway, Fc epsilon RI signaling pathway, GnRH signaling pathway, and MAPK signaling pathway (P?<?0.05). The regulatory role of miR-34a was more significant in the liver than in the brain of M. amblycephala. MiR-34a regulates glucose lipid homeostasis induced by high glucose diets by upregulating hepatic PI3K/Akt, FOXO, and TOR signaling pathways.     

Accumulation and tolerance characteristics of lead in <Emphasis Type="Italic">Althaea rosea</Emphasis> Cav. and <Emphasis Type="Italic">Malva crispa</Emphasis> L.

Two ornamental plants of Althaea rosea Cav. and Malva crispa L. were exposed to various concentrations of lead (Pb) (0, 50, 100, 200 and 500 mg·kg^?1) for 70 days to evaluate the accumulating potential and the tolerance characteristics. The results showed that both plant species grown normally under Pb stress, and A. rosea had a higher tolerance than M. crispa, while M. crispa had a higher ability in Pb accumulation than A. rosea. Besides, lower Pb concentration (50 mg·kg^?1) stimulated the shoot biomass in both plant species. Pb accumulation in plants was consistent with the increase of Pb levels, and the main accumulation sites were the roots and the older leaves. In addition, the photosynthetic pigments content and chlorophyll fluorescence parameters were influenced by Pb stress. In such case, both of the plants could improve the activities of antioxidant enzymes of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT) and ascorbate peroxidase (APX), and the contents of the total soluble sugar and soluble protein, which reached the highest value at Pb 100 mg·kg^?1, as well as the accumulation of the total thiols (T-SH) and non-protein thiols (NP-SH) to adapt to Pb stress. Thus, it provides the theoretical basis and possibility for ornamental plants of A. rosea and M. crispa in phytoremediation of Pb contaminated areas.