Incorporation of Non-Canonical Amino Acids into Proteins by Global Reassignment of Sense Codons

Non-canonical amino acids are finding increasing use in basic and applied research. Proteins that evolved naturally for biological function did so by exploiting the chemistries of the canonical amino acids; however, when proteins are repurposed for biomedical and pharmacological applications, they are often subject to conditions different from those characteristic of their original biological environments. Non-canonical amino acids can impart properties that are inaccessible within canonical protein sequence space, and can thereby lead to improved or new functionality. We describe simple methods for global replacement of canonical amino acids by their non-canonical counterparts in recombinant proteins made in high yield in bacterial expression hosts. These methods can be used to engineer both chemical and physical properties of recombinant proteins.

FnCpf1-Mediated Targeted Mutagenesis in Plants

Sequence-specific nucleases (SSNs) are nowadays fundamental tools to generate mutants that impaired in genes of interest. The bioactive molecules screened in the chemical genomics studies affect specific physiological process by disrupting the function of its target protein(s). Mutation analysis of the gene(s) of target protein(s) of the screened chemical is necessary to resolve how the chemical works in plants. Clustered regularly interspersed short palindromic repeats (CRISPR) from Prevotella and Francisella 1 (Cpf1) are newly characterized RNA-directed endonuclease. Several papers have shown clearly that Cpf1 could be a versatile SSN in plant genome engineering. Cfp1 from Francisella novicida (FnCpf1) recognizes TTN as its protospacer adjacent motif (PAM). FnCpf1 utilizes a shorter PAM compared to other known Cpf1s such as AsCpf1 or LbCpf1, which use TTTN as PAM. Since PAM length can be a limiting factor in target selection, this feature of FnCpf1 is practical for targeted mutagenesis experiments. The application of FnCpf1-mediated targeted mutagenesis to the chemical genomics could accelerate to figure out the mechanism of action of screened chemicals. Here, we describe procedures for targeted mutagenesis in rice and tobacco using FnCpf1.

Generation of Stable Amber Suppression Cell Lines

Noncanonical amino acid mutagenesis via amber suppression provides the means to tailor proteins inside living cells. A wide range of noncanonical amino acids have been incorporated using the Methanococcus pyrrolysyl-tRNA synthetase/tRNACUA (PylRS/PylT) in mammalian cell systems in proof of principle experiments, for (1) minimal genetically encoded fluorescence or affinity tagging, (2) photo-control of enzymes, (3) genetically encoded posttranslational protein modifications. We have developed a general and efficient method to genomically integrate the PylRS/PylT amber suppression machinery using PiggyBac-mediated transposition. A general protocol for the generation of stable amber suppression cell lines is described here. Using the modular plasmid system, homogenous and highly efficient amber suppression in a wide range of cell lines can be achieved.

A high-resolution genetic map of the cereal crown rot pathogen Fusarium pseudograminearum provides a near-complete genome assembly

Fusarium pseudograminearum is an important pathogen of wheat and barley, particularly in semi-arid environments. Previous genome assemblies for this organism were based entirely on short read data and are highly fragmented. In this work, a genetic map of F. pseudograminearum has been constructed for the first time based on a mapping population of 178 individuals. The genetic map, together with long read scaffolding of a short read-based genome assembly, was used to give a near-complete assembly of the four F. pseudograminearum chromosomes. Large regions of synteny between F. pseudograminearum and F. graminearum, the related pathogen that is the primary causal agent of cereal head blight disease, were previously proposed in the core conserved genome, but the construction of a genetic map to order and orient contigs is critical to the validation of synteny and the placing of species-specific regions. Indeed, our comparative analyses of the genomes of these two related pathogens suggest that rearrangements in the F. pseudograminearum genome have occurred in the chromosome ends. One of these rearrangements includes the transposition of an entire gene cluster involved in the detoxification of the benzoxazolinone (BOA) class of plant phytoalexins. This work provides an important genomic and genetic resource for F. pseudograminearum, which is less well characterized than F. graminearum. In addition, this study provides new insights into a better understanding of the sexual reproduction process in F. pseudograminearum, which informs us of the potential of this pathogen to evolve.

Genetic Incorporation of Unnatural Amino Acids into Proteins of Interest in Streptomyces venezuelae ATCC 15439

Site-specific, genetic incorporation of unnatural amino acids (UAAs) into proteins in living cells using engineered orthogonal aminoacyl-tRNA synthetase (aaRS)/tRNA pairs is a powerful tool for studying and manipulating protein structure and function. To date, UAA incorporation systems have been developed for several bacterial and eukaryotic model hosts. Due to the importance of Streptomyces as prolific producers of bioactive natural products and as model hosts for natural product biosynthesis and bioengineering studies, we have developed systems for the incorporation of the UAAs p-iodo-L-phenylalanine (pIPhe) and p-azido-L-phenylalanine (pAzPhe) into green fluorescent protein (GFP) in Streptomyces venezuelae ATCC 15439. Here, we describe the procedure for using this system to site-specifically incorporate pIPhe or pAzPhe into proteins of interest in S. venezuelae. The modular design of plasmids harboring UAA incorporation systems enables use of other aaRS or aaRS/tRNA pairs for the incorporation of other UAAs; and the vector backbone used allows the system to be transferred to diverse Streptomyces species via both protoplast transformation and conjugation.

Genomic Tagging of AGO1 Using CRISPR/Cas9-Mediated Homologous Recombination

Tagging of genes at the endogenous loci is a powerful strategy for the analysis of protein function. We have developed a homologous recombination-based approach for inserting epitope tag into Drosophila AGO1 locus by employing the CRISPR/Cas9 technology. The methodology involves co-expression of sgRNA (containing 20-nucleotide AGO1 targeting sequence) and Cas9 protein, together with a donor template that has HA-AGO1 cassette flanked by sequences homologous to the AGO1 locus. The integration is efficient and readily monitored by immunostaining of the transgenic cell line. This method facilitates rapid generation of stable cell lines and allows insertion of any tag sequence into endogenous loci, thus accelerating characterization of the tagged proteins.

Genetic Encoding of Unnatural Amino Acids in C. elegans

Site-specific incorporation of unnatural amino acids (UAAs) has greatly expanded the toolkit available to study biological phenomena in single cells. However, to address questions involving complex cellular interactions such as development, ageing, and the functions of the nervous system it is often necessary to use multicellular model organisms. The nematode Caenorhabditis elegans was the first organism to have its genetic code expanded. Due to its small size, ease of cultivation, and excellent UAA incorporation efficiency, C. elegans makes an ideal model organism to apply UAAs as tools to investigate the functioning of multicellular systems.Here, we describe methods to generate transgenic C. elegans capable of UAA incorporation, as well as how to deliver unnatural amino acids and test incorporation. Furthermore, we describe methods to uncage photosensitive unnatural amino acid derivatives.

Application of the Cre/lox System to Construct Auxotrophic Markers for Quantitative Genetic Analyses in Fusarium graminearum

The bacteriophage P1 Cre/lox system has been utilized in diverse fungi for marker recycling and exchange, generation of targeted chromosome translocations, and targeted deletion of interstitial chromosome segments. Here we show the application of this tool in the wheat and maize pathogen, Fusarium graminearum. We explored three different ways to introduce Cre into strains with floxed genes, namely transformation with an episomal or integrative plasmid (pLC28), fusion of protoplasts of strains carrying floxed genes with strains expressing Cre by forcing heterokaryons, and crosses between strains with floxed genes and strains expressing Cre to isolate progeny in which the target genes had been deleted during the cross. We used this system for the construction of strains bearing auxotrophic markers that were generated by gene replacement with positively selectable markers followed by Cre-mediated marker excision. In addition, updated protocols for transformation and crosses for F. graminearum are provided. In combination, strains and tools developed here add to the arsenal of methods that can be used to carry out molecular genetics with F. graminearum.

