Parallelization of MAFFT for large-scale multiple sequence alignments

Large-Scale Comparison of Toxin and Antitoxins in Listeria monocytogenes

Toxin–antitoxin systems (TASs) are widely distributed in prokaryotes and encode pairs of genes involved in many bacterial biological processes and mechanisms, including pathogenesis. The TASs have not been extensively studied in Listeria monocytogenes (Lm), a pathogenic bacterium of the Firmicutes phylum causing infections in animals and humans. Using our recently published TASmania database, we focused on the known and new putative TASs in 352 Listeria monocytogenes genomes and identified the putative core gene TASs (cgTASs) with the Pasteur BIGSdb-Lm database and, by complementarity, the putative accessory gene TAS (acTASs). We combined the cgTASs with those of an additional 227 L. monocytogenes isolates from our previous studies containing metadata information. We discovered that the differences in 14 cgTAS alleles are sufficient to separate the four main lineages of Listeria monocytogenes. Analyzing these differences in more details, we uncovered potentially co-evolving residues in some pairs of proteins in cgTASs, probably essential for protein–protein interactions within the TAS complex.

Bat IFITM3 restriction depends on S-palmitoylation and a polymorphic site within the CD225 domain

Host interferon-induced transmembrane proteins (IFITMs) are broad-spectrum antiviral restriction factors. Of these, IFITM3 potently inhibits viruses that enter cells through acidic endosomes, many of which are zoonotic and emerging viruses with bats (order Chiroptera) as their natural hosts. We previously demonstrated that microbat IFITM3 is antiviral. Here, we show that bat IFITMs are characterized by strong adaptive evolution and identify a highly variable and functionally important site-codon 70-within the conserved CD225 domain of IFITMs. Mutation of this residue in microbat IFITM3 impairs restriction of representatives of four different virus families that enter cells via endosomes. This mutant shows altered subcellular localization and reduced S-palmitoylation, a phenotype copied by mutation of conserved cysteine residues in microbat IFITM3. Furthermore, we show that microbat IFITM3 is S-palmitoylated on cysteine residues C71, C72, and C105, mutation of each cysteine individually impairs virus restriction, and a triple C71A-C72A-C105A mutant loses all restriction activity, concomitant with subcellular re-localization of microbat IFITM3 to Golgi-associated sites. Thus, we propose that S-palmitoylation is critical for Chiropteran IFITM3 function and identify a key molecular determinant of IFITM3 S-palmitoylation.

Subtypes of Neurohypophyseal Nonapeptide Receptors and Their Functions in Rat Kidneys

The nonapeptides of neurohypophysis, vasotocin and mesotocin, detected in most vertebrates, are replaced by vasopressin and oxytocin in mammals. Using bioinformatics methods, we determined the spectrum of receptor subtypes for these hormones in mammals and their physiological effects in the kidneys of rats. A search for sequences similar to the vertebrate vasotocin receptor by proteomes and transcriptomas of nine mammalian species and the rat genome revealed three subtypes of vasopressin receptors (V1a, V1b, and V2) and one type of oxytocin receptors. In the kidneys of non-anesthetized rats, which received a water load of 2 ml per 100 g of body weight, three effects of vasopressin were revealed: 1) increased reabsorption of water and sodium, 2) increased excretion of potassium ions, and 3) increased excretion of sodium ions. It has been suggested that each of the effects on the kidney is associated with selective stimulation of the vasopressin receptor subtypes V2, V1b, and V1a depending on the concentration of nonapeptide. In experiments on non-anaesthetized rats with a water load, the injection of oxytocin reduces the reabsorption of solute-free water in the kidneys and increases the excretion of sodium ions. The possible physiological mechanisms behind the realization of both effects with the participation of a single type of oxytocin receptors are being analyzed. Thus, the spectrum of activated receptor subtypes varies depending on the current concentration of neurohypophyseal hormones, as a result of which the predominant effect on renal function changes, which ensures precise regulation of water-salt homeostasis.

Long-term virus evolution in nature

Viruses spread to give rise to epidemics and pandemics, and some key parameters that include virus and host population numbers determine virus persistence or extinction in nature. Viruses evolve at different rates depending on the polymerase copying fidelity during genome replication and a number of environmental influences. Calculated rates of evolution in nature vary depending on the time interval between virus isolations. In particular, intrahost evolution is generally more rapid that interhost evolution, and several possible mechanisms for this difference are considered. The mechanisms by which the error-prone viruses evolve are very unlikely to render the operation of a molecular clock (constant rate of incorporation of mutations in the evolving genomes), although a clock is assumed in many calculations. Several computational tools permit the alignment of viral sequences and the establishment of phylogenetic relationships among viruses. The evolution of the virus in the form of dynamic mutant clouds in each infected individual, together with multiple environmental parameters renders the emergence and reemergence of viral pathogens an unpredictable event, another facet of biological complexity.

Efficient synthesis of L-phosphinothricin using a novel aminoacylase mined from Stenotrophomonas maltophilia

L-phosphinothricin (L-PPT) is a competitive and environmentally friendly herbicide. To develop an efficient approach for synthesis of l-PPT, a kinetic resolution route with a novel aminoacylase (SmAcy) mined from Stenotrophomonas maltophilia using N-acetyl-PPT as a substrate was first constructed. This SmAcy exhibited high hydrolytic activity and excellent enantioselectivity (E > 200) toward N-acetyl-PPT. Optically pure l-PPT (> 99.9 % ee_p) was acquired with high conversion (> 49 %) within 4 h by the whole cells. On the basis of the docking analysis, a main reason for high enantioselectivity (E > 200) of SmAcy towards l-enantiomer would be that the D-N-acetyl-PPT cannot interact with the key general acid-base residue and the metal ions. A low-cost and simple preparation process of the substrate from commercially available racemic PPT for production of L-PPT was provided. A chemical racemization method of the unreacted D-enantiomer of substrate was also provided to recycle the unwanted substrate enantiomer. This study provides a potential route for the industrial production of L-PPT.

The complete chloroplast genome of Tagetes erecta (Asteroideae), a wildly cultivated ornamental plant

Tagetes erecta (Asteraceae) has been wildly cultivated as ornamental and medicinal plant. Here, we reported the first chloroplast genome sequence of T. erecta. The chloroplast genome size is 152,065 bp with GC content of 37.4%, including a large single-copy (LSC) of 83,895 bp, a small single-copy (SSC) of 18,065 bp, and a pair of 25,048 bp IR (inverted repeat) regions. A total of 132 genes were annotated including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogenetic analysis revealed that T. erecta belongs to the subfamily Asteroideae.

The complete chloroplast genome sequence of a traditional Chinese medicine plant Bulbophyllum disciflorum Rolfe (Orchidaceae)

