ALBERT

Description: Background Long non-coding RNAs (lncRNAs) are involved in several immune processes, including the immune response to vaccination, but most of them remain uncharacterised in livestock species. The mechanism of action of aluminium adjuvants as vaccine components is neither not fully understood. Results We built a transcriptome from sheep PBMCs RNA-seq data in order to identify unannotated lncRNAs and analysed their expression patterns along protein coding genes. We found 2284 novel lncRNAs and assessed their conservation in terms of sequence and synteny. Differential expression analysis performed between animals inoculated with commercial vaccines or aluminium adjuvant alone and the co-expression analysis revealed lncRNAs related to the immune response to vaccines and adjuvants. A group of co-expressed genes enriched in cytokine signalling and production highlighted the differences between different treatments. A number of differentially expressed lncRNAs were correlated with a divergently located protein-coding gene, such as the OSM cytokine. Other lncRNAs were predicted to act as sponges of miRNAs involved in immune response regulation. Conclusions This work enlarges the lncRNA catalogue in sheep and puts an accent on their involvement in the immune response to repetitive vaccination, providing a basis for further characterisation of the non-coding sheep transcriptome within different immune cells.

Description: Background Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown. Results The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14–29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D. Conclusion Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

Description: Background There is a long-term interest in investigating the genetic basis of the horned/polled phenotype in domestic goats. Here, we report a genome-wide association study (GWAS) to detect the genetic loci affecting the polled phenotype in goats. Results We obtained a total of 13,980,209 biallelic SNPs, using the genotyping-by-sequencing data from 45 Jintang Black (JT) goats, which included 32 female and nine male goats, and four individuals with the polled intersex syndrome (PIS). Using a mixed-model based GWAS, we identified two association signals, which were located at 150,334,857–150,817,260 bp (P = 5.15 × 10− 119) and 128,286,704–131,306,537 bp (P = 2.74 × 10− 15) on chromosome 1. The genotype distributions of the 14 most significantly associated SNPs were completely correlated with horn status in goats, based on the whole-genome sequencing (WGS) data from JT and two other Chinese horned breeds. However, variant annotation suggested that none of the detected SNPs within the associated regions were plausible causal mutations. Via additional read-depth analyses and visual inspections of WGS data, we found a 10.1-kb deletion (CHI1:g. 129424781_129434939del) and a 480-kb duplication (CHI1:150,334,286–150,818,098 bp) encompassing two genes KCNJ15 and ERG in the associated regions of polled and PIS-affected goats. Notably, the 10.1-kb deletion also served as the insertion site for the 480-kb duplication, as validated by PCR and Sanger sequencing. Our WGS genotyping showed that all horned goats were homozygous for the reference alleles without either the structural variants (SVs), whereas the PIS-affected goats were homozygous for both the SVs. We also demonstrated that horned, polled, and PIS-affected individuals among 333 goats from JT and three other Chinese horned breeds can be accurately classified via PCR amplification and agarose gel electrophoresis of two fragments in both SVs. Conclusion Our results revealed that two genomic regions on chromosome 1 are major loci affecting the polled phenotypes in goats. We provided a diagnostic PCR to accurately classify horned, polled, and PIS-affected goats, which will enable a reliable genetic test for the early-in-life prediction of horn status in goats.

Description: Objective Colon cancer (CC) is one of the most common cancers whose progression is regulated by a number of factors, including circular RNAs (circRNAs). Nonetheless, circ_0038718 is a novel circRNA, and its regulatory mechanism in CC remains unclear. Methods Real-time quantitative PCR (qRT-PCR) was performed to detect the expression of circ_0038718, miR-195-5p and Axin2. Western blot was conducted to determine the protein expression of Axin2 and the key proteins on Wnt/β-catenin signaling pathway. Oligo (dT) 18 primers and RNase R were employed to identify the circular features of circ_0038718, and the location of circ_0038718 in cells was detected via nucleocytoplasmic separation. Dual-luciferase reporter assay and RNA binding protein immunoprecipitation experiment were carried out to investigate the molecular mechanism of circ_0038718/miR-195-5p/Axin2. Additionally, MTT assay was conducted to assess cell proliferation; Transwell assay was performed to evaluate cell migration and invasion, respectively. The effect of circ_0038718 on CC tumor growth was tested through tumor formation in nude mice. Results circ_0038718 was highly expressed in CC and could sponge miR-195-5p in cytoplasm. Silencing circ_0038718 suppressed the proliferative, migratory and invasive abilities of CC cells, while the promoting effect of high circ_0038718 expression on CC cells was reversed upon miR-195-5p over-expression. Axin2 was a downstream target of miR-195-5p and could regulate the Wnt/β-catenin signaling pathway. Axin2 expression was modulated by circ_0038718/miR-195-5p. Knockdown of Axin2 could also attenuate the promoting effect of high circ_0038718 expression on CC cell malignant progression, thus inhibiting tumor growth. Conclusion circ_0038718 is able to facilitate CC cell malignant progression via the miR-195-5p/Axin2 axis, which will provide a new idea for finding a novel targeted treatment of CC.

