Genomic Organization Analysis Research Articles

Molecular characterization of a new botybirnavirus that infects Alternaria sp. from tobacco.

The genus Alternaria comprises many important fungal pathogens that infect a wide variety of organisms. In this report, we present the discovery of a new double-stranded RNA (dsRNA) mycovirus called Alternaria botybirnavirus 2 (ABRV2) from a phytopathogenic strain, XC21-21C, of Alternaria sp. isolated from diseased tobacco leaves in China. The ABRV2 genome consists of two dsRNA components, namely dsRNA1 and dsRNA2, with lengths of 6,162 and 5,865 base pairs (bp), respectively. Each of these genomic dsRNAs is monocistronic, encoding hypothetical proteins of 201.6 kDa (P1) and 2193.3 kDa (P2). ABRV2 P1 and P2 share 50.54% and 63.13% amino acid sequence identity with the corresponding proteins encoded by dsRNA1 of Alternaria botybirnavirus 1 (ABRV1). Analysis of its genome organization and phylogenetic analysis revealed that ABRV2 is a new member of the genus Botybirnavirus.

The complete plastome sequences of invasive weed Parthenium hysterophorus: genome organization, evolutionary significance, structural features, and comparative analysis

Parthenium hysterophorus, a globally widespread weed, poses a significant threat to agricultural ecosystems due to its invasive nature. We investigated the chloroplast genome of P. hysterophorus in this study. Our analysis revealed that the chloroplast genome of P. hysterophorus spans a length of 151,881 base pairs (bp). It exhibits typical quadripartite structure commonly found in chloroplast genomes, including inverted repeat regions (IR) of 25,085 bp, a small single copy (SSC) region of 18,052 bp, and a large single copy (LSC) region of 83,588 bp. A total of 129 unique genes were identified in P. hysterophorus chloroplast genomes, including 85 protein-coding genes, 36 tRNAs, and eight rRNAs genes. Comparative analysis of the P. hysterophorus plastome with those of related species from the tribe Heliantheae revealed both conserved structures and intriguing variations. While many structural elements were shared among the species, we identified a rearrangement in the large single-copy region of P. hysterophorus. Moreover, our study highlighted notable gene divergence in several specific genes, namely matK, ndhF, clpP, rps16, ndhA, rps3, and ndhD. Phylogenetic analysis based on the 72 shared genes placed P. hysterophorus in a distinct clade alongside another species, P. argentatum. Additionally, the estimated divergence time between the Parthenium genus and Helianthus (sunflowers) was approximately 15.1 million years ago (Mya). These findings provide valuable insights into the evolutionary history and genetic relationships of P. hysterophorus, shedding light on its divergence and adaptation over time.

Co-infection with two novel mycoviruses affects the biocontrol activity of Trichoderma polysporum

This study reports the whole genome sequence of two double-stranded RNA viruses (mycoviruses) that co-infected an isolate of Trichoderma polysporum NFCF205 and the interactions between fungus and virus. Based on genome size, organization, sequence similarity, and phylogenic analysis, one of these is best characterized as a new member of the proposed family “Fusagraviridae”, and accordingly suggest that it be referred to as Trichoderma polysporum fusagravirus 1 (TpFGV1). The other is best classified as a new strain of the genus Alphapartitivirus in the family Partitiviridae, and should be designated Trichoderma polysporum partitivirus 1 from T. polysporum NFCF205 (TpPV1). By single spore isolation, virus-free or singly-infected isogenic strains were obtained and compared to determine the biological functions of the viruses. The colony morphology of the co-infected parental strain showed the characteristic spiraled hairy and fluffy mycelia, while both the singly-infected and virus-free strain rarely showed the conspicuous aerial mycelia. The growth rate of the co-infected strain was slightly, but nevertheless statistically significantly, higher on PDA. Interestingly, conidia production was observed to have significantly increased in the co-infected strain while no difference was seen in the others. Growth inhibition measured by dual-culturing with Rhizoctonia solani revealed a significant difference between the co-infected strain and a TpPV1-infected strain. However, the secreted soluble compounds in the co-infected strain showed the least antimicrobial activity. Moreover, antimicrobial activities via volatile organic compounds (VOCs), which signify strong biocontrol activity by Trichoderma spp., was the lowest in the co-infected strain. The profiles and amount of representative VOCs in the assessed strains, as measured by gas chromatography-mass spectrometry (GC/MS/MS), were affected by virus-infection. Virus infection affected colony morphology, conidia production, and biocontrol activities of T. polysporum. Most interestingly, virus-virus interaction was observed and measured in the co-infected strain, as their symptoms differed from those of the singly-infected strains.

