Search

The impacts of invertebrate RNA virus population dynamics on virulence and infection outcomes are poorly understood. Deformed wing virus (DWV), the main viral pathogen of honey bees, negatively impacts bee health, which can lead to colony death. Despite previous reports on the reduction of DWV diversity following the arrival of the parasitic mite Varroa destructor, the key DWV vector, we found high genetic diversity of DWV in infested United States honey bee colonies. Phylogenetic analysis showed that divergent US DWV genotypes are of monophyletic origin and were likely generated as a result of diversification after a genetic bottleneck. To investigate the population dynamics of this divergent DWV, we designed a series of novel infectious cDNA clones corresponding to coexisting DWV genotypes, thereby devising a reverse-genetics system for an invertebrate RNA virus quasispecies. Equal replication rates were observed for all clone-derived DWV variants in single infections. Surprisingly, individual clones replicated to the same high levels as their mixtures and even the parental highly diverse natural DWV population, suggesting that complementation between genotypes was not required to replicate to high levels. Mixed clone–derived infections showed a lack of strong competitive exclusion, suggesting that the DWV genotypes were adapted to coexist. Mutational and recombination events were observed across clone progeny, providing new insights into the forces that drive and constrain virus diversification. Accordingly, our results suggest that Varroa influences DWV dynamics by causing an initial selective sweep, which is followed by virus diversification fueled by negative frequency-dependent selection for new genotypes. We suggest that this selection might reflect the ability of rare lineages to evade host defenses, specifically antiviral RNA interference (RNAi). In support of this hypothesis, we show that RNAi induced against one DWV strain is less effective against an alternate strain from the same population. [ABSTRACT FROM AUTHOR]

Cultivated watermelon (Citrullus lanatus) is one of the most important food crops in the Cucurbitaceae family. Diversification after domestication has led cultivated watermelons to exhibit diverse fruit flesh colors, including red, yellow, and orange. Recently, there has been increased interest in red-fleshed watermelons because they contain the antioxidant cis-isomeric lycopene. We performed whole genome resequencing (WGRS) of 24 watermelons with different flesh colors to identify single-nucleotide polymorphisms (SNPs) related to high lycopene content. The resequencing data revealed 203,894–279,412 SNPs from read mapping between inbred lines and the 97103 reference genome. In total, 295,065 filtered SNPs were identified, which had an average polymorphism information content of 0.297. Most of these SNPs were intergenic (90.1%) and possessed a transversion (Tv) rate of 31.64%. Overall, 2,369 SNPs were chosen at 0.5 Mb physical intervals to analyze genetic diversity across the 24 inbred lines. A neighbor-joining dendrogram and principal coordinate analysis (PCA) based on the 2,369 SNPs revealed that the 24 inbred lines could be grouped into high and low lycopene-type watermelons. In addition, we analyzed SNPs that could discriminate high lycopene content, red-fleshed watermelon from low lycopene, yellow or orange watermelon inbred lines. For validation, 19 SNPs (designated as WMHL1–19) were chosen randomly, and cleavage amplified polymorphic sequence (CAPS) markers were designed. Genotyping of the above 24 lines and 12 additional commercial cultivars using WMHL1–19 CAPS markers resulted in match rates of over 0.92 for most validated markers in correlation with the flesh color phenotypes. Our results provide valuable genomic information regarding the high lycopene content phenotype of red-fleshed cultivated watermelons, and the identified SNPs will be useful for the development of molecular markers in the marker-assisted breeding of watermelons with high lycopene content. [ABSTRACT FROM AUTHOR]

The hoverfly, Eupeodes corollae, is a worldwide natural enemy of aphids and a plant pollinator. To provide insights into the biology of this species, we examined its population genetic structure by obtaining 1.15-GB random genomic sequences using next-generation sequencing and developing genome-wide microsatellite markers. A total of 79,138 microsatellite loci were initially isolated from the genomic sequences; after strict selection and further testing of 40 primer pairs in eight individuals, 24 polymorphic microsatellites with high amplification rates were developed. These microsatellites were used to examine the population genetic structure of 96 individuals from four field populations collected across southern to northern China. The number of alleles per locus ranged from 5 to 13 with an average of 8.75; the observed and expected heterozygosity varied from 0.235 to 0.768 and from 0.333 to 0.785, respectively. Population genetic structure analysis showed weak genetic differentiation among the four geographical populations of E. corollae, suggesting a high rate of gene flow reflecting likely widespread migration of E. corollae in China. [ABSTRACT FROM AUTHOR]

The origins of new genes are among the most fundamental questions in evolutionary biology. Our understanding of the ways that new genetic material appears and how that genetic material shapes population variation remains incomplete. De novo genes and duplicate genes are a key source of new genetic material on which selection acts. To better understand the origins of these new gene sequences, we explored the ways that structural variation might alter expression patterns and form novel transcripts. We provide evidence that chromosomal rearrangements are a source of novel genetic variation that facilitates the formation of de novo exons in Drosophila. We identify 51 cases of de novo exon formation created by chromosomal rearrangements in 14 strains of D. yakuba. These new genes inherit transcription start signals and open reading frames when the 5’ end of existing genes are combined with previously untranscribed regions. Such new genes would appear with novel peptide sequences, without the necessity for secondary transitions from non-coding RNA to protein. This mechanism of new peptide formations contrasts with canonical theory of de novo gene progression requiring non-coding intermediaries that must acquire new mutations prior to loss via pseudogenization. Hence, these mutations offer a means to de novo gene creation and protein sequence formation in a single mutational step, answering a long standing open question concerning new gene formation. We further identify gene expression changes to 134 existing genes, indicating that these mutations can alter gene regulation. Population variability for chromosomal rearrangements is considerable, with 2368 rearrangements observed across 14 inbred lines. More rearrangements were identified on the X chromosome than any of the autosomes, suggesting the X is more susceptible to chromosome alterations. Together, these results suggest that chromosomal rearrangements are a source of variation in populations that is likely to be important to explain genetic and therefore phenotypic diversity. [ABSTRACT FROM AUTHOR]

Whole mitogenome sequences (mtDNA) have been exploited for insect ecology studies, using them as molecular markers to reconstruct phylogenies, or to infer phylogeographic relationships and gene flow. Recent Anopheles phylogenomic studies have provided information regarding the time of deep lineage divergences within the genus. Here we report the complete 15,393 bp mtDNA sequences of Anopheles aquasalis, a Neotropical human malaria vector. When comparing its structure and base composition with other relevant and available anopheline mitogenomes, high similarity and conserved genomic features were observed. Furthermore, 22 mtDNA sequences comprising anopheline and Dipteran sibling species were analyzed to reconstruct phylogenies and estimate dates of divergence between taxa. Phylogenetic analysis using complete mtDNA sequences suggests that A. aquasalis diverged from the Anopheles albitarsis complex ~28 million years ago (MYA), and ~38 MYA from Anopheles darlingi. Bayesian analysis suggests that the most recent ancestor of Nyssorhynchus and Anopheles + Cellia was extant ~83 MYA, corroborating current estimates of ~79–100 MYA. Additional sampling and publication of African, Asian, and North American anopheline mitogenomes would improve the resolution of the Anopheles phylogeny and clarify early continental dispersal routes. [ABSTRACT FROM AUTHOR]

Wolbachia are alpha-proteobacteria known to infect arthropods, which are of interest for disease control since they have been associated with improved resistance to viral infection. Although several genomes for different strains have been sequenced, there is little knowledge regarding the relationship between this bacterium and their hosts, particularly on their dependency for survival. Motivated by the potential applications on disease control, we developed genome-scale models of four Wolbachia strains known to infect arthropods: wAlbB (Aedes albopictus), wVitA (Nasonia vitripennis), wMel and wMelPop (Drosophila melanogaster). The obtained metabolic reconstructions exhibit a metabolism relying mainly on amino acids for energy production and biomass synthesis. A gap analysis was performed to detect metabolic candidates which could explain the endosymbiotic nature of this bacterium, finding that amino acids, requirements for ubiquinone precursors and provisioning of metabolites such as riboflavin could play a crucial role in this relationship. This work provides a systems biology perspective for studying the relationship of Wolbachia with its host and the development of new approaches for control of the spread of arboviral diseases. This approach, where metabolic gaps are key objects of study instead of just additions to complete a model, could be applied to other endosymbiotic bacteria of interest. [ABSTRACT FROM AUTHOR]