Family-Free Genome Comparison

The comparison of genome structures across distinct species offers valuable insights into the species' phylogeny, genome organization, and gene associations. In this chapter, we review the family-free genome comparison tool FFGC which provides several methods for gene order analyses that do not require prior knowledge of evolutionary relationships between the genes across the studied genomes. Moreover, the tool features a complete workflow for genome comparison, requiring nothing but annotated genome sequences as input.

Protocol of Phytophthora capsici Transformation Using the CRISPR-Cas9 System

Phytophthora capsici is an important plant pathogen, which causes significant economic losses on multiple vegetable crops worldwide. It is an ideal model pathogen to study the role of important genes, plant-pathogen interactions, and fungicide resistance mechanisms etc. due to its wide range of hosts and genetic diversity. A more efficient gene editing tool is required to do these studies. Here, we describe a detailed experimental procedure using the CRISPR-Cas9 system to edit genes of interest in P. capsici, which has been proven to be an accurate and efficient gene editing method in P. capsici.

Impact of plasmid architecture on stability and yEGFP3 reporter gene expression in a set of isomeric multicopy vectors in yeast

Multicopy episomal plasmids in yeast, used whenever elevated levels of foreign or homologous gene expression are necessary, are known to be less stable compared to the endogenous 2-μm plasmid they are based on, at least without selective pressure. Considering that rich medium favors growth rate and, simultaneously, is less expensive than selective medium, enhancing stability in non-selective medium is extremely desirable. In this study, we changed the architecture of a multicopy model expression plasmid, creating six isoforms (same size, same DNA content but different positions and orientations of the expression block) and studied mitotic stability, copy number, as well as reporter yEGFP3 expression between isoforms. With one isoform being significantly more stable than the others and another one exhibiting elevated plasmid copy numbers in rich medium, we show that consideration of the arrangement of the plasmid elements might be crucial for productivity employing Saccharomyces cerevisiae as a host. We strongly believe that the ideal architecture has to be assessed for each case and assembly strategy has to begin by evaluating the stability of the vector backbone before insertion of the desired gene. For the plasmid set studied, yEGFP3 reporter production depends more on mitotic stability than on elevated plasmid copy numbers in a small number of cells retaining the plasmid under non-selective conditions.

The mitochondrial genome of the wolfberry fruit fly, Neoceratitis asiatica (Becker) (Diptera: Tephritidae) and the phylogeny of Neoceratitis Hendel genus

Neoceratitis asiatica (Becker) (Diptera: Tephritidae) is one of the most important fruit pestsof wolfberry which is a traditional Chinese medicinal herb. We characterized the complete mitochondrial genome of N. asiatica and described its organization in this study. This mitogenome had a total length of 15,481 bp, consisting of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and a non-coding region (A + T-rich control region). The overall base composition of N. asiatica in descending order was 40.6% A, 8.5% G, 38.4% T and 12.6% C. The phylogenetic relationships shows that Ceratitis capitata and N. asiatica may be sister taxa. This is the first report of the complete mitochondrial genome of a member of the Neoceratitis Genus and the complete mitochondrial genome sequence may provide useful information for phylogenetic analysis and studies between the genera Ceratitis and Neoceratitis.

miR-25 Tough Decoy Enhances Cardiac Function in Heart Failure

MicroRNAs are promising therapeutic targets, because their inhibition has the potential to normalize gene expression in diseased states. Recently, our group found that miR-25 is a key SERCA2a regulating microRNA, and we showed that multiple injections of antagomirs against miR-25 enhance cardiac contractility and function through SERCA2a restoration in a murine heart failure model. However, for clinical application, a more stable suppressor of miR-25 would be desirable. Tough Decoy (TuD) inhibitors are emerging as a highly effective method for microRNA inhibition due to their resistance to endonucleolytic degradation, high miRNA binding affinity, and efficient delivery. We generated a miR-25 TuD inhibitor and subcloned it into a cardiotropic AAV9 vector to evaluate its efficacy. The AAV9 TuD showed selective inhibition of miR-25 in vitro cardiomyoblast culture. In vivo, AAV9-miR-25 TuD delivered to the murine pressure-overload heart failure model selectively decreased expression of miR-25, increased levels of SERCA2a protein, and ameliorated cardiac dysfunction and fibrosis. Our data indicate that miR-25 TuD is an effective long-term suppressor of miR-25 and a promising therapeutic candidate to treat heart failure.

Characterization of the complete mitochondrial genomes of Maiestas dorsalis and Japananus hyalinus (Hemiptera: Cicadellidae) and comparison with other Membracoidea

Only six mitochondrial genomes (mitogenomes) have been previously published for Cicadellidae, the largest family of Hemiptera. This study provides complete, annotated mitogenomes of two additional cicadellid, species Maiestas dorsalis and Japananus hyalinus, and the first comparative mitogenome analysis across the superfamily Membracoidea. The mitogenomes of both sequenced species are similar to those of other studied hemipteran mitogenomes in organization and the lengths are 15,352 and 15,364 bp with an A + T content of 78.7% and 76.6%, respectively. In M. dorsalis, all sequenced genes are arranged in the putative ancestral insect gene arrangement, while the tRNA cluster trnW-trnC-trnY is rearranged to trnY-trnW-trnC in J. hyalinus, the first reported gene rearrangement in Membracoidea. Phylogenetic analyses of the 11 available membracoid mitogenomes and outgroups representing the other two cicadomorphan superfamilies supported the monophyly of Membracoidea, and indicated that treehoppers are a derived lineage of leafhoppers. ML and BI analyses yielded topologies that were congruent except for relationships among included representatives of subfamily Deltocephalinae. Exclusion of third codon positions of PCGs improved some node support values in ML analyses.

The proposed new species, cacao red vein virus, and three previously recognized badnavirus species are associated with cacao swollen shoot disease

BACKGROUND

Cacao swollen shoot virus (CSSV), Cacao swollen shoot CD virus (CSSCDV), and Cacao swollen shoot Togo A virus (CSSTAV) cause cacao swollen shoot disease (CSSD) in West Africa. During 2000-2003, leaf and shoot-swelling symptoms and rapid tree death were observed in cacao in Cote d'Ivoire and Ghana. Molecular tests showed positive infection in only ~50-60% of symptomatic trees, suggesting the possible emergence of an unknown badnavirus.

METHODS

The DNA virome was determined from symptomatic cacao samples using Illumina-Hi Seq, and sequence accuracy was verified by Sanger sequencing. The resultant 14, and seven previously known, full-length badnaviral genomic and RT-RNase H sequences were analyzed by pairwise distance analysis to resolve species relationships, and by Maximum likelihood (ML) to reconstruct phylogenetic relationships. The viral coding and non-coding sequences, genome organization, and predicted conserved protein domains (CPDs) were identified and characterized at the species level.

RESULTS

The 21 CSSD-badnaviral genomes and RT-RNase H sequences shared 70-100% and 72-100% identity, respectively. The RT-RNase H analysis predicted four species, based on an ≥80% species cutoff. The ML genome sequence tree resolved three well-supported clades, with ≥70% bootstrap, whereas, the RT-RNase H phylogeny was poorly resolved, however, both trees grouped CSSD isolates within one large clade, including the newly discovered Cacao red vein virus (CRVV) proposed species. The genome arrangement of the four species consists of four, five, or six predicted open reading frames (ORFs), and the CPDs have similar architectures. By comparison, two New World cacao-infecting badnaviruses encode four ORFs, and harbor CPDs like the West African species.

CONCLUSIONS

