Search

Microbial endosymbiosis is widespread in animals, with major ecological and evolutionary implications. Successful symbiosis relies on efficient vertical transmission through host generations. However, when symbionts negatively affect host fitness, hosts are expected to evolve suppression of symbiont effects or transmission. Here, we show that sex chromosomes control vertical transmission of feminizing Wolbachia endosymbionts in the isopod Armadillidium nasatum. Theory predicts that the invasion of an XY/XX species by cytoplasmic sex ratio distorters is unlikely because it leads to fixation of the unusual (and often lethal or infertile) YY genotype. We demonstrate that A. nasatum X and Y sex chromosomes are genetically highly similar and that YY individuals are viable and fertile, thereby enabling Wolbachia spread in this XY-XX species. Nevertheless, we show that Wolbachia cannot drive fixation of YY individuals, because infected YY females do not transmit Wolbachia to their offspring, unlike XX and XY females. The genetic basis fits the model of a Y-linked recessive allele (associated with an X-linked dominant allele), in which the homozygous state suppresses Wolbachia transmission. Moreover, production of all-male progenies by infected YY females restores a balanced sex ratio at the host population level. This suggests that blocking of Wolbachia transmission by YY females may have evolved to suppress feminization, thereby offering a whole new perspective on the evolutionary interplay between microbial symbionts and host sex chromosomes. [ABSTRACT FROM AUTHOR]

Lutzomyia longipalpis is the main vector of Leishmania infantum, the causative agent of visceral leishmaniasis in the Neotropical realm. Its taxonomic status has been widely discussed once it encompasses a complex of species. The knowledge about the genetic structure of insect vector populations helps the elucidation of components and interactions of the disease ecoepidemiology. Thus, the objective of this study was to genotypically analyze populations of the Lu. longipalpis complex from a macrogeographic perspective using Next Generation Sequencing. Polymorphism analysis of three molecular markers was used to access the levels of population genetic structure among nine different populations of sand flies. Illumina Amplicon Sequencing Protocol® was used to identify possible polymorphic sites. The library was sequenced on paired-end Illumina MiSeq platform. Significant macrogeographical population differentiation was observed among Lu. longipalpis populations via PCA and DAPC analyses. Our results revealed that populations of Lu. longipalpis from the nine municipalities were grouped into three clusters. In addition, it was observed that the levels of Lu. longipalpis population structure could be associated with distance isolation. This new sequencing method allowed us to study different molecular markers after a single sequencing run, and to evaluate population and inter-species differences on a macrogeographic scale. [ABSTRACT FROM AUTHOR]

Many disease resistance genes that have been transferred from wild relatives to cultivated wheat have played a significant role in wheat production worldwide. Ae. umbellulata is one of the species within the genus Aegilops that have been successfully used as sources of resistance genes to leaf rust, stem rust and powdery mildew. The objectives of the current work was to validate the map position of a major QTL that confers resistance to the stem rust pathogen races Ug99 (TTKSK) and TTTTF with an independent bi-parental mapping population and to refine the QTL region with a bulk segregant analysis approach. Two F2 bi-parental mapping populations were developed from stem rust resistant Ae. umbellulata accessions (PI 298905 and PI 5422375) and stem rust susceptible accessions (PI 542369 and PI 554395). Firstly, one of the two populations was used to map the chromosome location of the resistance gene. Later on, the 2nd population was used to validate the chromosome location in combination with a bulk segregant analysis approach. For the bulk segregant analysis, RNA was extracted from a bulk of leaf tissues of 12 homozygous resistant F3 families, and a separate bulk of 11 susceptible homozygous F3 families derived from the PI 5422375 and PI 554395 cross. The RNA samples of the two bulks and the two parents were sequenced for SNPs identification. Stem rust resistance QTL was validated on chromosome 2U of Ae. umbellulata in the same region in both populations. With bulk segregant analysis, the QTL position was delimited within 3.2 Mbp. Although there were a large number of genes in the orthologous region of the detected QTL on chromosome 2D of Ae. tauschii, we detected only two Ae. umbellulata NLR genes which can be considered as a potential candidate genes. [ABSTRACT FROM AUTHOR]

Newly emerged wheat blast disease is a serious threat to global wheat production. Wheat blast is caused by a distinct, exceptionally diverse lineage of the fungus causing rice blast disease. Through sequencing a recent field isolate, we report a reference genome that includes seven core chromosomes and mini-chromosome sequences that harbor effector genes normally found on ends of core chromosomes in other strains. No mini-chromosomes were observed in an early field strain, and at least two from another isolate each contain different effector genes and core chromosome end sequences. The mini-chromosome is enriched in transposons occurring most frequently at core chromosome ends. Additionally, transposons in mini-chromosomes lack the characteristic signature for inactivation by repeat-induced point (RIP) mutation genome defenses. Our results, collectively, indicate that dispensable mini-chromosomes and core chromosomes undergo divergent evolutionary trajectories, and mini-chromosomes and core chromosome ends are coupled as a mobile, fast-evolving effector compartment in the wheat pathogen genome. [ABSTRACT FROM AUTHOR]

Terpenes are responsible for most or all of the odor and flavor properties of Cannabis sativa, and may also impact effects users experience either directly or indirectly. We report the diversity of terpene profiles across samples bound for the Washington dispensary market. The remarkable degree of variation in terpene profiles ultimately results from action of a family of terpene synthase genes, only some of which have been described. Using a recently available genome assembly we describe 55 terpene synthases with genomic context, and tissue specific expression. The family is quite diverse from a protein similarity perspective, and subsets of the family are expressed in all tissues in the plant, including a set of root specific monoterpene synthases that could well have agronomic importance. Ultimately understanding and breeding for specific terpene profiles will require a good understanding of the gene family that underlies it. We intend for this work to serve as a foundation for that. [ABSTRACT FROM AUTHOR]

Mixed community or environmental DNA marker gene sequencing has become a commonly used technique for biodiversity analyses in freshwater systems. Many cytochrome c oxidase subunit I (COI) primer sets are now available for such work. The purpose of this study is to test whether COI primer choice affects the recovery of arthropod richness, beta diversity, and recovery of target assemblages in the benthos kick-net samples typically used in freshwater biomonitoring. We examine six commonly used COI primer sets on samples collected from six freshwater sites. Biodiversity analyses show that richness is sensitive to primer choice and the combined use of multiple COI amplicons recovers higher richness. Thus, to recover maximum richness, multiple primer sets should be used with COI metabarcoding. In ordination analyses based on community dissimilarity, samples consistently cluster by site regardless of amplicon choice or PCR replicate. Thus, for broadscale community analyses, overall beta diversity patterns are robust to COI marker choice. Recovery of traditional freshwater bioindicator assemblages such as Ephemeroptera, Trichoptera, Plectoptera, and Chironomidae as well as Arthropoda site indicators were differentially detected by each amplicon tested. This work will help future biodiversity and biomonitoring studies develop not just standardized, but optimized workflows that either maximize taxon-detection or the selection of amplicons for water quality or Arthropoda site indicators. [ABSTRACT FROM AUTHOR]

Advances in genomics have made whole genome studies increasingly feasible across the life sciences. However, new technologies and algorithmic advances do not guarantee flawless genomic sequences or annotation. Bias, errors, and artifacts can enter at any stage of the process from library preparation to annotation. When planning an experiment that utilizes a genome sequence as the basis for the design, there are a few basic checks that, if performed, may better inform the experimental design and ideally help avoid a failed experiment or inconclusive result. [ABSTRACT FROM AUTHOR]