Conserved Noncoding Elements (CNEs) are elements exhibiting extreme noncoding conservation in Metazoan genomes. They cluster around developmental genes and act as long-range enhancers, yet nothing that we know about their function explains the observed conservation levels. Clusters of CNEs coincide with topologically associating domains (TADs), indicating ancient origins and stability of TAD locations. This has suggested further hypotheses about the still elusive origin of CNEs, and has provided a comparative genomics-based method of estimating the position of TADs around developmentally regulated genes in genomes where chromatin conformation capture data is missing. To enable researchers in gene regulation and chromatin biology to start deciphering this phenomenon, we developed CNEr, a R/Bioconductor toolkit for large-scale identification of CNEs and for studying their genomic properties. We apply CNEr to two novel genome comparisons—fruit fly vs tsetse fly, and two sea urchin genomes—and report novel insights gained from their analysis. We also show how to reveal interesting characteristics of CNEs by coupling CNEr with existing Bioconductor packages. CNEr is available at Bioconductor () and maintained at github (). [ABSTRACT FROM AUTHOR]

Patterned expression of many developmental genes is specified by transcription factor gene expression, but is thought to be refined by chromatin-mediated repression. Regulatory DNA sequences called Polycomb Response Elements (PREs) are required to repress some developmental target genes, and are widespread in genomes, suggesting that they broadly affect developmental programs. While PREs in transgenes can nucleate trimethylation on lysine 27 of the histone H3 tail (H3K27me3), none have been demonstrated to be necessary at endogenous chromatin domains. This failure is thought to be due to the fact that most endogenous H3K27me3 domains contain many PREs, and individual PREs may be redundant. In contrast to these ideas, we show here that PREs near the wing selector gene vestigial have distinctive roles at their endogenous locus, even though both PREs are repressors in transgenes. First, a PRE near the promoter is required for vestigial activation and not for repression. Second, only the distal PRE contributes to H3K27me3, but even removal of both PREs does not eliminate H3K27me3 across the vestigial domain. Thus, endogenous chromatin domains appear to be intrinsically marked by H3K27me3, and PREs appear required to enhance this chromatin modification to high levels at inactive genes. [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]

Essential tremor (ET) is one of the most common movement disorders. The etiology of ET remains largely unexplained. Whole genome sequencing (WGS) is likely to be of value in understanding a large proportion of ET with Mendelian and complex disease inheritance patterns. In ET families with Mendelian inheritance patterns, WGS may lead to gene identification where WES analysis failed to identify the causative single nucleotide variant (SNV) or indel due to incomplete coverage of the entire coding region of the genome, in addition to accurate detection of larger structural variants (SVs) and copy number variants (CNVs). Alternatively, in ET families with complex disease inheritance patterns with gene x gene and gene x environment interactions enrichment of functional rare coding and non-coding variants may explain the heritability of ET. We performed WGS in eight ET families (n = 40 individuals) enrolled in the Family Study of Essential Tremor. The analysis included filtering WGS data based on allele frequency in population databases, rare SNV and indel classification and association testing using the Mixed-Model Kernel Based Adaptive Cluster (MM-KBAC) test. A separate analysis of rare SV and CNVs segregating within ET families was also performed. Prioritization of candidate genes identified within families was performed using phenolyzer. WGS analysis identified candidate genes for ET in 5/8 (62.5%) of the families analyzed. WES analysis in a subset of these families in our previously published study failed to identify candidate genes. In one family, we identified a deleterious and damaging variant (c.1367G>A, p.(Arg456Gln)) in the candidate gene, CACNA1G, which encodes the pore forming subunit of T-type Ca(2+) channels, CaV3.1, and is expressed in various motor pathways and has been previously implicated in neuronal autorhythmicity and ET. Other candidate genes identified include SLIT3 which encodes an axon guidance molecule and in three families, phenolyzer prioritized genes that are associated with hereditary neuropathies (family A, KARS, family B, KIF5A and family F, NTRK1). Functional studies of CACNA1G and SLIT3 suggest a role for these genes in ET disease pathogenesis. [ABSTRACT FROM AUTHOR]

Background: Schistosoma japonicum is a parasitic flatworm that causes human schistosomiasis, which is a significant cause of morbidity in China and the Philippines. A single draft genome was available for S. japonicum, yet this assembly is very fragmented and only covers 90% of the genome, which make it difficult to be applied as a reference in functional genome analysis and genes discovery. Findings: In this study, we present a high-quality assembly of the fluke S. japonicum genome by combining 20 G (~53X) long single molecule real time sequencing reads with 80 G (~ 213X) Illumina paired-end reads. This improved genome assembly is approximately 370.5 Mb, with contig and scaffold N50 length of 871.9 kb and 1.09 Mb, representing 142.4-fold and 6.2-fold improvement over the released WGS-based assembly, respectively. Additionally, our assembly captured 85.2% complete and 4.6% partial eukaryotic Benchmarking Universal Single-Copy Orthologs. Repetitive elements account for 46.80% of the genome, and 10,089 of the protein-coding genes were predicted from the improved genome, of which 96.5% have been functionally annotated. Lastly, using the improved assembly, we identified 20 significantly expanded gene families in S. japonicum, and those genes were primarily enriched in functions of proteolysis and protein glycosylation. Conclusions: Using the combination of PacBio and Illumina Sequencing technologies, we provided an improved high-quality genome of S. japonicum. This improved genome assembly, as well as the annotation, will be useful for the comparative genomics of the flukes and more importantly facilitate the molecular studies of this important parasite in the future. [ABSTRACT FROM AUTHOR]

Ageing populations pose one of the main public health crises of our time. Reprogramming gene expression by altering the activities of sequence-specific transcription factors (TFs) can ameliorate deleterious effects of age. Here we explore how a circuit of TFs coordinates pro-longevity transcriptional outcomes, which reveals a multi-tissue and multi-species role for an entire protein family: the E-twenty-six (ETS) TFs. In Drosophila, reduced insulin/IGF signalling (IIS) extends lifespan by coordinating activation of Aop, an ETS transcriptional repressor, and Foxo, a Forkhead transcriptional activator. Aop and Foxo bind the same genomic loci, and we show that, individually, they effect similar transcriptional programmes in vivo. In combination, Aop can both moderate or synergise with Foxo, dependent on promoter context. Moreover, Foxo and Aop oppose the gene-regulatory activity of Pnt, an ETS transcriptional activator. Directly knocking down Pnt recapitulates aspects of the Aop/Foxo transcriptional programme and is sufficient to extend lifespan. The lifespan-limiting role of Pnt appears to be balanced by a requirement for metabolic regulation in young flies, in which the Aop-Pnt-Foxo circuit determines expression of metabolic genes, and Pnt regulates lipolysis and responses to nutrient stress. Molecular functions are often conserved amongst ETS TFs, prompting us to examine whether other Drosophila ETS-coding genes may also affect ageing. We show that five out of eight Drosophila ETS TFs play a role in fly ageing, acting from a range of organs and cells including the intestine, adipose and neurons. We expand the repertoire of lifespan-limiting ETS TFs in C. elegans, confirming their conserved function in ageing and revealing that the roles of ETS TFs in physiology and lifespan are conserved throughout the family, both within and between species. [ABSTRACT FROM AUTHOR]

