Differential microRNA expression, microRNA arm switching, and microRNA:long noncoding RNA interaction in response to salinity stress in soybean

BACKGROUND

Soybean is a major legume crop with high nutritional and environmental values suitable for sustainable agriculture. Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are important regulators of gene functions in eukaryotes. However, the interactions between these two types of ncRNAs in the context of plant physiology, especially in response to salinity stress, are poorly understood.

RESULTS

Here, we challenged a cultivated soybean accession (C08) and a wild one (W05) with salt treatment and obtained their small RNA transcriptomes at six time points from both root and leaf tissues. In addition to thoroughly analyzing the differentially expressed miRNAs, we also documented the first case of miRNA arm-switching (miR166m), the swapping of dominant miRNA arm expression, in soybean in different tissues. Two arms of miR166m target different genes related to salinity stress (chloroplastic beta-amylase 1 targeted by miR166m-5p and calcium-dependent protein kinase 1 targeted by miR166m-3p), suggesting arm-switching of miR166m play roles in soybean in response to salinity stress. Furthermore, two pairs of miRNA:lncRNA interacting partners (miR166i-5p and lncRNA Gmax_MSTRG.35921.1; and miR394a-3p and lncRNA Gmax_MSTRG.18616.1) were also discovered in reaction to salinity stress.

CONCLUSIONS

This study demonstrates how ncRNA involves in salinity stress responses in soybean by miRNA arm switching and miRNA:lncRNA interactions. The behaviors of ncRNAs revealed in this study will shed new light on molecular regulatory mechanisms of stress responses in plants, and hence provide potential new strategies for crop improvement.

Horseshoe crab genomes reveal the evolution of genes and microRNAs after three rounds of whole genome duplication

Whole genome duplication (WGD) has occurred in relatively few sexually reproducing invertebrates. Consequently, the WGD that occurred in the common ancestor of horseshoe crabs ~135 million years ago provides a rare opportunity to decipher the evolutionary consequences of a duplicated invertebrate genome. Here, we present a high-quality genome assembly for the mangrove horseshoe crab Carcinoscorpius rotundicauda (1.7 Gb, N50 = 90.2 Mb, with 89.8% sequences anchored to 16 pseudomolecules, 2n = 32), and a resequenced genome of the tri-spine horseshoe crab Tachypleus tridentatus (1.7 Gb, N50 = 109.7 Mb). Analyses of gene families, microRNAs, and synteny show that horseshoe crabs have undergone three rounds (3R) of WGD. Comparison of C. rotundicauda and T. tridentatus genomes from populations from several geographic locations further elucidates the diverse fates of both coding and noncoding genes. Together, the present study represents a cornerstone for improving our understanding of invertebrate WGD events on the evolutionary fates of genes and microRNAs, at both the individual and population level. We also provide improved genomic resources for horseshoe crabs, of applied value for breeding programs and conservation of this fascinating and unusual invertebrate lineage.

Micro-RNA Clusters Integrate Evolutionary Constraints on Expression and Target Affinities: The miR-6/5/4/286/3/309 Cluster in Drosophila

A striking feature of micro-RNAs is that they are often clustered in the genomes of animals. The functional and evolutionary consequences of this clustering remain obscure. Here, we investigated a micro-RNA cluster miR-6/5/4/286/3/309 that is conserved across drosophilid lineages. Small RNA sequencing revealed expression of this micro-RNA cluster in Drosophila melanogaster leg discs, and conditional overexpression of the whole cluster resulted in leg appendage shortening. Transgenic overexpression lines expressing different combinations of micro-RNA cluster members were also constructed. Expression of individual micro-RNAs from the cluster resulted in a normal wild-type phenotype, but either the expression of several ancient micro-RNAs together (miR-5/4/286/3/309) or more recently evolved clustered micro-RNAs (miR-6-1/2/3) can recapitulate the phenotypes generated by the whole-cluster overexpression. Screening of transgenic fly lines revealed downregulation of leg-patterning gene cassettes in generation of the leg-shortening phenotype. Furthermore, cell transfection with different combinations of micro-RNA cluster members revealed a suite of downstream genes targeted by all cluster members, as well as complements of targets that are unique for distinct micro-RNAs. Considered together, the micro-RNA targets and the evolutionary ages of each micro-RNA in the cluster demonstrate the importance of micro-RNA clustering, where new members can reinforce and modify the selection forces on both the cluster regulation and the gene regulatory network of existing micro-RNAs. Key words: micro-RNA, cluster, evolution.