Description: Background SQUAMOSA promoter binding proteins (SBPs) genes encode a family of plant-specific transcription factors involved in various growth and development processes, including flower and fruit development, leaf initiation, phase transition, and embryonic development. The SBP gene family has been identified and characterized in many species, but no systematic analysis of the SBP gene family has been carried out in sugarcane. Results In the present study, a total of 50 sequences for 30 SBP genes were identified by the genome-wide analysis and designated SsSBP1 to SsSBP30 based on their chromosomal distribution. According to the phylogenetic tree, gene structure and motif features, the SsSBP genes were classified into eight groups (I to VIII). By synteny analysis, 27 homologous gene pairs existed in SsSBP genes, and 37 orthologous gene pairs between sugarcane and sorghum were found. Expression analysis in different tissues, including vegetative and reproductive organs, showed differential expression patterns of SsSBP genes, indicating their functional diversity in the various developmental processes. Additionally, 22 SsSBP genes were predicted as the potential targets of miR156. The differential expression pattern of miR156 exhibited a negative correlation of transcription levels between miR156 and the SsSBP gene in different tissues. Conclusions The sugarcane genome possesses 30 SsSBP genes, and they shared similar gene structures and motif features in their subfamily. Based on the transcriptional and qRT-PCR analysis, most SsSBP genes were found to regulate the leaf initial and female reproductive development. The present study comprehensively and systematically analyzed SBP genes in sugarcane and provided a foundation for further studies on the functional characteristics of SsSBP genes during different development processes.

Description: Background PIWI-interacting RNAs (piRNAs) are an abundant single-stranded type of small non-coding RNAs (sncRNAs), which initially were discovered in gonadal cells, with a role as defenders of genomic integrity in the germline, acting against the transposable elements. With a regular size range of 21-35 nt, piRNAs are associated with a PIWI-clade of Argonaute family proteins. The most widely accepted mechanisms of biogenesis for piRNAs involve the transcription of longer precursors of RNAs to be processed, by complexes of proteins, to functional size, preferentially accommodating uridine residues at the 5’ end and 3’ methylation to increase the stability of these molecules. piRNAs have also been detected in somatic cells, with diverse potential functions, indicating their high plasticity and pleiotropic activity. Discovered about two decades ago, piRNAs are a large and versatile type of sncRNA and that remain insufficiently identified and analyzed, through next-generation sequencing (NGS), to evaluate their processing, functions, and biogenesis in different cell types and during development. piRNAs’ distinction from other sncRNAs has led to controversial results and interpretation difficulties when using existing databases because of the heterogeneity of the criteria used in making the distinction. Description We present “piRNA-IPdb”, a database based uniquely on datasets obtaining after the defining characteristic of piRNAs: those small RNAs bound to PIWI proteins. We selected and analyzed sequences from piRBase that exclusively cover the binding to PIWI. We pooled a total of 18,821,815 sequences from RNA-seq after immunoprecipitation that included the binding to any of the mouse PIWI proteins (MILI, MIWI, or MIWI2). Conclusions In summary, we present the characteristics and potential use of piRNA-IPdb database for the analysis of bona fide piRNAs.

Description: Background Alternative splicing (AS) is an important mechanism of posttranscriptional modification and dynamically regulates multiple physiological processes in plants, including fruit ripening. However, little is known about alternative splicing during fruit development in fleshy fruits. Results We studied the alternative splicing at the immature and ripe stages during fruit development in cucumber, melon, papaya and peach. We found that 14.96–17.48% of multiexon genes exhibited alternative splicing. Intron retention was not always the most frequent event, indicating that the alternative splicing pattern during different developmental process differs. Alternative splicing was significantly more prevalent at the ripe stage than at the immature stage in cucumber and melon, while the opposite trend was shown in papaya and peach, implying that developmental stages adopt different alternative splicing strategies for their specific functions. Some genes involved in fruit ripening underwent stage-specific alternative splicing, indicating that alternative splicing regulates fruits ripening. Conserved alternative splicing events did not appear to be stage-specific. Clustering fruit developmental stages across the four species based on alternative splicing profiles resulted in species-specific clustering, suggesting that diversification of alternative splicing contributes to lineage-specific evolution in fleshy fruits. Conclusions We obtained high quality transcriptomes and alternative splicing events during fruit development across the four species. Dynamics and nonconserved alternative splicing were discovered. The candidate stage-specific AS genes involved in fruit ripening will provide valuable insight into the roles of alternative splicing during the developmental processes of fleshy fruits.