Bulbophyllum disciflorum is one of Orchidaceae species, which has important ornamental and economic value. Here, we reported the first chloroplast genome sequence of Bulbophyllum. The genome of B. disciflorum is 148,554 bp in length, including a large single-copy (LSC) of 79,001 bp, a small single-copy (SSC) of 16,797 bp, and a pair of inverted repeat regions (IR_a and IR_b) of 26,378bp. It contains 108 unique genes consisting of 74 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. The phylogenetic analysis revealed that Bulbophyllum is sister to the genus of Dendrobium.

Expression patterns of three JAK-STAT pathway genes in feather follicle development during chicken embryogenesis

The Janus kinase (JAK)-signal transducer and activator of transcription (STAT) (JAK-STAT) pathway is shown to restrain the hair follicles in catagen and telogen and prevent anagen reentry in murine hair follicle cycling. The early roles of JAK-STAT pathway genes in skin development remain uncharacterized in mouse and chicken models. Here, we revealed the expression patterns of three JAK-STAT pathway genes (JAK1, JAK2, and TYK2) in chicken embryonic skin at E6-E10 stages which are key to feather follicle morphogenesis. Multiple sequence alignment of the three genes from chicken and other species all showed a closely related homology with birds like quail and goose. Whole mount in situ hybridization (WISH) revealed weak expression of JAK1, JAK2, and TYK2 in chicken skin at E6 and E7, and followed with the focally restricted signals in the feather follicles of neck and body skin located dorsally at E8 for JAK1, E9 for TYK2 and E10 for JAK2 gene. All three genes displayed stronger expression in feather follicles of neck skin than that of body skin. The expression levels of JAK1 and TYK2 were much stronger than those of JAK2. Quantitative real-time PCR (qRT-PCR) analysis revealed the increased expression tendency for JAK2 both in the neck and body skin from E6 to E10, and the much stronger expression in neck and body skin at later stages (E8-E10) than earlier stages (E6 and E7) for JAK1 and TYK2. Overall, these findings suggest that JAK1 and TYK2, not JAK2 are important to specify the feather follicle primordia, and to arrange the proximal-distal axis of feather follicles, respectively, during the morphogenesis of feather follicles in embryonic chicken skin.

Characterization and mutagenesis of Chinese hamster ovary cells endogenous retroviruses to inactivate viral particle release

The Chinese hamster ovary (CHO) cells used to produce biopharmaceutical proteins are known to contain type‐C endogenous retrovirus (ERV) sequences in their genome and to release retroviral‐like particles. Although evidence for their infectivity is missing, this has raised safety concerns. As the genomic origin of these particles remained unclear, we characterized type‐C ERV elements at the genome, transcriptome, and viral particle RNA levels. We identified 173 type‐C ERV sequences clustering into three functionally conserved groups. Transcripts from one type‐C ERV group were full‐length, with intact open reading frames, and cognate viral genome RNA was loaded into retroviral‐like particles, suggesting that this ERV group may produce functional viruses. CRISPR‐Cas9 genome editing was used to disrupt the gag gene of the expressed type‐C ERV group. Comparison of CRISPR‐derived mutations at the DNA and RNA level led to the identification of a single ERV as the main source of the release of RNA‐loaded viral particles. Clones bearing a Gag loss‐of‐function mutation in this ERV showed a reduction of RNA‐containing viral particle release down to detection limits, without compromising cell growth or therapeutic protein production. Overall, our study provides a strategy to mitigate potential viral particle contaminations resulting from ERVs during biopharmaceutical manufacturing.

The complete chloroplast genome sequence of a Chinese endemic plant Dendrobium hancockii Rolfe (Orchidaceae)

Dendrobium hancockii Rolfe is a rare and endangered species endemic to China, with great medicinal value. Here, the first complete chloroplast genome sequence of D. hancockii was reported and characterized. The cpDNA exhibited the typical quadripartite structure of four parts: a long single-copy region, a short single-copy region and two inverted repeats. It encodes 106 genes, consisting of 72 unique protein-coding genes, 30 unique tRNA gene, and 4 unique rRNA genes. The phylogenetic analysis indicated that D. hancockii is basal-most species for the sect. Dendrobium.

"Unicorn from Hades", a new genus of Mysidae (Malacostraca: Mysida) from the Mariana Trench, with a systematic analysis of the deep-sea mysids

Since the 19th century, oceanic explorations have confirmed that the hadal zone (water depth > 6000 m) is not lifeless, but contains many fascinating organisms. Amongst them are the Mysida, which is a group of crustaceans found in many deep-sea trenches. Based on morphological observations and molecular phylogenetic analyses of an undescribed taxon within the subfamily Erythropinae, a new genus of deep-sea mysids, Xenomysis gen. n., is described from the Challenger Deep in the Mariana Trench. The new genus is not assigned to any of the tribes within the Erythropinae, as our analyses do not support the current classification scheme of Erythropinae. The result of a molecular-clock analysis with fossil calibration reveals that several groups of Mysida have independently colonized deep water habitats in different geological periods, from Triassic to Cretaceous. In addition, ancestral state reconstruction analyses indicate the degenerate eyes in both Mysidae and Petalophthalmidae is a result of parallel evolution, and the reduction of compound eyes to both "single fused eyeplate" and "two separate eyeplates" occurred multiple times independently in the evolution of Erythropinae.

Bayesian Estimation of Past Population Dynamics in BEAST 1.10 Using the Skygrid Coalescent Model