Functional genomic studies in diverse populations to characterize risk loci for Alzheimer Disease

AbstractBackgroundAs part of the ADSP Functional Genomics consortium, we aim to unravel ancestry specific effects of genomic regulation of AD loci. African Americans (AA) and Hispanics (HI) are at higher risk for Alzheimer’s disease (AD) compared to non‐Hispanic whites (NHW). Even though AA, HI and other diverse groups make up 40% of the US population, over 90% of participants in large‐scale genomic research studies have been NHW, representing European (EU) ancestry. Recent AD research including diverse groups has shown that genes and variants contributing to disease can have different risk effect dependent on the ancestral background they reside on (i.e. APOE). Together these clearly indicate a need to expand genome‐wide and gene‐specific analyses of genomic organization and regulation for other ancestries contributing to diverse, admixed populations, such as African (AF) or Amerindian (AI) ancestries.MethodWe analyzed interaction of variants in the APOE region with AD risk in WGS of diverse populations to identify modifiers of the APOEɛ4 risk effect. Additionally, to assess overall differences in chromatin regulation genome‐wide across ancestries, we performed bulk and single cell chromatin accessibility (ATAC‐seq) and bulk chromatin interaction (Hi‐C) assays in brains from admixed individuals with high EU, AF (>95%) or AI (>75%) ancestry contribution.ResultWe identified both close‐by and distant modifier variants (in the downstream PSG gene family region) inferring protection against the APOEɛ4 risk on the AF ancestry. Single cell ATAC‐seq data in brain showed higher accessibility of the APOE region and large regions of chr19 in astrocytes, but not microglia, of EU versus AF ancestry. Lastly, bulk Hi‐C analyses in brain strongly supported the hypothesis of ancestry dependent regulation by genome‐wide detection of differential chromatin interactions comparing EU to AF ancestry.ConclusionOur analyses using complementary functional genomic assay methods show clear differences of genome regulation between EU and AF ancestry. Overall, these data support genome regulation to be strongly influenced by ancestral origin. The generated data on samples with high single ancestral contributions will help elucidate risk effects seen for AD related genes in different population groups, as we describe for the APOE region.

Complete genome sequence of a novel virus isolated from the phytopathogenic fungus Ceratobasidium sp. AG-A strain SHX-YJLC-1.

A novel mycovirus, Ceratobasidium bipartite virus 1 (CBV1), was identified in Ceratobasidium sp. AG-A strain SHX-YJLC-1 isolated from diseased potato stems. The complete genome of CBV1 consists of two double-stranded RNA (dsRNA) segments: dsRNA1 (2311 bp) and dsRNA2 (1761 bp). dsRNA1 contains a single open reading frame (ORF1) encoding an RNA-dependent RNA polymerase (RdRp), while dsRNA2 contains a single ORF (ORF2) encoding a hypothetical protein (HP) with unknown function. BLASTp analysis revealed that RdRp (75.04%) and HP (61.86%) encoded by the two ORFs have the highest sequence similarity to their counterparts in Rhizoctonia solani dsRNA virus 11 (RsRV11). The genome organization and phylogenetic analysis indicated that the closest relatives to CBV1 are members of the proposed family "Bipartitiviridae". Based on the collective results, CBV1 is inferred to be a new member of the proposed family "Bipartitiviridae". This is the first report on the complete genome sequence of the novel bipartitivirus CBV1, which infects Ceratobasidium sp. AG-A strain SHX-YJLC-1.

Diversity of Mycoviruses Present in Strains of Binucleate Rhizoctonia and Multinucleate Rhizoctonia, Causal Agents for Potato Stem Canker or Black Scurf.