Chinese cabbage (Brassica rapa ssp. pekinensis) is a major crop that is widely cultivated, especially in Korea, Japan, and China. With the advent of next generation sequencing technology, the cost and time required for sequencing have decreased and the development of genome research accelerated. Genome sequencing of Chinese cabbage was completed in 2011 using the variety Chiifu-401-42, and since then the genome has been continuously updated. In the present study, we conducted whole-genome sequencing of Chinese cabbage inbred line CT001, a line widely used in traditional or molecular breeding, to improve the accuracy of genetic polymorphism analysis. The constructed CT001 pseudomolecule represented 85.4% (219.8 Mb) of the Chiifu reference genome, and a total of 38,567 gene models were annotated using RNA-Seq analysis. In addition, the spontaneous mutation rate of CT001 was estimated by resequencing DNA obtained from individual plants after sexual propagation for six generations to estimate the naturally occurring variations. The CT001 pseudomolecule constructed in this study will provide valuable resources for genomic studies on Chinese cabbage. [ABSTRACT FROM AUTHOR]

Unknown sequences, or gaps, are present in many published genomes across public databases. Gap filling is an important finishing step in de novo genome assembly, especially in large genomes. The gap filling problem is nontrivial and while there are many computational tools partially solving the problem, several have shortcomings as to the reliability and correctness of the output, i.e. the gap filled draft genome. SSPACE-LongRead is a scaffolding tool that utilizes long reads from multiple third-generation sequencing platforms in finding links between contigs and combining them. The long reads potentially contain sequence information to fill the gaps created in the scaffolding, but SSPACE-LongRead currently lacks this functionality. We present an automated pipeline called gapFinisher to process SSPACE-LongRead output to fill gaps after the scaffolding. gapFinisher is based on the controlled use of a previously published gap filling tool FGAP and works on all standard Linux/UNIX command lines. We compare the performance of gapFinisher against two other published gap filling tools PBJelly and GMcloser. We conclude that gapFinisher can fill gaps in draft genomes quickly and reliably. In addition, the serial design of gapFinisher makes it scale well from prokaryote genomes to larger genomes with no increase in the computational footprint. [ABSTRACT FROM AUTHOR]

Plant mitochondrial genomes are usually assembled and displayed as circular maps based on the widely-held view across the broad community of life scientists that circular genome-sized molecules are the primary form of plant mitochondrial DNA, despite the understanding by plant mitochondrial researchers that this is an inaccurate and outdated concept. Many plant mitochondrial genomes have one or more pairs of large repeats that can act as sites for inter- or intramolecular recombination, leading to multiple alternative arrangements (isoforms). Most mitochondrial genomes have been assembled using methods unable to capture the complete spectrum of isoforms within a species, leading to an incomplete inference of their structure and recombinational activity. To document and investigate underlying reasons for structural diversity in plant mitochondrial DNA, we used long-read (PacBio) and short-read (Illumina) sequencing data to assemble and compare mitochondrial genomes of domesticated (Lactuca sativa) and wild (L. saligna and L. serriola) lettuce species. We characterized a comprehensive, complex set of isoforms within each species and compared genome structures between species. Physical analysis of L. sativa mtDNA molecules by fluorescence microscopy revealed a variety of linear, branched, and circular structures. The mitochondrial genomes for L. sativa and L. serriola were identical in sequence and arrangement and differed substantially from L. saligna, indicating that the mitochondrial genome structure did not change during domestication. From the isoforms in our data, we infer that recombination occurs at repeats of all sizes at variable frequencies. The differences in genome structure between L. saligna and the two other Lactuca species can be largely explained by rare recombination events that rearranged the structure. Our data demonstrate that representations of plant mitochondrial genomes as simple, circular molecules are not accurate descriptions of their true nature and that in reality plant mitochondrial DNA is a complex, dynamic mixture of forms. [ABSTRACT FROM AUTHOR]

Background: Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data. Methods and results: In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1. Conclusions: We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest. [ABSTRACT FROM AUTHOR]

Teleost fishes, thanks to their rapid evolution of sex determination mechanisms, provide remarkable opportunities to study the formation of sex chromosomes and the mechanisms driving the birth of new master sex determining (MSD) genes. However, the evolutionary interplay between the sex chromosomes and the MSD genes they harbor is rather unexplored. We characterized a male-specific duplicate of the anti-Müllerian hormone (amh) as the MSD gene in Northern Pike (Esox lucius), using genomic and expression evidence as well as by loss-of-function and gain-of-function experiments. Using RAD-Sequencing from a family panel, we identified Linkage Group (LG) 24 as the sex chromosome and positioned the sex locus in its sub-telomeric region. Furthermore, we demonstrated that this MSD originated from an ancient duplication of the autosomal amh gene, which was subsequently translocated to LG24. Using sex-specific pooled genome sequencing and a new male genome sequence assembled using Nanopore long reads, we also characterized the differentiation of the X and Y chromosomes, revealing a small male-specific insertion containing the MSD gene and a limited region with reduced recombination. Our study depicts reveals an unexpectedly low level of differentiation between a pair of sex chromosomes harboring an old MSD gene in a wild teleost fish population, and highlights both the pivotal role of genes from the amh pathway in sex determination, as well as the importance of gene duplication as a mechanism driving the turnover of sex chromosomes in this clade. [ABSTRACT FROM AUTHOR]

Long-read sequencing and novel long-range assays have revolutionized de novo genome assembly by automating the reconstruction of reference-quality genomes. In particular, Hi-C sequencing is becoming an economical method for generating chromosome-scale scaffolds. Despite its increasing popularity, there are limited open-source tools available. Errors, particularly inversions and fusions across chromosomes, remain higher than alternate scaffolding technologies. We present a novel open-source Hi-C scaffolder that does not require an a priori estimate of chromosome number and minimizes errors by scaffolding with the assistance of an assembly graph. We demonstrate higher accuracy than the state-of-the-art methods across a variety of Hi-C library preparations and input assembly sizes. The Python and C++ code for our method is openly available at . [ABSTRACT FROM AUTHOR]

Dengue fever is caused by dengue viruses (DENV) from the Flavivirus genus and is the most prevalent arboviral disease. DENV exists in four immunogenically distinct and genetically-related serotypes (DENV-1 to 4), each subdivided in genotypes. Despite the endemicity of all four DENV serotypes in Thailand, no prior study has characterized the circulation of DENV in the southern provinces of the country. To determine the genetic diversity of DENV circulating in Southern Thailand in 2015 and 2016, we investigated 46 viruses from 182 patients’ sera confirmed positive for DENV by serological and Nested RT-PCR tests. Our dataset included 2 DENV-1, 20 DENV-2, 9 DENV-3 and 15 DENV-4. Phylogenetic analysis was performed on viral envelop sequences. This revealed that part of the identified genotypes from DENV-1 and DENV-4 had been predominant in Asia (genotype I for both serotypes), while genotype II for DENV-4 and the Cosmopolitan genotype DENV-2 were also circulating. Whereas DENV-3 genotype II had been predominantly detected in South East Asia during the previous decades, we found genotype III and genotype I in Southern Thailand. All DENV genotype identified in this study were closely related to contemporary strains circulating in Southeast Asian countries, emphasizing the regional circulation of DENV. These results provide new insights into the co-circulation of all four DENV serotypes in Southern Thailand, confirming the hyperendemicity of DENV in the region. These findings also suggest a new trend of dissemination for some DENV serotypes with a possible shift in genotype distribution; as recently observed in other Asian countries. [ABSTRACT FROM AUTHOR]