Males of the Drosophila melanogaster mutant croaker (cro) generate a polycyclic pulse song dissimilar to the monocyclic songs typical of wild-type males during courtship. However, cro has not been molecularly mapped to any gene in the genome. We demonstrate that cro is a mutation in the gene encoding the Calmodulin-binding transcription factor (Camta) by genetic complementation tests with chromosomal deficiencies, molecular cloning of genomic fragments that flank the cro-mutagenic P-insertion, and phenotypic rescue of the cro mutant phenotype by Camta+-encoding cDNA as well as a BAC clone containing the gene for Camta. We further show that knockdown of the Camta-encoding gene phenocopies cro mutant songs when targeted to a subset of fruitless-positive neurons that include the mcALa and AL1 clusters in the brain. cro-GAL4 and an anti-Camta antibody labeled a large number of brain neurons including mcALa. We conclude that the Camta-encoding gene represents the cro locus, which has been implicated in a species-specific difference in courtship songs between D. sechellia and simulans. [ABSTRACT FROM AUTHOR]

Y chromosomes often contain amplified genes which can increase dosage of male fertility genes and counteract degeneration via gene conversion. Here we identify genes with increased copy number on both X and Y chromosomes in various species of Drosophila, a pattern that has previously been associated with sex chromosome drive involving the Slx and Sly gene families in mice. We show that recurrent X/Y co-amplification appears to be an important evolutionary force that has shaped gene content evolution of sex chromosomes in Drosophila. We also demonstrate that convergent acquisition and amplification of testis expressed gene families are common on Drosophila sex chromosomes, and especially on recently formed ones, and we carefully characterize one putative novel X/Y co-amplification system. We find that co-amplification of the S-Lap1/GAPsec gene pair on both the X and the Y chromosome occurred independently several times in members of the D. obscura group, where this normally autosomal gene pair is sex-linked due to a sex chromosome—autosome fusion. We explore several evolutionary scenarios that would explain this pattern of co-amplification. Investigation of gene expression and short RNA profiles at the S-Lap1/GAPsec system suggest that, like Slx/Sly in mice, these genes may be remnants of a cryptic sex chromosome drive system, however additional transgenic experiments will be necessary to validate this model. Regardless of whether sex chromosome drive is responsible for this co-amplification, our findings suggest that recurrent gene duplications between X and Y sex chromosomes could have a widespread effect on genomic and evolutionary patterns, including the epigenetic regulation of sex chromosomes, the distribution of sex-biased genes, and the evolution of hybrid sterility. [ABSTRACT FROM AUTHOR]

Background: With the expansion of soil transmitted helminth (STH) intervention efforts and the corresponding decline in infection prevalence, there is an increased need for sensitive and specific STH diagnostic assays. Previously, through next generation sequencing (NGS)-based identification and targeting of non-coding, high copy-number repetitive DNA sequences, we described the development of a panel of improved quantitative real-time PCR (qPCR)-based assays for the detection of Necator americanus, Ancylostoma duodenale, Ancylostoma ceylanicum, Trichuris trichiura, and Strongyloides stercoralis. However, due to the phenomenon of chromosome diminution, a similar assay based on high copy-number repetitive DNA was not developed for the detection of Ascaris lumbricoides. Recently, the publication of a reference-level germline genome sequence for A. lumbricoides has facilitated our development of an improved assay for this human pathogen of vast global importance. Methodology/Principal findings: Repurposing raw DNA sequence reads from a previously published Illumina-generated, NGS-based A. lumbricoides germline genome sequencing project, we performed a cluster-based repeat analysis utilizing RepeatExplorer2 software. This analysis identified the most prevalent repetitive DNA element of the A. lumbricoides germline genome (AGR, Ascaris germline repeat), which was then used to develop an improved qPCR assay. During experimental validation, this assay demonstrated a fold increase in sensitivity of ~3,100, as determined by relative Cq values, when compared with an assay utilizing a previously published, frequently employed, ribosomal ITS DNA target. A comparative analysis of 2,784 field-collected samples was then performed, successfully verifying this improved sensitivity. Conclusions/Significance: Through analysis of the germline genome sequence of A. lumbricoides, a vastly improved qPCR assay has been developed. This assay, utilizing a high copy number repeat target found in eggs and embryos (the AGR repeat), will improve prevalence estimates that are fundamental to the programmatic decision-making process, while simultaneously strengthening mathematical models used to examine STH infection rates. Furthermore, through the identification of an optimal target for PCR, future assay development efforts will also benefit, as the identity of the optimized repeat DNA target is likely to remain unchanged despite continued improvement in PCR-based diagnostic technologies. [ABSTRACT FROM AUTHOR]

The bertha armyworm (BAW), Mamestra configurata, is a significant pest of canola (Brassica napus L. and B. rapa L.) in western North America that undergoes cyclical outbreaks every 6–8 years. During peak outbreaks millions of dollars are spent on insecticidal control and, even with control efforts, subsequent damage can result in losses worth millions of dollars. Despite the importance of this pest insect, information is lacking on the dispersal ability of BAW and the genetic variation of populations from across its geographic range which may underlie potential differences in their susceptibility to insecticides or pathogens. Here, we examined the genetic diversity of BAW populations during an outbreak across its geographic range in western North America. First, mitochondrial cytochrome oxidase 1 (CO1) barcode sequences were used to confirm species identification of insects captured in a network of pheromone traps across the range, followed by haplotype analyses. We then sequenced the BAW genome and used double-digest restriction site associated DNA sequencing, mapped to the genome, to identify 1000s of single nucleotide polymorphisms (SNP) markers. CO1 haplotype analysis identified 9 haplotypes distributed across 28 sample locations and three laboratory-reared colonies. Analysis of genotypic data from both the CO1 and SNP markers revealed little population structure across BAW’s vast range. The CO1 haplotype pattern showed a star-like phylogeny which is often associated with species whose population abundance and range has recently expanded and combined with pheromone trap data, indicates the outbreak may have originated from a single focal point in central Saskatchewan. The relatively recent introduction of canola and rapid expansion of the canola growing region across western North America, combined with the cyclical outbreaks of BAW caused by precipitous population crashes, has likely selected for a genetically homogenous BAW population adapted to this crop. [ABSTRACT FROM AUTHOR]

The decline of many bumblebee species (Bombus spp.) has been linked to an increased prevalence of pathogens caused by spillover from managed bees. Although poorly understood, RNA viruses are suspected of moving from managed honeybees (Apis mellifera) into wild bumblebees through shared floral resources. We examined if RNA viruses spillover from managed honeybees, the extent to which viruses are replicating within bumblebees, and the role of flowers in transmission. Prevalence and active infections of deformed wing virus (DWV) were higher in bumblebees collected near apiaries and when neighboring honeybees had high infection levels. We found no DWV in bumblebees where honeybee foragers and honeybee apiaries were absent. The prevalence of black queen cell virus (BQCV) was also higher in bumblebees collected near apiaries. Furthermore, we detected viruses on 19% of flowers, all of which were collected within apiaries. Our results corroborate the hypothesis that viruses are spilling over from managed honeybees to wild bumblebees and that flowers may be an important route for transmission. [ABSTRACT FROM AUTHOR]

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei—the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])—within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp–exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes. [ABSTRACT FROM AUTHOR]

Stop-codon read-through refers to the phenomenon that a ribosome goes past the stop codon and continues translating into the otherwise untranslated region (UTR) of a transcript. Recent ribosome-profiling experiments in eukaryotes uncovered widespread stop-codon read-through that also varies among tissues, prompting the adaptive hypothesis that stop-codon read-through is an important, regulated mechanism for generating proteome diversity. Here we propose and test a competing hypothesis that stop-codon read-through arises mostly from molecular errors and is largely nonadaptive. The error hypothesis makes distinct predictions about the probability of read-through, frequency of sequence motifs for read-through, and conservation of the read-through region, each of which is supported by genome-scale data from yeasts and fruit flies. Thus, except for the few cases with demonstrated functions, stop-codon read-through is generally nonadaptive. This finding, along with other molecular errors recently quantified, reveals a much less precise or orderly cellular life than is commonly thought. [ABSTRACT FROM AUTHOR]