Inferring past population dynamics over time from heterochronous molecular sequence data is often achieved using the Bayesian Skygrid model, a nonparametric coalescent model that estimates the effective population size over time. Available in BEAST, a cross-platform program for Bayesian analysis of molecular sequences using Markov chain Monte Carlo, this coalescent model is often estimated in conjunction with a molecular clock model to produce time-stamped phylogenetic trees. We here provide a practical guide to using BEAST and its accompanying applications for the purpose of drawing inference under these models. We focus on best practices, potential pitfalls, and recommendations that can be generalized to other software packages for Bayesian inference. This protocol shows how to use TempEst, BEAUti, and BEAST 1.10 (http://beast.community/; last accessed July 29, 2019), LogCombiner as well as Tracer in a complete workflow.

Comparison of chloroplast genomes of Gynura species: sequence variation, genome rearrangement and divergence studies

Background

Some Gynura species have been reported to be natural anti-diabetic plants. Improvement of their traits towards application relies on hybridization. Clearly, phylogenetic relationships could optimize compatible hybridizations. For flowerings plants, chloroplast genomes have been used to solve many phylogenetic relationships. To date, the chloroplast genome sequences of 4 genera of the tribe Senecioneae have been uploaded to GenBank. The internal relationships within the genus Gynura and the relationship of the genus Gynura with other genera in the tribe Senecioneae need further research.

Results

The chloroplast genomes of 4 Gynura species were sequenced, assembled and annotated. In comparison with those of 12 other Senecioneae species, the Gynura chloroplast genome features were analysed in detail. Subsequently, differences in the microsatellite and repeat types in the tribe were found. From the comparison, it was found that IR expansion and contraction are conserved in the genera Gynura, Dendrosenecio and Ligularia. Compared to other regions on the chloroplast genome, the region from 25,000 to 50,000 bp was not conserved. Seven ndh genes in this region are under purifying selection, with small changes in amino acids. The whole chloroplast genome sequences of 16 Senecioneae species were used to build a phylogenetic tree. Based on the oldest Artemisia pollen fossil, the divergence time was estimated.

Conclusions

Sequencing the chloroplast genomes of 4 Gynura species helps us to solve many problems. The phylogenetic relationships and divergence time among 4 Gynura and 16 Senecioneae species were evaluated by comparing their chloroplast genomes. The phylogenetic relationship of the genera Gynura and Ligularia was different from that observed previous work. In a previous phylogenetic tree, the genus Ligularia belonged to the Tussilagininae subtribe, which was in a lineage that diverged earlier than other genera. Further morphology and genome-wide analyses are needed to clarify the genus relationships.

Pan-genome analysis of Riemerella anatipestifer reveals its genomic diversity and acquired antibiotic resistance associated with genomic islands

Riemerella anatipestifer is a gram-negative bacterium that leads to severe contagious septicemia in ducks, turkeys, chickens, and wild waterfowl. Here, a pan-genome with 32 R. anatipestifer genomes is re-established, and the mathematical model is calculated to evaluate the expansion of R. anatipestifer genomes, which were determined to be open. Average nucleotide identity (ANI) and phylogenetic analysis preliminarily clarify intraspecies variation and distance. Comparative genomic analysis of R. anatipestifer found that horizontal gene transfer events, which provide an expressway for the recruitment of novel functionalities and facilitate genetic diversity in microbial genomes, play a key role in the process of acquiring and transmitting antibiotic-resistance genes in R. anatipestifer. Furthermore, a new antibiotic-resistance gene cluster was identified in the same loci in 14 genomes. The uneven distribution of virulence factors was also confirmed by our results. Our study suggests that the ability to acquire foreign genes (such as antibiotic-resistance genes) increases the adaptability of R. anatipestifer, and the virulence genes with little mobility are highly conserved in R. anatipestifer.

TRIM69 Inhibits Vesicular Stomatitis Indiana Virus

Vesicular stomatitis Indiana virus (VSIV), formerly known as vesicular stomatitis virus (VSV) Indiana (VSVIND), is a model virus that is exceptionally sensitive to the inhibitory action of interferons (IFNs). Interferons induce an antiviral state by stimulating the expression of hundreds of interferon-stimulated genes (ISGs). These ISGs can constrain viral replication, limit tissue tropism, reduce pathogenicity, and inhibit viral transmission. Since VSIV is used as a backbone for multiple oncolytic and vaccine strategies, understanding how ISGs restrict VSIV not only helps in understanding VSIV-induced pathogenesis but also helps us evaluate and understand the safety and efficacy of VSIV-based therapies. Thus, there is a need to identify and characterize the ISGs that possess anti-VSIV activity. Using arrayed ISG expression screening, we identified TRIM69 as an ISG that potently inhibits VSIV. This inhibition was highly specific as multiple viruses, including influenza A virus, HIV-1, Rift Valley fever virus, and dengue virus, were unaffected by TRIM69. Indeed, just one amino acid substitution in VSIV can govern sensitivity/resistance to TRIM69. Furthermore, TRIM69 is highly divergent in human populations and exhibits signatures of positive selection that are consistent with this gene playing a key role in antiviral immunity. We propose that TRIM69 is an IFN-induced inhibitor of VSIV and speculate that TRIM69 could be important in limiting VSIV pathogenesis and might influence the specificity and/or efficacy of vesiculovirus-based therapies.IMPORTANCE Vesicular stomatitis Indiana virus (VSIV) is a veterinary pathogen that is also used as a backbone for many oncolytic and vaccine strategies. In natural and therapeutic settings, viral infections like VSIV are sensed by the host, and as a result the host cells make proteins that can protect them from viruses. In the case of VSIV, these antiviral proteins constrain viral replication and protect most healthy tissues from virus infection. In order to understand how VSIV causes disease and how healthy tissues are protected from VSIV-based therapies, it is crucial that we identify the proteins that inhibit VSIV. Here, we show that TRIM69 is an antiviral defense that can potently and specifically block VSIV infection.

Draft Genome Sequences of Enterococcus faecalis Strains Isolated from Healthy Japanese Individuals

Enterococcus faecalis is a common commensal of the intestines of humans and other mammals but is also a frequent cause of serious ailments. Here, we report 14 draft genome sequences of strains of Enterococcus faecalis, a normal inhabitant and Gram-positive bacterium that was isolated from 7 healthy Japanese volunteers.

Draft genome sequence of Solanum aethiopicum provides insights into disease resistance, drought tolerance, and the evolution of the genome

BACKGROUND

The African eggplant (Solanum aethiopicum) is a nutritious traditional vegetable used in many African countries, including Uganda and Nigeria. It is thought to have been domesticated in Africa from its wild relative, Solanum anguivi. S. aethiopicum has been routinely used as a source of disease resistance genes for several Solanaceae crops, including Solanum melongena. A lack of genomic resources has meant that breeding of S. aethiopicum has lagged behind other vegetable crops.

RESULTS

We assembled a 1.02-Gb draft genome of S. aethiopicum, which contained predominantly repetitive sequences (78.9%). We annotated 37,681 gene models, including 34,906 protein-coding genes. Expansion of disease resistance genes was observed via 2 rounds of amplification of long terminal repeat retrotransposons, which may have occurred ∼1.25 and 3.5 million years ago, respectively. By resequencing 65 S. aethiopicum and S. anguivi genotypes, 18,614,838 single-nucleotide polymorphisms were identified, of which 34,171 were located within disease resistance genes. Analysis of domestication and demographic history revealed active selection for genes involved in drought tolerance in both "Gilo" and "Shum" groups. A pan-genome of S. aethiopicum was assembled, containing 51,351 protein-coding genes; 7,069 of these genes were missing from the reference genome.

CONCLUSIONS

The genome sequence of S. aethiopicum enhances our understanding of its biotic and abiotic resistance. The single-nucleotide polymorphisms identified are immediately available for use by breeders. The information provided here will accelerate selection and breeding of the African eggplant, as well as other crops within the Solanaceae family.

Deianiraea, an extracellular bacterium associated with the ciliate Paramecium, suggests an alternative scenario for the evolution of Rickettsiales

Rickettsiales are a lineage of obligate intracellular Alphaproteobacteria, encompassing important human pathogens, manipulators of host reproduction, and mutualists. Here we report the discovery of a novel Rickettsiales bacterium associated with Paramecium, displaying a unique extracellular lifestyle, including the ability to replicate outside host cells. Genomic analyses show that the bacterium possesses a higher capability to synthesise amino acids, compared to all investigated Rickettsiales. Considering these observations, phylogenetic and phylogenomic reconstructions, and re-evaluating the different means of interaction of Rickettsiales bacteria with eukaryotic cells, we propose an alternative scenario for the evolution of intracellularity in Rickettsiales. According to our reconstruction, the Rickettsiales ancestor would have been an extracellular and metabolically versatile bacterium, while obligate intracellularity would have evolved later, in parallel and independently, in different sub-lineages. The proposed new scenario could impact on the open debate on the lifestyle of the last common ancestor of mitochondria within Alphaproteobacteria.

Using MARRVEL v1.2 for Bioinformatics Analysis of Human Genes and Variant Pathogenicity

One of the greatest challenges in the bioinformatic analysis of human sequencing data is identifying which variants are pathogenic. Numerous databases and tools have been generated to address this difficulty. However, these many useful data and tools are broadly dispersed, requiring users to search for their variants of interest through human genetic databases, variant function prediction tools, and model organism databases. To solve this problem, we collected data and observed workflows of human geneticists, clinicians, and model organism researchers to carefully select and display valuable information that facilitates the evaluation of whether a variant is likely to be pathogenic. This program, Model organism Aggregated Resources for Rare Variant ExpLoration (MARRVEL) v1.2, allows users to collect relevant data from 27 public sources for further efficient bioinformatic analysis of the pathogenicity of human variants. © 2019 by John Wiley & Sons, Inc.

On the Catalytic Activity of a GT1 Family Glycosyltransferase from Streptomyces venezuelae ISP5230

GT1 family glycosyltansferase, Sv0189, from Streptomyces venezuelae ISP5230 (ATCC 10721) was characterized. The recombinantly produced protein Sv0189 possessed UDP-glycosyltransferase activity. Screening, using an assay employing unnatural nitrophenyl glycosides as activated donors, resulted in the discovery of a broad substrate scope with respect to both acceptor molecules and donor sugars. In addition to polyphenols, including anthraquinones, simple aromatics containing primary or secondary alcohols, a variety of complex natural products and synthetic drugs were glucosylated or xylosylated by Sv0189. Regioselectivity was established through the isolation and characterization of glucosylated products. Sv0189 and homologous proteins are widely distributed among Streptomyces species, and their apparent substrate promiscuity reveals potential for their development as biocatalysts for glycodiversification.

Complete chloroplast genome sequences of four Allium species: comparative and phylogenetic analyses

The genus Allium is one of the largest monocotyledonous genera, containing over 850 species, and most of these species are found in temperate climates of the Northern Hemisphere. Furthermore, as a large number of new Allium species continue to be identified, phylogenetic classification based on morphological characteristics and a few genetic markers will gradually exhibit extremely low discriminatory power. In this study, we present the use of complete chloroplast genome sequences in genome-scale phylogenetic studies of Allium. We sequenced and assembled four Allium chloroplast genomes and retrieved five published chloroplast genomes from GenBank. All nine chloroplast genomes were used for genomic comparison and phylogenetic inference. The chloroplast genomes, ranging from 152,387 bp to 154,482 bp in length, exhibited conservation of genomic structure, and gene organization and order. Subsequently, we observed the expansion of IRs from the basal monocot Acorus americanus to Allium, identified 814 simple sequence repeats, 131 tandem repeats, 154 dispersed repeats and 109 palindromic repeats, and found six highly variable regions. The phylogenetic relationships of the Allium species inferred from the chloroplast genomes obtained high support, indicating that chloroplast genome data will be useful for further resolution of the phylogeny of the genus Allium.

Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species

Amomum villosum is an important medicinal and edible plant with several pharmacologically active volatile oils. However, identifying A. villosum from A. villosum var. xanthioides and A. longiligulare which exhibit similar morphological characteristics to A. villosum, is difficult. The main goal of this study, therefore, is to mine genetic resources and improve molecular methods that could be used to distinguish these species. A total of eight complete chloroplasts (cp) genomes of these Amomum species which were collected from the main producing areas in China were determined to be 163,608–164,069 bp in size. All genomes displayed a typical quadripartite structure with a pair of inverted repeat (IR) regions (29,820–29,959 bp) that separated a large single copy (LSC) region (88,680–88,857 bp) from a small single copy (SSC) region (15,288–15,369 bp). Each genome encodes 113 different genes with 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. More than 150 SSRs were identified in the entire cp genomes of these three species. The Sanger sequencing results based on 32 Amomum samples indicated that five highly divergent regions screened from cp genomes could not be used to distinguish Amomum species. Phylogenetic analysis showed that the cp genomes could not only accurately identify Amomum species, but also provide a solid foundation for the establishment of phylogenetic relationships of Amomum species. The availability of cp genome resources and the comparative analysis is beneficial for species authentication and phylogenetic analysis in Amomum.

The complete chloroplast genome of copper-tolerance plant Elsholtzia splendens

Elsholtzia splendens is a copper-tolerance plant colonized in copper mines in southern China. In this study, we sequenced and de novo assembled the complete chloroplast genome of E. splendens. The complete chloroplast genome is 150,761 bp (37.8% of GC) in length and contains 87 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Phylogenetic analysis revealed that among the 11 Lamiaceae species, Perilla citriodor is the closest relative of E. splendens. The complete chloroplast genome of E. splendens provides a valuable resource for comparative and evolutionary analysis among Lamiaceae species and may be helpful in understanding the molecular mechanism of copper tolerance in E. splendens.

Characterization of the complete chloroplast genome of Dalbergia cultrata (Leguminosae)

Dalbergia cultrata is a Near Threatened species with high ecological and economic values. In this study, its chloroplast genome was assembled using Illumina pair-end sequencing dataset. The chloroplast genome has a quadripartite structure with 156,385 bp in length and contains a pair of 16,392 bp inverted repeat (IR) regions, which were separated by large single copy (LSC: 86,040 bp) region and small single copy (SSC: 37,561 bp) region. A total of 121 genes were annotated, including 77 protein-coding genes (PCGs), 36 tRNAs, and 8 rRNAs. The overall GC content was 36.1%. The phylogenetic analysis revealed that D. cultrata has close relationship to D. hainanensis and D. odorifera. This complete chloroplast genome can be readily used for population genetic studies of D. cultrata.

The complete chloroplast genome sequence of the Siraitia Grosvenorii (Cucurbitaceae)

Siraitia grosvenorii is a famous Chinese plant used in traditional food and medicine with pharmacological effects. The complete chloroplast genome sequence of S. grosvenorii has been determined in this study. The total genome size is 158,834 bp in length and contains a pair of inverted repeats (IRs) of 26,288 bp, which were separated by large single-copy (LSC) and small single-copy (SSC) of 87,702 bp and 18,556 bp length, respectively. A total of 131 genes were predicted including 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis showed that S. grosvenorii belongs to the family Cucurbitaceae. The complete chloroplast genome of S. grosvenorii would play a significant role in the development of molecular markers in plant phylogenetic and population genetic studies.

Sequencing and comparative analysis of three Chlorella genomes provide insights into strain-specific adaptation to wastewater

Microalgal Chlorella has been demonstrated to process wastewater efficiently from piggery industry, yet optimization through genetic engineering of such a bio-treatment is currently challenging, largely due to the limited data and knowledge in genomics. In this study, we first investigated the differential growth rates among three wastewater-processing Chlorella strains: Chlorella sorokiniana BD09, Chlorella sorokiniana BD08 and Chlorella sp. Dachan, and the previously published Chlorella sorokiniana UTEX 1602, showing us that BD09 maintains the best tolerance in synthetic wastewater. We then performed genome sequencing and analysis, resulting in a high-quality assembly for each genome with scaffold N50 > 2 Mb and genomic completeness ≥91%, as well as genome annotation with 9,668, 10,240, 9,821 high-confidence gene models predicted for BD09, BD08, and Dachan, respectively. Comparative genomics study unravels that metabolic pathways, which are involved in nitrogen and phosphorus assimilation, were enriched in the faster-growing strains. We found that gene structural variation and genomic rearrangement might contribute to differential capabilities in wastewater tolerance among the strains, as indicated by gene copy number variation, domain reshuffling of orthologs involved, as well as a ~1 Mb-length chromosomal inversion we observed in BD08 and Dachan. In addition, we speculated that an associated bacterium, Microbacterium chocolatum, which was identified within Dachan, play a possible role in synergizing nutrient removal. Our three newly sequenced Chlorella genomes provide a fundamental foundation to understand the molecular basis of abiotic stress tolerance in wastewater treatment, which is essential for future genetic engineering and strain improvement.

The evolutionarily conserved indolergic receptors of the non-hematophagous elephant mosquito Toxorhynchites amboinensis

The conservation of the mosquito indolergic receptors across the Culicinae and Anophelinae mosquito lineages, which spans 200 million years of evolution, is a testament to the central role of indolic compounds in the biology of these insects. Indole and skatole have been associated with the detection of oviposition sites and animal hosts. To evaluate the potential ecological role of these two compounds, we have used a pharmacological approach to characterize homologs of the indolergic receptors Or2 and Or10 in the non-hematophagous elephant mosquito Toxorhynchites amboinensis. We provide evidence that both receptors are narrowly tuned to indole and skatole like their counterparts from hematophagous mosquitoes. These findings indicate that Toxorhynchites detects indole and skatole in an ecological context to be determined and underscore the importance of understanding the role of these compounds in mosquitoes.

SMCHD1 mutation spectrum for facioscapulohumeral muscular dystrophy type 2 (FSHD2) and Bosma arhinia microphthalmia syndrome (BAMS) reveals disease-specific localisation of variants in the ATPase domain

BACKGROUND

Variants in the Structural Maintenance of Chromosomes flexible Hinge Domain-containing protein 1 (SMCHD1) can cause facioscapulohumeral muscular dystrophy type 2 (FSHD2) and the unrelated Bosma arhinia microphthalmia syndrome (BAMS). In FSHD2, pathogenic variants are found anywhere in SMCHD1 while in BAMS, pathogenic variants are restricted to the extended ATPase domain. Irrespective of the phenotypic outcome, both FSHD2-associated and BAMS-associated SMCHD1 variants result in quantifiable local DNA hypomethylation. We compared FSHD2, BAMS and non-pathogenic SMCHD1 variants to derive genotype-phenotype relationships.

METHODS

Examination of SMCHD1 variants and methylation of the SMCHD1-sensitive FSHD locus DUX4 in 187 FSHD2 families, 41 patients with BAMS and in control individuals. Analysis of variants in a three-dimensional model of the ATPase domain of SMCHD1.

RESULTS

DUX4 methylation analysis is essential to establish pathogenicity of SMCHD1 variants. Although the FSHD2 mutation spectrum includes all types of variants covering the entire SMCHD1 locus, missense variants are significantly enriched in the extended ATPase domain. Identification of recurrent variants suggests disease-specific residues for FSHD2 and in BAMS, consistent with a largely disease-specific localisation of variants in SMCHD1.

CONCLUSIONS

The localisation of missense variants within the ATPase domain of SMCHD1 may contribute to the differences in phenotypic outcome.

Genome wide identification and comparative analysis of glutathione transferases (GST) family genes in Brassica napus

Glutathione transferases (GSTs) are multifunctional enzymes that play important roles in plant development and responses to biotic and abiotic stress. However, a systematic analysis of GST family members in Brassica napus has not yet been reported. In this study, we identified 179 full-length GST genes in B. napus, 44.2% of which are clustered on various chromosomes. In addition, we identified 141 duplicated GST gene pairs in B. napus. Molecular evolutionary analysis showed that speciation and whole-genome triplication played important roles in the divergence of the B. napus GST duplicated genes. Transcriptome analysis of 21 tissues at different developmental stages showed that 47.6% of duplicated GST gene pairs have divergent expression patterns, perhaps due to structural divergence. We constructed a GST gene coexpression network with genes encoding various transcription factors (NAC, MYB, WRKY and bZIP) and identified six modules, including genes expressed during late seed development (after 40 days; BnGSTU19, BnGSTU20 and BnGSTZ1) and in the seed coat (BnGSTF6 and BnGSTF12), stamen and anther (BnGSTF8), root and stem (BnGSTU21), leaves and funiculus, as well as during the late stage of pericarp development (after 40 days; BnGSTU12 and BnGSTF2) and in the radicle during seed germination (BnGSTF14, BnGSTU1, BnGSTU28, and BnGSTZ1). These findings lay the foundation for elucidating the roles of GSTs in B. napus.

Overexpression of Osmosensitive Ca2+-Permeable Channel TMEM63B Promotes Migration in HEK293T Cells

The recent discovery of the osmosensitive calcium (Ca²⁺) channel OSCA has revealed the potential mechanism by which plant cells sense diverse stimuli. Osmosensory transporters and mechanosensitive channels can detect and respond to osmotic shifts that play an important role in active cell homeostasis. Members of the TMEM63 family of proteins are described as the closest homologues of OSCAs. Here, we characterize TMEM63B, a mammalian homologue of OSCAs, recently classified as mechanosensitive. In HEK293T cells, TMEM63B localizes to the plasma membrane and is associated with F-actin. This Ca²⁺-permeable channel specifically induces Ca²⁺ influx across the membrane in response to extracellular Ca²⁺ concentration and hyperosmolarity. In addition, overexpression of TMEM63B in HEK293T cells significantly enhanced cell migration and wound healing. The link between Ca²⁺ osmosensitivity and cell migration might help to establish TMEM63B's pathogenesis, for example, in cancer in which it is frequently overexpressed.

Phylogenomics Provides New Insights into Gains and Losses of Selenoproteins among Archaeplastida

Selenoproteins that contain selenocysteine (Sec) are found in all kingdoms of life. Although they constitute a small proportion of the proteome, selenoproteins play essential roles in many organisms. In photosynthetic eukaryotes, selenoproteins have been found in algae but are missing in land plants (embryophytes). In this study, we explored the evolutionary dynamics of Sec incorporation by conveying a genomic search for the Sec machinery and selenoproteins across Archaeplastida. We identified a complete Sec machinery and variable sizes of selenoproteomes in the main algal lineages. However, the entire Sec machinery was missing in the Bangiophyceae-Florideophyceae clade (BV) of Rhodoplantae (red algae) and only partial machinery was found in three species of Archaeplastida, indicating parallel loss of Sec incorporation in different groups of algae. Further analysis of genome and transcriptome data suggests that all major lineages of streptophyte algae display a complete Sec machinery, although the number of selenoproteins is low in this group, especially in subaerial taxa. We conclude that selenoproteins tend to be lost in Archaeplastida upon adaptation to a subaerial or acidic environment. The high number of redox-active selenoproteins found in some bloom-forming marine microalgae may be related to defense against viral infections. Some of the selenoproteins in these organisms may have been gained by horizontal gene transfer from bacteria.

Genetic Characterization of a Recombinant Myxoma Virus in the Iberian Hare (Lepus granatensis)

Myxomatosis is a lethal disease in wild European and domestic rabbits (Oryctolagus cuniculus), which is caused by a Myxoma virus (MYXV) infection-a leporipoxvirus that is found naturally in some Sylvilagus rabbit species in South America and California. The introduction of MYXV into feral European rabbit populations of Australia and Europe, in the early 1950s, demonstrated the best-documented field example of host-virus coevolution, following a cross-species transmission. Recently, a new cross-species jump of MYXV has been suggested in both Great Britain and Spain, where European brown hares (Lepus europaeus) and Iberian hares (Lepus granatensis) were found dead with lesions consistent with those observed in myxomatosis. To investigate the possibility of a new cross-species transmission event by MYXV, tissue samples collected from a wild Iberian hare found dead in Spain (Toledo region) were analyzed and deep sequenced. Our results reported a new MYXV isolate (MYXV Toledo) in the tissues of this species. The genome of this new virus was found to encode three disruptive genes (M009L, M036L, and M152R) and a novel ~2.8 kb recombinant region, which resulted from an insertion of four novel poxviral genes towards the 3' end of the negative strand of its genome. From the open reading frames inserted into the MYXV Toledo virus, a new orthologue of a poxvirus host range gene family member was identified, which was related to the MYXV gene M064R. Overall, we confirmed the identity of a new MYXV isolate in Iberian hares, which, we hypothesized, was able to more effectively counteract the host defenses in hares and start an infectious process in this new host.

Frischella japonica sp. nov., an anaerobic member of the Orbales in the Gammaproteobacteria, isolated from the gut of the eastern honey bee, Apis cerana japonica Fabricius

The gut of honey bees is characterized by a stable and relatively simple community of bacteria, consisting of seven to ten phylotypes. Two closely related honey bees, Apis mellifera (western honey bee) and Apis cerana (eastern honey bee), show a largely comparable occurrence of those phylotypes, but a distinct set of bacterial species and strains within each bee species. Here, we describe the isolation and characterization of Ac13T, a new species within the rare proteobacterial genus Frischella from A. cerana japonica Fabricius. Description of Ac13T as a new species is supported by low identity of the 16S rRNA gene sequence (97.2 %), of the average nucleotide identity based on orthologous genes (77.5 %) and digital DNA-DNA hybridization relatedness (24.7 %) to the next but far related type strain Frischella perrara PEB0191T, isolated from A. mellifera. Cells of Ac13T are mesophilic and have a mean length of 2-4 µm and a width of 0.5 µm. Optimal growth was achieved in anoxic conditions, whereas growth was not observed in oxic conditions and strongly reduced in microaerophilic environment. Strain Ac13T shares several features with other members of the Orbaceae, such as the major fatty acid profile, the respiratory quinone type and relatively low DNA G+C content, in accordance with its evolutionary relationship. Unlike F. perrara, strain Ac13T is susceptible to a broad range of antibiotics, which could be indicative for an antibiotic-free A. cerana bee keeping. In conclusion, we propose strain Ac13T as a novel species for which we propose the name Frischella japonica sp. nov. with the type strain Ac13T (=NCIMB 15259=JCM 34075).

Vagococcus bubulae sp. nov., isolated from ground beef, and Vagococcus vulneris sp. nov., isolated from a human foot wound

Two unusual catalase-negative, Gram-stain-positive, Vagococcus-like isolates that were referred to the CDC Streptococcus Laboratory for identification are described. Strain SS1994^T was isolated from ground beef and strain SS1995^T was isolated from a human foot wound. Comparative 16S rRNA gene sequence analysis of isolates SS1994^T and SS1995^T against Vagococcus type strain sequences supported their inclusion in the genus Vagococcus. Strain SS1994^T showed high sequence similarity (>97.0 %) to the two most recently proposed species, Vagococcus martis (99.2 %) and Vagococcus teuberi (99.0 %) followed by Vagococcus penaei (98.8 %), strain SS1995^T (98.6 %), Vagococcus carniphilus (98.0 %), Vagococcus acidifermentans (98.0 %) and Vagococcus fluvialis (97.9 %). The 16S rRNA gene sequence of strain SS1995^T was most similar to V. penaei (99.1 %), followed by SS1994^T (98.6 %), V. martis (98.4 %), V. teuberi (98.1 %), V. acidifermentans (97.8 %), and both V. carniphilus and V. fluvialis (97.5 %). A polyphasic taxonomic study using conventional biochemical and the rapid ID 32 STREP system, MALDI-TOF MS, cell fatty acid analysis, pairwise sequence comparisons of the 16S rRNA, rpoA, rpoB, pheS and groL genes, and comparative core and whole genome sequence analyses revealed that strains SS1994^T and SS1995^T were two novel Vagococcus species. The novel taxonomic status of the two isolates was confirmed with core genome phylogeny, average nucleotide identity <84 % and in silico DNA-DNA hybridization <28 % to any other Vagococcus species. The names Vagococcusbubulae SS1994^T=(CCUG 70831^T=LMG 30164^T) and Vagococcusvulneris SS1995^T=(CCUG 70832^T=LMG 30165^T) are proposed.

Profiling of chromatin accessibility and identification of general cis-regulatory mechanisms that control two ocular lens differentiation pathways

BACKGROUND

Promoters and enhancers are cis-regulatory DNA sequences that control specificity and quantity of transcription. Both are rich on clusters of cis-acting sites that interact with sequence-specific DNA-binding transcription factors (TFs). At the level of chromatin, these regions display increased nuclease sensitivity, reduced nucleosome density, including nucleosome-free regions, and specific combinations of posttranslational modifications of core histone proteins. Together, "open" and "closed" chromatins represent transcriptionally active and repressed states of individual genes, respectively. Cellular differentiation is marked by changes in local chromatin structure. Lens morphogenesis, regulated by TF Pax6, includes differentiation of epithelial precursor cells into lens fibers in parallel with differentiation of epithelial precursors into the mature lens epithelium.

RESULTS

Using ATAC-seq, we investigated dynamics of chromatin changes during mouse lens fibers and epithelium differentiation. Tissue-specific features of these processes are demonstrated via comparative studies of embryonic stem cells, forebrain, and liver chromatins. Unbiased analysis reveals cis-regulatory logic of lens differentiation through known (e.g., AP-1, Ets, Hsf4, Maf, and Pax6 sites) and novel (e.g., CTCF, Tead, and NF1) motifs. Twenty-six DNA-binding TFs, recognizing these cis-motifs, are markedly up-regulated in differentiating lens fibers. As specific examples, our ATAC-seq data uncovered both the regulatory regions and TF binding motifs in Foxe3, Prox1, and Mip loci that are consistent with previous, though incomplete, experimental data. A cross-examination of Pax6 binding with ATAC-seq data demonstrated that Pax6 bound to both open (H3K27ac and P300-enriched) and closed chromatin domains in lens and forebrain.

CONCLUSIONS

Our study has generated the first lens chromatin accessibility maps that support a general model of stage-specific chromatin changes associated with transcriptional activities of batteries of genes required for lens fiber cell formation. Analysis of active (or open) promoters and enhancers reveals important cis-DNA motifs that establish the molecular foundation for temporally and spatially regulated gene expression in lens. Together, our data and models open new avenues for the field to conduct mechanistic studies of transcriptional control regions, reconstruction of gene regulatory networks that govern lens morphogenesis, and identification of cataract-causing mutations in noncoding sequences.

Effect of CO2 Concentration on Uptake and Assimilation of Inorganic Carbon in the Extreme Acidophile Acidithiobacillus ferrooxidans

This study was motivated by surprising gaps in the current knowledge of microbial inorganic carbon (Ci) uptake and assimilation at acidic pH values (pH < 3). Particularly striking is the limited understanding of the differences between Ci uptake mechanisms in acidic versus circumneutral environments where the Ci predominantly occurs either as a dissolved gas (CO2) or as bicarbonate (HCO3 -), respectively. In order to gain initial traction on the problem, the relative abundance of transcripts encoding proteins involved in Ci uptake and assimilation was studied in the autotrophic, polyextreme acidophile Acidithiobacillus ferrooxidans whose optimum pH for growth is 2.5 using ferrous iron as an energy source, although they are able to grow at pH 5 when using sulfur as an energy source. The relative abundance of transcripts of five operons (cbb1-5) and one gene cluster (can-sulP) was monitored by RT-qPCR and, in selected cases, at the protein level by Western blotting, when cells were grown under different regimens of CO2 concentration in elemental sulfur. Of particular note was the absence of a classical bicarbonate uptake system in A. ferrooxidans. However, bioinformatic approaches predict that sulP, previously annotated as a sulfate transporter, is a novel type of bicarbonate transporter. A conceptual model of CO2 fixation was constructed from combined bioinformatic and experimental approaches that suggests strategies for providing ecological flexibility under changing concentrations of CO2 and provides a portal to elucidating Ci uptake and regulation in acidic conditions. The results could advance the understanding of industrial bioleaching processes to recover metals such as copper at acidic pH. In addition, they may also shed light on how chemolithoautotrophic acidophiles influence the nutrient and energy balance in naturally occurring low pH environments.

Characterization of Two Pseudomonas aeruginosa Viruses vB_PaeM_SCUT-S1 and vB_PaeM_SCUT-S2

The sophisticated antibiotic resistance mechanism of Pseudomonas aeruginosa has urged the development of alternative antibacterial strategies. Phage therapy has been proven successful for the treatment of multidrug-resistant infections. In this study, we reported two virulent P. aeruginosa phages, vB_PaeM_SCUT-S1 (S1) and vB_PaeM_SCUT-S2 (S2), which were characterized at morphological, genomic, and proteomic levels. Phages S1 and S2 were assigned to the Myoviridae family. The genome sequencing showed that the genome size of Phage S1 was 66,046 bp and that of Phage S2 was 94,434 bp. The phylogenetic tree indicated that the two phages were distantly related to each other and were classified in the genera Pbunavirus and Pakpunavirus respectively. Thirty-one proteins were identified for each phage by mass spectrometry and were used to substantiate the function of the predicted coding genes. The two phages inhibited the growth of P. aeruginosa strain PAO1 at low multiplicity of infection levels and had good performance both on preventing biofilm formation and eradicating preformed biofilms. They were also stable over a wide range of temperature and pH values, supporting their potential use in the treatment of P. aeruginosa infections.

In-depth transcriptome characterization uncovers distinct gene family expansions for Cupressus gigantea important to this long-lived species' adaptability to environmental cues

Background

Cupressus gigantea, a rare and endangered tree species with remarkable medicinal value, is endemic to the Tibetan Plateau. Yet, little is known about the underlying genetics of the unique ecological adaptability of this extremely long-lived conifer with a large genome size. Here, we present its first de novo and multi-tissue transcriptome in-depth characterization.

Results

We performed Illumina paired-end sequencing and RNA libraries assembly derived from terminal buds, male and female strobili, biennial leaves, and cambium tissues taken from adult C. gigantea. In total, large-scale high-quality reads were assembled into 101,092 unigenes, with an average sequence length of 1029 bp, and 6848 unigenes (6.77%) were mapped against the KEGG databases to identify 292 pathways. A core set of 41,373 genes belonging to 2412 orthologous gene families shared between C. gigantea and nine other plants was revealed. In addition, we identified 2515 small to larger-size gene families containing in total 9223 genes specific to C. gigantea, and enriched for gene ontologies relating to biotic interactions. We identified an important terpene synthases gene family expansion with its 121 putative members.

Conclusions

This study presents the first comprehensive transcriptome characterization of C. gigantea. Our results will facilitate functional genomic studies to support genetic improvement and conservation programs for this endangered conifer.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5584-6) contains supplementary material, which is available to authorized users.