In this study, the diversity of putative mycoviruses present in 66 strains of binucleate Rhizoctonia (BNR, including anastomosis group (AG)-A, AG-Fa, AG-K, and AG-W) and 192 strains of multinucleate Rhizoctonia (MNR, including AG-1-IA, AG-2-1, AG-3 PT, AG-4HGI, AG-4HGII, AG-4HGIII, and AG-5), which are the causal agents of potato stem canker or black scurf, was studied using metatranscriptome sequencing. The number of contigs related to mycoviruses identified from BNR and MNR was 173 and 485, respectively. On average, each strain of BNR accommodated 2.62 putative mycoviruses, while each strain of MNR accommodated 2.53 putative mycoviruses. Putative mycoviruses detected in both BNR and MNR contained positive single-stranded RNA (+ssRNA), double-stranded RNA (dsRNA), and negative single-stranded RNA (-ssRNA) genomes, with +ssRNA genome being the prevalent nucleic acid type (82.08% in BNR and 75.46% in MNR). Except for 3 unclassified, 170 putative mycoviruses found in BNR belonged to 13 families; excluding 33 unclassified, 452 putative mycoviruses found in MNR belonged to 19 families. Through genome organization, multiple alignments, and phylogenetic analyses, 4 new parititviruses, 39 novel mitoviruses, and 4 new hypoviruses with nearly whole genome were detected in the 258 strains of BNR and MNR.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains.

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary Tcell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ Tcells, highlighting the importance of host genome organization in HIV-1 infection.

Lineage Tracing and Single-Nucleus Multiomics Reveal Novel Features of Adaptive and Maladaptive Repair after Acute Kidney Injury.

Understanding the mechanisms underlying adaptive and maladaptive renal repair after AKI and their long-term consequences is critical to kidney health. The authors used lineage tracing of cycling cells and single-nucleus multiomics (profiling transcriptome and chromatin accessibility) after AKI. They demonstrated that AKI triggers a cell-cycle response in most epithelial and nonepithelial kidney cell types. They also showed that maladaptive proinflammatory proximal tubule cells (PTCs) persist until 6 months post-AKI, although they decreased in abundance over time, in part, through cell death. Single-nucleus multiomics of lineage-traced cells revealed regulatory features of adaptive and maladaptive repair. These included activation of cell state-specific transcription factors and cis-regulatory elements, and effects in PTCs even after adaptive repair, weeks after the injury event. AKI triggers a proliferative response as part of an intrinsic cellular repair program, which can lead to adaptive renal repair, restoring kidney structure and function, or maladaptive repair with the persistence of injured proximal tubule cells (PTCs) and an altered kidney structure. However, the cellular and molecular understanding of these repair programs is limited. To examine chromatin and transcriptional responses in the same cell upon ischemia-reperfusion injury (IRI), we combined genetic fate mapping of cycling ( Ki67+ ) cells labeled early after IRI with single-nucleus multiomics-profiling transcriptome and chromatin accessibility in the same nucleus-and generated a dataset of 83,315 nuclei. AKI triggered a broad cell cycle response preceded by cell type-specific and global transcriptional changes in the nephron, the collecting and vascular systems, and stromal and immune cell types. We observed a heterogeneous population of maladaptive PTCs throughout proximal tubule segments 6 months post-AKI, with a marked loss of maladaptive cells from 4 weeks to 6 months. Gene expression and chromatin accessibility profiling in the same nuclei highlighted differences between adaptive and maladaptive PTCs in the activity of cis-regulatory elements and transcription factors, accompanied by corresponding changes in target gene expression. Adaptive repair was associated with reduced expression of genes encoding transmembrane transport proteins essential to kidney function. Analysis of genome organization and gene activity with single-cell resolution using lineage tracing and single-nucleus multiomics offers new insight into the regulation of renal injury repair. Weeks to months after mild-to-moderate IRI, maladaptive PTCs persist with an aberrant epigenetic landscape, and PTCs exhibit an altered transcriptional profile even following adaptive repair.

Accurate Chromosome Identification in the Prunus Subgenus Cerasus (Prunus pseudocerasus) and its Relatives by Oligo-FISH.

A precise, rapid and straightforward approach to chromosome identification is fundamental for cytogenetics studies. However, the identification of individual chromosomes was not previously possible for Chinese cherry or other Prunus species due to the small size and similar morphology of their chromosomes. To address this issue, we designed a pool of oligonucleotides distributed across specific pseudochromosome regions of Chinese cherry. This oligonucleotide pool was amplified through multiplex PCR with specific internal primers to produce probes that could recognize specific chromosomes. External primers modified with red and green fluorescence tags could produce unique signal barcoding patterns to identify each chromosome concomitantly. The same oligonucleotide pool could also discriminate all chromosomes in other Prunus species. Additionally, the 5S/45S rDNA probes and the oligo pool were applied in two sequential rounds of fluorescence in situ hybridization (FISH) localized to chromosomes and showed different distribution patterns among Prunus species. At the same time, comparative karyotype analysis revealed high conservation among P. pseudocerasus, P. avium, and P. persica. Together, these findings establish this oligonucleotide pool as the most effective tool for chromosome identification and the analysis of genome organization and evolution in the genus Prunus.