Three previously recognized West African cacao-infecting badnaviral species were identified, and a fourth, previously unidentified species, CRVV, is described for the first time. The CRVV is a suspect causal agent of the rapid decline phenotype, however Koch's Postulates have not been proven. To reconcile viral evolutionary with epidemiology considerations, more detailed information about CSSD-genomic variability is essential. Also, the functional basis for the multiple genome arrangements and subtly distinct CPD architectures among cacao-infecting badnaviruses is poorly understood. New knowledge about functional relationships may help explain the diverse symptomatologies observed in affected cacao trees.

Phylogeny analysis from gene-order data with massive duplications

BACKGROUND

Gene order changes, under rearrangements, insertions, deletions and duplications, have been used as a new type of data source for phylogenetic reconstruction. Because these changes are rare compared to sequence mutations, they allow the inference of phylogeny further back in evolutionary time. There exist many computational methods for the reconstruction of gene-order phylogenies, including widely used maximum parsimonious methods and maximum likelihood methods. However, both methods face challenges in handling large genomes with many duplicated genes, especially in the presence of whole genome duplication.

METHODS

In this paper, we present three simple yet powerful methods based on maximum-likelihood (ML) approaches that encode multiplicities of both gene adjacency and gene content information for phylogenetic reconstruction.

RESULTS

Extensive experiments on simulated data sets show that our new method achieves the most accurate phylogenies compared to existing approaches. We also evaluate our method on real whole-genome data from eleven mammals. The package is publicly accessible at http://www.geneorder.org .

CONCLUSIONS

Our new encoding schemes successfully incorporate the multiplicity information of gene adjacencies and gene content into an ML framework, and show promising results in reconstruct phylogenies for whole-genome data in the presence of massive duplications.

Sequencing of the complete mitochondrial genome of a fish-parasitic flatworm Paratetraonchoides inermis (Platyhelminthes: Monogenea): tRNA gene arrangement reshuffling and implications for phylogeny

BACKGROUND

Paratetraonchoides inermis (Monogenea: Tetraonchoididae) is a flatworm parasitising the gills of uranoscopid fishes. Its morphological characteristics are ambiguous, and molecular data have never been used to study its phylogenetic relationships, which makes its taxonomic classification controversial. Also, several decades of unsuccessful attempts to resolve the relationships within the Monogenea present a strong indication that morphological datasets may not be robust enough to be used to infer evolutionary histories. As the use of molecular data is currently severely limited by their scarcity, we have sequenced and characterized the complete mitochondrial (mt) genome of P. inermis. To investigate its phylogenetic position, we performed phylogenetic analyses using Bayesian inference and maximum likelihood approaches using concatenated amino acid sequences of all 12 protein-coding genes on a dataset containing all available monogenean mt genomes.

RESULTS

The circular mt genome of P. inermis (14,654 bp) contains the standard 36 genes: 22 tRNAs, two rRNAs, 12 protein-encoding genes (PCGs; Atp8 is missing) and a major non-coding region (mNCR). All genes are transcribed from the same strand. The A + T content of the whole genome (82.6%), as well as its elements, is the highest reported among the monogeneans thus far. Three tRNA-like cloverleaf structures were found in mNCR. Several results of the phylogenomic analysis are in disagreement with previously proposed relationships: instead of being closely related to the Gyrodactylidea, Tetraonchidea exhibit a phylogenetic affinity with the Dactylogyridea + Capsalidea clade; and the order Capsalidea is neither basal within the subclass Monopisthocotylea, nor groups with the Gyrodactylidea, but instead forms a sister clade with the Dactylogyridea. The mt genome of P. inermis exhibits a unique gene order, with an extensive reorganization of tRNAs. Monogenea exhibit exceptional gene order plasticity within the Neodermata.

CONCLUSIONS

This study shows that gene order within monopisthocotylid mt genomes is evolving at uneven rates, which creates misleading evolutionary signals. Furthermore, our results indicate that all previous attempts to resolve the evolutionary history of the Monogenea may have produced at least partially erroneous relationships. This further corroborates the necessity to generate more molecular data for this group of parasitic animals.

CASAAV: A CRISPR-Based Platform for Rapid Dissection of Gene Function In Vivo

In vivo loss-of-function studies are currently limited by the need for appropriate conditional knockout alleles. CRISPR/Cas9 is a powerful tool commonly used to induce loss-of-function mutations in vitro. However, CRISPR components have been difficult to deploy in vivo. To address this problem, we developed the CASAAV (CRISPR/Cas9/AAV-based somatic mutagenesis) platform, in which recombinant adeno-associated virus (AAV) is used to deliver tandem guide RNAs and Cre recombinase to Cre-dependent Cas9-P2A-GFP mice. Because Cre is under the control of a tissue-specific promoter, this system allows temporally controlled, cell type-selective knockout of virtually any gene to be obtained within a month using only one mouse line. Here, we focus on gene disruption in cardiomyocytes, but the system could easily be adapted to inactivate genes in other cell types transduced by AAV. © 2017 by John Wiley & Sons, Inc.

Scalable Production of iPSC-Derived Human Neurons to Identify Tau-Lowering Compounds by High-Content Screening

Lowering total tau levels is an attractive therapeutic strategy for Alzheimer's disease and other tauopathies. High-throughput screening in neurons derived from human induced pluripotent stem cells (iPSCs) is a powerful tool to identify tau-targeted therapeutics. However, such screens have been hampered by heterogeneous neuronal production, high cost and low yield, and multi-step differentiation procedures. We engineered an isogenic iPSC line that harbors an inducible neurogenin 2 transgene, a transcription factor that rapidly converts iPSCs to neurons, integrated at the AAVS1 locus. Using a simplified two-step protocol, we differentiated these iPSCs into cortical glutamatergic neurons with minimal well-to-well variability. We developed a robust high-content screening assay to identify tau-lowering compounds in LOPAC and identified adrenergic receptors agonists as a class of compounds that reduce endogenous human tau. These techniques enable the use of human neurons for high-throughput screening of drugs to treat neurodegenerative disease.

Identification, chromosomal arrangements and expression analyses of the evolutionarily conserved prmt1 gene in chicken in comparison with its vertebrate paralogue prmt8

Nine protein arginine methyltransferases (PRMTs) are conserved in mammals and fish. Among these, PRMT1 is the major type I PRMT for asymmetric dimethylarginine (ADMA) formation and is the most conserved and widely distributed one. Two chicken prmt1 splicing variants were assembled and confirmed by RT-PCR experiments. However, only two scaffolds containing single separate prmt1 exon with high GC contents are present in the current chicken genome assembly. Besides, prmt1 exons are scattered in separate small scaffolds in most avian species. Complete prmt1 gene has only been predicted from two falcon species with few neighboring genes. Crocodilians are considered close to the common ancestor shared by crocodilians and birds. The gene arrangements around prmt1 in American alligator are different from that in birds but are largely conserved in human. Orthologues of genes in a large segment of human chromosomal 19 around PRMT1 are missing or not assigned to the current chicken chromosomes. In comparison, prmt8, the prmt1 paralogue, is on chicken chromosome 1 with the gene arrangements downstream of prmt8 highly conserved in birds, crocodilians, and human. However, the ones upstream vary greatly in birds. Biochemically, we found that though prmt1 transcripts were detected, limited or none PRMT1 protein was present in chicken tissues. Moreover, a much higher level of PRMT8 protein was detected in chicken brain than in mouse brain. While PRMT8 is brain specific in other vertebrate species studied, low level of PRMT8 was present in chicken but not mouse liver and muscle. We also showed that the ADMA level in chicken was similar to that in mouse. This study provides the critical information of chicken PRMT1 and PRMT8 for future analyses of the function of protein arginine methyltransferases in birds.

Genomic structure and evolution of the mating type locus in the green seaweed Ulva partita