The pepper weevil, Anthonomus eugenii, is a major pest on Capsicum species. Apart from natural spread, there is a risk of spread via international pepper trade. In the Netherlands, a pepper weevil outbreak occurred in 2012 and affected six greenhouses producing different sweet pepper varieties. The following year, a pepper weevil outbreak occurred in Italy. To trace the origin of the Dutch outbreak and to establish if the Dutch and Italian outbreaks were linked, we determined the mitogenomes of A. eugenii specimens collected at outbreak locations, and compared these with specimens from the native area, and other areas where the pest was introduced either by natural dispersal or via trade. The circular 17,257 bp A. eugenii mitogenome comprises thirteen mitochondrial genes typically found in insect species. Intra-species variation of these mitochondrial genes revealed four main mitochondrial lineages encompassing 41 haplotypes. The highest diversity was observed for specimens from its presumed native area (i.e. Mexico). The Dutch outbreak specimens represented three highly similar haplotypes, suggesting a single introduction of the pest. The major Dutch haplotype was also found in two specimens from New Jersey. As the Netherlands does not have pepper trade with New Jersey, it is likely that the specimens sampled in New Jersey and those sampled in the Netherlands originate from a shared source that was not included in this study. In addition, our analysis shows that the Italian and Dutch outbreaks were not linked. The mitochondrial genome is a useful tool to trace outbreak populations and the methodology presented in this paper could prove valuable for other invasive pest species, such as the African fruit moth Thaumatotibia leucotreta and emerald ash borer Agrilus planipennis. [ABSTRACT FROM AUTHOR]

Babesia bovis establishes persistent infections of long duration in cattle, despite the development of effective anti-disease immunity. One mechanism used by the parasite to achieve persistence is rapid antigenic variation of the VESA1 cytoadhesion ligand through segmental gene conversion (SGC), a phenomenon thought to be a form of homologous recombination (HR). To begin investigation of the enzymatic basis for SGC we initially identified and knocked out the Bbrad51 gene encoding the B. bovis Rad51 ortholog. BbRad51 was found to be non-essential for in vitro growth of asexual-stage parasites. However, its loss resulted in hypersensitivity to methylmethane sulfonate (MMS) and an apparent defect in HR. This defect rendered attempts to complement the knockout phenotype by reinsertion of the Bbrad51 gene into the genome unsuccessful. To circumvent this difficulty, we constructed an artificial chromosome, BbACc3, into which the complete Bbrad51 locus was inserted, for expression of BbRad51 under regulation by autologous elements. Maintenance of BbACc3 makes use of centromeric sequences from chromosome 3 and telomeric ends from chromosome 1 of the B. bovis C9.1 line. A selection cassette employing human dihydrofolate reductase enables recovery of transformants by selection with pyrimethamine. We demonstrate that the BbACc3 platform is stably maintained once established, assembles nucleosomes to form native chromatin, and expands in telomere length over time. Significantly, the MMS-sensitivity phenotype observed in the absence of Bbrad51 was successfully complemented at essentially normal levels. We provide cautionary evidence, however, that in HR-competent parasites BbACc3 can recombine with native chromosomes, potentially resulting in crossover. We propose that, under certain circumstances this platform can provide a useful alternative for the genetic manipulation of this group of parasites, particularly when regulated gene expression under the control of autologous elements may be important. [ABSTRACT FROM AUTHOR]

Castanopsis sclerophylla (Lindl.) Schott is an important species of evergreen broad-leaved tree in subtropical areas and has high ecological and economic value. However, there are few studies on its chloroplast genome. In this study, the complete chloroplast genome sequence of C. sclerophylla was determined using the Illumina HiSeq 2500 platform. The complete chloroplast genome of C. sclerophylla is 160,497 bp long, including a pair of inverted repeat (IR) regions (25,675 bp) separated by a large single-copy (LSC) region of 90,255 bp and a small single-copy (SSC) region of 18,892 bp. The overall GC content of the chloroplast genome is 36.82%. A total of 131 genes were found; of these, 111 genes are unique and annotated, including 79 protein-coding genes, 27 transfer RNA genes (tRNAs), and four ribosomal RNA genes (rRNAs). Twenty-one genes were found to be duplicated in the IR regions. Comparative analysis indicated that IR contraction might be the reason for the smaller chloroplast genome of C. sclerophylla compared to three congeneric species. Sequence analysis indicated that the LSC and SSC regions are more divergent than IR regions within Castanopsis; furthermore, greater divergence was found in noncoding regions than in coding regions. The maximum likelihood phylogenetic analysis showed that four species of the genus Castanopsis form a monophyletic clade and that C. sclerophylla is closely related to Castanopsis hainanensis with strong bootstrap values. These results not only provide a basic understanding of Castanopsis chloroplast genomes, but also illuminate Castanopsis species evolution within the Fagaceae family. Furthermore, these findings will be valuable for future studies of genetic diversity and enhance our understanding of the phylogenetic evolution of Castanopsis. [ABSTRACT FROM AUTHOR]

Whole genome sequencing can provide essential public health information. However, it is now known that widely used short-read methods have the potential to miss some randomly-distributed segments of genomes. This can prevent phages, plasmids, and virulence factors from being detected or properly identified. Here, we compared assemblies of three complete Shiga toxin-producing Escherichia coli (STEC) O26:H11/H- genomes from two different sequence types (ST21 and 29), each acquired using the Nextera XT MiSeq, MinION nanopore-based sequencing, and Pacific Biosciences (PacBio) sequencing. Each closed genome consisted of a single chromosome, approximately 5.7 Mb for CFSAN027343, 5.6 Mb for CFSAN027346, and 5.4 MB for CFSAN027350. However, short-read whole genome sequencing (WGS) using Nextera XT MiSeq failed to identify some virulence genes in plasmids and on the chromosome, both of which were detected using the long-read platforms. Results from long-read MinION and PacBio allowed us to identify differences in plasmid content: a single 88 kb plasmid in CFSAN027343; a 157kb plasmid in CFSAN027350; and two plasmids in CFSAN027346 (one 95 Kb, one 72 Kb). These data enabled rapid characterization of the virulome, detection of antimicrobial genes, and composition/location of Stx phages. Taken together, positive correlations between the two long-read methods for determining plasmids, virulome, antimicrobial resistance genes, and phage composition support MinION sequencing as one accurate and economical option for closing STEC genomes and identifying specific virulence markers. [ABSTRACT FROM AUTHOR]

Criniviruses accumulate in the phloem tissue and damage crops by reducing chlorophyll which is essential for plant growth and development. Tomato chlorosis crinivirus (ToCV) is vectored by several whitefly species that damage tomato crops throughout the world. In South Africa, ToCV is a poorly studied pathogen of global economic importance. Therefore, a national survey was initiated to investigate the occurrence and distribution of criniviruses infecting tomato crops in South Africa. Whitefly infested tomato crops exhibiting interveinal leaf chlorosis and chlorotic flecking symptoms were assayed for crinivirus infections using a multiplex reverse transcription polymerase reaction (RT-PCR) approach to assess for the presence of crinivirus species that are known to infect solanaceous hosts. Next-generation sequencing (NGS) was used to generate the complete genome of ToCV from South Africa. Results from the survey indicated that ToCV is presently the only crinivirus species infecting tomatoes in South Africa. Blast analysis showed that the RNA-1 segment of ToCV from South Africa (ToCR1-186) matched 99% to Spanish isolates. On the other hand, the RNA-2 (ToCR2-186) segment matched 98% to a South Korean isolate and three Spanish isolates. Although recombination events were not detected, phylogenetic studies showed inconsistencies in the grouping of RNA-1 and RNA-2 segments for some of the ToCV isolates analyzed in this study. Therefore, we suggest the possibility of intraspecific reassortment. This is the first comprehensive study and full genome sequence of ToCV from South Africa. The information generated from this study is intended to raise awareness of ToCV infections on tomato crops in South Africa. [ABSTRACT FROM AUTHOR]