Identification, Genomic Organization, and Comprehensive Expression Analysis Reveals the Implication of Cicer arietinum SKP1-like Genes in Abiotic Stress

Identification, Genomic Organization, and Comprehensive Expression Analysis Reveals the Implication of Cicer arietinum SKP1-like Genes in Abiotic Stress

Genome-Wide Identification, Evolutionary Analysis, and Expression Patterns of Cathepsin Superfamily in Black Rockfish (Sebastes schlegelii) following Aeromonas salmonicida Infection.

Cathepsins are lysosomal cysteine proteases belonging to the papain family and play crucial roles in intracellular protein degradation/turnover, hormone maturation, antigen processing, and immune responses. In the present study, 18 cathepsins were systematically identified from the fish S. schlegelii genome. Phylogenetic analysis indicated that cathepsin superfamilies are categorized into eleven major clusters. Synteny and genome organization analysis revealed that whole-genome duplication led to the expansion of S. schlegelii cathepsins. Evolutionary rate analyses indicated that the lowest Ka/Ks ratios were observed in CTSBa (0.13) and CTSBb (0.14), and the highest Ka/Ks ratios were observed in CTSZa (1.97) and CTSZb (1.75). In addition, cathepsins were ubiquitously expressed in all examined tissues, with high expression levels observed in the gill, intestine, head kidney, and spleen. Additionally, most cathepsins were differentially expressed in the head kidney, gill, spleen, and liver following Aeromonas salmonicida infection, and their expression signatures showed tissue-specific and time-dependent patterns. Finally, protein–protein interaction network (PPI) analyses revealed that cathepsins are closely related to a few immune-related genes, such as interleukins, chemokines, and TLR genes. These results are expected to be valuable for comparative immunological studies and provide insights for further functional characterization of cathepsins in fish species.

High Diversity of Novel Viruses in the Tree Pathogen Phytophthora castaneae Revealed by High-Throughput Sequencing of Total and Small RNA.

Phytophthora castaneae, an oomycete pathogen causing root and trunk rot of different tree species in Asia, was shown to harbor a rich diversity of novel viruses from different families. Four P. castaneae isolates collected from Chamaecyparis hodginsii in a semi-natural montane forest site in Vietnam were investigated for viral presence by traditional and next-generation sequencing (NGS) techniques, i.e., double-stranded RNA (dsRNA) extraction and high-throughput sequencing (HTS) of small RNAs (sRNAs) and total RNA. Genome organization, sequence similarity, and phylogenetic analyses indicated that the viruses were related to members of the order Bunyavirales and families Endornaviridae, Megabirnaviridae, Narnaviridae, Totiviridae, and the proposed family “Fusagraviridae.” The study describes six novel viruses: Phytophthora castaneae RNA virus 1–5 (PcaRV1-5) and Phytophthora castaneae negative-stranded RNA virus 1 (PcaNSRV1). All six viruses were detected by sRNA sequencing, which demonstrates an active RNA interference (RNAi) system targeting viruses in P. castaneae. To our knowledge, this is the first report of viruses in P. castaneae and the whole Phytophthora major Clade 5, as well as of the activity of an RNAi mechanism targeting viral genomes among Clade 5 species. PcaRV1 is the first megabirnavirus described in oomycetes and the genus Phytophthora.

Characterization and Complete Genomic Analysis of Vibrio Parahaemolyticus-Infecting Phage KIT05.

Vibrio parahaemolyticus is a bacterial pathogen in marine aquaculture systems and a major cause of food-borne illnesses worldwide. In the present study, Vibrio phage KIT05 was isolated from water collected from a shrimp farm in the Mekong Delta, Vietnam. It was characterized based on its morphology, growth curve, lytic properties, and genome sequence. Under the electron microscope, KIT05 particles had an icosahedral head with a diameter of 62.3nm and a short tail of 24.1nm. The one-step growth curve of KIT05 showed that its latency time was approximately 40min and burst size was 18 plaque-forming units/cell. The genome of KIT05 comprises 50,628bp with a GC content of 41.63%. It contains 60 open reading frames that are encoded within both strands and four tRNAs. The presence of direct terminal repeats of 130bp at both ends of the KIT05 DNA was determined. According to phage morphology, genomic organization, and phylogeny analysis, Vibrio phage KIT05 was classified into the family Podoviridae. The genome annotation revealed that KIT05 had no virulent or lysogenic genes. This study may help identify a novel candidate for developing biocontrol agents for Vibrio parahaemolyticus.

