Selection of housekeeping genes as internal controls for quantitative RT-PCR analysis of the veined rapa whelk (Rapana venosa)

The MYST Family Histone Acetyltransferase SasC Governs Diverse Biological Processes in Aspergillus fumigatus

The conserved MYST proteins form the largest family of histone acetyltransferases (HATs) that acetylate lysines within the N-terminal tails of histone, enabling active gene transcription. Here, we have investigated the biological and regulatory functions of the MYST family HAT SasC in the opportunistic human pathogenic fungus Aspergillus fumigatus using a series of genetic, biochemical, pathogenic, and transcriptomic analyses. The deletion (Δ) of sasC results in a drastically reduced colony growth, asexual development, spore germination, response to stresses, and the fungal virulence. Genome-wide expression analyses have revealed that the ΔsasC mutant showed 2402 significant differentially expressed genes: 1147 upregulated and 1255 downregulated. The representative upregulated gene resulting from ΔsasC is hacA, predicted to encode a bZIP transcription factor, whereas the UV-endonuclease UVE-1 was significantly downregulated by ΔsasC. Furthermore, our Western blot analyses suggest that SasC likely catalyzes the acetylation of H3K9, K3K14, and H3K29 in A. fumigatus. In conclusion, SasC is associated with diverse biological processes and can be a potential target for controlling pathogenic fungi.

Functional Characterization of the GNAT Family Histone Acetyltransferase Elp3 and GcnE in Aspergillus fumigatus

Post-translational modifications of chromatin structure by histone acetyltransferase (HATs) play a pivotal role in the regulation of gene expression and diverse biological processes. However, the function of GNAT family HATs, especially Elp3, in the opportunistic human pathogenic fungus Aspergillus fumigatus is largely unknown. To investigate the roles of the GNAT family HATs Elp3 and GcnE in the A. fumigatus, we have generated and characterized individual null Δelp3 and ΔgcnE mutants. The radial growth of fungal colonies was significantly decreased by the loss of elp3 or gcnE, and the number of asexual spores (conidia) in the ΔgcnE mutant was significantly reduced. Moreover, the mRNA levels of the key asexual development regulators were also significantly low in the ΔgcnE mutant compared to wild type (WT). Whereas both the Δelp3 and ΔgcnE mutants were markedly impaired in the formation of adherent biofilms, the ΔgcnE mutant showed a complete loss of surface structure and of intercellular matrix. The ΔgcnE mutant responded differently to oxidative stressors and showed significant susceptibility to triazole antifungal agents. Furthermore, Elp3 and GcnE function oppositely in the production of secondary metabolites, and the ΔgcnE mutant showed attenuated virulence. In conclusion, Elp3 and GcnE are associated with diverse biological processes and can be potential targets for controlling the pathogenic fungus.

Integrated mRNA and miRNA transcriptomic analysis reveals the response of Rapana venosa to the metamorphic inducer (juvenile oysters)

Metamorphosis, as a critical developmental event, controls the population dynamics of most marine invertebrates, especially some carnivorous gastropods that feed on bivalves, whose population dynamics not only affect the maintenance of the ecological balance but also impact the protection of bivalve resources; therefore, the metamorphosis of carnivorous gastropods deserve attention. Here, we investigated the mechanism underlying the response of the carnivorous gastropod Rapana venosa to its metamorphic inducer juvenile oysters through integrated analysis of miRNA and mRNA profiles. According to the results, we speculated that the AMPK signaling pathway may be the critical regulator in the response to juvenile oysters in R. venosa competent larvae. The NF-kB and JAK-STAT signaling pathways that regulated apoptosis were also activated by the metamorphic inducer, which may result in the degeneration of the velum. Additionally, the significant changes in the expression of the SARP-19 precursor gene and protein cibby homolog 1-like gene may indicate that these signaling pathways also regulate growth and development during metamorphosis. This study provides further evidence that juvenile oysters can induce metamorphosis of R. venosa at the transcriptional level, which expands our understanding of the metamorphosis mechanism in carnivorous gastropods.

Systems Drug Discovery for Diffuse Large B Cell Lymphoma Based on Pathogenic Molecular Mechanism via Big Data Mining and Deep Learning Method