Whole Genome Sequencing of Mycobacterium tuberculosis Clinical Isolates From India Reveals Genetic Heterogeneity and Region-Specific Variations That Might Affect Drug Susceptibility

Whole genome sequencing (WGS) of Mycobacterium tuberculosis has been constructive in understanding its evolution, genetic diversity and the mechanisms involved in drug resistance. A large number of sequencing efforts from across the globe have revealed genetic diversity among clinical isolates and the genetic determinants for their resistance to anti-tubercular drugs. Considering the high TB burden in India, the availability of WGS studies is limited. Here we present, WGS results of 200 clinical isolates of M. tuberculosis from North India which are categorized as sensitive to first-line drugs, mono-resistant, multi-drug resistant and pre-extensively drug resistant isolates. WGS revealed that 20% of the isolates were co-infected with M. tuberculosis and non-tuberculous mycobacteria species. We identified 12,802 novel genetic variations in M. tuberculosis isolates including 343 novel SNVs in 38 genes which are known to be associated with drug resistance and are not currently used in the diagnostic kits for detection of drug resistant TB. We also identified M. tuberculosis lineage 3 to be predominant in the northern region of India. Additionally, several novel SNVs, which may potentially confer drug resistance were found to be enriched in the drug resistant isolates sampled. This study highlights the significance of employing WGS in diagnosis and for monitoring further development of MDR-TB strains.