The Pol32 protein is one of the universal subunits of DNA polymerase δ (Pol δ), which is responsible for genome replication in eukaryotic cells. Although the role of Pol32 in DNA repair has been well-characterized, its exact function in genome replication remains obscure as studies in single cell systems have not established an essential role for Pol32 in the process. Here we characterize Pol32 in the context of Drosophila melanogaster development. In the rapidly dividing embryonic cells, loss of Pol32 halts genome replication as it specifically disrupts Pol δ localization to the nucleus. This function of Pol32 in facilitating the nuclear import of Pol δ would be similar to that of accessory subunits of DNA polymerases from mammalian Herpes viruses. In post-embryonic cells, loss of Pol32 reveals mitotic fragile sites in the Drosophila genome, a defect more consistent with Pol32’s role as a polymerase processivity factor. Interestingly, these fragile sites do not favor repetitive sequences in heterochromatin, with the rDNA locus being a striking exception. Our study uncovers a possibly universal function for DNA polymerase ancillary factors and establishes a powerful system for the study of chromosomal fragile sites in a non-mammalian organism. [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]

For many genes, proper gene expression requires coordinated and dynamic interactions between multiple regulatory elements, each of which can either promote or silence transcription. In Drosophila, the complexity of the regulatory landscape is further complicated by the tight physical pairing of homologous chromosomes, which can permit regulatory elements to interact in trans, a phenomenon known as transvection. To better understand how gene expression can be programmed through cis- and trans-regulatory interactions, we analyzed transvection effects for a collection of alleles of the eyes absent (eya) gene. We find that trans-activation of a promoter by the eya eye-specific enhancers is broadly supported in many allelic backgrounds, and that the availability of an enhancer to act in trans can be predicted based on the molecular lesion of an eya allele. Furthermore, by manipulating promoter availability in cis and in trans, we demonstrate that the eye-specific enhancers of eya show plasticity in their promoter preference between two different transcriptional start sites, which depends on promoter competition between the two potential targets. Finally, we show that certain alleles of eya demonstrate pairing-sensitive silencing resulting from trans-interactions between Polycomb Response Elements (PREs), and genetic and genomic data support a general role for PcG proteins in mediating transcriptional silencing at eya. Overall, our data highlight how eya gene regulation relies upon a complex but plastic interplay between multiple enhancers, promoters, and PREs. [ABSTRACT FROM AUTHOR]

Pythium guiyangense, an oomycete from a genus of mostly plant pathogens, is an effective biological control agent that has wide potential to manage diverse mosquitoes. However, its mosquito-killing mechanisms are almost unknown. In this study, we observed that P. guiyangense could utilize cuticle penetration and ingestion of mycelia into the digestive system to infect mosquito larvae. To explore pathogenic mechanisms, a high-quality genome sequence with 239 contigs and an N50 contig length of 1,009 kb was generated. The genome assembly is approximately 110 Mb, which is almost twice the size of other sequenced Pythium genomes. Further genome analysis suggests that P. guiyangense may arise from a hybridization of two related but distinct parental species. Phylogenetic analysis demonstrated that P. guiyangense likely evolved from common ancestors shared with plant pathogens. Comparative genome analysis coupled with transcriptome sequencing data suggested that P. guiyangense may employ multiple virulence mechanisms to infect mosquitoes, including secreted proteases and kazal-type protease inhibitors. It also shares intracellular Crinkler (CRN) effectors used by plant pathogenic oomycetes to facilitate the colonization of plant hosts. Our experimental evidence demonstrates that CRN effectors of P. guiyangense can be toxic to insect cells. The infection mechanisms and putative virulence effectors of P. guiyangense uncovered by this study provide the basis to develop improved mosquito control strategies. These data also provide useful knowledge on host adaptation and evolution of the entomopathogenic lifestyle within the oomycete lineage. A deeper understanding of the biology of P. guiyangense effectors might also be useful for management of other important agricultural pests. [ABSTRACT FROM AUTHOR]

The evolution of sexual dimorphism is constrained by a shared genome, leading to ‘sexual antagonism’, in which different alleles at given loci are favoured by selection in males and females. Despite its wide taxonomic incidence, we know little about the identity, genomic location, and evolutionary dynamics of antagonistic genetic variants. To address these deficits, we use sex-specific fitness data from 202 fully sequenced hemiclonal Drosophila melanogaster fly lines to perform a genome-wide association study (GWAS) of sexual antagonism. We identify approximately 230 chromosomal clusters of candidate antagonistic single nucleotide polymorphisms (SNPs). In contradiction to classic theory, we find no clear evidence that the X chromosome is a hot spot for sexually antagonistic variation. Characterising antagonistic SNPs functionally, we find a large excess of missense variants but little enrichment in terms of gene function. We also assess the evolutionary persistence of antagonistic variants by examining extant polymorphism in wild D. melanogaster populations and closely related species. Remarkably, antagonistic variants are associated with multiple signatures of balancing selection across the D. melanogaster distribution range and in their sister species D. simulans, indicating widespread and evolutionarily persistent (about 1 million years) genomic constraints on the evolution of sexual dimorphism. Based on our results, we propose that antagonistic variation accumulates because of constraints on the resolution of sexual conflict over protein coding sequences, thus contributing to the long-term maintenance of heritable fitness variation. [ABSTRACT FROM AUTHOR]

Bats are implicated as natural reservoirs for a wide range of zoonotic viruses including SARS and MERS coronaviruses, Ebola, Marburg, Nipah, Hendra, Rabies and other lyssaviruses. Accordingly, many One Health surveillance and viral discovery programs have focused on bats. In this report we present viral metagenomic data from bats collected in the Kingdom of Saudi Arabia [KSA]. Unbiased high throughput sequencing of fecal samples from 72 bat individuals comprising four species; lesser mouse-tailed bat (Rhinopoma hardwickii), Egyptian tomb bat (Taphozous perforatus), straw-colored fruit bat (Eidolon helvum), and Egyptian fruit bat (Rousettus aegyptiacus) revealed molecular evidence of a diverse set of viral families: Picornaviridae (hepatovirus, teschovirus, parechovirus), Reoviridae (rotavirus), Polyomaviridae (polyomavirus), Papillomaviridae (papillomavirus), Astroviridae (astrovirus), Caliciviridae (sapovirus), Coronaviridae (coronavirus), Adenoviridae (adenovirus), Paramyxoviridae (paramyxovirus), and unassigned mononegavirales (chuvirus). Additionally, we discovered a bastro-like virus (Middle East Hepe-Astrovirus), with a genomic organization similar to Hepeviridae. However, since it shared homology with Hepeviridae and Astroviridae at ORF1 and in ORF2, respectively, the newly discovered Hepe-Astrovirus may represent a phylogenetic bridge between Hepeviridae and Astroviridae. [ABSTRACT FROM AUTHOR]

The innate immune system in insects is regulated by specific signalling pathways. Most immune related pathways were identified and characterized in holometabolous insects such as Drosophila melanogaster, and it was assumed they would be highly conserved in all insects. The hemimetabolous insect, Rhodnius prolixus, has served as a model to study basic insect physiology, but also is a major vector of the human parasite, Trypanosoma cruzi, that causes 10,000 deaths annually. The publication of the R. prolixus genome revealed that one of the main immune pathways, the Immune-deficiency pathway (IMD), was incomplete and probably non-functional, an observation shared with other hemimetabolous insects including the pea aphid (Acyrthosiphon pisum) and the bedbug (Cimex lectularius). It was proposed that the IMD pathway is inactive in R. prolixus as an adaptation to prevent eliminating beneficial symbiont gut bacteria. We used bioinformatic analyses based on reciprocal BLAST and HMM-profile searches to find orthologs for most of the “missing” elements of the IMD pathway and provide data that these are regulated in response to infection with Gram-negative bacteria. We used RNAi strategies to demonstrate the role of the IMD pathway in regulating the expression of specific antimicrobial peptides (AMPs) in the fat body of R. prolixus. The data indicate that the IMD pathway is present and active in R. prolixus, which opens up new avenues of research on R. prolixus-T. cruzi interactions. [ABSTRACT FROM AUTHOR]