Millipede genomes reveal unique adaptations during myriapod evolution

The Myriapoda, composed of millipedes and centipedes, is a fascinating but poorly understood branch of life, including species with a highly unusual body plan and a range of unique adaptations to their environment. Here, we sequenced and assembled 2 chromosomal-level genomes of the millipedes Helicorthomorpha holstii (assembly size = 182 Mb; shortest scaffold/contig length needed to cover 50% of the genome [N50] = 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (assembly size = 449 Mb, N50 = 26.78 Mb mainly on 17 pseudomolecules). Unique genomic features, patterns of gene regulation, and defence systems in millipedes, not observed in other arthropods, are revealed. Both repeat content and intron size are major contributors to the observed differences in millipede genome size. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute (AGO) proteins for loading small RNAs are duplicated in both millipedes, but unlike in insects, an AGO duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs—including species-specific microRNA arm switching—providing differential gene regulation is also obtained. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules found in centipedes. We identify sets of genes associated with the ozadene that play roles in chemical defence as well as antimicrobial activity. Macro-synteny analyses revealed highly conserved genomic blocks between the 2 millipedes and deuterostomes. Collectively, our analyses of millipede genomes reveal that a series of unique adaptations have occurred in this major lineage of arthropod diversity. The 2 high-quality millipede genomes provided here shed new light on the conserved and lineage-specific features of millipedes and centipedes. These findings demonstrate the importance of the consideration of both centipede and millipede genomes—and in particular the reconstruction of the myriapod ancestral situation—for future research to improve understanding of arthropod evolution, and animal evolutionary genomics more widely.

Myriapods were among the first arthropods to invade the land over 400 million years ago, and survive today as the herbivorous millipedes and venomous centipedes. This study describes the genome sequences of two millipedes, Helicorthomorpha holstii and Trigoniulus corallinus, revealing unique adaptations not found in other arthropods.

Infection patterns of dengue, Zika and endosymbiont Wolbachia in the mosquito Aedes albopictus in Hong Kong

BACKGROUND

The mosquito Aedes albopictus is a vector of dengue and Zika viruses. Insecticide-resistant mosquito populations have evolved in recent decades, suggesting that new control strategies are needed. Hong Kong has a monsoon-influenced humid subtropical climate, which favours the spread of mosquitoes. However, baseline information on the composition and dynamics of the occurrence of endosymbiont Wolbachia in local Ae. albopictus is lacking, hindering the development of scientifically-informed control measures. This study identifies the presence and absence of dengue and Zika viruses, and Wolbachia infection in Aedes albopictus in Hong Kong.

METHODS

Oviposition traps were set at 57 areas in Hong Kong, and both immature and adult mosquitoes were collected on a monthly basis between April 2018 and April 2019 as the study sample. Each individual mosquito in this sample was processed and screened for the presence of the dengue and Zika viruses and the endosymbionts Wolbachia wAlbA and wAlbB with PCR.

RESULTS

Totals of 967 and 984 mosquitoes were tested respectively for the presence of dengue and Zika viruses, and no trace of either infection was found in these samples. The presence of wAlbA and wAlbB was also tested in 1582 individuals. Over 80% of these individuals were found to be stably infected with Wolbachia throughout the thirteen-month collection period (~ 47% singly-infected; ~ 36.8% doubly infected with both wAlbA and wAlbB).

CONCLUSIONS

The high degree of Wolbachia wAlbA and wAlbB infection in Ae. albopictus mosquitoes in Hong Kong, coupled with the absence of any signs of infection by dengue and Zika viruses, contrasts significantly with the pattern of mosquito infection in other parts of Asia. Further studies of the infection pattern in local mosquitoes are warranted before mosquito control strategies used in other regions are implemented in Hong Kong.

Chromosomal-level reference genome of the incense tree Aquilaria sinensis

Trees in the genus Aquilaria (Thymelaeaceae) are known as lign aloes, and are native to the forests of southeast Asia. Lign aloes produce agarwood as an antimicrobial defence. Agarwood has a long history of cultural and medicinal use, and is of considerable commercial value. However, due to habitat destruction and over collection, lign aloes are threatened in the wild. We present a chromosomal‐level assembly for Aquilaria sinensis, a lign aloe endemic to China known as the incense tree, based on Illumina short‐read, 10X Genomics linked‐read, and Hi‐C sequencing data. Our 783.8 Mbp A. sinensis genome assembly is of high physical contiguity, with a scaffold N50 of 87.6 Mbp, and high completeness, with a 95.8% BUSCO score for eudicotyledon genes. We include 17 transcriptomes from various plant tissues, providing a total of 35,965 gene models. We reveal the first complete set of genes involved in sesquiterpenoid production, plant defence, and agarwood production for the genus Aquilaria, including genes involved in the biosynthesis of sesquiterpenoids via the mevalonic acid (MVA), 1‐deoxy‐D‐xylulose‐5‐phosphate (DXP), and methylerythritol phosphate (MEP) pathways. We perform a detailed repeat content analysis, revealing that transposable elements account for ~61% of the genome, with major contributions from gypsy‐like and copia‐like LTR retroelements. We also provide a comparative analysis of repeat content across sequenced species in the order Malvales. Our study reveals the first chromosomal‐level genome assembly for a tree in the genus Aquilaria and provides an unprecedented opportunity to address a variety of applied, genomic and evolutionary questions in the Thymelaeaceae more widely.

miRNA-Mediated Interactions in and between Plants and Insects

Our understanding of microRNA (miRNA) regulation of gene expression and protein translation, as a critical area of cellular regulation, has blossomed in the last two decades. Recently, it has become apparent that in plant-insect interactions, both plants and insects use miRNAs to regulate their biological processes, as well as co-opting each others' miRNA systems. In this review article, we discuss the current paradigms of miRNA-mediated cellular regulation and provide examples of plant-insect interactions that utilize this regulation. Lastly, we discuss the potential biotechnological applications of utilizing miRNAs in agriculture.

High Diversity of Genogroup I Picobirnaviruses in Mammals

In a molecular epidemiology study using 791 fecal samples collected from different terrestrial and marine mammals in Hong Kong, genogroup I picobirnaviruses (PBVs) were positive by RT-PCR targeting the partial RdRp gene in specimens from five cattle, six monkeys, 17 horses, nine pigs, one rabbit, one dog, and 12 California sea lions, with 11, 9, 23, 17, 1, 1, and 15 sequence types in the positive specimens from the corresponding animals, respectively. Phylogenetic analysis showed that the PBV sequences from each kind of animal were widely distributed in the whole tree with high diversity, sharing 47.4–89.0% nucleotide identities with other genogroup I PBV strains based on the partial RdRp gene. Nine complete segment 1 (viral loads 1.7 × 10⁴ to 5.9 × 10⁶/ml) and 15 segment 2 (viral loads 4.1 × 10³ to 1.3 × 10⁶/ml) of otarine PBVs from fecal samples serially collected from California sea lions were sequenced. In the two phylogenetic trees constructed using ORF2 and ORF3 of segment 1, the nine segment 1 sequences were clustered into four distinct clades (C1–C4). In the tree constructed using RdRp gene of segment 2, the 15 segment 2 sequences were clustered into nine distinct clades (R1–R9). In four sea lions, PBVs were detected in two different years, with the same segment 1 clade (C3) present in two consecutive years from one sea lion and different clades present in different years from three sea lions. A high diversity of PBVs was observed in a variety of terrestrial and marine mammals. Multiple sequence types with significant differences, representing multiple strains of PBV, were present in the majority of PBV-positive samples from different kinds of animals.

Evaluation of 16SpathDB 2.0, an automated 16S rRNA gene sequence database, using 689 complete bacterial genomes

Interpretation of 16S rRNA sequences is a difficult problem faced by clinical microbiologists and technicians. In this study, we evaluated the updated 16SpathDB 2.0 database, using 689 16S rRNA sequences from 689 complete genomes of medically important bacteria. Among these 689 16S rRNA sequences, none was wrongly identified, with 35.8% reported as a single bacterial species having >98% identity with the query sequence (category 1), 63.9% reported as more than 1 bacterial species having >98% identity with the query sequence (category 2), 0.3% reported to the genus level (category 3), and none reported as no match (category 4). For the 16S rRNA sequences of non-duplicated bacterial species reported as category 1 or 2, the percentage of bacterial species reported as category 1 was significantly higher for anaerobic Gram-positive/Gram-negative bacteria than aerobic/facultative anaerobic Gram-positive/Gram-negative bacteria. 16SpathDB 2.0 is a user-friendly and accurate database for 16S rRNA sequence interpretation in clinical laboratories.

Identification of microRNA-like RNAs in mycelial and yeast phases of the thermal dimorphic fungus Penicillium marneffei

Background

Penicillium marneffei is the most important thermal dimorphic fungus causing systemic mycosis in China and Southeast Asia. While miRNAs are increasingly recognized for their roles in post-transcriptional regulation of gene expression in animals and plants, miRNAs in fungi were less well studied and their potential roles in fungal dimorphism were largely unknown. Based on P. marneffei genome sequence, we hypothesize that miRNA-like RNAs (milRNAs) may be expressed in the dimorphic fungus.

Methodology/Principal Findings

We attempted to identify milRNAs in P. marneffei in both mycelial and yeast phase using high-throughput sequencing technology. Small RNAs were more abundantly expressed in mycelial than yeast phase. Sequence analysis revealed 24 potential milRNA candidates, including 17 candidates in mycelial and seven in yeast phase. Two genes, dcl-1 and dcl-2, encoding putative Dicer-like proteins and the gene, qde-2, encoding Argonaute-like protein, were identified in P. marneffei. Phylogenetic analysis showed that dcl-2 of P. marneffei was more closely related to the homologues in other thermal dimorphic pathogenic fungi than to Penicillium chrysogenum and Aspergillus spp., suggesting the co-evolution of dcl-2 among the thermal dimorphic fungi. Moreover, dcl-2 demonstrated higher mRNA expression levels in mycelial than yeast phase by 7 folds (P<0.001). Northern blot analysis confirmed the expression of two milRNAs, PM-milR-M1 and PM-milR-M2, only in mycelial phase. Using dcl-1^KO, dcl-2^KO, dcl^DKO and qde-2^KO deletion mutants, we showed that the biogenesis of both milRNAs were dependent on dcl-2 but not dcl-1 or qde-2. The mRNA expression levels of three predicted targets of PM-milR-M1 were upregulated in knockdown strain PM-milR-M1^KD, supporting regulatory function of milRNAs.

Conclusions/Significance

Our findings provided the first evidence for differential expression of milRNAs in different growth phases of thermal dimorphic fungi and shed light on the evolution of fungal proteins involved in milRNA biogenesis and possible role of post-transcriptional control in governing thermal dimorphism.

Penicillium marneffei is the most important thermal dimorphic pathogenic fungus in Southeast Asia. Despite findings on diverse genes and mechanisms involved in dimorphic switching, the key to signally pathways governing the switch is still unknown. Since miRNAs are important regulatory molecules in eukaryotes, we attempt to define if miRNAs are expressed in different growth phases of P. marneffei. Using high-throughput sequencing, we identified 24 potential milRNA candidates in P. marneffei, which were more abundantly expressed in mycelial than yeast phase. Two genes, dcl-1 and dcl-2, encoding Dicer-like proteins and the gene, qde-2, encoding Argonaute-like protein, were also identified. Phylogenetic analysis showed that dcl-2 of P. marneffei was more closely related to the homologues in other thermal dimorphic pathogenic fungi than to Penicillium chrysogenum and Aspergillus spp.. dcl-2 demonstrated higher mRNA levels in mycelial than yeast phase. Northern blot analysis confirmed expression of two milRNAs, PM-milR-M1 and PM-milR-M2, only in mycelial phase, whose expression was dependent on dcl-2 but not dcl-1 or qde-2. The mRNA levels of three predicted targets of PM-milR-M1 were upregulated in knockdown strain PM-milR-M1^KD, supporting its regulatory function. This study represents the first discovery of milRNAs in thermal dimorphic fungi, with differential expression in different growth phases.

Identification and characterization of a novel paramyxovirus, porcine parainfluenza virus 1, from deceased pigs

We describe the discovery and characterization of a novel paramyxovirus, porcine parainfluenza virus 1 (PPIV-1), from swine. The virus was detected in 12 (3.1 %) of 386 nasopharyngeal and two (0.7 %) of 303 rectal swab samples from 386 deceased pigs by reverse transcription-PCR, with viral loads of up to 10(6) copies ml(-1). Complete genome sequencing and phylogenetic analysis showed that PPIV-1 represented a novel paramyxovirus within the genus Respirovirus, being most closely related to human parainfluenza virus 1 (HPIV-1) and Sendai virus (SeV). In contrast to HPIV-1, PPIV-1 possessed a mRNA editing function in the phosphoprotein gene. Moreover, PPIV-1 was unique among respiroviruses in having two G residues instead of three to five G residues following the A6 run at the editing site. Nevertheless, PPIV-1, HPIV-1 and SeV share common genomic features and may belong to a separate group within the genus Respirovirus. The presence of PPIV-1 in mainly respiratory samples suggests a possible association with respiratory disease, similar to HPIV-1 and SeV.

Analysis of multilocus sequence typing schemes for 35 different bacteria revealed that gene loci of 10 bacteria could be replaced to improve cost-effectiveness

Although multilocus sequence typing (MLST) has been widely used for bacterial typing, the contribution of the gene loci to the discriminatory power of each MLST scheme is unknown. We analyzed the discriminatory powers of 36 MLST schemes using all combinations of the 7 loci and contributions of each locus to the schemes. In 10 schemes, sequencing 6 loci can achieve the discriminatory powers of 7 loci. For the other 26 schemes, the median marginal increase in discriminatory power when 7 instead of 6 loci were used is 0.0004. Sequencing the 7 loci of 50 strains each of Pseudomonas aeruginosa and Acinetobacter baumannii revealed that the discriminatory power for P. aeruginosa was 0.9861 when either 6 (without trp) or 7 loci were used and that for A. baumannii was 0.9363 when 5, 6, or 7 loci were used. Genes that have no additional or minimal contribution to the overall discriminatory powers should be replaced.

Identification and complete genome analysis of three novel paramyxoviruses, Tuhoko virus 1, 2 and 3, in fruit bats from China

Among 489 bats of 11 species in China, three novel paramyxoviruses [Tuhokovirus 1, 2 and 3 (ThkPV-1, ThkPV-2 and ThkPV-3)] were discovered in 15 Leschenault's rousettes. Phylogenetically, the three viruses are most closely related to Menangle and Tioman virus. Genome analysis showed that their 3'-leader sequences are unique by possessing GA instead of AG at the 5th and 6th positions. Unlike Menangle and Tioman virus, key amino acids for neuraminidase activity characteristic of rubulavirus attachment proteins are present. The genome of ThkPV-1 represents the largest rubulavirus genome. Unique features between the three viruses include perfect complementary 5'-trailer and 3'-leader sequence and a unique cysteine pair in attachment protein of ThkPV-1, G at + 1 position in all predicted mRNA sequences of ThkPV-2, and amino acid substitutions in the conserved N-terminal motif of nucleocapsid of ThkPV-3. Analysis of phosphoprotein gene mRNA products confirmed mRNA editing. Antibodies to the viruses were detected in 48–60% of Leschenault's rousettes.