Porcine epidemic diarrhea virus strain CH/HLJ/18 isolated in China: characterization and phylogenetic analysis

BACKGROUND

Porcine epidemic diarrhea (PED) is an infectious disease of the digestive tract caused by the porcine epidemic diarrhea virus (PEDV), characterized by vomiting, severe diarrhea, and high mortality rates in piglets. In recent years, the distribution of this disease in China has remarkably increased, and its pathogenicity has also increased. PEDV has been identified as the main cause of viral diarrhea in piglets. This study aimed to understand the genetic evolution and diversity of PEDV to provide a theoretical basis for the development of new vaccines and the prevention and treatment of PED.

METHODS

A PEDV strain was isolated from the small intestine of a diarrheal piglet using Vero cells. The virus was identified using reverse transcription-polymerase chain reaction (RT-PCR), indirect immunofluorescence assay (IFA), and transmission electron microscopy. The whole genome sequence was sequenced, phylogenetic analysis was conducted using MEGA (version 7.0), and recombination analysis was performed using RDP4 and SimPlot. The S protein amino acid sequence was aligned using Cluster X (version 2.0), and the S protein was modeled using SWISS-MODEL to compare differences in structure and antigenicity. Finally, the piglets were inoculated with PEDV to evaluate its pathogenicity in newborn piglets.

RESULT

PEDV strain CH/HLJ/18 was isolated. CH/HLJ/18 shared 89.4-99.2% homology with 52 reference strains of PEDV belonging to the GII-a subgroup. It was a recombinant strain of PEDV BJ-2011-1 and PEDV CH_hubei_2016 with a breakpoint located in ORF1b. Unique amino acid deletions and mutations were observed in the CH/HLJ/18 S protein. The piglets then developed severe watery diarrhea and died within 7 d of inoculation with CH/HLJ/18, suggesting that CH/HLJ/18 was highly pathogenic to newborn piglets.

CONCLUSION

A highly pathogenic recombinant PEDV GII-a strain, CH/HLJ/18, was identified in China, with unique deletion and mutation of amino acids in the S protein that may lead to changes in protein structure and antigenicity. These results will be crucial for understanding the prevalence and variation of PEDV and for preventing and controlling PED.

Assembly and analysis of the complete mitochondrial genome of Forsythia suspensa (Thunb.) Vahl

BACKGROUND

Forsythia suspensa (Thunb.) Vahl is a valuable ornamental and medicinal plant. Although the nuclear and chloroplast genomes of F. suspensa have been published, its complete mitochondrial genome sequence has yet to be reported. In this study, the genomic DNA of F. suspensa yellowish leaf material was extracted, sequenced by using a mixture of Illumina Novaseq6000 short reads and Oxford Nanopore PromethION long reads, and the sequencing data were assembled and annotated.

RESULT

The F. suspensa mitochondrial genome was obtained in the length of 535,692 bp with a circular structure, and the GC content was 44.90%. The genome contains 60 genes, including 36 protein-coding genes, 21 tRNA genes, and three rRNA genes. We further analyzed RNA editing of the protein-coding genes, relative synonymous codon usage, and sequence repeats based on the genomic data. There were 25 homologous sequences between F. suspensa mitochondria and chloroplast genome, which involved the transfer of 8 mitochondrial genes, and 9473 homologous sequences between mitochondrial and nuclear genomes. Analysis of the nucleic acid substitution rate, nucleic acid diversity, and collinearity of protein-coding genes of the F. suspensa mitochondrial genome revealed that the majority of genes may have undergone purifying selection, exhibiting a slower rate of evolution and a relatively conserved structure. Analysis of the phylogenetic relationships among different species revealed that F. suspensa was most closely related to Olea europaea subsp. Europaea.

CONCLUSION

In this study, we sequenced, assembled, and annotated a high-quality F. suspensa mitochondrial genome. The results of this study will enrich the mitochondrial genome data of Forsythia, lay a foundation for the phylogenetic development of Forsythia, and promote the evolutionary analysis of Oleaceae species.

Genome-wide identification of BBX gene family and their expression patterns under salt stress in soybean

BACKGROUND

BBX genes are key players in the regulation of various developmental processes and stress responses, which have been identified and functionally characterized in many plant species. However, our understanding of BBX family was greatly limited in soybean.

RESULTS

In this study, 59 BBX genes were identified and characterized in soybean, which can be phylogenetically classified into 5 groups. GmBBXs showed diverse gene structures and motif compositions among the groups and similar within each group. Noticeably, synteny analysis suggested that segmental duplication contributed to the expansion of GmBBX family. Moreover, our RNA-Seq data indicated that 59 GmBBXs showed different transcript profiling under salt stress, and qRT-PCR analysis confirmed their expression patterns. Among them, 22 GmBBXs were transcriptionally altered with more than two-fold changes by salt stress, supporting that GmBBXs play important roles in soybean tolerance to salt stress. Additionally, Computational assay suggested that GmBBXs might potentially interact with GmGI3, GmTOE1b, GmCOP1, GmCHI and GmCRY, while eight types of transcription factors showed potentials to bind the promoter regions of GmBBX genes.