Germ cell immortality, or transgenerational maintenance of the germ line, could be promoted by mechanisms that could occur in either mitotic or meiotic germ cells. Here we report for the first time that the GSP-2 PP1/Glc7 phosphatase promotes germ cell immortality. Small RNA-induced genome silencing is known to promote germ cell immortality, and we identified a separation-of-function allele of C. elegans gsp-2 that is compromised for germ cell immortality and is also defective for small RNA-induced genome silencing and meiotic but not mitotic chromosome segregation. Previous work has shown that GSP-2 is recruited to meiotic chromosomes by LAB-1, which also promoted germ cell immortality. At the generation of sterility, gsp-2 and lab-1 mutant adults displayed germline degeneration, univalents, histone methylation and histone phosphorylation defects in oocytes, phenotypes that mirror those observed in sterile small RNA-mediated genome silencing mutants. Our data suggest that a meiosis-specific function of GSP-2 ties small RNA-mediated silencing of the epigenome to germ cell immortality. We also show that transgenerational epigenomic silencing at hemizygous genetic elements requires the GSP-2 phosphatase, suggesting a functional link to small RNAs. Given that LAB-1 localizes to the interface between homologous chromosomes during pachytene, we hypothesize that small localized discontinuities at this interface could promote genomic silencing in a manner that depends on small RNAs and the GSP-2 phosphatase. [ABSTRACT FROM AUTHOR]

The whitefly Bemisia tabaci MEAM1 is a notorious vector capable of transmitting many plant viruses, resulting in serious crop loss and food shortage around the world. To investigate potential sRNA-mediated regulatory mechanisms in whiteflies that are affected by virus acquisition and transmission, we conducted small RNA (sRNA) deep sequencing and performed genome-wide profiling of piwi-interacting RNAs (piRNAs) in whiteflies that were fed on tomato yellow leaf curl virus (TYLCV)-infected or non-infected tomato plants for 24, 48, and 72 h. In the present study, piRNA reads ranging from 564,395 to 1,715,652 per library were identified and shown to distribute unevenly in clusters (57 to 96 per library) on the whitefly (B. tabaci MEAM1) genome. Among them, 53 piRNA clusters were common for all treatments. Comparative analysis between libraries generated from viruliferous and non-viruliferous whiteflies identified five TYLCV-induced and 24 TYLCV-suppressed piRNA clusters. Approximately 62% of piRNAs were derived from non-coding sequences including intergenic regions, introns, and untranslated regions (UTRs). The remaining 38% were derived from coding sequences (CDS) or repeat elements. Interestingly, six protein coding genes were targeted by the TYLCV-induced piRNAs. We identified a large number of piRNAs that were distributed in clusters across the whitefly genome, with 60% being derived from non-coding regions. Comparative analysis revealed that feeding on a virus-infected host caused induction and suppression of only a small number of piRNA clusters in whiteflies. Although piRNAs primarily regulate the activity of transposable elements, our results suggest that they may have additional functions in regulating protein coding genes and in insect-virus interactions. [ABSTRACT FROM AUTHOR]

We evaluated the impact of bacterial rhabduscin synthesis on bacterial virulence and phenoloxidase inhibition in a Spodoptera model. We first showed that the rhabduscin cluster of the entomopathogenic bacterium Xenorhabdus nematophila was not necessary for virulence in the larvae of Spodoptera littoralis and Spodoptera frugiperda. Bacteria with mutations affecting the rhabduscin synthesis cluster (ΔisnAB and ΔGT mutants) were as virulent as the wild-type strain. We then developed an assay for measuring phenoloxidase activity in S. frugiperda and assessed the ability of bacterial culture supernatants to inhibit the insect phenoloxidase. Our findings confirm that the X. nematophila rhabduscin cluster is required for the inhibition of S. frugiperda phenoloxidase activity. The X. nematophila ΔisnAB mutant was unable to inhibit phenoloxidase, whereas ΔGT mutants displayed intermediate levels of phenoloxidase inhibition relative to the wild-type strain. The culture supernatants of Escherichia coli and of two entomopathogenic bacteria, Serratia entomophila and Xenorhabdus poinarii, were unable to inhibit S. frugiperda phenoloxidase activity. Heterologous expression of the X. nematophila rhabduscin cluster in these three strains was sufficient to restore inhibition. Interestingly, we observed pseudogenization of the X. poinarii rhabduscin gene cluster via the insertion of a 120 bp element into the isnA promoter. The inhibition of phenoloxidase activity by X. poinarii culture supernatants was restored by expression of the X. poinarii rhabduscin cluster under the control of an inducible Ptet promoter, consistent with recent pseudogenization. This study paves the way for advances in our understanding of the virulence of several entomopathogenic bacteria in non-model insects, such as the new invasive S. frugiperda species in Africa. [ABSTRACT FROM AUTHOR]

Graph-based representations are considered to be the future for reference genomes, as they allow integrated representation of the steadily increasing data on individual variation. Currently available tools allow de novo assembly of graph-based reference genomes, alignment of new read sets to the graph representation as well as certain analyses like variant calling and haplotyping. We here present a first method for calling ChIP-Seq peaks on read data aligned to a graph-based reference genome. The method is a graph generalization of the peak caller MACS2, and is implemented in an open source tool, Graph Peak Caller. By using the existing tool vg to build a pan-genome of Arabidopsis thaliana, we validate our approach by showing that Graph Peak Caller with a pan-genome reference graph can trace variants within peaks that are not part of the linear reference genome, and find peaks that in general are more motif-enriched than those found by MACS2. [ABSTRACT FROM AUTHOR]

Endogenous small RNAs (sRNAs) and Argonaute proteins are ubiquitous regulators of gene expression in germline and somatic tissues. sRNA-Argonaute complexes are often expressed in gametes and are consequently inherited by the next generation upon fertilization. In Caenorhabditis elegans, 26G-RNAs are primary endogenous sRNAs that trigger the expression of downstream secondary sRNAs. Two subpopulations of 26G-RNAs exist, each of which displaying strongly compartmentalized expression: one is expressed in the spermatogenic gonad and associates with the Argonautes ALG-3/4; plus another expressed in oocytes and in embryos, which associates with the Argonaute ERGO-1. The determinants and dynamics of gene silencing elicited by 26G-RNAs are largely unknown. Here, we provide diverse new insights into these endogenous sRNA pathways of C. elegans. Using genetics and deep sequencing, we dissect a maternal effect of the ERGO-1 branch of the 26G-RNA pathway. We find that maternal primary sRNAs can trigger the production of zygotic secondary sRNAs that are able to silence targets, even in the absence of zygotic primary triggers. Thus, the interaction of maternal and zygotic sRNA populations, assures target gene silencing throughout animal development. Furthermore, we explore other facets of 26G-RNA biology related to the ALG-3/4 branch. We find that sRNA abundance, sRNA pattern of origin and the 3’ UTR length of target transcripts are predictors of the regulatory outcome by the Argonautes ALG-3/4. Lastly, we provide evidence suggesting that ALG-3 and ALG-4 regulate their own mRNAs in a negative feedback loop. Altogether, we provide several new regulatory insights on the dynamics, target regulation and self-regulation of the endogenous RNAi pathways of C. elegans. [ABSTRACT FROM AUTHOR]

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well. [ABSTRACT FROM AUTHOR]