Speciation genomic studies aim to interpret patterns of genome-wide variation in light of the processes that give rise to new species. However, interpreting the genomic “landscape” of speciation is difficult, because many evolutionary processes can impact levels of variation. Facilitated by the first chromosome-level assembly for the group, we use whole-genome sequencing and simulations to shed light on the processes that have shaped the genomic landscape during a radiation of monkeyflowers. After inferring the phylogenetic relationships among the 9 taxa in this radiation, we show that highly similar diversity (π) and differentiation (FST) landscapes have emerged across the group. Variation in these landscapes was strongly predicted by the local density of functional elements and the recombination rate, suggesting that the landscapes have been shaped by widespread natural selection. Using the varying divergence times between pairs of taxa, we show that the correlations between FST and genome features arose almost immediately after a population split and have become stronger over time. Simulations of genomic landscape evolution suggest that background selection (BGS; i.e., selection against deleterious mutations) alone is too subtle to generate the observed patterns, but scenarios that involve positive selection and genetic incompatibilities are plausible alternative explanations. Finally, tests for introgression among these taxa reveal widespread evidence of heterogeneous selection against gene flow during this radiation. Combined with previous evidence for adaptation in this system, we conclude that the correlation in FST among these taxa informs us about the processes contributing to adaptation and speciation during a rapid radiation. [ABSTRACT FROM AUTHOR]

Garden asparagus (Asparagus officinalis L.) is a popular vegetable cultivated worldwide. The secondary metabolites in its shoot are helpful for human health. We analyzed A. officinalis transcriptomes and identified differentially expressed genes (DEGs) involved in the biosynthesis of rutin and protodioscin, which are health-promoting functional compounds, and determined their association with stem color. We sequenced the complete mRNA transcriptome using the Illumina high-throughput sequencing platform in one white, three green, and one purple asparagus cultivars. A gene set was generated by de novo assembly of the transcriptome sequences and annotated using a BLASTx search. To investigate the relationship between the contents of rutin and protodioscin and their gene expression levels, rutin and protodioscin were analyzed using high-performance liquid chromatography. A secondary metabolite analysis using high-performance liquid chromatography showed that the rutin content was higher in green asparagus, while the protodioscin content was higher in white asparagus. We studied the genes associated with the biosynthesis of the rutin and protodioscin. The transcriptomes of the five cultivars generated 336 599 498 high-quality clean reads, which were assembled into 239 873 contigs with an average length of 694 bp, using the Trinity v2.4.0 program. The green and white asparagus cultivars showed 58 932 DEGs. A comparison of rutin and protodioscin biosynthesis genes revealed that 12 of the 57 genes associated with rutin and two of the 50 genes associated with protodioscin showed more than four-fold differences in expression. These DEGs might have caused a variation in the contents of these two metabolites between green and white asparagus. The present study is possibly the first to report transcriptomic gene sets in asparagus. The DEGs putatively involved in rutin and protodioscin biosynthesis might be useful for molecular engineering in asparagus. [ABSTRACT FROM AUTHOR]

Human Parainfluenza viruses (HPIV) type 1 and 3 are important causes of respiratory tract infections in young children globally. HPIV infections do not confer complete protective immunity so reinfections occur throughout life. Since no effective vaccine is available for the two virus subtypes, comprehensive understanding of HPIV-1 and HPIV-3 genetic and epidemic features is important for diagnosis, prevention, and treatment of HPIV-1 and HPIV-3 infections. Relatively few whole genome sequences are available for both HPIV-1 and HPIV-3 viruses, so our study sought to provide whole genome sequences from multiple countries to further the understanding of the global diversity of HPIV at a whole-genome level. We collected HPIV-1 and HPIV-3 samples and isolates from Argentina, Australia, France, Mexico, South Africa, Switzerland, and USA from the years 2003–2011 and sequenced the genomes of 40 HPIV-1 and 75 HPIV-3 viruses with Sanger and next-generation sequencing with the Ion Torrent, Illumina, and 454 platforms. Phylogenetic analysis showed that the HPIV-1 genome is evolving at an estimated rate of 4.97 × 10−4 mutations/site/year (95% highest posterior density 4.55 × 10−4 to 5.38 × 10−4) and the HPIV-3 genome is evolving at a similar rate (3.59 × 10−4 mutations/site/year, 95% highest posterior density 3.26 × 10−4 to 3.94 × 10−4). There were multiple genetically distinct lineages of both HPIV-1 and 3 circulating on a global scale. Further surveillance and whole-genome sequencing are greatly needed to better understand the spatial dynamics of these important respiratory viruses in humans. [ABSTRACT FROM AUTHOR]

Gene manipulation and epitope tagging are essential tools for understanding the molecular function of specific genes. The opportunistic human pathogen Candida albicans is a diploid fungus that utilizes a non-canonical genetic code. Since selection markers available in this organism are scarce, several tools based on recyclable markers have been developed for gene disruption, such as the Clox system. This system relies on the Cre recombinase, which recycles selection markers flanked by loxP sites with high efficiency, facilitating single marker or multi-marker recycling. However, PCR-based modules for epitope tagging, such the pFA-modules, mainly use limited non-recyclable auxotrophic markers. To solve this problem, we have used a Gibson assembly strategy to construct a set of new plasmids where the auxotrophic markers of the pFA vectors were swapped with five recyclable marker modules of the Clox system, enhancing the versatility of the pFA plasmids. This new toolkit, named pFA-Clox, is composed of 36 new vectors for gene disruption and epitope tagging (GFP, 3xGFP, mCherry, 3xHA, 5xmyc and TAP). These plasmids contain the dominant NAT1 marker, as well as URA3, HIS1 and ARG4 cassettes, thereby permitting functional analysis of laboratory strains as well as clinical isolates of C. albicans. In summary, we have adapted the Clox system to the pFA-backbone vectors. Thus, the set of primers used for the amplification of previously published pFA modules can also be utilized in this new pFA-Clox system. Therefore, this new toolkit harbors the advantages of both systems, allowing accelerated gene modification strategies that could reduce time and costs in strain construction for C. albicans. [ABSTRACT FROM AUTHOR]

Larix populations at the tundra-taiga ecotone in northern Siberia are highly under-represented in population genetic studies, possibly due to the remoteness of these regions that can only be accessed at extraordinary expense. The genetic signatures of populations in these boundary regions are therefore largely unknown. We aim to generate organelle reference genomes for the detection of single nucleotide polymorphisms (SNPs) that can be used for paleogenetic studies. We present 19 complete chloroplast genomes and mitochondrial genomic sequences of larches from the southern lowlands of the Taymyr Peninsula (northernmost range of Larix gmelinii (Rupr.) Kuzen.), the lower Omoloy River, and the lower Kolyma River (both in the range of Larix cajanderi Mayr). The genomic data reveal 84 chloroplast SNPs and 213 putatively mitochondrial SNPs. Parsimony-based chloroplast haplotype networks show no spatial structure of individuals from different geographic origins, while the mitochondrial haplotype network shows at least a slight spatial structure with haplotypes from the Omoloy and Kolyma populations being more closely related to each other than to most of the haplotypes from the Taymyr populations. Whole genome alignments with publicly available complete chloroplast genomes of different Larix species show that among official plant barcodes only the rcbL gene contains sufficient polymorphisms, but has to be sequenced completely to distinguish the different provenances. We provide 8 novel mitochondrial SNPs that are putatively diagnostic for the separation of L. gmelinii and L. cajanderi, while 4 chloroplast SNPs have the potential to distinguish the L. gmelinii/L. cajanderi group from other Larix species. Our organelle references can be used for a targeted primer and probe design allowing the generation of short amplicons. This is particularly important with regard to future investigations of, for example, the biogeographic history of Larix by screening ancient sedimentary DNA of Larix. [ABSTRACT FROM AUTHOR]