Diffuse large B cell lymphoma (DLBCL) is an aggressive heterogeneous disease. The most common subtypes of DLBCL include germinal center b-cell (GCB) type and activated b-cell (ABC) type. To learn more about the pathogenesis of two DLBCL subtypes (i.e., DLBCL ABC and DLBCL GCB), we firstly construct a candidate genome-wide genetic and epigenetic network (GWGEN) by big database mining. With the help of two DLBCL subtypes’ genome-wide microarray data, we identify their real GWGENs via system identification and model order selection approaches. Afterword, the core GWGENs of two DLBCL subtypes could be extracted from real GWGENs by principal network projection (PNP) method. By comparing core signaling pathways and investigating pathogenic mechanisms, we are able to identify pathogenic biomarkers as drug targets for DLBCL ABC and DLBCL GCD, respectively. Furthermore, we do drug discovery considering drug-target interaction ability, drug regulation ability, and drug toxicity. Among them, a deep neural network (DNN)-based drug-target interaction (DTI) model is trained in advance to predict potential drug candidates holding higher probability to interact with identified biomarkers. Consequently, two drug combinations are proposed to alleviate DLBCL ABC and DLBCL GCB, respectively.

The Lysine Demethylases KdmA and KdmB Differently Regulate Asexual Development, Stress Response, and Virulence in Aspergillus fumigatus

Histone demethylases govern diverse cellular processes, including growth, development, and secondary metabolism. In the present study, we investigated the functions of two lysine demethylases, KdmA and KdmB, in the opportunistic human pathogenic fungus Aspergillus fumigatus. Experiments with mutants harboring deletions of genes encoding KdmA (ΔkdmA) and KdmB (ΔkdmB) showed that KdmA is necessary for normal growth and proper conidiation, whereas KdmB negatively regulates vegetative growth and conidiation. In both mutant strains, tolerance to H2O2 was significantly decreased, and the activities of both conidia-specific catalase (CatA) and mycelia-specific catalase (Cat1) were decreased. Both mutants had significantly increased sensitivity to the guanine nucleotide synthesis inhibitor 6-azauracil (6AU). The ΔkdmA mutant produced more gliotoxin (GT), but the virulence was not changed significantly in immunocompromised mice. In contrast, the production of GT and virulence were markedly reduced by the loss of kdmB. Comparative transcriptomic analyses revealed that the expression levels of developmental process-related genes and antioxidant activity-related genes were downregulated in both mutants. Taken together, we concluded that KdmA and KdmB have opposite roles in vegetative growth, asexual sporulation, and GT production. However, the two proteins were equally important for the development of resistance to 6AU.

Symbiotic microbiome and metabolism profiles reveal the effects of induction by oysters on the metamorphosis of the carnivorous gastropod Rapana venosa

Most marine mollusks have a pelagic larval phase, and they need to undergo metamorphosis to develop into adults. Metamorphosis is affected by many factors, including abiotic factors such as temperature, salinity and illumination as well as biological factors such as food and microorganisms. In our previous study, we found that the metamorphosis of Rapana venosa requires induction by juvenile oysters, which are the food source of R. venosa. However, the regulatory mechanism of this induction is largely unknown. In the present study, we evaluated the impacts of induction by juvenile oysters on competent larvae of R. venosa. Competent larvae were experimentally divided into two pools, and scallop shells without juvenile oysters and scallop shells with juvenile oysters were added for 2 h and 12 h to monitor alterations in critical gene expression, symbiotic microbiota and metabolomic responses. The carboxypeptidase gene was increased while the cellulase gene was decreased, which may mean that the food habit transition was induced by juvenile oysters. Meanwhile, critical genes in the neuroendocrine system were also significantly altered in juvenile oysters. Furthermore, dramatic changes in the symbiotic microbiota and metabolism profiles were observed, with many of them associated with the digestive system and neuroendocrine system. In conclusion, juveniles as food resources may induce metamorphosis in R. venosa by regulating the neuroendocrine system and promoting the development of the digestive system and changes in digestive enzymes. This study may provide evidence that induction by juvenile oysters can promote food habit transition and metamorphosis in R. venosa by regulating the metabolome and microbiome and further altering the digestive and neuroendocrine systems of R. venosa, which expands our understanding of the regulatory mechanism of metamorphosis in R. venosa. However, further studies are needed to explore the specific substance inducing metamorphosis released by juvenile oysters.