Within family Baculoviridae, members of the Betabaculovirus genus are employed as biocontrol agents against lepidopteran pests, either alone or in combination with selected members of the Alphabaculovirus genus. Epinotia aporema granulovirus (EpapGV) is a fast killing betabaculovirus that infects the bean shoot borer (E. aporema) and is a promising biopesticide. Because occlusion bodies (OBs) play a key role in baculovirus horizontal transmission, we investigated the composition of EpapGV OBs. Using mass spectrometry-based proteomics we could identify 56 proteins that are included in the OBs during the final stages of larval infection. Our data provides experimental validation of several annotated hypothetical coding sequences. Proteogenomic mapping against genomic sequence detected a previously unannotated ac110-like core gene and a putative translation fusion product of ORFs epap48 and epap49. Comparative studies of the proteomes available for the family Baculoviridae highlight the conservation of core gene products as parts of the occluded virion. Two proteins specific for betabaculoviruses (Epap48 and Epap95) are incorporated into OBs. Moreover, quantification based on emPAI values showed that Epap95 is one of the most abundant components of EpapGV OBs. [ABSTRACT FROM AUTHOR]

The repeatability or predictability of evolution is a central question in evolutionary biology and most often addressed in experimental evolution studies. Here, we infer how genetically heterogeneous natural systems acquire the same molecular changes to address how genomic background affects adaptation in natural populations. In particular, we take advantage of independently formed neo-sex chromosomes in Drosophila species that have evolved dosage compensation by co-opting the dosage-compensation male-specific lethal (MSL) complex to study the mutational paths that have led to the acquisition of hundreds of novel binding sites for the MSL complex in different species. This complex recognizes a conserved 21-bp GA-rich sequence motif that is enriched on the X chromosome, and newly formed X chromosomes recruit the MSL complex by de novo acquisition of this binding motif. We identify recently formed sex chromosomes in the D. melanica and D. robusta species groups by genome sequencing and generate genomic occupancy maps of the MSL complex to infer the location of novel binding sites. We find that diverse mutational paths were utilized in each species to evolve hundreds of de novo binding motifs along the neo-X, including expansions of microsatellites and transposable element (TE) insertions. However, the propensity to utilize a particular mutational path differs between independently formed X chromosomes and appears to be contingent on genomic properties of that species, such as simple repeat or TE density. This establishes the “genomic environment” as an important determinant in predicting the outcome of evolutionary adaptations. [ABSTRACT FROM AUTHOR]

Genome annotation is the process of identifying the location and function of a genome's encoded features. Improving the biological accuracy of annotation is a complex and iterative process requiring researchers to review and incorporate multiple sources of information such as transcriptome alignments, predictive models based on sequence profiles, and comparisons to features found in related organisms. Because rapidly decreasing costs are enabling an ever-growing number of scientists to incorporate sequencing as a routine laboratory technique, there is widespread demand for tools that can assist in the deliberative analytical review of genomic information. To this end, we present Apollo, an open source software package that enables researchers to efficiently inspect and refine the precise structure and role of genomic features in a graphical browser-based platform. Some of Apollo’s newer user interface features include support for real-time collaboration, allowing distributed users to simultaneously edit the same encoded features while also instantly seeing the updates made by other researchers on the same region in a manner similar to Google Docs. Its technical architecture enables Apollo to be integrated into multiple existing genomic analysis pipelines and heterogeneous laboratory workflow platforms. Finally, we consider the implications that Apollo and related applications may have on how the results of genome research are published and made accessible. [ABSTRACT FROM AUTHOR]

Palmitoylation is an important posttranslational modification regulating diverse cellular functions. Consequently, aberrant palmitoylation can lead to diseases such as neuronal disorders or cancer. In humans there are roughly one hundred times more palmitoylated proteins than enzymes catalyzing palmitoylation (palmitoyltransferases). Therefore, it is an important challenge to establish the links between palmitoyltransferases and their targets. From publicly available data, we find that expression of human ZDHHC8 correlates significantly with cancer survival. To elucidate the organismal function of ZDHHC8, we study the Drosophila ortholog of hZDHHC8, CG34449/dZDHHC8. Knockdown of dZDHHC8 causes tissue overgrowth while dZDHHC8 mutants are larval lethal. We provide a list of 159 palmitoylated proteins in Drosophila and present data suggesting that scribble and Ras64B are targets of dZDHHC8. [ABSTRACT FROM AUTHOR]

Spodoptera exempta nucleopolyhedrovirus (SpexNPV) is a viral pathogen of the African armyworm, Spodoptera exempta (Lepidoptera: Noctuidae), a significant agricultural pest of cereal crops in Africa. SpexNPV has been evaluated as a potential insecticide for control of this pest and has served as the subject of research on baculovirus pathology and transmission. Occlusion bodies (OBs) of SpexNPV isolate 244.1 were examined, and the nucleotide sequence of the genome was determined and characterized. SpexNPV-244.1 OBs consisted of irregular polyhedra with a size and appearance typical for alphabaculoviruses. Virions within the polyhedra contained 1–8 nucleocapsids per unit envelope. The SpexNPV-244.1 genome was comprised of a 129,528 bp circular sequence, in which 139 ORFs were annotated. Five homologous regions (hrs) consisting of a variable number of 28-bp imperfect palindromes were identified in the genome. The genome sequence contained the 38 core genes of family Baculoviridae, as well as three ORFs unique to the SpexNPV sequence and one ORF that was apparently acquired by horizontal gene transfer with a betabaculovirus ancestor. Phylogenetic inference with core gene amino acid sequence alignments placed SpexNPV-244.1 in a lineage containing alphabaculoviruses of Spodoptera frugiperda and Spodopotera exigua which in turn is part of a larger group of alphabaculoviruses from the subfamily Noctuinae in the lepidopteran family Noctuidae. Kimura-2-parameter pairwise nucleotide distances indicated that SpexNPV-244.1 represented a different and previously unlisted species in the genus Alphabaculovirus. Gene parity plots indicated that the gene order of SpexNPV-244.l was extensively collinear with that of Spodoptera exigua NPV (SeMNPV). These plots also revealed a group of 17 core genes whose order was conserved in other alpha- and betabaculoviruses. [ABSTRACT FROM AUTHOR]

The impact of pesticides on the health of bee pollinators is determined in part by the capacity of bee detoxification systems to convert these compounds to less toxic forms. For example, recent work has shown that cytochrome P450s of the CYP9Q subfamily are critically important in defining the sensitivity of honey bees and bumblebees to pesticides, including neonicotinoid insecticides. However, it is currently unclear if solitary bees have functional equivalents of these enzymes with potentially serious implications in relation to their capacity to metabolise certain insecticides. To address this question, we sequenced the genome of the red mason bee, Osmia bicornis, the most abundant and economically important solitary bee species in Central Europe. We show that O. bicornis lacks the CYP9Q subfamily of P450s but, despite this, exhibits low acute toxicity to the N-cyanoamidine neonicotinoid thiacloprid. Functional studies revealed that variation in the sensitivity of O. bicornis to N-cyanoamidine and N-nitroguanidine neonicotinoids does not reside in differences in their affinity for the nicotinic acetylcholine receptor or speed of cuticular penetration. Rather, a P450 within the CYP9BU subfamily, with recent shared ancestry to the Apidae CYP9Q subfamily, metabolises thiacloprid in vitro and confers tolerance in vivo. Our data reveal conserved detoxification pathways in model solitary and eusocial bees despite key differences in the evolution of specific pesticide-metabolising enzymes in the two species groups. The discovery that P450 enzymes of solitary bees can act as metabolic defence systems against certain pesticides can be leveraged to avoid negative pesticide impacts on these important pollinators. [ABSTRACT FROM AUTHOR]