A new method is introduced allowing seamless assembly of independent, functionally tested, blunt-end double strand nucleic acid (DNA fragments not supplied in vectors such as plasmids) into more complex biological (e.g. protein expression vectors) and higher order multi-device (e.g. biochemical pathways). Individual parts include bacterial selection markers and origins of replication, promoters useful in a variety of species, transcription terminators, shuttle sequences and a variety of “N” and “C” terminal solubility/affinity protein tags. Parts are not subjected to pre-assembly manipulation with nucleic acid modifying enzymes. Instead, they are simply mixed in appropriate pre-defined combinations and concentrations and then seamlessly linked into devices employing a specialized thermostable enzyme blend. Combinatorial assembly of parts is an inherent time-saving feature of the new method, in contrast to hierarchical binary assembly (“one part at a time”) methods. This feature substantially simplifies and speeds optimization of device and system development. The versatility and functionality of the new method was shown by combinatorial assembly of parts into vector devices, one of which optimally expressed protein from a model gene. Also, a four-enzyme biosynthetic pathway system was re-created by combinatorial construction from parts and devices. Concepts discussed in this paper provide synthetic biologists, chemists and bio-engineers with improved and expanded capability to create novel biological molecules and systems. [ABSTRACT FROM AUTHOR]

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of these shifts remain largely elusive. Here, we used classic quantitative trait analysis, combined with genomic and transcriptomic information to dissect the genetic basis of asexual, parthenogenetic reproduction in the brown alga Ectocarpus. We found that parthenogenesis is controlled by the sex locus, together with two additional autosomal loci, highlighting the key role of the sex chromosome as a major regulator of asexual reproduction. We identify several negative effects of parthenogenesis on male fitness, and different fitness effects of parthenogenetic capacity depending on the life cycle generation. Although allele frequencies in natural populations are currently unknown, we discuss the possibility that parthenogenesis may be under both sex-specific selection and generation/ploidally-antagonistic selection, and/or that the action of fluctuating selection on this trait may contribute to the maintenance of polymorphisms in populations. Importantly, our data provide the first empirical illustration, to our knowledge, of a trade-off between the haploid and diploid stages of the life cycle, where distinct parthenogenesis alleles have opposing effects on sexual and asexual reproduction and may help maintain genetic variation. These types of fitness trade-offs have profound evolutionary implications in natural populations and may structure life history evolution in organisms with haploid-diploid life cycles. [ABSTRACT FROM AUTHOR]

In 2010, the rare OP354-like P[8]b rotavirus subtype was detected in children less than 2 years old in Ghana. In this follow-up study, to provide insight into the evolutionary history of the genome of Ghanaian P[8]b strains RVA/Human-wt/GHA/GHDC949/2010/G9P[8] and RVA/Human-wt/GHA/GHM0094/2010/G9P[8] detected in an infant and a 7-month old child hospitalised for acute gastroenteritis, we sequenced the complete genome using both Sanger sequencing and Illumina MiSeq technology followed by phylogenetic analysis of the near-full length sequences. Both strains possessed the Wa-like/genotype 1 constellation G9P[8]b-I1-R1-C1-M1-A1-N1-T1-E1-H1. Sequence comparison and phylogenetic inference showed that both strains were identical at the lineage level throughout the 11 genome segments. Their VP7 sequences belonged to the major sub-lineage of the G9-lineage III whereas their VP4 sequences belonged to P[8]b cluster I. The VP7 and VP4 genes of the study strains were closely related to a Senegalese G9P[8]b strain detected in 2009. In the remaining nine genome segments, both strains consistently clustered together with Wa-like RVA strains possessing either P[8]a or P[8]b mostly of African RVA origin. The introduction of a P[8]b subtype VP4 gene into the stable Wa-like strain backbone may result in strains that might propagate easily in the human population, with a potential to become an important public health concern, especially because it is not certain if the monovalent rotavirus vaccine (Rotarix) used in Ghana will be efficacious against such strains. Our analysis of the full genomes of GHM0094 and GHDC949 adds to knowledge of the genetic make-up and evolutionary dynamics of P[8]b rotavirus strains. [ABSTRACT FROM AUTHOR]

Background: Global estimates for cholera annually approximate 4 million cases worldwide with 95,000 deaths. Recent outbreaks, including Haiti and Yemen, are reminders that cholera is still a global health concern. Cholera outbreaks can rapidly induce high death tolls by overwhelming the capacity of health facilities, especially in remote areas or areas of civil unrest. Recent studies demonstrated that stool specimens preserved on filter paper facilitate molecular analysis of Vibrio cholerae in resource limited settings. Specimens preserved in a rapid, low-cost, safe and sustainable manner for sequencing provides previously unavailable data about circulating cholera strains. This may ultimately provide new information to shape public policy response on cholera control and elimination. Methodology/Principal findings: Whole genome sequencing (WGS) recovered close to a complete sequence of the V. cholerae O1 genome with satisfactory genome coverage from stool specimens enriched in alkaline peptone water (APW) and V. cholerae culture isolates, both spotted on filter paper. The minimum concentration of V. cholerae DNA sufficient to produce quality genomic information was 0.02 ng/μL. The genomic data confirmed the presence or absence of genes of epidemiological interest, including cholera toxin and pilus loci. WGS identified a variety of diarrheal pathogens from APW-enriched specimen spotted filter paper, highlighting the potential for this technique to explore the gut microbiome, potentially identifying co-infections, which may impact the severity of disease. WGS demonstrated that these specimens fit within the current global cholera phylogenetic tree, identifying the strains as the 7th pandemic El Tor. Conclusions: WGS results allowed for mapping of short reads from APW-enriched specimen and culture isolate spotted filter papers this provided valuable molecular epidemiological sequence information on V. cholerae strains from remote, low-resource settings. These results identified the presence of co-infecting pathogens while providing rare insight into the specific V. cholerae strains causing outbreaks in cholera-endemic areas. [ABSTRACT FROM AUTHOR]

COLOSTRONONI is a dietary supplement consisting of bovine colostrum and Morinda citrifolia fruit (Noni). In this study, we tested the capability of COLOSTRONONI to influence gut microbiota composition using an in vivo evaluation in rats. Furthermore, we analyzed the effect of COLOSTRONONI on the systemic inflammatory responses as well as on the gut permeability of the animals. Altogether, our analyses supported the concept of COLOSTRONONI as a natural food supplement that doesn't affect (neither negatively nor positively) gut microbiota homeostasis in healthy conditions. Moreover, COLOSTRONONI highlighted a lower effect in the expression of genes coding for IL-10, Il-12 and TNF-α response allowing us to hypothesize an immunomodulatory activity of this dietary supplement. [ABSTRACT FROM AUTHOR]

The Oxford Nanopore MinION DNA sequencing device can produce large amounts of long sequences, typically several kilobases, within a few hours. This long read capacity was exploited to detect antimicrobial resistance genes (ARGs) in a large veterinary teaching hospital environment, and to assess their taxonomic origin, genetic organisation and association with mobilisation markers concurrently. Samples were collected on eight occasions between November 2016 and May 2017 (inclusive) in a longitudinal study. Nanopore sequencing was performed on total DNA extracted from the samples after a minimal enrichment step in broth. Many ARGs present in the veterinary hospital environment could potentially confer resistance to antimicrobials widely used in treating infections of companion animals, including aminoglycosides, extended-spectrum beta-lactams, sulphonamides, macrolides, and tetracyclines. High-risk ARGs, defined here as single or multiple ARGs associated with pathogenic bacterial species or with mobile genetic elements, were shared between the intensive care unit (ICU) patient cages, a dedicated laundry trolley and a floor cleaning mop-bucket. By contrast, a floor surface from an office corridor without animal contact and located outside the veterinary hospital did not contain such high-risk ARGs. Relative abundances of high-risk ARGs and co-localisation of these genes on the same sequence read were higher in the laundry trolley and mop bucket samples, compared to the ICU cages, suggesting that amplification of ARGs is likely to occur in the collection points for hospital waste. These findings have prompted the implementation of targeted intervention measures in the veterinary hospital to mitigate the risks of transferring clinically important ARGs between sites and to improve biosecurity practices in the facility. [ABSTRACT FROM AUTHOR]