The evolution of sex chromosomes and mating loci in organisms with UV systems of sex/mating type determination in haploid phases via genes on UV chromosomes is not well understood. We report the structure of the mating type (MT) locus and its evolutionary history in the green seaweed Ulva partita, which is a multicellular organism with an isomorphic haploid-diploid life cycle and mating type determination in the haploid phase. Comprehensive comparison of a total of 12.0 and 16.6 Gb of genomic next-generation sequencing data for mt- and mt+ strains identified highly rearranged MT loci of 1.0 and 1.5 Mb in size and containing 46 and 67 genes, respectively, including 23 gametologs. Molecular evolutionary analyses suggested that the MT loci diverged over a prolonged period in the individual mating types after their establishment in an ancestor. A gene encoding an RWP-RK domain-containing protein was found in the mt- MT locus but was not an ortholog of the chlorophycean mating type determination gene MID. Taken together, our results suggest that the genomic structure and its evolutionary history in the U. partita MT locus are similar to those on other UV chromosomes and that the MT locus genes are quite different from those of Chlorophyceae.

The impact of the protein interactome on the syntenic structure of mammalian genomes

Conserved synteny denotes evolutionary preserved gene order across species. It is not well understood to which degree functional relationships between genes are preserved in syntenic blocks. Here we investigate whether protein-coding genes conserved in mammalian syntenic blocks encode gene products that serve the common functional purpose of interacting at protein level, i.e. connectivity. High connectivity among protein-protein interactions (PPIs) was only moderately associated with conserved synteny on a genome-wide scale. However, we observed a smaller subset of 3.6% of all syntenic blocks with high-confidence PPIs that had significantly higher connectivity than expected by random. Additionally, syntenic blocks with high-confidence PPIs contained significantly more chromatin loops than the remaining blocks, indicating functional preservation among these syntenic blocks. Conserved synteny is typically defined by sequence similarity. In this study, we also examined whether a functional relationship, here PPI connectivity, can identify syntenic blocks independently of orthology. While orthology-based syntenic blocks with high-confident PPIs and the connectivity-based syntenic blocks largely overlapped, the connectivity-based approach identified additional syntenic blocks that were not found by conventional sequence-based methods alone. Additionally, the connectivity-based approach enabled identification of potential orthologous genes between species. Our analyses demonstrate that subsets of syntenic blocks are associated with highly connected proteins, and that PPI connectivity can be used to detect conserved synteny even if sequence conservation drifts beyond what orthology algorithms normally can identify.

YeeJ is an inverse autotransporter from Escherichia coli that binds to peptidoglycan and promotes biofilm formation

Escherichia coli is a commensal or pathogenic bacterium that can survive in diverse environments. Adhesion to surfaces is essential for E. coli colonization, and thus it is important to understand the molecular mechanisms that promote this process in different niches. Autotransporter proteins are a class of cell-surface factor used by E. coli for adherence. Here we characterized the regulation and function of YeeJ, a poorly studied but widespread representative from an emerging class of autotransporter proteins, the inverse autotransporters (IAT). We showed that the yeeJ gene is present in ~40% of 96 completely sequenced E. coli genomes and that YeeJ exists as two length variants, albeit with no detectable functional differences. We demonstrated that YeeJ promotes biofilm formation in different settings through exposition at the cell-surface. We also showed that YeeJ contains a LysM domain that interacts with peptidoglycan and thus assists its localization into the outer membrane. Additionally, we identified the Polynucleotide Phosphorylase PNPase as a repressor of yeeJ transcription. Overall, our work provides new insight into YeeJ as a member of the recently defined IAT class, and contributes to our understanding of how commensal and pathogenic E. coli colonise their environments.

Generation of a new mouse model of glaucoma characterized by reduced expression of the AP-2β and AP-2δ proteins

We generated 6 transgenic lines with insertion of an expression plasmid for the R883/M xanthine dehydrogenase (XDH) mutant protein. Approximately 20% of the animals deriving from one of the transgenic lines show ocular abnormalities and an increase in intra-ocular pressure which are consistent with glaucoma. The observed pathologic phenotype is not due to expression of the transgene, but rather the consequence of the transgene insertion site, which has been defined by genome sequencing. The insertion site maps to chromosome 1qA3 in close proximity to the loci encoding AP-2β and AP-2δ, two proteins expressed in the eye. The insertion leads to a reduction in AP-2β and AP-2δ levels. Down-regulation of AP-2β expression is likely to be responsible for the pathologic phenotype, as conditional deletion of the Tfap2b gene in the neural crest has recently been shown to cause defective development of the eye anterior segment and early-onset glaucoma. In these conditional knock-out and our transgenic mice, the morphological/histological features of the glaucomatous pathology are surprisingly similar. Our transgenic mouse represents a model of angle-closure glaucoma and a useful tool for the study of the pathogenesis and the development of innovative therapeutic strategies.

Identification and genomic comparison of temperate bacteriophages derived from emetic Bacillus cereus

Cereulide-producing Bacillus cereus isolates can cause serious emetic (vomiting) syndrome and even acute lethality. As mobile genetic elements, the exploration of prophages derived from emetic B. cereus isolates will help in our understanding of the genetic diversity and evolution of these pathogens. In this study, five temperate phages derived from cereulide-producing B. cereus strains were induced, with four of them undergoing genomic sequencing. Sequencing revealed that they all belong to the Siphoviridae family, but presented in different forms in their hosts. PfNC7401 and PfIS075 have typical icosahedral heads, probably existing alone as phagemids in the host with self-replicating capability in the lysogenic state. PfEFR-4, PfEFR-5, and PfATCC7953 have elongated heads, with the genomes of the former two identified as linear dsDNA, which could be integrated into the host genome during the lysogenic state. Genomic comparison of the four phages with others also derived from emetic B. cereus isolates showed similar genome structures and core genes, thus displaying host spectrum specificity. In addition, phylogenic analysis based on the complete genome and conserved tail fiber proteins of 36 Bacillus species-derived phages confirmed that the phages derived from emetic B. cereus strains were highly similar. Furthermore, one endolysin LysPfEFR-4 was cloned and showed lytic activity against all tested emetic B. cereus strains and cross-lytic activity against some other pathogenic bacteria, implying a potential to control bacterial contamination in the food supply.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

Assigning ecological roles to the populations belonging to a phenanthrene-degrading bacterial consortium using omic approaches

The present study describes the behavior of a natural phenanthrene-degrading consortium (CON), a synthetic consortium (constructed with isolated strains from CON) and an isolated strain form CON (Sphingobium sp. AM) in phenanthrene cultures to understand the interactions among the microorganisms present in the natural consortium during phenanthrene degradation as a sole carbon and energy source in liquid cultures. In the contaminant degradation assay, the defined consortium not only achieved a major phenanthrene degradation percentage (> 95%) but also showed a more efficient elimination of the intermediate metabolite. The opposite behavior occurred in the CON culture where the lowest phenanthrene degradation and the highest HNA accumulation were observed, which suggests the presence of positive and also negative interaction in CON. To consider the uncultured bacteria present in CON, a metagenomic library was constructed with total CON DNA. One of the resulting scaffolds (S1P3) was affiliated with the Betaproteobacteria class and resulted in a significant similarity with a genome fragment from Burkholderia sp. HB1 chromosome 1. A complete gene cluster, which is related to one of the lower pathways (meta-cleavage of catechol) involved in PAH degradation (ORF 31–43), mobile genetic elements and associated proteins, was found. These results suggest the presence of at least one other microorganism in CON besides Sphingobium sp. AM, which is capable of degrading PAH through the meta-cleavage pathway. Burkholderiales order was further found, along with Sphingomonadales order, by a metaproteomic approach, which indicated that both orders were metabolically active in CON. Our results show the presence of negative interactions between bacterial populations found in a natural consortium selected by enrichment techniques; moreover, the synthetic syntrophic processing chain with only one microorganism with the capability of degrading phenanthrene was more efficient in contaminant and intermediate metabolite degradation than a generalist strain (Sphingobium sp. AM).