CONCLUSIONS

Fifty-nine BBX genes were identified and characterized in soybean, and their expression patterns under salt stress and computational assays suggested their functional roles in response to salt stress. These findings will contribute to future research in regard to functions and regulatory mechanisms of soybean BBX genes in response to salt stress.

Identification and responding to exogenous hormone of HB-KNOX family based on transcriptome data of Caucasian clover

Caucasian clover (Trifolium ambiguum M. Bieb.) is a strongly rhizomatous, low-crowned perennial leguminous and ground-covering grass. The species is resistant to cold, arid temperatures and grazing due to a well-developed underground rhizome system and a strong clonal reproduction capacity. KNOTTED1-LIKE HOMEOBOX (KNOX) genes are a family of plant-specific homeobox transcription factors with important roles in plant development. Preliminary transcriptome analysis enabled us to understand the gene expression in five different tissues, which helped us to screen the predetermined genes of the HB-KNOX family genes for the rhizome growth and development of Caucasian clover. The study identified 41 TaKNOX genes from the Caucasian clover transcriptome database. Gene length, MW and pl of TaKNOX family transcription factors varied, but the gene structure and motifs were relatively conserved in bioinformatics analysis. Phylogenetic analyses of Arabidopsis thaliana, soybean, Medicago truncatula and Caucasian clover were performed to study the evolutionary and functional relationships in various species. Prediction and verification of the subcellular localizations revealed the diverse subcellular localization of these 41 TaKNOX proteins. The expression profile of exogenous hormones showed that the TaKNOX gene showed multiple expression regulation patterns, and was involved in 6-BA, IAA and KT signaling pathways. Our results reveal the characteristics of the TaKNOX gene family, thus laying a foundation for further functional analysis of the KNOX family in Caucasian clover.

Genome wide analysis of the heavy-metal-associated (HMA) gene family in tomato and expression profiles under different stresses

The heavy-metal-associated (HMA) family plays a major role in the transportation of metals. Despite having the genome sequence of the tomato (Solanum lycopersicum), the HMA gene family has not been studied yet. In this study, we identified 48 HMA genes and categorized them into Cu/Ag P1B-ATPase and Zn/Co/Cd/Pb P1BATPase sub-families according to their phylogenic relationship with Arabidopsis and rice. The SlHMA genes were distributed throughout the 12 chromosomes. Analysis of gene structure, chromosomal position, and synteny, revealed that segmental duplications bestowed their evolution. The high numbers of stress-related cis-elements were found to be present in the putative promoter regions indicate the involvement of SlHMAs in stress modulation pathways. RNA-seq data revealed that SlHMAs had divergent expression in different tissues and developmental stages, where members of Cu/Ag P1B-ATPase subfamily were strongly expressed in the roots. RT-qPCR analysis of nine selected SlHMAs showed that most of the genes were up-regulated in response to heavy metals and moderately regulated in response to different abiotic stresses such as salt, drought, and cold.

A detailed comparative analysis of codon usage bias in Alongshan virus

Alongshan virus (ALSV) is an emerging tick-borne pathogen that infects humans, causing febrile disease. ALSV uses Ixodes Persulcatus ticks to infect humans with a wide range of signs, from asymptomatic to encephalitis-like syndrome. There is an increasing public health concern about the ALSV infection. To get insight into the impacts of viral relations with their hosts on viral ability, survival, and evasion from hosts immune systems remain unknown. The codon usage is a driving force in viral genome evolution; therefore, we enrolled 41 ALSV strains in codon usage analysis to elucidate the molecular evolutionary dynamics of ALSV. The results indicate that the overall codon usage among ALSV isolates is relatively similar and slightly biased. Base compositions for the cds were in order of G >A >C >U and in the third position of codons G3 >A3 >C3 >T3. The RSCU values revealed that the more frequently used codons were mostly GC ended. Different codon preferences in ALSV genes in relation to codon usage of H. sapiens and Ixodes Persulcatus genes were found. Neutrality plot was determined to reveal the superiority of natural selection over directional mutation pressure in causing CUB based on GC12 versus GC3 contents. The results of these studies suggest that the emergence of ALSV in China, Russia and Finland may also be reflected in ALSV codon usage. Altogether, the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of ALSV.

The characteristic of the synonymous codon usage and phylogenetic analysis of hepatitis B virus