HIV-1 strain diversity in Bulgaria is extensive and includes contributions from nearly all major subtypes and the Circulating Recombinant Forms (CRF): 01_AE, 02_AG, and 05_DF. Prior to this study, HIV-1 sequence information from Bulgaria has been based solely on the pro-RT gene, which represent less than 15% of the viral genome. To further characterize HIV-1 in Bulgaria, assess participant risk behaviors, and strengthen knowledge of circulating strains in the region, the study “Genetic Subtypes of HIV-1 in Bulgaria (RV240)” was conducted. This study employed the real time-PCR based Multi-region Hybridization Assay (MHA) B/non-B and HIV-1 sequencing to survey 215 of the approximately 1100 known HIV-1 infected Bulgarian adults (2008–2009) and determine if they were infected with subtype B HIV-1. The results indicated a subtype B prevalence of 40% and demonstrate the application of the MHA B/non-B in an area containing broad HIV-1 strain diversity. Within the assessed risk behaviors, the proportion of subtype B infection was greatest in men who have sex with men and lowest among those with drug use risk factors. During this study, 15 near full-length genomes and 22 envelope sequences were isolated from study participants. Phylogenetic analysis shows the presence of subtypes A1, B, C, F1, and G, CRF01_AE, CRF02_AG, CRF05_DF, and one unique recombinant form (URF). These sequences also show the presence of two strain groups containing participants with similar risk factors. Previous studies in African and Asian cohorts have shown that co-circulation of multiple subtypes can lead to viral recombination within super-infected individuals and the emergence of new URFs. The low prevalence of URFs in the presence of high subtype diversity in this study, may be the result of successful infection prevention and control programs. Continued epidemiological monitoring and support of infection prevention programs will help maintain control of the HIV-1 epidemic in Bulgaria. [ABSTRACT FROM AUTHOR]

In addition to the BAC-based reference sequence of the accession Columbia-0 from the year 2000, several short read assemblies of THE plant model organism Arabidopsis thaliana were published during the last years. Also, a SMRT-based assembly of Landsberg erecta has been generated that identified translocation and inversion polymorphisms between two genotypes of the species. Here we provide a chromosome-arm level assembly of the A. thaliana accession Niederzenz-1 (AthNd-1_v2c) based on SMRT sequencing data. The best assembly comprises 69 nucleome sequences and displays a contig length of up to 16 Mbp. Compared to an earlier Illumina short read-based NGS assembly (AthNd-1_v1), a 75 fold increase in contiguity was observed for AthNd-1_v2c. To assign contig locations independent from the Col-0 gold standard reference sequence, we used genetic anchoring to generate a de novo assembly. In addition, we assembled the chondrome and plastome sequences. Detailed analyses of AthNd-1_v2c allowed reliable identification of large genomic rearrangements between A. thaliana accessions contributing to differences in the gene sets that distinguish the genotypes. One of the differences detected identified a gene that is lacking from the Col-0 gold standard sequence. This de novo assembly extends the known proportion of the A. thaliana pan-genome. [ABSTRACT FROM AUTHOR]

Recent data support that the vaginal microbiota may alter mucosal pharmacokinetics (PK) of topically delivered microbicides. Our team developed an intravaginal ring (IVR) that delivers tenofovir (TFV) (8–10 mg/day) alone or with levonorgestrel (LNG) (20 ug/day). We evaluated the effect of IVRs on the vaginal microbiota, and describe how the vaginal microbiota impacts mucosal PK of TFV. CONRAD A13-128 was a randomized, placebo controlled phase I study. We randomized 51 women to TFV, TFV/LNG or placebo IVR. We assessed the vaginal microbiota by sequencing the V3-V4 regions of 16S rRNA genes prior to IVR insertion and after approximately 15 days of use. We measured the concentration of TFV in the cervicovaginal (CV) aspirate, and TFV and TFV-diphosphate (TFV-DP) in vaginal tissue at the end of IVR use. The change in relative or absolute abundance of vaginal bacterial phylotypes was similar among active and placebo IVR users (all q values >0.13). TFV concentrations in CV aspirate and vaginal tissue, and TFV-DP concentrations in vaginal tissue were not significantly different among users with community state type (CST) 4 versus those with Lactobacillus dominated microbiota (all p values >0.07). The proportions of participants with CV aspirate concentrations of TFV >200,000 ng/mL and those with tissue TFV-DP concentrations >1,000 fmol/mg were similar among women with anaerobe versus Lactobacillus dominated microbiota (p = 0.43, 0.95 respectively). There were no significant correlations between the CV aspirate concentration of TFV and the relative abundances of Gardnerella vaginalis or Prevotella species. Tissue concentrations of TFV-DP did not correlate with any the relative abundances of any species, including Gardnerella vaginalis. In conclusion, active IVRs did not differ from the placebo IVR on the effect on the vaginal microbiota. Local TFV and TFV-DP concentrations were high and similar among IVR users with Lactobacillus dominated microbiota versus CST IV vaginal microbiota. Trial registration: ClinicalTrials.gov . [ABSTRACT FROM AUTHOR]

Microbial colonization can be detrimental to the integrity of metal surfaces and lead to microbiologically influenced corrosion (MIC). Biocorrosion is a serious problem for aquatic and marine industries in the world. In Minnesota (USA), where this study was conducted, biocorrosion severely affects the maritime transportation industry. The anticorrosion activity of a variety of compounds, including chemical (magnesium peroxide) and biological (surfactin, capsaicin, and gramicidin) molecules were investigated as coating additives. We also evaluated a previously engineered, extremely stable, non-biocidal enzyme known to interfere in bacterial signaling, SsoPox (a quorum quenching lactonase). Experimental steel coupons were submerged in water from the Duluth Superior Harbor (DSH) for 8 weeks in the laboratory. Biocorrosion was evaluated by counting the number and the coverage of corrosion tubercles on coupons and also by ESEM imaging of the coupon surface. Three experimental coating additives significantly reduced the formation of corrosion tubercles: surfactin, magnesium peroxide and the quorum quenching lactonase by 31%, 36% and 50%, respectively. DNA sequence analysis of the V4 region of the bacterial 16S rRNA gene revealed that these decreases in corrosion were associated with significant changes in the composition of bacterial communities on the steel surfaces. These results demonstrate the potential of highly stable quorum quenching lactonases to provide a reliable, cost-effective method to treat steel structures and prevent biocorrosion. [ABSTRACT FROM AUTHOR]

Centromeres are essential chromosomal regions that mediate kinetochore assembly and spindle attachments during cell division. Despite their functional conservation, centromeres are among the most rapidly evolving genomic regions and can shape karyotype evolution and speciation across taxa. Although significant progress has been made in identifying centromere-associated proteins, the highly repetitive centromeres of metazoans have been refractory to DNA sequencing and assembly, leaving large gaps in our understanding of their functional organization and evolution. Here, we identify the sequence composition and organization of the centromeres of Drosophila melanogaster by combining long-read sequencing, chromatin immunoprecipitation for the centromeric histone CENP-A, and high-resolution chromatin fiber imaging. Contrary to previous models that heralded satellite repeats as the major functional components, we demonstrate that functional centromeres form on islands of complex DNA sequences enriched in retroelements that are flanked by large arrays of satellite repeats. Each centromere displays distinct size and arrangement of its DNA elements but is similar in composition overall. We discover that a specific retroelement, G2/Jockey-3, is the most highly enriched sequence in CENP-A chromatin and is the only element shared among all centromeres. G2/Jockey-3 is also associated with CENP-A in the sister species D. simulans, revealing an unexpected conservation despite the reported turnover of centromeric satellite DNA. Our work reveals the DNA sequence identity of the active centromeres of a premier model organism and implicates retroelements as conserved features of centromeric DNA. [ABSTRACT FROM AUTHOR]