CC and CXC chemokines in turbot (Scophthalmus maximus L.): Identification, evolutionary analyses, and expression profiling after Aeromonas salmonicida infection

Chemokines are a superfamily of structurally related cytokines, which exert essential roles in guiding cell migration in development, homeostasis, and immunity. CC and CXC chemokines are the two major subfamilies in teleost species. In this study, a total of seventeen CC and CXC chemokines, with inclusion of twelve CC and five CXC chemokines, were systematically identified from the turbot genome, making turbot the teleost harboring the least number of CC and CXC chemokines among all teleost species ever reported. Phylogeny, synteny, and genomic organization analyses were performed to annotate these genes, and multiple chemokine genes were identified in the turbot genome, due to the tandem duplications (CCL19 and CCL20), the whole genome duplications (CCL20, CCL25, and CXCL12), and the teleost-specific members (CCL34-36, CCL44, and CXCL18). In addition, chemokines were ubiquitously expressed in nine examined healthy tissues, with high expression levels observed in liver, gill, and spleen. Moreover, most chemokines were significantly differentially expressed in gill and spleen after Aeromonas salmonicida infection, and exhibited tissue-specific and time-dependent manner. Finally, protein-protein interaction network (PPI) analysis indicated that turbot chemokines interacted with a few immune-related genes such as interleukins, cathepsins, stats, and TLRs. These results should be valuable for comparative immunological studies and provide insights for further functional characterization of chemokines in teleost.

Isolation of High Molecular Weight DNA from the Model Beetle Tribolium for Nanopore Sequencing.

The long-read Nanopore sequencing has been recently applied for assembly of complex genomes and analysis of linear genome organization. The most critical factor for successful long-read sequencing is extraction of high molecular weight (HMW) DNA of sufficient purity and quantity. The challenges associated with input DNA quality are further amplified when working with extremely small insects with hard exoskeletons. Here, we optimized the isolation of HMW DNA from the model beetle Tribolium and tested for use in Nanopore sequencing. We succeeded in overcoming all the difficulties in HMW handling and library preparation that were encountered when using published protocols and commercial kits. Isolation of nuclei and subsequent purification of DNA on an anion-exchange chromatography column resulted in genomic HMW DNA that was efficiently relaxed, of optimal quality and in sufficient quantity for Nanopore MinION sequencing. DNA shearing increased average N50 read values up to 26 kb and allowed us to use a single flow cell in multiple library loads for a total output of more than 13 Gb. Although our focus was on T. castaneum and closely related species, we expect that this protocol, with appropriate modifications, could be extended to other insects, particularly beetles.

Genome-guided investigation of anti-inflammatory sesterterpenoids with 5-15 trans-fused ring system from phytopathogenic fungi.

Fungal terpenoids catalyzed by bifunctional terpene synthases (BFTSs) possess interesting bioactive and chemical properties. In this study, an integrated approach of genome mining, heterologous expression, and in vitro enzymatic activity assay was used, and these identified a unique BFTS sub-clade critical to the formation of a 5-15 trans-fused bicyclic sesterterpene preterpestacin I (1). The 5-15 bicyclic BFTS gene clusters were highly conserved but showed relatively wide phylogenetic distribution across several species of the diverged fungal classes Dothideomycetes and Sordariomycetes. Further genomic organization analysis of these homologous biosynthetic gene clusters from this clade revealed a glycosyltransferase from the graminaceous pathogen Bipolaris sorokiniana isolate BS11134, which was absent in other 5-15 bicyclic BFTS gene clusters. Targeted isolation guided by BFTS gene deletion led to the identification of two new sesterterpenoids (4, and 6) from BS11134. Compounds 2 and 4 showed moderate effects on LPS-induced nitrous oxide production in the murine macrophage-like cell line RAW264.7 with in vitro inhibition rates of 36.6 ± 2.4% and 24.9 ± 2.1% at 10 μM, respectively. The plausible biosynthetic pathway of these identified compounds was proposed as well. This work revealed that phytopathogenic fungi can serve as important sources of active terpenoids via systematic analysis of the genomic organization of BFTS biosynthetic gene clusters, their phylogenetic distribution in fungi, and cyclization properties of their metabolic products. KEY POINTS: • Genome mining of the first BFTS BGC harboring a glycosyltransferase. • Gene-deletion guided isolation revealed three novel 5-15 bicyclic sesterterpenoids. • Biosynthetic pathway of isolated sesterterpenoids was proposed.