Description: Background The phytopatogen Claviceps paspali is the causal agent of Ergot disease in Paspalum spp., which includes highly productive forage grasses such as P. dilatatum. This disease impacts dairy and beef production by affecting seed quality and producing mycotoxins that can affect performance in feeding animals. The molecular basis of pathogenicity of C. paspali remains unknown, which makes it more difficult to find solutions for this problem. Secreted proteins are related to fungi virulence and can manipulate plant immunity acting on different subcellular localizations. Therefore, identifying and characterizing secreted proteins in phytopathogenic fungi will provide a better understanding of how they overcome host defense and cause disease. The aim of this work is to analyze the whole genome sequences of three C. paspali isolates to obtain a comparative genome characterization based on possible secreted proteins and pathogenicity factors present in their genome. In planta RNA-seq analysis at an early stage of the interaction of C. paspali with P. dilatatum stigmas was also conducted in order to determine possible secreted proteins expressed in the infection process. Results C. paspali isolates had compact genomes and secretome which accounted for 4.6–4.9% of the predicted proteomes. More than 50% of the predicted secretome had no homology to known proteins. RNA-Seq revealed that three protein-coding genes predicted as secreted have mayor expression changes during 1 dpi vs 4 dpi. Also, three of the first 10 highly expressed genes in both time points were predicted as effector-like. CAZyme-like proteins were found in the predicted secretome and the most abundant family could be associated to pectine degradation. Based on this, pectine could be a main component affected by the cell wall degrading enzymes of C. paspali. Conclusions Based on predictions from DNA sequence and RNA-seq, unique probable secreted proteins and probable pathogenicity factors were identified in C. paspali isolates. This information opens new avenues in the study of the biology of this fungus and how it modulates the interaction with its host. Knowledge of the diversity of the secretome and putative pathogenicity genes should facilitate future research in disease management of Claviceps spp.

Description: Background miRNAs regulate circadian patterns by modulating the biological clocks of animals. In our previous study, we found that the clock gene exhibited a cosine expression pattern in the fallopian tube of chicken uterus. Clock-controlled miRNAs are present in mammals and Drosophila; however, whether there are clock-controlled miRNAs in the chicken uterus and, if so, how they regulate egg-laying rhythms is unclear. In this study, we selected 18 layer hens with similar ovipositional rhythmicity (each of three birds were sacrificed for study per 4 h throughout 24 h); their transcriptomes were scanned to identify the circadian miRNAs and to explore regulatory mechanisms within the uterus of chickens. Results We identified six circadian miRNAs that are mainly associated with several biological processes including ion trans-membrane transportation, response to calcium ion, and enrichment of calcium signaling pathways. Verification of the experimental results revealed that miR-449c-5p exhibited a cosine expression pattern in the chicken uterus. Ca2+-transporting ATPase 4 (ATP2B4) in the plasma membrane is the predicted target gene of circadian miR-449c-5p and is highly enriched in the calcium signaling pathway. We speculated that clock-controlled miR-449c-5p regulated Ca2+ transportation during eggshell calcification in the chicken uterus by targeting ATP2B4. ATP2B4 mRNA and protein were rhythmically expressed in the chicken uterus, and dual-luciferase reporter gene assays confirmed that ATP2B4 was directly targeted by miR-449c-5p. The expression of miR-449c-5p showed an opposite trend to that of ATP2B4 within a 24 h cycle in the chicken uterus; it inhibited mRNA and protein expression of ATP2B4 in the uterine tubular gland cells. In addition, overexpression of ATP2B4 significantly decreased intracellular Ca2+ concentration (P 

Description: Background Mammals have wide variations in testicular position, with scrotal testes in some species and ascrotal testes in others. Although cryptorchidism is hazardous to human health, some mammalian taxa are natural cryptorchids. However, the evolution of testicular position and the molecular mechanisms underlying the maintenance of health, including reproductive health, in ascrotal mammals are not clear. Results In the present study, comparative genomics and evolutionary analyses revealed that genes associated with the extracellular matrix and muscle, contributing to the development of the gubernaculum, were involved in the evolution of testicular position in mammals. Moreover, genes related to testicular position were significantly associated with spermatogenesis and sperm fertility. These genes showed rapid evolution and the signature of positive selection, with specific substitutions in ascrotal mammals. Genes associated with testicular position were significantly enriched in functions and pathways related to cancer, DNA repair, DNA replication, and autophagy. Conclusions Our results revealed that alterations in gubernaculum development contributed to the evolution of testicular position in mammals and provided the first support for two hypotheses for variation in testicular position in mammals, the “cooling hypothesis”, which proposes that the scrotum provides a cool environment for acutely heat-sensitive sperm and the “training hypothesis”, which proposes that the scrotum develops the sperm by exposing them to an exterior environment. Further, we identified cancer resistance and DNA repair as potential protective mechanisms in natural cryptorchids. These findings provide general insights into cryptorchidism and have implications for health and infertility both in humans and domestic mammals.