Screening of reference genes in tiger puffer (Takifugu rubripes) across tissues and under different nutritional conditions

The present study was aimed at screening suitable reference genes for quantitative real-time polymerase chain reaction (qRT-PCR) in tiger puffer (Takifugu rubripes), an important aquaculture species in Asia and also a good model species for lipid research. Specifically, this reference gene screening was targeted at standardization of gene expression in different tissues (liver, muscle, brain, intestine, heart, eye, skin, and spleen) or under different nutritional conditions (starvation and different dietary lipid levels). Eight candidate reference genes (ribosomal protein L19 and L13 (RPL19 and RPL13), elongation factor-1 alpha (EF1α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine guanine phosphoribosyl transferase1 (HPRT1), beta-2-Microglobulin (B2M), 18S ribosomal RNA (18SrRNA), and beta actin (ACTB)) were evaluated with four algorithms (geNorm, NormFinder, BestKeeper, and comparative ΔCt method). The results showed that different algorithms generated inconsistent results. Based on these findings, RPL19, EF1α, 18SrRNA, and RPL13 were relatively stable in different tissues of tiger puffer. During starvation conditions, ACTB/RPL19 was the best reference gene combination. Under different dietary lipid levels, ACTB/RPL13 was the most suitable reference gene combination. The present results will help researchers to obtain more accurate results in future qRT-PCR analysis in tiger puffer.

Evaluation of candidate reference genes for quantitative RTqPCR analysis in goldfish (Carassius auratus L.) in healthy and CyHV-2 infected fish

The accuracy of quantitative real time PCR (RTqPCR) can be attained only when a suitable reference gene is used. The gene expression for a particular gene may vary within different cells at different conditions. Hence, the suitability and stability of various potential reference genes have to be determined for expression studies. In this study, we have examined the potential of four different reference genes including β-Actin (ACTB), 18S ribosomal RNA (18S), glyceraldehyde-3P-dehydrogenase (GAPDH), and elongation factor 1 alpha (EF1A_A) in seven different tissues including gill, liver, kidney, spleen, heart, muscle and intestine of goldfish (Carassius auratus). The housekeeping genes were analyzed from healthy fish and in CyHV-2 challenged fish. Based upon the real time PCR results the gene expression varied among the genes and in tissues. The expression levels of the housekeeping genes were then compared and evaluated with the RefFinder web tool which analyses results using four different algorithms - BestKeeper, delta Ct, geNorm and NormFinder. EF1A_A was ranked to be the best gene in healthy fish by BestKeeper and geNorm analysis. The delta Ct and NormFinder algorithm have found 18S to be a stable gene in healthy fish but 18S was given to be least expressed in challenged fish. ACTB was also given as a stable gene by geNorm analysis in both healthy and challenged fish. Also, in CyHV-2 challenged fish, EF1A_A was identified as the best gene by all the three analysis except by BestKeeper analysis, where it has ranked GADPH as the best housekeeping gene. Expression of the four candidate reference genes differed across all tissue types tested, inferring that a thorough study of the reference genes is necessary for cross tissue comparison. These results can be further used in the immune gene response study of goldfish infected with any viral pathogen to develop better health strategies in the disease management of goldfish aquaculture.

Characterization of the mbsA Gene Encoding a Putative APSES Transcription Factor in Aspergillus fumigatus

The APSES family proteins are transcription factors (TFs) with a basic helix-loop-helix domain, known to regulate growth, development, secondary metabolism, and other biological processes in Aspergillus species. In the genome of the human opportunistic pathogenic fungus Aspergillus fumigatus, five genes predicted to encode APSES TFs are present. Here, we report the characterization of one of these genes, called mbsA (Afu7g05620). The deletion (Δ) of mbsA resulted in significantly decreased hyphal growth and asexual sporulation (conidiation), and lowered mRNA levels of the key conidiation genes abaA, brlA, and wetA. Moreover, ΔmbsA resulted in reduced spore germination rates, elevated sensitivity toward Nikkomycin Z, and significantly lowered transcripts levels of genes associated with chitin synthesis. The mbsA deletion also resulted in significantly reduced levels of proteins and transcripts of genes associated with the SakA MAP kinase pathway. Importantly, the cell wall hydrophobicity and architecture of the ΔmbsA asexual spores (conidia) were altered, notably lacking the rodlet layer on the surface of the ΔmbsA conidium. Comparative transcriptomic analyses revealed that the ΔmbsA mutant showed higher mRNA levels of gliotoxin (GT) biosynthetic genes, which was corroborated by elevated levels of GT production in the mutant. While the ΔmbsA mutant produced higher amount of GT, ΔmbsA strains showed reduced virulence in the murine model, likely due to the defective spore integrity. In summary, the putative APSES TF MbsA plays a multiple role in governing growth, development, spore wall architecture, GT production, and virulence, which may be associated with the attenuated SakA signaling pathway.

Changes in Symbiotic Microbiota and Immune Responses in Early Development Stages of Rapana venosa (Valenciennes, 1846) Provide Insights Into Immune System Development in Gastropods

The symbiotic microbiota can stimulate modulation of immune system, which also can promote immune system mature in critical developmental periods. In this study, we have investigated the symbiotic microbiota in Rapana venosa at five early development stages using Illumina high-throughput sequencing, and detected immune responses in larvae. Analysis of the symbiotic microbiota sequences identified that the most abundant phylum was Proteobacteria. Beta diversity analysis indicated that the structure of the symbiotic microbiota dramatically shifted in early development stages. The abundance of immune-related KEGG Orthologs (KOs) also increased in competent larval (J4, 30-day post-hatching) and postlarval after 3 days of metamorphosis (Y5, 33-day post-hatching) stages. Acid phosphatase activity decreased significantly in the Y5 stage, and alkaline phosphatase activity also at a lower level in Y5 stage, whereas lysozyme activities exhibited no remarkable change. Also, the activities of catalase and superoxide dismutase activities decreased dramatically during early development stages of R. venosa. Dramatic changes in the symbiotic microbiota and the immune response mainly occurred in the initially hatched veliger (C1), competent larval (J4) and postlarval (Y5) stages, during which the hosts might experience substantial environmental changes or changes in physiological structure and function. These findings expand our understanding of the stage-specific symbiotic microbiota in R. venosa and the close association between immune system and symbiotic microbiota in mollusks, however, the specific relationship may need more researches are needed to investigated in the future.

Screening and evaluation of the stability of expression of reference genes in Lymantria dispar (Lepidoptera: Erebidae) using qRT-PCR

The quantitative real-time polymerase chain reaction (qRT-PCR) has rapidly become the most sensitive and accurate method for the quantitative analysis of gene expression. Normalization of gene expression to that of relatively stably expressed housekeeping genes is required to facilitate the study of gene expression and to obtain more accurate RT-PCR data. However, no studies of the stability of expression of housekeeping genes in Lymantria dispar have been reported. In the present study, BestKeeper, GeNorm and NormFinder statistical software was used to evaluate the expression of thirteen candidate reference genes in L. dispar under different conditions. The expression levels of candidate reference genes were determined for two biological factors (developmental stages and tissues) and four abiotic treatments (temperature, insecticide, CO₂ and starvation). The results showed that the best candidate reference genes in L. dispar were TUB, AK, RPS15 for developmental stages, RPL32 and GAPDH for tissues, ACTB and EF1-α for CO₂ stress, GAPDH and RPL32 for temperature stress, RPS3 and GAPDH for insecticide stress, and GAPDH and RPS3 for starvation stress. In summary, EF1-α and TUB are preferential housekeeping genes in L. dispar under various conditions. These results provide a basis for the further study of functional genes of L. dispar.

Identification of candidate reference genes for qRT-PCR normalization studies of salinity stress and injury in Onchidium reevesii

Real-time quantitative reverse transcription-PCR (qRT-PCR) is an undeniably effective tool for measuring levels of gene expression, but the accuracy and reliability of the statistical data obtained depend mainly on the basal expression of selected housekeeping genes in many samples. To date, there have been few analyses of stable housekeeping genes in Onchidium reevesii under salinity stress and injury. In this study, the gene expression stabilities of seven commonly used housekeeping genes, CYC, RPL28S, ACTB, TUBB, EF1a, Ubiq and 18S RNA, were investigated using BestKeeper, geNorm, NormFinder and RefFinfer. Although the results of the four programs varied to some extent, in general, RPL28S, TUBB, ACTB and EF1a were ranked highly. ACTB and TUBB were found to be the most stable housekeeping genes under salinity stress, and EF1a plus TUBB was the most stable combination under injury stress. When analysing target gene expression in different tissues, RPL28S or EF1a should be selected as the reference gene according to the level of target gene expression. Under extreme environmental stress (salinity) conditions, ACTB (0 ppt, 5 ppt, 15 ppt, 25 ppt) and TUBB (35 ppt) are reasonable reference gene choices when expression stability and abundance are considered. Under conditions of 15 ppt salinity and injury stress, our results showed that the best two-gene combination was TUBB plus EF1a. Therefore, we suggest that RPL28S, ACTB and TUBB are suitable reference genes for evaluating mRNA transcript levels. Based on candidate gene expression analysis, the tolerance of O. reevesii to low salinity (low osmotic pressure) is reduced compared to its tolerance to high salinity (high osmotic pressure). These findings will help researchers obtain accurate results in future quantitative gene expression analyses of O. reevesii under other stress conditions.

Transcriptome-based identification of the optimal reference genes as internal controls for quantitative RT-PCR in razor clam (Sinonovacula constricta)

Quantitative real-time PCR (qRT-PCR) is a standard method to measure gene expression in function exploring. Accurate and reproducible data of qRT-PCR requires appropriate reference genes, which are stably expressed under different experimental conditions. However, no housekeeping genes were validated as internal controls for qRT-PCR in Sinonovacula constricta. In this study, we classified the transcriptome data of two tissues for Vibrio infection and Cd²⁺ stress into ten clusters based on the gene expression patterns. Among them, cluster 5 had the most stable gene expression patterns regardless of tissues and treatments as the database for candidate reference genes. A total of 55 orthologs of classical housekeeping genes in the clam transcriptome were annotated. Combined the expression profiles and housekeeping genes in S. constricta, we chose eight candidate reference genes and validated their expression in Vibrio-infected samples and different tissues by qRT-PCR. Their expression stability was analyzed by three different algorithms geNorm, NormFinder and BestKeeper. Although the rank of the eight candidate reference genes is different in different treatments using different software, RS9 could be the best reference genes for normalization of qRT-PCR expression data in S. constricta under various treatments considering the above analysis. Meanwhile, the ranking of genes based on the CV values of transcriptomic data was similar to the validation results. This study provides for the first time a list of suitable reference genes for S. constricta and a valuable resource for further studies of clam immune defense systems.

Understanding microRNA Regulation Involved in the Metamorphosis of the Veined Rapa Whelk (Rapana venosa)

The veined rapa whelk (Rapana venosa) is widely consumed in China. Nevertheless, it preys on oceanic bivalves, thereby reducing this resource worldwide. Its larval metamorphosis comprises a transition from pelagic to benthic form, which involves considerable physiological and structural changes and has vital roles in its natural populations and commercial breeding. Thus, understanding the endogenous microRNAs (miRNAs) that drive metamorphosis is of great interest. This is the first study to use high-throughput sequencing to examine the alterations in miRNA expression that occur during metamorphosis in a marine gastropod. A total of 195 differentially expressed miRNAs were obtained. Sixty-five of these were expressed during the transition from precompetent to competent larvae. Thirty-three of these were upregulated and the others were downregulated. Another 123 miRNAs were expressed during the transition from competent to postlarvae. Ninety-six of these were upregulated and the remaining 27 were downregulated. The expression of miR-276-y, miR-100-x, miR-183-x, and miR-263-x showed a >100-fold change during development, while the miR-242-x and novel-m0052-3p expression levels changed over 3000-fold. Putative target gene coexpression, gene ontology, and pathway analyses suggest that these miRNAs play important parts in cell proliferation, migration, apoptosis, metabolic regulation, and energy absorption. Twenty miRNAs and their target genes involved in ingestion, digestion, cytoskeleton, cell adhesion, and apoptosis were identified. Nine of them were analyzed with real-time polymerase chain reaction (PCR), which showed an inverse correlation between the miRNAs and their relative expression levels. Our data elucidate the role of miRNAs in R. venosa metamorphic transition and serve as a solid basis for further investigations into regulatory mechanisms of gastropod metamorphosis.