Complete Nucleotide Sequence, Genome Organization, and Comparative Genomic Analyses of Citrus Yellow-Vein Associated Virus (CYVaV)

Citrus yellow-vein disease (CYVD) was first reported in California in 1957. We now report that CYVD is associated with a virus-like agent, provisionally named citrus yellow-vein associated virus (CYVaV). The CYVaV RNA genome has 2,692 nucleotides and codes for two discernable open reading frames (ORFs). ORF1 encodes a protein of 190 amino acid (aa) whereas ORF2 is presumably generated by a −1 ribosomal frameshifting event just upstream of the ORF1 termination signal. The frameshift product (717 aa) encodes the RNA-dependent RNA polymerase (RdRp). Phylogenetic analyses suggest that CYVaV is closely related to unclassified virus-like RNAs in the family Tombusviridae. Bio-indexing and RNA-seq experiments indicate that CYVaV can induce yellow vein symptoms independently of known citrus viruses or viroids.

Understanding 3D genome organization by multidisciplinary methods.

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

Comparative characterization of two putative duplicated sodium channel genes in the red flour beetle, Tribolium castaneum

The insect voltage-gated sodium channel is the primary target of pyrethroids and novel efficient insecticides such as indoxacarb and metaflumizone. In this study, we cloned and characterized two putative sodium channel genes, TcNav1 and TcNav2, in Tribolium castaneum. The composite TcNav1 and TcNav2 encode a protein of 2045 and 2037 amino acid residues, sharing 76.1% and 75.5% amino acid identity with Drosophila para, respectively. Comparative analysis of genomic organization showed that TcNav1 and TcNav2 contain 26 and 27 exons, respectively. Analysis of the expression patterns showed that the mRNA levels of TcNav1 and TcNav2 were predominantly expressed in the head. RNAi-mediated knockdown of both TcNav1 and TcNav2 adversely affected adult emergence and significantly decreased sensitivity to deltamethrin. Significantly reduced pupation rate and sensitivity to beta-cypermethrin were observed after injection of siRNA targeting TcNav1 but not TcNav2. Taken together, we provide evidence that sodium channel gene has undergone duplication in T. castaneum, resulting in diversified developmental and toxicological functions.

Molecular and biological characterization of a new mulberry idaeovirus

A virus-like disease with symptoms including mosaic structure, deformation, vein clearing and necrosis on the leaves and deformation, crumbling, and scab on the fruits was detected in black mulberry trees (Morus nigra L.) in Kayseri province of Turkey. A novel positive single‐stranded RNA virus with a bipartite genome and the mulberry badnavirus 1 (MBV-1) were detected in the black mulberry trees by high throughput sequencing and bioinformatic analyses. The novel virus RNA1 (5,796/7 nt) encodes a polyprotein (1,808 aa, 204.31 kDa) with three conserved domains, [MTR (aa 294–705), Hel (aa 971–1,226) and RdRp (aa 1,348–1,788)], whereas RNA2 (2,243 nt) encodes two putative proteins, MP (374 aa, 40.98 kDa), and CP (272 aa, 30.59 kDa), separated by an intergenic region of 97 nt. The highest amino acids identities were 70, 57 and 70 % with raspberry bushy dwarf virus (RBDV) for ORF1, MP and CP genes, respectively. The genome organization and phylogenetic analyses suggested that the novel virus is likely a putative new member of the genus Idaeovirus and it has been tentatively named black mulberry idaeovirus (BMIV). Virus survey showed both the BMIV and MBV-1 are likely prevalent in the region. Seven complete (six Turkish and one Iranian) and 41 partial genome sequences of the BMIV isolates revealed moderate genetic diversity (0.033 ± 0.001 %, 0.020 ± 0.002 % and 0.016 ± 0.002 % for RNA1, RNA2, and partial genomes, respectively). Both the BMIV and MBV-1 were detected in all tested pollens (n = 24, 100 %), in seed-borne balck mulberry saplings (n = 96, 100 %).This situation clearly revealed the potential spread risk of both viruses in black mulberry plantations and the necessity of taking precautions.