Genome-Wide Screen for Genes Involved in Caenorhabditis elegans Developmentally Timed Sleep

In Caenorhabditis elegans, Notch signaling regulates developmentally timed sleep during the transition from L4 larval stage to adulthood (L4/A) . To identify core sleep pathways and to find genes acting downstream of Notch signaling, we undertook the first genome-wide, classical genetic screen focused on C. elegans developmentally timed sleep. To increase screen efficiency, we first looked for mutations that suppressed inappropriate anachronistic sleep in adult hsp::osm-11 animals overexpressing the Notch coligand OSM-11 after heat shock. We retained suppressor lines that also had defects in L4/A developmentally timed sleep, without heat shock overexpression of the Notch coligand. Sixteen suppressor lines with defects in developmentally timed sleep were identified. One line carried a new allele of goa-1; loss of GOA-1 Gαo decreased C. elegans sleep. Another line carried a new allele of gpb-2, encoding a Gβ5 protein; Gβ5 proteins have not been previously implicated in sleep. In other scenarios, Gβ5 GPB-2 acts with regulators of G protein signaling (RGS proteins) EAT-16 and EGL-10 to terminate either EGL-30 Gαq signaling or GOA-1 Gαo signaling, respectively. We found that loss of Gβ5 GPB-2 or RGS EAT-16 decreased L4/A sleep. By contrast, EGL-10 loss had no impact. Instead, loss of RGS-1 and RGS-2 increased sleep. Combined, our results suggest that, in the context of L4/A sleep, GPB-2 predominantly acts with EAT-16 RGS to inhibit EGL-30 Gαq signaling. These results confirm the importance of G protein signaling in sleep and demonstrate that these core sleep pathways function genetically downstream of the Notch signaling events promoting sleep.

Using CRISPR-Cas9 to Generate Gene-Corrected Autologous iPSCs for the Treatment of Inherited Retinal Degeneration

Patient-derived induced pluripotent stem cells (iPSCs) hold great promise for autologous cell replacement. However, for many inherited diseases, treatment will likely require genetic repair pre-transplantation. Genome editing technologies are useful for this application. The purpose of this study was to develop CRISPR-Cas9-mediated genome editing strategies to target and correct the three most common types of disease-causing variants in patient-derived iPSCs: (1) exonic, (2) deep intronic, and (3) dominant gain of function. We developed a homology-directed repair strategy targeting a homozygous Alu insertion in exon 9 of male germ cell-associated kinase (MAK) and demonstrated restoration of the retinal transcript and protein in patient cells. We generated a CRISPR-Cas9-mediated non-homologous end joining (NHEJ) approach to excise a major contributor to Leber congenital amaurosis, the IVS26 cryptic-splice mutation in CEP290, and demonstrated correction of the transcript and protein in patient iPSCs. Lastly, we designed allele-specific CRISPR guides that selectively target the mutant Pro23His rhodopsin (RHO) allele, which, following delivery to both patient iPSCs in vitro and pig retina in vivo, created a frameshift and premature stop that would prevent transcription of the disease-causing variant. The strategies developed in this study will prove useful for correcting a wide range of genetic variants in genes that cause inherited retinal degeneration.

Functional expression of the entire adhesiome of Salmonella enterica serotype Typhimurium

Adhesins are crucial virulence factors of pathogenic bacteria involved in colonization, transmission and pathogenesis. Many bacterial genomes contain the information for a surprisingly large number of diverse adhesive structures. One prominent example is the invasive and facultative intracellular pathogen Salmonella enterica with an adhesiome of up to 20 adhesins. Such large repertoire of adhesins contributes to colonization of a broad range of host species and may allow adaptation to various environments within the host, as well as in non-host environments. For S. enterica, only few members of the adhesiome are functionally expressed under laboratory conditions, and accordingly the structural and functional understanding of the majority of adhesins is sparse. We have devised a simple and versatile approach to functionally express all adhesins of S. enterica serotype Typhimurium, either within Salmonella or within heterologous hosts such as Escherichia coli. We demonstrate the surface expression of various so far cryptic adhesins and show ultrastructural features using atomic force microscopy and transmission electron microscopy. In summary, we report for the first time the expression of the entire adhesiome of S. enterica serotype Typhimurium.

The complete chloroplast genome sequence of Aconitum coreanum and Aconitum carmichaelii and comparative analysis with other Aconitum species

Aconitum species (belonging to the Ranunculaceae) are well known herbaceous medicinal ingredients and have great economic value in Asian countries. However, there are still limited genomic resources available for Aconitum species. In this study, we sequenced the chloroplast (cp) genomes of two Aconitum species, A. coreanum and A. carmichaelii, using the MiSeq platform. The two Aconitum chloroplast genomes were 155,880 and 157,040 bp in length, respectively, and exhibited LSC and SSC regions separated by a pair of inverted repeat regions. Both cp genomes had 38% GC content and contained 131 unique functional genes including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. The gene order, content, and orientation of the two Aconitum cp genomes exhibited the general structure of angiosperms, and were similar to those of other Aconitum species. Comparison of the cp genome structure and gene order with that of other Aconitum species revealed general contraction and expansion of the inverted repeat regions and single copy boundary regions. Divergent regions were also identified. In phylogenetic analysis, Aconitum species positon among the Ranunculaceae was determined with other family cp genomes in the Ranunculales. We obtained a barcoding target sequence in a divergent region, ndhC–trnV, and successfully developed a SCAR (sequence characterized amplified region) marker for discrimination of A. coreanum. Our results provide useful genetic information and a specific barcode for discrimination of Aconitum species.

Molecular Structure and Phylogenetic Analyses of Complete Chloroplast Genomes of Two Aristolochia Medicinal Species

The family Aristolochiaceae, comprising about 600 species of eight genera, is a unique plant family containing aristolochic acids (AAs). The complete chloroplast genome sequences of Aristolochia debilis and Aristolochia contorta are reported here. The results show that the complete chloroplast genomes of A. debilis and A. contorta comprise circular 159,793 and 160,576 bp-long molecules, respectively and have typical quadripartite structures. The GC contents of both species were 38.3% each. A total of 131 genes were identified in each genome including 85 protein-coding genes, 37 tRNA genes, eight rRNA genes and one pseudogene (ycf1). The simple-sequence repeat sequences mainly comprise A/T mononucletide repeats. Phylogenetic analyses using maximum parsimony (MP) revealed that A. debilis and A. contorta had a close phylogenetic relationship with species of the family Piperaceae, as well as Laurales and Magnoliales. The data obtained in this study will be beneficial for further investigations on A. debilis and A. contorta from the aspect of evolution, and chloroplast genetic engineering.

An enhanced hTERT promoter-driven CRISPR/Cas9 system selectively inhibits the progression of bladder cancer cells