Background

Hepatitis B virus (HBV) infection is a crucial medical issue worldwide. The dependence of HBV replication on host cell machineries and their co-evolutionary interactions prompt the codon usage pattern of viral genes to translation selection and mutation pressure.

Objective

The evolutionary characteristics of HBV and the natural selection effects of the human genome on the codon usage characteristics were analyzed to provide a basis for medication development for HBV infection.

Methods

The codon usage pattern of sequences from different HBV genotypes of our isolates and reference HBV genome sequences downloaded from the National Center for Biotechnology Information (NCBI) database were analyzed by computing the relative synonymous codon usage (RSCU), nucleotide content, codon adaptation index (CAI) and the effective number of codons (ENC).

Results

The highest ENC values were observed in the C genotypes, followed by the B genotypes. The ENC values indicated a weak codon usage bias (CUB) in HBV genome. The number of codons differentially used between the three genotypes was markedly higher than that of similarly used codons. High CAI values indicated a good adaptability of HBV to its host. The ENC plot indicated the occurrence of mutational pressure in the three genotypes. The mean Ka/Ks ratios in the three genotypes were lower than 1, which indicated a negative selection pressure. The CAI and GC3% plot indicated the existence of CUB in the HBV genome.

Conclusions

Nucleotide composition, mutation bias, negative selection and mutational pressure are key factors influencing the CUB and phylogenetic diversity in HBV genotypes. The data provided here could be useful for developing drugs for HBV infection.

Electronic supplementary material

The online version of this article (10.1007/s13258-020-00932-w) contains supplementary material, which is available to authorized users.

Identification of eukaryotic translation initiation factors and the temperature-dependent nature of Turnip mosaic virus epidemics in allopolyploid Brassica juncea

Eukaryotic translation initiation factors (eIFs) are essential protein complexes involved in the translation of mRNA into proteins. These initiation factors are generally used as targets in the control of plant RNA virus infections. In the present study, we identified a total 190 eIFs, clustered phylogenetically into 40 distinct subfamilies in the allopolyploid Brassica juncea. Extensive evolutionary duplications of the eIFs in B. juncea suggest their increased genetic diversity and wide adaptability. The induction of expressions in some of the eIFs after inoculation against Turnip mosaic virus (TuMV) provided candidate targets to be used in the control of viral infections. In addition, the expression profiles of eIFs under different temperatures suggested that the TuMV epidemic was temperature dependent. The eIFs expressions suggested that the systemic viral infections were more acute in plants grown between 20 °C and 28 °C. In addition, our results revealed that new subgroups of eIFs, eIF2β, eIF2α, eIF2Bβ, EF1A, and PABP could be represented as targets for antiviral strategies in B. juncea. In summary, our findings would be helpful in studying the complex mechanisms of eIF-mediated, temperature-dependent RNA virus control in B. juncea.

Genome-wide analysis of the MATE gene family in potato

The multidrug and toxic compound extrusion (MATE) protein family is a newly discovered family of secondary transporters that extrude metabolic waste and a variety of antibiotics out of the cell using an electrochemical gradient of H⁺ or Na⁺ across the membrane. The main function of MATE gene family is to participate in the process of plant detoxification and morphogenesis. The genome-wide analysis of the MATE genes in potato genome was conducted. At least 48 genes were initially identified and classified into six subfamilies. The chromosomal localization of MATE gene family showed that they could be distributed on 11 chromosomes except chromosome 9. The number of amino acids is 145-616, the molecular weight of proteins is 15.96-66.13 KD, the isoelectric point is 4.97-9.17, and they were located on the endoplasmic reticulum with having 4-13 transmembrane segments. They contain only two parts of the exons and UTR without introns. Some members of the first subfamily of potato MATE gene family are clustered with At2g04070 and they may be related to the transport of toxic compounds such as alkaloids and heavy metal. The function of the members of the second subfamily may be similar to that of At3g23560, which is related to tetramethylammonium transport. Some members of the third subfamily are clustered with At3g59030 and they may be involved in the transport of flavonoids. The fifth subfamily may be related to the transport of iron ions. The function of the sixth subfamily may be similar to that of At4g39030, which is related to salicylic acid transport. There are three kinds of conserved motifs in potato MATE genes, including the motif 1, motif 2, and motif 3. Each motif has 50 amino acids. The number of each motif is different in the gene sequence, of which 45 MATE genes contain at least a motif, but there is no motif in ST0015301, ST0045283, and ST0082336. These results provide a reference for further research on the function of potato MATE genes.

Comparative genomics to explore phylogenetic relationship, cryptic sexual potential and host specificity of Rhynchosporium species on grasses

BACKGROUND

The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts.

RESULTS

The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type.

CONCLUSION

Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.