Optimization of qRT-PCR assay for zika virus detection in human serum and urine

Zika Virus (ZIKV) is a mosquito-borne flavivirus that the World Health Organization (WHO) declared a global concern due to the severity of infection. This study focuses on determining the level of detection of ZIKV RNA in human serum and urine. Known amounts of Zika virus were added to uninfected human serum and urine samples. Different reverse transcriptases were compared to select the optimal enzyme for this application. Zika RNA in these samples was then quantified with qRT-PCR to determine the lower limit of detection in these fluids and to construct a standard curve. Student's t-test of paired samples was used in order to identify statistical differences. The SuperScript III enzyme was able to produce more ZIKV cDNA when compared to PrimeScript. Zika virus RNA was found to be detectable at lower levels (2.5 PFU/mL) in urine than in serum (250 PFU/mL) when using SuperScript III. This study demonstrates how the selection of both the human clinical specimen, and the reverse transcriptase enzyme involved in the molecular detection of ZIKV by quantitative real-time polymerase chain reaction (qRT-PCR), play an important role in enabling improved detection of the virus.

Culturable endophytic fungal communities associated with plants in organic and conventional farming systems and their effects on plant growth

As compared to organic farming system, conventional farming system relies on higher inputs of synthetic agrochemicals, which may reduce the abundance, diversity, and beneficial effects of plant endophytic fungal communities. This study compares the diversity and abundance of culturable endophytic fungal communities associated with four plant species –corn, tomato, pepper, and watermelon grown in separate organic and conventional fields. In all, 740 fungal isolates were identified, of which 550 were from the organic fields and 190 from the conventional ones. These fungal isolates were grouped into eight orders and 22 species, with the two most abundant species being Trichoderma sp. and Pichia guilliermondi. The fungal species diversity and abundance were both significantly higher in the organic than in the conventional fields. All the isolated endophytic fungi improved tomato plants’ shoot growth and biomass significantly, as compared with the water control. Six fungal isolates also exhibited activity that enhanced tomato fruit yields. These results suggest that these endophytic fungi might be a considerable boost to sustainable agricultural production, while also reducing the agricultural application of chemicals and thus benefiting the environment and human health.