The virulence plasmid pSLT as exemplified by the 94 Kb plasmid in Salmonella Typhimurium strain LT2 is only found in isolates of serovar Typhimurium. While it occurs commonly among such isolates recent genotyping methods have shown that it is mostly confined to certain genotypes. Although pSLT plasmids are capable of self-transmissibility under experimental conditions their confinement to certain host genotypes suggests that in practice they are maintained by vertical rather than by horizontal transmission. This would imply that evolution of the pSLT plasmid proceeds in parallel with evolution of its host. The development of a phylogenetic evolutionary framework for genotypes of S. Typhimurium based on single-nucleotide-polymorphism (SNPs) typing provided an opportunity to test whether the pSLT plasmid coevolves with its host genotype. Accordingly SNPs analysis was applied to the pSLT plasmids from 71 strains S. Typhimurium of Australian and international origins representing most of the genotypes which commonly have a pSLT. The phylogenetic tree showed that pSLT sequences clustered into almost the same groups as the host chromosomes so that each pSLT genotype was associated with a single host genotype. A search for tandem repeats in pSLT plasmids showed that a 9 bp VNTR in the traD gene occurred in the pSLT from all isolates belonging to Clade II but not from isolates belonging to Clade I. Another 9 bp repeat occurred only in three Clade I genotypes with a recent common ancestor. The evidence relating to both of these VNTRs supports the proposition that the pSLT plasmid is only transmitted vertically. Some isolates belonging to one S. Typhimurium genotype were found to have pSLTs which have lost a large block of genes when a resistance gene cassette has been acquired. Examples were found of pSLT plasmids which have recombined with other plasmids to form fusion plasmids sometimes with loss of some pSLT genes. In all cases the underlying genotype of the modified pSLT was the same as the genotype of regular pSLTs with the same host genotype implying that these changes have occurred within the host cell of the pSLT plasmid. [ABSTRACT FROM AUTHOR]

Larvae of the pest Protaetia brevitarsis are used to treat infections in traditional Chinese medicine. However, genomic information about this non-model species is currently lacking. To better understand the fundamental biology of this non-model species, its transcriptome was obtained using next generation sequencing and then analyzed. A total of 7.62 Gb of clean reads were obtained, which were assembled into 169,087 transcripts corresponding to 142,000 annotated unigenes. These unigenes were functionally classified according to Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. A total of 41,921 unigenes were assigned to 56 GO terms, 21,454 unigenes were divided among 26 KOG categories, and 16,368 unigenes were assigned to 32 KEGG pathways. In addition, 19,144 simple sequence repeats (SSRs) were identified. Furthermore, several kinds of natural antimicrobial peptides and proteins, 4 histones with potential antimicrobial activity, and 41 potential antimicrobial peptide sequences were identified. These data are the first reported whole transcriptome sequence of P. brevitarsis larvae, which represents a valuable genomic resource for studying this species, thus promoting the utilization of its medical potential. [ABSTRACT FROM AUTHOR]

Mussels of the genus Bathymodiolus are among the most widespread colonizers of hydrothermal vent and cold seep environments, sustained by endosymbiosis with chemosynthetic bacteria. Presumed species of Bathymodiolus are abundant at newly discovered cold seeps on the Mid-Atlantic continental slope, however morphological taxonomy is challenging, and their phylogenetic affinities remain unestablished. Here we used mitochondrial sequence to classify species found at three seep sites (Baltimore Canyon seep (BCS; ~400m); Norfolk Canyon seep (NCS; ~1520m); and Chincoteague Island seep (CTS; ~1000m)). Mitochondrial COI (N = 162) and ND4 (N = 39) data suggest that Bathymodiolus childressi predominates at these sites, although single B. mauritanicus and B. heckerae individuals were detected. As previous work had suggested that methanotrophic and thiotrophic interactions can both occur at a site, and within an individual mussel, we investigated the symbiont communities in gill tissues of a subset of mussels from BCS and NCS. We constructed metabarcode libraries with four different primer sets spanning the 16S gene. A methanotrophic phylotype dominated all gill microbial samples from BCS, but sulfur-oxidizing Campylobacterota were represented by a notable minority of sequences from NCS. The methanotroph phylotype shared a clade with globally distributed Bathymodiolus spp. symbionts from methane seeps and hydrothermal vents. Two distinct Campylobacterota phylotypes were prevalent in NCS samples, one of which shares a clade with Campylobacterota associated with B. childressi from the Gulf of Mexico and the other with Campylobacterota associated with other deep-sea fauna. Variation in chemosynthetic symbiont communities among sites and individuals has important ecological and geochemical implications and suggests shifting reliance on methanotrophy. Continued characterization of symbionts from cold seeps will provide a greater understanding of the ecology of these unique environments as well and their geochemical footprint in elemental cycling and energy flux. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) are key players in regulation of gene expression at post-transcription level in eukaryotic cells. MiRNAs have been intensively studied in plants, animals and viruses. The investigations of bacterial miRNAs have gained less attention, except for the recent studies on miRNAs derived from Streptococcus mutans ATCC 25175 and Escherichia coli DH10B. In this study, high-throughput sequencing approach was employed to investigate the miRNA population in bacteria of the genus Thalassospira using both the miRDeep2 and CID-miRNA methods. A total of 984 putative miRNAs were identified in 9 species of the genus Thalassospira using both miRDeep and CID-miRNA analyses. Fifty seven conserved putative miRNAs were found in different species of the genus Thalassospira, and up to 6 miRNAs were found to be present at different locations in the T. alkalitolerans JCM 18968T, T. lucentensis QMT2T and T. xianhensis P-4T. None of the putative miRNAs was found to share sequence to the reported miRNAs in E. coli DH10B and S. mutans ATCC 25175. The findings provide a comprehensive list of computationally identified miRNAs in 9 bacterial species of the genus Thalassospira and addressed the existing knowledge gap on the presence of miRNAs in the Thalassospira genomes. [ABSTRACT FROM AUTHOR]

Amoebophrya is part of an enigmatic, diverse, and ubiquitous marine alveolate lineage known almost entirely from anonymous environmental sequencing. Two cultured Amoebophrya strains grown on core dinoflagellate hosts were used for transcriptome sequencing. BLASTx using different genetic codes suggests that Amoebophyra sp. ex Karlodinium veneficum uses the three typical stop codons (UAA, UAG, and UGA) to encode amino acids. When UAA and UAG are translated as glutamine about half of the alignments have better BLASTx scores, and when UGA is translated as tryptophan one fifth have better scores. However, the sole stop codon appears to be UGA based on conserved genes, suggesting contingent translation of UGA. Neither host sequences, nor sequences from the second strain, Amoebophrya sp. ex Akashiwo sanguinea had similar results in BLASTx searches. A genome survey of Amoebophyra sp. ex K. veneficum showed no evidence for transcript editing aside from mitochondrial transcripts. The dynein heavy chain (DHC) gene family was surveyed and of 14 transcripts only two did not use UAA, UAG, or UGA in a coding context. Overall the transcriptome displayed strong bias for A or U in third codon positions, while the tRNA genome survey showed bias against codons ending in U, particularly for amino acids with two codons ending in either C or U. Together these clues suggest contingent translation mechanisms in Amoebophyra sp. ex K. veneficum and a phylogenetically distinct instance of genetic code modification. [ABSTRACT FROM AUTHOR]

