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.

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.

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.

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.

Feline morbillivirus, a previously undescribed paramyxovirus associated with tubulointerstitial nephritis in domestic cats.

We describe the discovery and isolation of a paramyxovirus, feline morbillivirus (FmoPV), from domestic cat (Felis catus). FmoPV RNA was detected in 56 (12.3%) of 457 stray cats (53 urine, four rectal swabs, and one blood sample) by RT-PCR. Complete genome sequencing of three FmoPV strains showed genome sizes of 16,050 bases, the largest among morbilliviruses, because of unusually long 5' trailer sequences of 400 nt. FmoPV possesses identical gene contents (3'-N-P/V/C-M-F-H-L-5') and is phylogenetically clustered with other morbilliviruses. IgG against FmoPV N protein was positive in 49 sera (76.7%) of 56 RT-PCR-positive cats, but 78 (19.4%) of 401 RT-PCR-negative cats (P < 0.0001) by Western blot. FmoPV was isolated from CRFK feline kidney cells, causing cytopathic effects with cell rounding, detachment, lysis, and syncytia formation. FmoPV could also replicate in subsequent passages in primate Vero E6 cells. Infected cell lines exhibited finely granular and diffuse cytoplasmic fluorescence on immunostaining for FmoPV N protein. Electron microscopy showed enveloped virus with typical "herringbone" appearance of helical N in paramyxoviruses. Histological examination of necropsy tissues in two FmoPV-positive cats revealed interstitial inflammatory infiltrate and tubular degeneration/necrosis in kidneys, with decreased cauxin expression in degenerated tubular epithelial cells, compatible with tubulointerstitial nephritis (TIN). Immunohistochemical staining revealed FmoPV N protein-positive renal tubular cells and mononuclear cells in lymph nodes. A case-control study showed the presence of TIN in seven of 12 cats with FmoPV infection, but only two of 15 cats without FmoPV infection (P < 0.05), suggesting an association between FmoPV and TIN.

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.

Complete genome sequence of a novel paramyxovirus, Tailam virus, discovered in Sikkim rats.

We discovered a novel paramyxovirus, Tailam virus, of subfamily Paramyxovirinae, in the kidneys and spleens of Sikkim rats. The coding potential of its genome (3'-N-P/V/C-M-F-SH-TM-G-L-5') is similar to those of Beilong virus and J virus, with putative proteins having 59.1 to 94.4% and 23.8 to 80.1% amino acid identities to those of Beilong virus and J virus, respectively.

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.

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.

Novel Picobirnaviruses in Respiratory and Alimentary Tracts of Cattle and Monkeys with Large Intra- and Inter-Host Diversity.

Picobirnaviruses (PBVs) are mostly found in animal alimentary samples. In this study, among 576 respiratory specimens from 476 mammals and 100 chickens, genogroup I PBVs were detected in three cattle and three monkeys, and a genogroup II PBV-positive sample was collected from one cattle specimen. More than one PBV sequence type was observed in two and one genogroup I PBV-positive samples from cattle and monkeys, respectively. Twenty-four complete/near-complete segments 2 (nine from respiratory and 15 from alimentary samples) from the cattle and monkey genogroup I PBVs and one complete segment 2 from the cattle genogroup II PBV were sequenced. Similar to other studies, the cattle PBVs also showed a high diversity. In contrast, the monkey PBVs observed in this study were clustered into three distinct clades. Within each clade, all the sequences showed >99% amino acid identities. This unique phenomenon is probably due to the fact that monkeys in our locality reside in separated troops with minimal inter-troop contact.

Comparative genome and evolutionary analysis of naturally occurring Beilong virus in brown and black rats.

Recently, we reported the presence of Beilong virus in spleen and kidney samples of brown rats and black rats, suggesting that these rodents could be natural reservoirs of Beilong virus. In this study, four genomes of Beilong virus from brown rats and black rats were sequenced. Similar to the Beilong virus genome sequenced from kidney mesangial cell line culture, those of J-virus from house mouse and Tailam virus from Sikkim rats, these four genomes from naturally occurring Beilong virus also contain the eight genes (3'-N-P/V/C-M-F-SH-TM-G-L-5'). In these four genomes, the attachment glycoprotein encoded by the G gene consists of 1046 amino acids; but for the original Beilong virus genome sequenced from kidney mesangial cell line, the G CDS was predicted to be prematurely terminated at position 2205 (TGG→TAG), resulting in a 734-amino-acid truncated G protein. This phenomenon of a lack of nonsense mutation in naturally occurring Beilong viruses was confirmed by sequencing this region of 15 additional rodent samples. Phylogenetic analyses showed that the cell line and naturally occurring Beilong viruses were closely clustered, without separation into subgroups. In addition, these viruses were further clustered with J-virus and Tailam virus, with high bootstrap supports of >90%, forming a distinct group in Paramyxoviridae. Brown rats and black rats are natural reservoirs of Beilong virus. Our results also supports that the recently proposed genus, Jeilongvirus, should encompass Beilong virus, J-virus and Tailam virus as members.

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 × 104 to 5.9 × 106/ml) and 15 segment 2 (viral loads 4.1 × 103 to 1.3 × 106/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.

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.