The genetic architecture of adaptive traits is of key importance to predict evolutionary responses. Most adaptive traits are polygenic—i.e., result from selection on a large number of genetic loci—but most molecularly characterized traits have a simple genetic basis. This discrepancy is best explained by the difficulty in detecting small allele frequency changes (AFCs) across many contributing loci. To resolve this, we use laboratory natural selection to detect signatures for selective sweeps and polygenic adaptation. We exposed 10 replicates of a Drosophila simulans population to a new temperature regime and uncovered a polygenic architecture of an adaptive trait with high genetic redundancy among beneficial alleles. We observed convergent responses for several phenotypes—e.g., fitness, metabolic rate, and fat content—and a strong polygenic response (99 selected alleles; mean s = 0.059). However, each of these selected alleles increased in frequency only in a subset of the evolving replicates. We discerned different evolutionary paradigms based on the heterogeneous genomic patterns among replicates. Redundancy and quantitative trait (QT) paradigms fitted the experimental data better than simulations assuming independent selective sweeps. Our results show that natural D. simulans populations harbor a vast reservoir of adaptive variation facilitating rapid evolutionary responses using multiple alternative genetic pathways converging at a new phenotypic optimum. This key property of beneficial alleles requires the modification of testing strategies in natural populations beyond the search for convergence on the molecular level. [ABSTRACT FROM AUTHOR]

More than any other genome components, Transposable Elements (TEs) have the capacity to move across species barriers through Horizontal Transfer (HT), with substantial evolutionary consequences. Previous large-scale surveys, based on full-genomes comparisons, have revealed the transposition mode as an important predictor of HT rates variation across TE superfamilies. However, host biology could represent another major explanatory factor, one that needs to be investigated through extensive taxonomic sampling. Here we test this hypothesis using a field collection of 460 Arthropod species from Tahiti and surrounding islands. Through targeted massive parallel sequencing, we uncover patterns of HT in three widely-distributed TE superfamilies with contrasted modes of transposition. In line with earlier findings, the DNA transposons under study (TC1-Mariner) were found to transfer horizontally at the highest frequency, closely followed by the LTR superfamily (Copia), in contrast with the non-LTR superfamily (Jockey), that mostly diversifies through vertical inheritance and persists longer within genomes. Strikingly, across all superfamilies, we observe a marked excess of HTs in Lepidoptera, an insect order that also commonly hosts baculoviruses, known for their ability to transport host TEs. These results turn the spotlight on baculoviruses as major potential vectors of TEs in Arthropods, and further emphasize the importance of non-vertical TE inheritance in genome evolution. [ABSTRACT FROM AUTHOR]

Intestinal infection by the parasitic blood fluke Schistosoma bovis is a common veterinary problem in Africa and the Middle East and occasionally in the Mediterranean Region. The species also has the ability to form interspecific hybrids with the human parasite S. haematobium with natural hybridisation observed in West Africa, presenting possible zoonotic transmission. Additionally, this exchange of alleles between species may dramatically influence disease dynamics and parasite evolution. We have generated a 374 Mb assembly of the S. bovis genome using Illumina and PacBio-based technologies. Despite infecting different hosts and organs, the genome sequences of S. bovis and S. haematobium appeared strikingly similar with 97% sequence identity. The two species share 98% of protein-coding genes, with an average sequence identity of 97.3% at the amino acid level. Genome comparison identified large continuous parts of the genome (up to several 100 kb) showing almost 100% sequence identity between S. bovis and S. haematobium. It is unlikely that this is a result of genome conservation and provides further evidence of natural interspecific hybridization between S. bovis and S. haematobium. Our results suggest that foreign DNA obtained by interspecific hybridization was maintained in the population through multiple meiosis cycles and that hybrids were sexually reproductive, producing viable offspring. The S. bovis genome assembly forms a highly valuable resource for studying schistosome evolution and exploring genetic regions that are associated with species-specific phenotypic traits. [ABSTRACT FROM AUTHOR]

Like other domesticates, the efficient utilization of nitrogen resources is also important for the only fully domesticated insect, the silkworm. Deciphering the way in which artificial selection acts on the silkworm genome to improve the utilization of nitrogen resources and to advance human-favored domestication traits, will provide clues from a unique insect model for understanding the general rules of Darwin's evolutionary theory on domestication. Storage proteins (SPs), which belong to a hemocyanin superfamily, basically serve as a source of amino acids and nitrogen during metamorphosis and reproduction in insects. In this study, through blast searching on the silkworm genome and further screening of the artificial selection signature on silkworm SPs, we discovered a candidate domestication gene, i.e., the methionine-rich storage protein 1 (SP1), which is clearly divergent from other storage proteins and exhibits increased expression in the ova of domestic silkworms. Knockout of SP1 via the CRISPR/Cas9 technique resulted in a dramatic decrease in egg hatchability, without obvious impact on egg production, which was similar to the effect in the wild silkworm compared with the domestic type. Larval development and metamorphosis were not affected by SP1 knockout. Comprehensive ova comparative transcriptomes indicated significant higher expression of genes encoding vitellogenin, chorions, and structural components in the extracellular matrix (ECM)-interaction pathway, enzymes in folate biosynthesis, and notably hormone synthesis in the domestic silkworm, compared to both the SP1 mutant and the wild silkworm. Moreover, compared with the wild silkworms, the domestic one also showed generally up-regulated expression of genes enriched in the structural constituent of ribosome and amide, as well as peptide biosynthesis. This study exemplified a novel case in which artificial selection could act directly on nitrogen resource proteins, further affecting egg nutrients and eggshell formation possibly through a hormone signaling mediated regulatory network and the activation of ribosomes, resulting in improved biosynthesis and increased hatchability during domestication. These findings shed new light on both the understanding of artificial selection and silkworm breeding from the perspective of nitrogen and amino acid resources. [ABSTRACT FROM AUTHOR]

Bats are newly identified reservoirs of hantaviruses (HVs) among which very divergent HVs have been discovered in recent years. However, their significance for public health remains unclear since their seroprevalence as well as antigenic relationship with human-infecting HVs have not been investigated. In the present study archived tissues of 1,419 bats of 22 species from 6 families collected in 5 south and southwest provinces in China were screened by pan-HV RT-PCR following viral metagenomic analysis. As a result nine HVs have been identified in two bat species in two provinces and phylogenetically classified into two species, Laibin virus (LAIV, ICTV approved species, 1 strain) and Xuan son virus (XSV, proposed species, 8 strains). Additionally, 709 serum samples of these bats were also analyzed by ELISA to investigate the seroprevalence and cross-reactivity between different HVs using expressed recombinant nucleocapsid proteins (rNPs) of LAIV, XSV and Seoul virus (SEOV). The cross-reactivity of some bat sera were further confirmed by western blot (WB) using three rNPs followed by fluorescent antibody virus neutralization test (FAVNT) against live SEOV. Results showed that the total HV seropositive rate of bat sera was 18.5% (131/709) with many cross reacting with two or all three rNPs and several able to neutralize SEOV. WB analysis using the three rNPs and their specific hyperimmune sera demonstrated cross-reactivity between XSV/SEOV and LAIV/XSV, but not LAIV/SEOV, indicating that XSV is antigenically closer to human-infecting HVs. In addition a study of the distribution of the viruses identified an area covering the region between Chinese Guangxi and North Vietnam, in which XSV and LAIV circulate within different bat colonies with a high seroprevalence. A circulation sphere of bat-borne HVs has therefore been proposed. [ABSTRACT FROM AUTHOR]