The current therapies for treating tumors are lacking in efficacy and specificity. Synthetic biology principles may bring some new possible methods for curing cancer. Here we present a synthetic logic circuit based on the CRISPR/Cas9 system. The CRISPR/Cas9 technology has been applied in many biological fields, including cancer research. In this study, the expression of Cas9 nuclease was controlled indirectly by an enhanced hTERT promoter using the GAL4/upstream activating sequence (UAS) binding system. Cas9 was driven by 5XUAS, single guide RNA (sgRNA) was used to target mutant or wild-type HRAS, and the fusion gene GAL4-P65 was driven by the enhanced hTERT promoter. The system was tested in bladder cancer cells (T24 and 5637) and the results showed that the enhanced hTERT promoter could drive the expression of GAL4-P65 in these bladder cancer cell lines. Then all these devices were packed into lentivirus and the results of quantitative real-time PCR showed that the mRNA expression level of HRAS was selectively inhibited in the T24 and 5637 cells. The results of functional experiments suggested that the proliferation, cell migration and invasion were selectively suppressed, and that the apoptosis rate was increased in bladder cancer cells but not in human foreskin fibroblasts (HFF). In conclusion, we successfully constructed an enhanced hTERT promoter-driven CRISPR/Cas9 system and data showed that it could selectively suppress the progression of bladder cancer cells.

Mitochondrial genome evolution in the Saccharomyces sensu stricto complex

Exploring the evolutionary patterns of mitochondrial genomes is important for our understanding of the Saccharomyces sensu stricto (SSS) group, which is a model system for genomic evolution and ecological analysis. In this study, we first obtained the complete mitochondrial sequences of two important species, Saccharomyces mikatae and Saccharomyces kudriavzevii. We then compared the mitochondrial genomes in the SSS group with those of close relatives, and found that the non-coding regions evolved rapidly, including dramatic expansion of intergenic regions, fast evolution of introns and almost 20-fold higher rearrangement rates than those of the nuclear genomes. However, the coding regions, and especially the protein-coding genes, are more conserved than those in the nuclear genomes of the SSS group. The different evolutionary patterns of coding and non-coding regions in the mitochondrial and nuclear genomes may be related to the origin of the aerobic fermentation lifestyle in this group. Our analysis thus provides novel insights into the evolution of mitochondrial genomes.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of the Medicinal Plant Artemisia annua

The complete chloroplast genome of Artemisia annua (Asteraceae), the primary source of artemisinin, was sequenced and analyzed. The A. annua cp genome is 150,995 bp, and harbors a pair of inverted repeat regions (IRa and IRb), of 24,850 bp each that separate large (LSC, 82,988 bp) and small (SSC, 18,267 bp) single-copy regions. Our annotation revealed that the A. annua cp genome contains 113 genes and 18 duplicated genes. The gene order in the SSC region of A. annua is inverted; this fact is consistent with the sequences of chloroplast genomes from three other Artemisia species. Fifteen (15) forward and seventeen (17) inverted repeats were detected in the genome. The existence of rich SSR loci in the genome suggests opportunities for future population genetics work on this anti-malarial medicinal plant. In A. annua cpDNA, the rps19 gene was found in the LSC region rather than the IR region, and the rps19 pseudogene was absent in the IR region. Sequence divergence analysis of five Asteraceae species indicated that the most highly divergent regions were found in the intergenic spacers, and that the differences between A. annua and A. fukudo were very slight. A phylogenetic analysis revealed a sister relationship between A. annua and A. fukudo. This study identified the unique characteristics of the A. annua cp genome. These results offer valuable information for future research on Artemisia species identification and for the selective breeding of A. annua with high pharmaceutical efficacy.

Complete mitochondrial genomes of Anopheles stephensi and An. dirus and comparative evolutionary mitochondriomics of 50 mosquitoes

To better understand the phylogeny and evolution of mosquitoes, the complete mitochondrial genome (mitogenome) of Anopheles stephensi and An. dirus were sequenced and annotated, and a total of 50 mosquito mitogenomes were comparatively analyzed. The complete mitogenome of An. stephensi and An. dirus is 1,5371 bp and 1,5406 bp long, respectively. The main features of the 50 mosquito mitogenomes are conservative: 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, positive AT-skew and negative GC-skew. The gene order trnA-trnR in ancestral insects is rearranged. All tRNA genes have the typical clover leaf secondary structure but tRNA Ser . The control regions are highly variable in size. PCGs show signals of purifying selection, but evidence for positive selection in ND2, ND4 and ND6 is found. Bayesian and Maximum Likelihood phylogenetic analyses based on all PCG nucleotides produce an identical tree topology and strongly support the monophyly of subgenera Cellia, Anopheles, Keterszia and Nyssorhynchus, the sister relationship of the subgenera Nyssorhynchus and Keterszia, and Cellia and Anopheles. The most recent ancestor of the genus Anopheles and Culicini + Aedini exited ~145 Mya ago. This is the first comprehensive study of mosquito mitogenomes, which are effective for mosquito phylogeny at various taxonomic levels.

Chloroplast genomes of Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea: Structures and comparative analysis

We investigated the complete chloroplast (cp) genomes of non-model Arabidopsis halleri ssp. gemmifera and Arabidopsis lyrata ssp. petraea using Illumina paired-end sequencing to understand their genetic organization and structure. Detailed bioinformatics analysis revealed genome sizes of both subspecies ranging between 154.4~154.5 kbp, with a large single-copy region (84,197~84,158 bp), a small single-copy region (17,738~17,813 bp) and pair of inverted repeats (IRa/IRb; 26,264~26,259 bp). Both cp genomes encode 130 genes, including 85 protein-coding genes, eight ribosomal RNA genes and 37 transfer RNA genes. Whole cp genome comparison of A. halleri ssp. gemmifera and A. lyrata ssp. petraea, along with ten other Arabidopsis species, showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. The location and distribution of repeat sequences were determined, and sequence divergences of shared genes were calculated among related species. Comparative phylogenetic analysis of the entire genomic data set and 70 shared genes between both cp genomes confirmed the previous phylogeny and generated phylogenetic trees with the same topologies. The sister species of A. halleri ssp. gemmifera is A. umezawana, whereas the closest relative of A. lyrata spp. petraea is A. arenicola.

Improved base excision repair inhibition and bacteriophage Mu Gam protein yields C:G-to-T:A base editors with higher efficiency and product purity

We recently developed base editing, the programmable conversion of target C:G base pairs to T:A without inducing double-stranded DNA breaks (DSBs) or requiring homology-directed repair using engineered fusions of Cas9 variants and cytidine deaminases. Over the past year, the third-generation base editor (BE3) and related technologies have been successfully used by many researchers in a wide range of organisms. The product distribution of base editing-the frequency with which the target C:G is converted to mixtures of undesired by-products, along with the desired T:A product-varies in a target site-dependent manner. We characterize determinants of base editing outcomes in human cells and establish that the formation of undesired products is dependent on uracil N-glycosylase (UNG) and is more likely to occur at target sites containing only a single C within the base editing activity window. We engineered CDA1-BE3 and AID-BE3, which use cytidine deaminase homologs that increase base editing efficiency for some sequences. On the basis of these observations, we engineered fourth-generation base editors (BE4 and SaBE4) that increase the efficiency of C:G to T:A base editing by approximately 50%, while halving the frequency of undesired by-products compared to BE3. Fusing BE3, BE4, SaBE3, or SaBE4 to Gam, a bacteriophage Mu protein that binds DSBs greatly reduces indel formation during base editing, in most cases to below 1.5%, and further improves product purity. BE4, SaBE4, BE4-Gam, and SaBE4-Gam represent the state of the art in C:G-to-T:A base editing, and we recommend their use in future efforts.

Transcriptional Control of the Lateral-Flagellar Genes of Bradyrhizobium diazoefficiens