Many plants have natural partnerships with microbes that can boost their nitrogen (N) and/or phosphorus (P) acquisition. To assess whether wheat may have undiscovered associations of these types, we tested if N/P-starved Triticum aestivum show microbiome profiles that are simultaneously different from those of N/P-amended plants and those of their own bulk soils. The bacterial and fungal communities of root, rhizosphere, and bulk soil samples from the Historical Dryland Plots (Lethbridge, Canada), which hold T. aestivum that is grown both under N/P fertilization and in conditions of extreme N/P-starvation, were taxonomically described and compared (bacterial 16S rRNA genes and fungal Internal Transcribed Spacers—ITS). As the list may include novel N- and/or P-providing wheat partners, we then identified all the operational taxonomic units (OTUs) that were proportionally enriched in one or more of the nutrient starvation- and plant-specific communities. These analyses revealed: a) distinct N-starvation root and rhizosphere bacterial communities that were proportionally enriched, among others, in OTUs belonging to families Enterobacteriaceae, Chitinophagaceae, Comamonadaceae, Caulobacteraceae, Cytophagaceae, Streptomycetaceae, b) distinct N-starvation root fungal communities that were proportionally enriched in OTUs belonging to taxa Lulworthia, Sordariomycetes, Apodus, Conocybe, Ascomycota, Crocicreas, c) a distinct P-starvation rhizosphere bacterial community that was proportionally enriched in an OTU belonging to genus Agrobacterium, and d) a distinct P-starvation root fungal community that was proportionally enriched in OTUs belonging to genera Parastagonospora and Phaeosphaeriopsis. Our study might have exposed wheat-microbe connections that can form the basis of novel complementary yield-boosting tools. [ABSTRACT FROM AUTHOR]

Retrotransposon expression during arbuscular mycorrhizal (AM) fungal colonisation of sunflower roots (Helianthus annuus) was analysed using Illumina RNA-Seq, in order to verify whether mycorrhizal symbiosis can activate retrotransposable elements. Illumina cDNA libraries were produced from RNAs isolated from the roots of sunflower plants at 4 and 16 days after inoculation with the AM fungus Rhizoglomus irregulare and from their respective control plants. Illumina reads were mapped to a library of reverse transcriptase-encoding sequences, putatively belonging to long terminal repeat retrotransposons of Gypsy and Copia superfamilies. Forty-six different reverse transcriptase sequences were transcribed, although at a low rate, in mycorrhizal or control roots and only four were significantly over-expressed at day 16, compared with control roots. Almost all expressed or over-expressed sequences belonged to low-copy elements, mostly, of the Copia superfamily. A meta-analysis, using publicly available Illumina cDNA libraries obtained from sunflower plants treated with different hormones and chemicals, mimicking stimuli produced by abiotic and biotic stresses, was also conducted. Such analyses indicated that the four reverse transcriptase sequences over-expressed in mycorrhizal roots were explicitly induced only by AM symbiosis, showing the specificity of AM stimuli compared to that of other fungal/plant interactions. [ABSTRACT FROM AUTHOR]

The diversity of viruses probably exceeds biodiversity of eukaryotes, but little is known about the origin and emergence of novel virus species. Experimentation and disease outbreak investigations have allowed the characterization of rapid molecular virus adaptation. However, the processes leading to the establishment of functionally distinct virus taxa in nature remain obscure. Here, we demonstrate that incipient speciation in a natural host species has generated distinct ecological niches leading to adaptive isolation in an RNA virus. We found a very strong association between the distributions of two major phylogenetic clades in Tula orthohantavirus (TULV) and the rodent host lineages in a natural hybrid zone of the European common vole (Microtus arvalis). The spatial transition between the virus clades in replicated geographic clines is at least eight times narrower than between the hybridizing host lineages. This suggests a strong barrier for effective virus transmission despite frequent dispersal and gene flow among local host populations, and translates to a complete turnover of the adaptive background of TULV within a few hundred meters in the open, unobstructed landscape. Genetic differences between TULV clades are homogenously distributed in the genomes and mostly synonymous (93.1%), except for a cluster of nonsynonymous changes in the 5′ region of the viral envelope glycoprotein gene, potentially involved in host-driven isolation. Evolutionary relationships between TULV clades indicate an emergence of these viruses through rapid differential adaptation to the previously diverged host lineages that resulted in levels of ecological isolation exceeding the progress of speciation in their vertebrate hosts. [ABSTRACT FROM AUTHOR]

Heterogeneity is prevalent in cancer both between and within individuals. Although a few studies have identified several circulating microRNAs (miRNAs) for cancer diagnosis, the complete plasma miRNA profile for hepatocellular carcinoma (HCC) remains undefined, and whether the plasma miRNA profiles are heterogeneous is unknown. Here, we obtained individualized plasma miRNA profiles of both healthy subjects and HCC patients via genome-wide deep sequencing. Compared with the highly stable miRNA profile of the healthy subjects, the profile of the HCC patients was highly variable. Seven miRNAs were optimized as potential plasma-based biomarkers for HCC diagnosis. Combined with the clinical data of The Cancer Genome Atlas (TCGA) cohort, three out of the seven miRNAs were correlated with the survival of the HCC patients. To investigate the effect of cancer cells on the plasma miRNAs profile, we compared the most differentially expressed miRNAs between plasma and tissues. Furthermore, miRNAseq data of HCC patients from TCGA were recruited for comparisons. We found that the differences between plasma and tissue were inconsistent, suggesting that other cells in addition to cancer cells also contribute to plasma miRNAs. Using two HCC cancer cell lines, we examined the levels of seven differentially expressed miRNAs. The reverse direction of certain miRNAs alterations between cancer cells and media further confirmed that miRNAs may be selectively pump out by cancer cells. [ABSTRACT FROM AUTHOR]

The bacterium Agrobacterium tumefaciens has been the workhorse in plant genome engineering. Customized replacement of native tumor-inducing (Ti) plasmid elements enabled insertion of a sequence of interest called Transfer-DNA (T-DNA) into any plant genome. Although these transfer mechanisms are well understood, detailed understanding of structure and epigenomic status of insertion events was limited by current technologies. Here we applied two single-molecule technologies and analyzed Arabidopsis thaliana lines from three widely used T-DNA insertion collections (SALK, SAIL and WISC). Optical maps for four randomly selected T-DNA lines revealed between one and seven insertions/rearrangements, and the length of individual insertions from 27 to 236 kilobases. De novo nanopore sequencing-based assemblies for two segregating lines partially resolved T-DNA structures and revealed multiple translocations and exchange of chromosome arm ends. For the current TAIR10 reference genome, nanopore contigs corrected 83% of non-centromeric misassemblies. The unprecedented contiguous nucleotide-level resolution enabled an in-depth study of the epigenome at T-DNA insertion sites. SALK_059379 line T-DNA insertions were enriched for 24nt small interfering RNAs (siRNA) and dense cytosine DNA methylation, resulting in transgene silencing via the RNA-directed DNA methylation pathway. In contrast, SAIL_232 line T-DNA insertions are predominantly targeted by 21/22nt siRNAs, with DNA methylation and silencing limited to a reporter, but not the resistance gene. Additionally, we profiled the H3K4me3, H3K27me3 and H2A.Z chromatin environments around T-DNA insertions using ChIP-seq in SALK_059379, SAIL_232 and five additional T-DNA lines. We discovered various effect s ranging from complete loss of chromatin marks to the de novo incorporation of H2A.Z and trimethylation of H3K4 and H3K27 around the T-DNA integration sites. This study provides new insights into the structural impact of inserting foreign fragments into plant genomes and demonstrates the utility of state-of-the-art long-range sequencing technologies to rapidly identify unanticipated genomic changes. [ABSTRACT FROM AUTHOR]