Since the initial description of the genomic patterns expected under models of positive selection acting on standing genetic variation and on multiple beneficial mutations—so-called soft selective sweeps—researchers have sought to identify these patterns in natural population data. Indeed, over the past two years, large-scale data analyses have argued that soft sweeps are pervasive across organisms of very different effective population size and mutation rate—humans, Drosophila, and HIV. Yet, others have evaluated the relevance of these models to natural populations, as well as the identifiability of the models relative to other known population-level processes, arguing that soft sweeps are likely to be rare. Here, we look to reconcile these opposing results by carefully evaluating three recent studies and their underlying methodologies. Using population genetic theory, as well as extensive simulation, we find that all three examples are prone to extremely high false-positive rates, incorrectly identifying soft sweeps under both hard sweep and neutral models. Furthermore, we demonstrate that well-fit demographic histories combined with rare hard sweeps serve as the more parsimonious explanation. These findings represent a necessary response to the growing tendency of invoking parameter-heavy, assumption-laden models of pervasive positive selection, and neglecting best practices regarding the construction of proper demographic null models. [ABSTRACT FROM AUTHOR]

Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method–integrative modeling of genomic regions (IMGR)–to increase the genomic resolution of our traces to 10 kb. [ABSTRACT FROM AUTHOR]

The earliest stage of animal development is controlled by maternally deposited mRNA transcripts and proteins. Once the zygote is able to transcribe its own genome, maternal transcripts are degraded, in a tightly regulated process known as the maternal to zygotic transition (MZT). While this process has been well-studied within model species, we have little knowledge of how the pools of maternal and zygotic transcripts evolve. To characterize the evolutionary dynamics and functional constraints on early embryonic expression, we created a transcriptomic dataset for 14 Drosophila species spanning over 50 million years of evolution, at developmental stages before and after the MZT, and compared our results with a previously published Aedes aegypti developmental time course. We found deep conservation over 250 million years of a core set of genes transcribed only by the zygote. This select group is highly enriched in transcription factors that play critical roles in early development. However, we also identify a surprisingly high level of change in the transcripts represented at both stages over the phylogeny. While mRNA levels of genes with maternally deposited transcripts are more highly conserved than zygotic genes, those maternal transcripts that are completely degraded at the MZT vary dramatically between species. We also show that hundreds of genes have different isoform usage between the maternal and zygotic genomes. Our work suggests that maternal transcript deposition and early zygotic transcription are remarkably dynamic over evolutionary time, despite the widespread conservation of early developmental processes. [ABSTRACT FROM AUTHOR]

Dragon’s blood collected from the genus Dracaena is used as a renowned traditional medicine in various cultures worldwide. However, the genetics of the genus Dracaena and the formation mechanism of dragon’s blood remain poorly understood. Here, we generate the first draft genome reference assembly of an elite Chinese Dracaena species, Dracaena cambodiana, from next-generation sequencing data with 89.46× coverage. The reads were assembled into 2,640,704 contigs with an N50 length of 1.87 kb, and a 1.05 Gb assembly was finally assembled with 2,379,659 scaffolds. Furthermore, 97.75% of the 267,243 simple sequence repeats identified from these scaffolds were mononucleotide, dinucleotide, and trinucleotide repeats. Among all 53,700 predicted genes, 158 genes involved in cell wall and plant hormone synthesis and reactive oxygen species scavenging showed altered regulation during the formation of dragon’s blood. This study provides a genomic characterization of D. cambodiana and improves understanding of the molecular mechanism of dragon’s blood formation. This report represents the first genome-wide characterization of a Dracaena species in the Asparagaceae. [ABSTRACT FROM AUTHOR]

In Drosophila melanogaster, the male-specific lethal (MSL) complex plays a key role in dosage compensation by stimulating expression of male X-chromosome genes. It consists of MSL proteins and two long noncoding RNAs, roX1 and roX2, that are required for spreading of the complex on the chromosome and are redundant in the sense that loss of either does not affect male viability. However, despite rapid evolution, both roX species are present in diverse Drosophilidae species, raising doubts about their full functional redundancy. Thus, we have investigated consequences of deleting roX1 and/or roX2 to probe their specific roles and redundancies in D. melanogaster. We have created a new mutant allele of roX2 and show that roX1 and roX2 have partly separable functions in dosage compensation. In larvae, roX1 is the most abundant variant and the only variant present in the MSL complex when the complex is transmitted (physically associated with the X-chromosome) in mitosis. Loss of roX1 results in reduced expression of the genes on the X-chromosome, while loss of roX2 leads to MSL-independent upregulation of genes with male-biased testis-specific transcription. In roX1 roX2 mutant, gene expression is strongly reduced in a manner that is not related to proximity to high-affinity sites. Our results suggest that high tolerance of mis-expression of the X-chromosome has evolved. We propose that this may be a common property of sex-chromosomes, that dosage compensation is a stochastic process and its precision for each individual gene is regulated by the density of high-affinity sites in the locus. [ABSTRACT FROM AUTHOR]

The fall armyworm (FAW), Spodoptera frugiperda (J.E. Smith), is an important pest of maize in the Americas and has recently been introduced into Africa. Fall armyworm populations have developed resistance to control strategies that depend on insecticides and transgenic plants expressing Bacillus thuringiensis toxins. The study of various resistance mechanisms at the molecular level and the development novel control strategies have been hampered by a lack of functional genomic tools such as gene editing in this pest. In the current study, we explored the possibility of using the CRISPR/Cas9 system to modify the genome of FAW. We first identified and characterized the abdominal-A (Sfabd-A) gene of FAW. Sfabd-A single guide RNA (sgRNA) and Cas9 protein were then injected into 244 embryos of FAW. Sixty-two embryos injected with Sfabd-A sgRNA hatched. Of these hatched embryos, twelve developed into larvae that displayed typical aba-A mutant phenotypes such as fused segments. Of the twelve mutant larvae, three and five eventually developed into female and male moths, respectively. Most mutant moths were sterile, and one female produced a few unviable eggs when it was outcrossed to a wild-type male. Genotyping of 20 unhatched Sfabd-A sgRNA-injected embryos and 42 moths that developed from Sfabd-A sgRNA-injected embryos showed that 100% of the unhatched embryos and 50% of the moths contained indel mutations at the Sfabd-A genomic locus near the guide RNA target site. These results suggest that the CRISPR/Cas9 system is highly efficient in editing FAW genome. Importantly, this gene editing technology can be used to validate gene function to facilitate an understanding of the resistance mechanism and lead to the development of novel pest management approaches. [ABSTRACT FROM AUTHOR]

The parasitic mite, Varroa destructor, in combination with the viruses it vectors, is the main cause for global colony losses of the European honeybee, Apis mellifera. However, an isolated honeybee population established in 1999 on the Island of Gotland, Sweden has naturally acquired resistance to the mite, and has survived without mite control treatment for more than 18 years. A recent study has shown that this mite resistant (MR) population also appears to be resistant to Black queen cell virus (BQCV) and Sacbrood virus (SBV) and tolerant to Deformed wing virus (DWV), relative to nearby mite susceptible (MS) honeybee populations. In this study, RNA sequencing was employed to corroborate these previous findings and identify other viral factors that may play a role in the enhanced survival of this mite resistant honeybee population. Two additional honeybee-infecting viruses, Apis rhabdovirus-1 (ARV-1) and Lake Sinai virus (LSV), were identified and near-complete genomes of these two viruses were obtained. Phylogenetic analyses of the assembled virus sequences revealed consistent separation between the MR and MS honeybee populations, although it is unclear whether this is due to pre-existing differences between the viruses in the two populations when they were established, and isolated, or due to virus genetic adaptation towards reduced virulence in the MR population, to promote colony survival. Reverse transcription quantitative polymerase chain reaction(RT-qPCR) analyses show higher ARV and LSV titres in MS colonies compared to MR colonies, gradually increasing from summer to autumn 2009, and reaching maximum titres in the following spring 2010. While the DWV and BQCV titres in MR colonies increased between autumn 2009 and spring 2010, the SBV practically disappeared entirely by spring 2010. Possible explanations for the apparent virus tolerance or resistance in the Gotland mite-resistant honeybee population are discussed. [ABSTRACT FROM AUTHOR]