Bradyrhizobium diazoefficiens, a soybean N2-fixing symbiont, possesses a dual flagellar system comprising a constitutive subpolar flagellum and inducible lateral flagella. Here, we analyzed the genomic organization and biosynthetic regulation of the lateral-flagellar genes. We found that these genes are located in a single genomic cluster, organized in two monocistronic transcriptional units and three operons, one possibly containing an internal transcription start site. Among the monocistronic units is blr6846, homologous to the class IB master regulators of flagellum synthesis in Brucella melitensis and Ensifer meliloti and required for the expression of all the lateral-flagellar genes except lafA2, whose locus encodes a single lateral flagellin. We therefore named blr6846 lafR (lateral-flagellar regulator). Despite its similarity to two-component response regulators and its possession of a phosphorylatable Asp residue, lafR behaved as an orphan response regulator by not requiring phosphorylation at this site. Among the genes induced by lafR is flbTL , a class III regulator. We observed different requirements for FlbTL in the synthesis of each flagellin subunit. Although the accumulation of lafA1, but not lafA2, transcripts required FlbTL, the production of both flagellin polypeptides required FlbTL Moreover, the regulation cascade of this lateral-flagellar regulon appeared to be not as strictly ordered as those found in other bacterial species.IMPORTANCE Bacterial motility seems essential for the free-living style in the environment, and therefore these microorganisms allocate a great deal of their energetic resources to the biosynthesis and functioning of flagella. Despite energetic costs, some bacterial species possess dual flagellar systems, one of which is a primary system normally polar or subpolar, and the other is a secondary, lateral system that is produced only under special circumstances. Bradyrhizobium diazoefficiens, an N2-fixing symbiont of soybean plants, possesses dual flagellar systems, including the lateral system that contributes to swimming in wet soil and competition for nodulation and is expressed under high energy availability, as well as under requirement for high torque by the flagella. The structural organization and transcriptional regulation of the 41 genes that comprise this secondary flagellar system seem adapted to adjust bacterial energy expenditures for motility to the soil's environmental dynamics.

Mitochondrial Genomes of Two Bombycoidea Insects and Implications for Their Phylogeny

The mitochondrial genome (mt genome) provides important information for understanding molecular evolution and phylogenetics. As such, the two complete mt genomes of Ampelophaga rubiginosa and Rondotia menciana were sequenced and annotated. The two circular genomes of A. rubiginosa and R. menciana are 15,282 and 15,636 bp long, respectively, including 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. The nucleotide composition of the A. rubiginosa mt genome is A + T rich (81.5%) but is lower than that of R. menciana (82.2%). The AT skew is slightly positive and the GC skew is negative in these two mt genomes. Except for cox1, which started with CGA, all other 12PCGs started with ATN codons. The A + T-rich regions of A. rubiginosa and R. menciana were 399 bp and 604 bp long and consist of several features common to Bombycoidea insects. The order and orientation of A. rubiginosa and R. menciana mitogenomes with the order trnM-trnI-trnQ-nad2 is different from the ancestral insects in which trnM is located between trnQ and nad2 (trnI-trnQ-trnM-nad2). Phylogenetic analyses indicate that A. rubiginosa belongs in the Sphingidae family, and R. menciana belongs in the Bombycidae family.

Phylogenetic relationships of diurnal, phytotelm-breeding Melanophryniscus (Anura: Bufonidae) based on mitogenomic data

Melanophryniscus is a bufonid frog genus with a broad geographic distribution over southeastern South America. In recent years, several new species of Melanophryniscus have been discovered in southern Brazil showing a distinctive life-history strategy for the genus - breeding in phytotelmata - as well as a strong association with high-altitude regions. In this study, we use mitogenomic data to infer the phylogenetic relationships among diurnal, phytotelm-breeding Melanophryniscus and to determine the timing of their divergence. We obtained the mitochondrial genomes (not including the control region) for eight individuals of Melanophryniscus representing all three described species (M. alipioi, M. milanoi, and M. xanthostomus), as well as some recently-discovered and potentially new species. Gene order was conserved in all species and corresponded to the general order found in bufonids. Although the phylogenetic relationships among the studied species was poorly supported, dating confirmed that they diverged during the Pleistocene, suggesting that phytotelm breeding could have arisen during drier periods in the glacial/interglacial cycles due to a decrease in the availability of permanent streams or ephemeral/temporary streams or ponds in which Melanophryniscus species commonly breed.

Deep sequencing analysis reveals a TMV mutant with a poly(A) tract reduces host defense responses in Nicotiana benthamiana

Tobacco mosaic virus (TMV) possesses an upstream pseudoknotted domain (UPD), which is important for replication. After substituting the UPD with an internal poly(A) tract (43 nt), a mutant TMV-43A was constructed. TMV-43A replicated slower than TMV and induced a non-lethal mosaic symptom in Nicotiana benthamiana. In this study, deep sequencing was performed to detect the differences of small RNA profiles between TMV- and TMV-43A-infected N. benthamiana. The results showed that TMV-43A produced lesser amount of virus-derived interfering RNAs (vsiRNAs) than that of TMV. However, the distributions of vsiRNAs generation hotspots between TMV and TMV-43A were similar. Expression of genes related to small RNA biogenesis in TMV-43A-infected N. benthamiana was significantly lower than that of TMV, which leads to generation of lesser vsiRNAs. The expressions of host defense response genes were up-regulated after TMV infection, as compared to TMV-43A-infected plants. Host defense response to TMV-43A infection was lower than that to TMV. The absence of UPD might contribute to the reduced host response to TMV-43A. Our study provides valuable information in the role of the UPD in eliciting host response genes after TMV infection in N. benthamiana. (187 words).

Whole-genome characterization of Uruguayan strains of avian infectious bronchitis virus reveals extensive recombination between the two major South American lineages

Infectious bronchitis virus (Gammacoronavirus, Coronaviridae) is a genetically variable RNA virus that causes one of the most persistent respiratory diseases in poultry. The virus is classified in genotypes and lineages with different epidemiological relevance. Two lineages of the GI genotype (11 and 16) have been widely circulating for decades in South America. GI-11 is an exclusive South American lineage while the GI-16 lineage is distributed in Asia, Europe and South America. Here, we obtained the whole genome of two Uruguayan strains of the GI-11 and GI-16 lineages using Illumina high-throughput sequencing. The strains here sequenced are the first obtained in South America for the infectious bronchitis virus and provide new insights into the origin, spreading and evolution of viral variants. The complete genome of the GI-11 and GI-16 strains have 27,621 and 27,638 nucleotides, respectively, and possess the same genomic organization. Phylogenetic incongruence analysis reveals that both strains have a mosaic genome that arose by recombination between Euro Asiatic strains of the GI-16 lineage and ancestral South American GI-11 viruses. The recombination occurred in South America and produced two viral variants that have retained the full-length S1 sequences of the parental lineages but are extremely similar in the rest of their genomes. These recombinant virus have been extraordinary successful, persisting in the continent for several years with a notorious wide geographic distribution. Our findings reveal a singular viral dynamics and emphasize the importance of complete genomic characterization to understand the emergence and evolutionary history of viral variants.

•

Genomic analysis was performed in two main lineages of Infectious bronchitis virus.

•

Lineages differ in their S1 sequences but are similar in the rest of the genome.

•

Genomic similarity between both lineages arise by inter-lineage recombination.

•

Inter-lineage recombination occurred in South America between European/Asiatic and local strain.

•

Recombinant forms have persisted in the continent for several years with wide geographic distribution.

Towards an Understanding of Enhanced Biomass Digestibility by In Planta Expression of a Family 5 Glycoside Hydrolase

In planta expression of a thermophilic endoglucanase (AcCel5A) reduces recalcitrance by creating voids and other irregularities in cell walls of Arabidopsis thaliana that increase enzyme accessibility without negative impacts on plant growth or cell wall composition. Our results suggest that cellulose β-1-4 linkages can be cut sparingly in the assembling wall and that these minimal changes, made at the proper time, have an impact on plant cell wall recalcitrance without negative effects on overall plant development.

The highly rearranged mitochondrial genomes of the crabs Maja crispata and Maja squinado (Majidae) and gene order evolution in Brachyura