Supersensitive Odorant Receptor Underscores Pleiotropic Roles of Indoles in Mosquito Ecology

Mosquitoes exhibit highly diverse and fast evolving odorant receptors (ORs). The indole-sensitive OR gene clade, comprised of Or2 and Or10 is a notable exception on account of its conservation in both mosquito subfamilies. This group of paralogous genes exhibits a complex developmental expression pattern in Aedes aegypti: AaegOr2 is expressed in both adults and larvae, AaegOr10 is adult-specific and a third member named AaegOr9 is larva-specific. OR2 and OR10 have been deorphanized and are selectively activated by indole and skatole, respectively. Using the two-electrode voltage clamp of Xenopus oocytes expressing Ae. aegypti ORs, we show that AaegOR9 is supersensitive and narrowly tuned to skatole. Our findings suggest that Ae. aegypti has evolved two distinct molecular strategies to detect skatole in aquatic and terrestrial environments, highlighting the central ecological roles of indolic compounds in the evolutionary and life histories of these insects.

Eriobotrya Belongs to Rhaphiolepis (Maleae, Rosaceae): Evidence From Chloroplast Genome and Nuclear Ribosomal DNA Data

The Eriobotrya-Rhaphiolepis (ER) clade consists of about 46 species distributed in East and Southeast Asia. Although Eriobotrya and Rhaphiolepis have been supported to form a clade, the monophyly of Eriobotrya and Rhaphiolepis at the genus level has never been well tested and their phylogenetic positions in Maleae still remain uncertain. This study aims to reconstruct a robust phylogeny of the ER clade in the framework of Maleae with a broad taxon sampling and clarify the phylogenetic relationship between Eriobotrya and Rhaphiolepis. This study employed sequences of the whole plastome (WP) and entire nuclear ribosomal DNA (nrDNA) repeats assembled from the genome skimming approach and included 83 samples representing 76 species in 32 genera of Rosaceae, especially Maleae. The Maximum Likelihood (ML) and Bayesian Analysis (BI) based on three datasets, i.e., WP, coding sequences of plastome (CDS), and nrDNA, strongly supported the paraphyly of Eriobotrya, within which Rhaphiolepis was nested. Our plastid tree supported the sister relationship between the ER clade and Heteromeles, and the nrDNA tree, however, did not resolve the phylogenetic placement of the ER clade in Maleae. Strong incongruence between the plastid and the nuclear trees is most likely explained by hybridization events, which may have played an important role in the evolutionary history of the ER clade. Molecular, morphological, and geographic evidence all supports the merge of Eriobotrya with Rhaphiolepis, which has the nomenclatural priority. We herein transferred 36 taxa of Eriobotrya to Rhaphiolepis. We also proposed a new name, Rhaphiolepis loquata B.B.Liu & J.Wen, for the economically important loquat, as the specific epithet "japonica" was pre-occupied in Rhaphiolepis.