We sequenced the mitochondrial genomes of the spider crabs Maja crispata and Maja squinado (Majidae, Brachyura). Both genomes contain the whole set of 37 genes characteristic of Bilaterian genomes, encoded on both α- and β-strands. Both species exhibit the same gene order, which is unique among known animal genomes. In particular, all the genes located on the β-strand form a single block. This gene order was analysed together with the other nine gene orders known for the Brachyura. Our study confirms that the most widespread gene order (BraGO) represents the plesiomorphic condition for Brachyura and was established at the onset of this clade. All other gene orders are the result of transformational pathways originating from BraGO. The different gene orders exhibit variable levels of genes rearrangements, which involve only tRNAs or all types of genes. Local homoplastic arrangements were identified, while complete gene orders remain unique and represent signatures that can have a diagnostic value. Brachyura appear to be a hot-spot of gene order diversity within the phylum Arthropoda. Our analysis, allowed to track, for the first time, the fully evolutionary pathways producing the Brachyuran gene orders. This goal was achieved by coupling sophisticated bioinformatic tools with phylogenetic analysis.

Complete chloroplast genome sequence of common bermudagrass (Cynodon dactylon (L.) Pers.) and comparative analysis within the family Poaceae

Common bermudagrass (Cynodon dactylon (L.) Pers.) belongs to the subfamily Chloridoideae of the Poaceae family, one of the most important plant families ecologically and economically. This grass has a long connection with human culture but its systematics is relatively understudied. In this study, we sequenced and investigated the chloroplast genome of common bermudagrass, which is 134,297 bp in length with two single copy regions (LSC: 79,732 bp; SSC: 12,521 bp) and a pair of inverted repeat (IR) regions (21,022 bp). The annotation contains a total of 128 predicted genes, including 82 protein-coding, 38 tRNA, and 8 rRNA genes. Additionally, our in silico analyses identified 10 sets of repeats longer than 20 bp and predicted the presence of 36 RNA editing sites. Overall, the chloroplast genome of common bermudagrass resembles those from other Poaceae lineages. Compared to most angiosperms, the accD gene and the introns of both clpP and rpoC1 genes are missing. Additionally, the ycf1, ycf2, ycf15, and ycf68 genes are pseudogenized and two genome rearrangements exist. Our phylogenetic analysis based on 47 chloroplast protein-coding genes supported the placement of common bermudagrass within Chloridoideae. Our phylogenetic character mapping based on the parsimony principle further indicated that the loss of the accD gene and clpP introns, the pseudogenization of four ycf genes, and the two rearrangements occurred only once after the most recent common ancestor of the Poaceae diverged from other monocots, which could explain the unusual long branch leading to the Poaceae when phylogeny is inferred based on chloroplast sequences.

A genome-scale metabolic reconstruction of Lysinibacillus sphaericus unveils unexploited biotechnological potentials

The toxic lineage (TL) of Lysinibacillus sphaericus has been extensively studied because of its potential biotechnological applications in biocontrol of mosquitoes and bioremediation of toxic metals. We previously proposed that L. sphaericus TL should be considered as a novel species based on a comparative genomic analysis. In the current work, we constructed the first manually curated metabolic reconstruction for this species on the basis of the available genomes. We elucidated the central metabolism of the proposed species and, beyond confirming the reported experimental evidence with genomic a support, we found insights to propose novel applications and traits to be considered in further studies. The strains belonging to this lineage exhibit a broad repertory of genes encoding insecticidal factors, some of them remain uncharacterized. These strains exhibit other unexploited biotechnological important traits, such as lactonases (quorum quenching), toxic metal resistance, and potential for aromatic compound degradation. In summary, this study provides a guideline for further research aimed to implement this organism in biocontrol and bioremediation. Similarly, we highlighted the unanswered questions to be responded in order to gain a deeper understanding of the L. sphaericus TL biology.

Short DNA Hairpins Compromise Recombinant Adeno-Associated Virus Genome Homogeneity

Short hairpin (sh)RNAs delivered by recombinant adeno-associated viruses (rAAVs) are valuable tools to study gene function in vivo and a promising gene therapy platform. Our data show that incorporation of shRNA transgenes into rAAV constructs reduces vector yield and produces a population of truncated and defective genomes. We demonstrate that sequences with hairpins or hairpin-like structures drive the generation of truncated AAV genomes through a polymerase redirection mechanism during viral genome replication. Our findings reveal the importance of genomic secondary structure when optimizing viral vector designs. We also discovered that shDNAs could be adapted to act as surrogate mutant inverted terminal repeats (mTRs), sequences that were previously thought to be required for functional self-complementary AAV vectors. The use of shDNAs as artificial mTRs opens the door to engineering a new generation of AAV vectors with improved potency, genetic stability, and safety for both preclinical studies and human gene therapy.

HIV-1 Conserved Mosaics Delivered by Regimens with Integration-Deficient DC-Targeting Lentiviral Vector Induce Robust T Cells

To be effective against HIV type 1 (HIV-1), vaccine-induced T cells must selectively target epitopes, which are functionally conserved (present in the majority of currently circulating and reactivated HIV-1 strains) and, at the same time, beneficial (responses to which are associated with better clinical status and control of HIV-1 replication), and rapidly reach protective frequencies upon exposure to the virus. Heterologous prime-boost regimens using virally vectored vaccines are currently the most promising vaccine strategies; nevertheless, induction of robust long-term memory remains challenging. To this end, lentiviral vectors induce high frequencies of memory cells due to their low-inflammatory nature, while typically inducing only low anti-vector immune responses. Here, we describe construction of novel candidate vaccines ZVex.tHIVconsv1 and ZVex.tHIVconsv2, which are based on an integration-deficient lentiviral vector platform with preferential transduction of human dendritic cells and express a bivalent mosaic of conserved-region T cell immunogens with a high global HIV-1 match. Each of the two mosaic vaccines was individually immunogenic. When administered together in heterologous prime-boost regimens with chimpanzee adenovirus and/or poxvirus modified vaccinia virus Ankara (MVA) vaccines to BALB/c and outbred CD1-Swiss mice, they induced a median frequency of over 6,000 T cells/10⁶ splenocytes, which were plurifunctional, broadly specific, and cross-reactive. These results support further development of this vaccine concept.

Expression and functional analysis of the lysine decarboxylase and copper amine oxidase genes from the endophytic fungus Colletotrichum gloeosporioides ES026

Huperzine A (HupA) isolated from Huperzia serrata is an important compound used to treat Alzheimer's disease (AD). Recently, HupA was reported in various endophytic fungi, with Colletotrichum gloeosporioides ES026 previously isolated from H. serrata shown to produce HupA. In this study, we performed next-generation sequencing and de novo RNA sequencing of C. gloeosporioides ES026 to elucidate the molecular functions, biological processes, and biochemical pathways of these unique sequences. Gene ontology and Kyoto Encyclopedia of Genes and Genomes assignments allowed annotation of lysine decarboxylase (LDC) and copper amine oxidase (CAO) for their conversion of L-lysine to 5-aminopentanal during HupA biosynthesis. Additionally, we constructed a stable, high-yielding HupA-expression system resulting from the overexpression of CgLDC and CgCAO from the HupA-producing endophytic fungus C. gloeosporioides ES026 in Escherichia coli. Quantitative reverse transcription polymerase chain reaction analysis confirmed CgLDC and CgCAO expression, and quantitative determination of HupA levels was assessed by liquid chromatography high-resolution mass spectrometry, which revealed that elevated expression of CgLDC and CgCAO produced higher yields of HupA than those derived from C. gloeosporioides ES026. These results revealed CgLDC and CgCAO involvement in HupA biosynthesis and their key role in regulating HupA content in C. gloeosporioides ES026.