Whole genome phylogeny for 21 Drosophila species using predicted 2b-RAD fragments

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BACKGROUND

Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.

RESULTS

Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.

CONCLUSIONS

Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

Genetic variation in disease resistance in Drosophila spp. is mitigated in Drosophila sechellia by specialization to a toxic host

We found that Drosophila species vary in their susceptibility to the broad-spectrum entomopathogen, Metarhizium anisopliae (strain Ma549). Generalist species were generally more resistant than dietary specialists, with the cactophilic Drosophila buzzatii and Drosophila sechellia, a specialist of the Morinda citrifolia (Morinda) fruit, being most susceptible. Morinda fruit is reported to be toxic to most herbivores because it contains Octanoic Acid (OA). We confirmed that OA is toxic to Drosophila spp., other than D. sechellia, and we also found that OA is highly toxic to entomopathogenic fungi including Ma549 and Beauveria bassiana. Drosophila sechellia fed a diet containing OA, even at levels much less than found in Morinda fruit, had greatly reduced susceptibility to Ma549. This suggests that specializing to Morinda may have provided an enemy-free space, reducing adaptive prioritization on a strong immune response. Our results demonstrate that M. anisopliae and Drosophila species with divergent lifestyles provide a versatile model system for understanding the mechanisms of host-pathogen interactions at different scales and in environmental context.

Comparison of Preference for Chemicals Associated with Fruit Fermentation between Drosophila melanogaster and Drosophila suzukii and between Virgin and Mated D. melanogaster

Two taxonomically similar Drosophila species, Drosophila melanogaster and Drosophila suzukii, are known to have distinct habitats: D. melanogaster is mostly found near overripe and fermented fruits, whereas D. suzukii is attracted to fresh fruits. Since chemical concentrations are typically higher in overripe and fermented fruits than in fresh fruits, D. melanogaster is hypothesized to be attracted to higher concentrations of volatiles than D. suzukii. Therefore, the chemical preferences of the two flies were compared via Y-tube olfactometer assays and electroantennogram (EAG) experiments using various concentrations of 2-phenylethanol, ethanol, and acetic acid. D. melanogaster exhibited a higher preference for high concentrations of all the chemicals than that of D. suzukii. In particular, since acetic acid is mostly produced at the late stage of fruit fermentation, the EAG signal distance to acetic acid between the two flies was higher than those to 2-phenylethanol and ethanol. This supports the hypothesis that D. melanogaster prefers fermented fruits compared to D. suzukii. When comparing virgin and mated female D. melanogaster, mated females showed a higher preference for high concentrations of chemicals than that of virgin females. In conclusion, high concentrations of volatiles are important attraction factors for mated females seeking appropriate sites for oviposition.

Telomeric retrotransposons show propensity to form G-quadruplexes in various eukaryotic species

BACKGROUND

Canonical telomeres (telomerase-synthetised) are readily forming G-quadruplexes (G4) on the G-rich strand. However, there are examples of non-canonical telomeres among eukaryotes where telomeric tandem repeats are invaded by specific retrotransposons. Drosophila melanogaster represents an extreme example with telomeres composed solely by three retrotransposons-Het-A, TAHRE and TART (HTT). Even though non-canonical telomeres often show strand biased G-distribution, the evidence for the G4-forming potential is limited.

RESULTS

Using circular dichroism spectroscopy and UV absorption melting assay we have verified in vitro G4-formation in the HTT elements of D. melanogaster. Namely 3 in Het-A, 8 in TART and 2 in TAHRE. All the G4s are asymmetrically distributed as in canonical telomeres. Bioinformatic analysis showed that asymmetric distribution of potential quadruplex sequences (PQS) is common in telomeric retrotransposons in other Drosophila species. Most of the PQS are located in the gag gene where PQS density correlates with higher DNA sequence conservation and codon selection favoring G4-forming potential. The importance of G4s in non-canonical telomeres is further supported by analysis of telomere-associated retrotransposons from various eukaryotic species including green algae, Diplomonadida, fungi, insects and vertebrates. Virtually all analyzed telomere-associated retrotransposons contained PQS, frequently with asymmetric strand distribution. Comparison with non-telomeric elements showed independent selection of PQS-rich elements from four distinct LINE clades.

CONCLUSION

Our findings of strand-biased G4-forming motifs in telomere-associated retrotransposons from various eukaryotic species support the G4-formation as one of the prerequisites for the recruitment of specific retrotransposons to chromosome ends and call for further experimental studies.

2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics

Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.

Expression of antimicrobial peptides associated with different susceptibilities to environmental chemicals in Drosophila suzukii and Drosophila melanogaster

Drosophila suzukii is a serious agricultural pest. The evolved morphology of the female D. suzukii assists in penetrating the surface of fresh fruit and spawns eggs with its unique ovipositor. Conversely, Drosophila melanogaster, a taxonomically close species with D. suzukii, largely inhabits decaying and fermenting fruits and is consistently exposed to extensive environmental chemicals, such as 2-phenylethanol, ethanol, and acetic acid, produced by microorganisms. Considering the distinct habitats of the two flies, D. suzukii is thought to be more susceptible to environmental chemicals than D. melanogaster. We investigated the significantly higher survival rate of D. melanogaster following exposure to 2-phenylethanol, ethanol, and acetic acid. A comparison of the expression of antimicrobial peptides (AMPs) between the two flies treated with chemicals established that AMPs were generally more abundantly induced in D. melanogaster than in D. suzukii, particularly in the gut and fat body. Among the AMPs, the induction of genes (Diptericin A, Diptericin B, and Metchnikowin), which are regulated by the immune deficiency (IMD) pathway, was significantly higher than that of Drosomycin, which belongs to the Toll pathway in chemical-treated D. melanogaster. A transgenic RNAi fly (D. melanogaster) with silenced expression of AMPs and Relish, a transcription factor of the IMD pathway, exhibited significantly reduced survival rates than the control fly. Our results suggest that AMPs regulated by the IMD pathway play an important role in the chemical tolerance of D. melanogaster, and these flies are adapted to their habitats by physiological response.

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.

Insights into the Divergence of Chinese Ips Bark Beetles during Evolutionary Adaptation

Simple Summary

Bark beetle species of the genus Ips are among the major pests of Chinese conifer forests. Based on mitochondrial genome and SNP, we investigated the phylogenetic relationships and evolutionary trends of 19 populations of six Ips species that had serious outbreaks in recent years. Our results demonstrated the relationships between Ips evolution and host plants, pheromones, and altitudinal differences, and provided new insights into the mechanism of adaptive evolution of Ips bark beetles.

Abstract

Many bark beetles of the genus Ips are economically important insect pests that cause severe damage to conifer forests worldwide. In this study, sequencing the mitochondrial genome and restriction site-associated DNA of Ips bark beetles helps us understand their phylogenetic relationships, biogeographic history, and evolution of ecological traits (e.g., pheromones and host plants). Our results show that the same topology in phylogenetic trees constructed in different ways (ML/MP/BI) and with different data (mtDNA/SNP) helps us to clarify the phylogenetic relationships between Chinese Ips bark beetle populations and Euramerican species and their higher order clades; Ips bark beetles are polyphyletic. The structure of the mitochondrial genome of Ips bark beetles is similar and conserved to some extent, especially in the sibling species Ips typographus and Ips nitidus. Genetic differences among Ips species are mainly related to their geographic distribution and different hosts. The evolutionary pattern of aggregation pheromones of Ips species reflects their adaptations to the environment and differences among hosts in their evolutionary process. The evolution of Ips species is closely related to the uplift of the Qinghai-Tibet Plateau and host switching. Our study addresses the evolutionary trend and phylogenetic relationships of Ips bark beetles in China, and also provides a new perspective on the evolution of bark beetles and their relationships with host plants and pheromones.

Genetic Linkage and Physical Mapping for an Oyster Mushroom (Pleurotus cornucopiae) and Quantitative Trait Locus Analysis for Cap Color

Oyster mushrooms are grown commercially worldwide, especially in many developing countries, for their easy cultivation and high biological efficiency. Pleurotus cornucopiae is one of the main oyster mushroom species because of its gastronomic value and nutraceutical properties. Cap color is an important trait, since consumers prefer dark mushrooms, which are now represented by only a small portion of the commercial varieties. Breeding efforts are required to improve quality-related traits to satisfy various demands of consumers. Here, we present a saturated genetic linkage map of P. cornucopiae constructed by using a segregating population of 122 monokaryons and 3,449 single nucleotide polymorphism (SNP) markers generated by the 2b-RAD approach. The map contains 11 linkage groups covering 961.6 centimorgans (cM), with an average marker spacing of 0.27 cM. The genome of P. cornucopiae was de novo sequenced, resulting in 425 scaffolds (>1,000 bp) with a total genome size of 35.1 Mb. The scaffolds were assembled to the pseudochromosome level with the assistance of the genetic linkage map. A total of 97% SNP markers (3,357) were physically localized on 140 scaffolds that were assigned to 11 pseudochromosomes, with a total of 32.5 Mb, representing 92.5% of the whole genome. Six quantitative trait loci (QTL) controlling cap color of P. cornucopiae were detected, accounting for a total phenotypic variation of 65.6%, with the highest value for the QTL on pseudochromosome 5 (18%). The results of our study provide a solid base for marker-assisted breeding for agronomic traits and especially for studies on biological mechanisms controlling cap color in oyster mushrooms. IMPORTANCE Oyster mushrooms are produced and consumed all over the world. Pleurotus cornucopiae is one of the main oyster mushroom species. Dark-cap oyster mushrooms are becoming more and more popular with consumers, but dark varieties are rare on the market. Prerequisites for efficient breeding programs are the availability of high-quality whole genomes and genetic linkage maps. Genetic studies to fulfill some of these prerequisites have hardly been done for P. cornucopiae. In this study, we de novo sequenced the genome and constructed a saturated genetic linkage map for P. cornucopiae. The genetic linkage map was effectively used to assist the genome assembly and identify QTL that genetically control the trait cap color. As well, the genome characteristics of P. cornucopiae were compared to the closely related species Pleurotus ostreatus. The results provided a basis for understanding the genetic background and marker-assisted breeding of this economically important mushroom species.

A framework for studying behavioral evolution by reconstructing ancestral repertoires

Although different animal species often exhibit extensive variation in many behaviors, typically scientists examine one or a small number of behaviors in any single study. Here, we propose a new framework to simultaneously study the evolution of many behaviors. We measured the behavioral repertoire of individuals from six species of fruit flies using unsupervised techniques and identified all stereotyped movements exhibited by each species. We then fit a Generalized Linear Mixed Model to estimate the intra- and inter-species behavioral covariances, and, by using the known phylogenetic relationships among species, we estimated the (unobserved) behaviors exhibited by ancestral species. We found that much of intra-specific behavioral variation has a similar covariance structure to previously described long-time scale variation in an individual's behavior, suggesting that much of the measured variation between individuals of a single species in our assay reflects differences in the status of neural networks, rather than genetic or developmental differences between individuals. We then propose a method to identify groups of behaviors that appear to have evolved in a correlated manner, illustrating how sets of behaviors, rather than individual behaviors, likely evolved. Our approach provides a new framework for identifying co-evolving behaviors and may provide new opportunities to study the mechanistic basis of behavioral evolution.

Mass spectrometric characterisation of the major peptides of the male ejaculatory duct, including a glycopeptide with an unusual zwitterionic glycosylation

Peptides present in the seminal fluid of Drosophila melanogaster can function as antimicrobial agents, enzyme inhibitors and as pheromones that elicit physiological and behavioural responses in the post-mated female. Understanding the molecular interactions by which these peptides influence reproduction requires detailed knowledge of their molecular structures. However, this information is often lacking and cannot be gleaned from just gene sequences and standard proteomic data. We now report the native structures of four seminal fluid peptides (andropin, CG42782, Met75C and Acp54A1) from the ejaculatory duct of male D. melanogaster. The mature CG42782, Met75C and Acp54A1 peptides each have a cyclic structure formed by a disulfide bond, which will reduce conformational freedom and enhance metabolic stability. In addition, the presence of a penultimate Pro in CG42782 and Met75C will help prevent degradation by carboxypeptidases. Met75C has undergone more extensive post-translational modifications with the formation of an N-terminal pyroglutamyl residue and the attachment of a mucin-like O-glycan to the side chain of Thr₄. Both of these modifications are expected to further enhance the stability of the secreted peptide. The glycan has a rare zwitterionic structure comprising an O-linked N-acetyl hexosamine, a hexose and, unusually, phosphoethanolamine. A survey of various genomes showed that andropin, CG42782, and Acp54A1 are relatively recent genes and are restricted to the melanogaster subgroup. Met75C, however, was also found in members of the obscura species groups and in Scaptodrosophila lebanonensis. Andropin is related to the cecropin gene family and probably arose by tandem gene duplication, whereas CG42782, Met75C and Acp54A1 possibly emerged de novo. We speculate that the post-translational modifications that we report for these gene products will be important not only for a biological function, but also for metabolic stability and might also facilitate transport across tissue barriers, such as the blood-brain barrier of the female insect. BIOLOGICAL SIGNIFICANCE: Seminal fluid peptides of D. melanogaster function as antimicrobials, enzyme inhibitors and as pheromones, eliciting physiological and behavioural responses in the post-mated female. A fuller understanding of how these peptides influence reproduction requires knowledge not only of their primary structure, but also of their post-translational modification. However, this information is often lacking and difficult to glean from standard proteomic data. The reported modifications, including the unusual glycosylation, adds much to our knowledge of this important class of peptides in this model organism, par excellence.

The UDP-Glycosyltransferase Family in Drosophila melanogaster: Nomenclature Update, Gene Expression and Phylogenetic Analysis

UDP-glycosyltransferases (UGTs) are important conjugation enzymes found in all kingdoms of life, catalyzing a sugar conjugation with small lipophilic compounds and playing a crucial role in detoxification and homeostasis. The UGT gene family is defined by a signature motif in the C-terminal domain where the uridine diphosphate (UDP)-sugar donor binds. UGTs have been identified in a number of insect genomes over the last decade and much progress has been achieved in characterizing their expression patterns and molecular functions. Here, we present an update of the complete repertoire of UGT genes in Drosophila melanogaster and provide a brief overview of the latest research in this model insect. A total of 35 UGT genes are found in the D. melanogaster genome, localized to chromosomes 2 and 3 with a high degree of gene duplications on the chromosome arm 3R. All D. melanogaster UGT genes have now been named in FlyBase according to the unified UGT nomenclature guidelines. A phylogenetic analysis of UGT genes shows lineage-specific gene duplications. Analysis of anatomical and induced gene expression patterns demonstrate that some UGT genes are differentially expressed in various tissues or after environmental treatments. Extended searches of UGT orthologs from 18 additional Drosophila species reveal a diversity of UGT gene numbers and composition. The roles of Drosophila UGTs identified to date are briefly reviewed, and include xenobiotic metabolism, nicotine resistance, olfaction, cold tolerance, sclerotization, pigmentation, and immunity. Together, the updated genomic information and research overview provided herein will aid further research in this developing field.

Identification of transcriptional responsive genes to acetic acid, ethanol, and 2-phenylethanol exposure in Drosophila melanogaster

The fruit fly, Drosophila melanogaster, is predominantly found in overripe, rotten, fermenting, or decaying fruits and is constantly exposed to chemical stressors such as acetic acid, ethanol, and 2-phenylethanol. D. melanogaster has been employed as a model system for studying the molecular bases of various types of chemical-induced tolerance. Expression profiling using Illumina sequencing has been performed for identifying changes in gene expression that may be associated with evolutionary adaptation to exposure of acetic acid, ethanol, and 2-phenylethanol. We identified a total of 457 differentially expressed genes that may affect sensitivity or tolerance to three chemicals in the chemical treatment group as opposed to the control group. Gene-set enrichment analysis revealed that the genes involved in metabolism, multicellular organism reproduction, olfaction, regulation of signal transduction, and stress tolerance were over-represented in response to chemical exposure. Furthermore, we also detected a coordinated upregulation of genes in the Toll- and Imd-signaling pathways after the chemical exposure. Quantitative reverse transcription PCR analysis revealed that the expression levels of nine genes within the set of genes identified by RNA sequencing were up- or downregulated owing to chemical exposure. Taken together, our data suggest that such differentially expressed genes are coordinately affected by chemical exposure. Transcriptional analyses after exposure of D. melanogaster with three chemicals provide unique insights into subsequent functional studies on the mechanisms underlying the evolutionary adaptation of insect species to environmental chemical stressors.

Genome Wide Analysis for Growth at Two Growth Stages in A New Fast-Growing Common Carp Strain (Cyprinus carpio L.)

In order to identify candidate genes or loci associated with growth performance of the newly established common carp strain, Xinlong, we conducted a genome-wide association analysis using 2b-RAD technology on 123 individuals. We constructed two sets of libraries associated with growth-related parameters (weight, length, width and depth) measured at two different grow-out stages. Among the 413,059 SNPs identified using SOAP SNP calling, 147,131 were tested for GWAS after quality filtering. Finally, 39 overlapping SNPs, assigned to four genomic locations, were associated with growth traits in two stages. These loci were assigned to functional classes related to immune response, response to stress, neurogenesis, cholesterol metabolism and development, and proliferation and differentiation of cells. By overlapping results of Plink and EMMAX analyses, we identified three genes: TOX, PLK2 and CD163 (both methods P < 0.05). Our study results could be used for marker-assisted selection to further improve the growth of the Xinlong strain, and illustrate that largely different sets of genes drive the growth of carp in the early and late grow-out stages.

Helping decision making for reliable and cost-effective 2b-RAD sequencing and genotyping analyses in non-model species

High-throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b-RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b-RAD protocols on non-model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b-RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade-off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective-base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b-RAD protocols on non-model organisms with different genome sizes, helping decision-making for a reliable and cost-effective genotyping.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

Molecular evolution of the sex peptide network in Drosophila

Successful reproduction depends on interactions between numerous proteins beyond those involved directly in gamete fusion. Although such reproductive proteins evolve in response to sexual selection pressures, how networks of interacting proteins arise and evolve as reproductive phenotypes change remains an open question. Here, we investigated the molecular evolution of the 'sex peptide network' of Drosophila melanogaster, a functionally well-characterized reproductive protein network. In this species, the peptide hormone sex peptide (SP) and its interacting proteins cause major changes in female physiology and behaviour after mating. In contrast, females of more distantly related Drosophila species do not respond to SP. In spite of these phenotypic differences, we detected orthologs of all network proteins across 22 diverse Drosophila species and found evidence that most orthologs likely function in reproduction throughout the genus. Within SP-responsive species, we detected the recurrent, adaptive evolution of several network proteins, consistent with sexual selection acting to continually refine network function. We also found some evidence for adaptive evolution of several proteins along two specific phylogenetic lineages that correspond with increased expression of the SP receptor in female reproductive tracts or increased sperm length, respectively. Finally, we used gene expression profiling to examine the likely degree of functional conservation of the paralogs of an SP network protein that arose via gene duplication. Our results suggest a dynamic history for the SP network in which network members arose before the onset of robust SP-mediated responses and then were shaped by both purifying and positive selection.

"What You Need, Baby, I Got It": Transposable Elements as Suppliers of Cis-Operating Sequences in Drosophila

Transposable elements (TEs) are constitutive components of both eukaryotic and prokaryotic genomes. The role of TEs in the evolution of genes and genomes has been widely assessed over the past years in a variety of model and non-model organisms. Drosophila is undoubtedly among the most powerful model organisms used for the purpose of studying the role of transposons and their effects on the stability and evolution of genes and genomes. Besides their most intuitive role as insertional mutagens, TEs can modify the transcriptional pattern of host genes by juxtaposing new cis-regulatory sequences. A key element of TE biology is that they carry transcriptional control elements that fine-tune the transcription of their own genes, but that can also perturb the transcriptional activity of neighboring host genes. From this perspective, the transposition-mediated modulation of gene expression is an important issue for the short-term adaptation of physiological functions to the environmental changes, and for long-term evolutionary changes. Here, we review the current literature concerning the regulatory and structural elements operating in cis provided by TEs in Drosophila. Furthermore, we highlight that, besides their influence on both TEs and host genes expression, they can affect the chromatin structure and epigenetic status as well as both the chromosome's structure and stability. It emerges that Drosophila is a good model organism to study the effect of TE-linked regulatory sequences, and it could help future studies on TE-host interactions in any complex eukaryotic genome.

Chemosensation and Evolution of Drosophila Host Plant Selection

The ability to respond to chemosensory cues is critical for survival of most organisms. Among insects, Drosophila melanogaster has the best characterized olfactory system, and the availability of genome sequences of 30 Drosophila species provides an ideal scenario for studies on evolution of chemosensation. Gene duplications of chemoreceptor genes allow for functional diversification of the rapidly evolving chemoreceptor repertoire. Although some species of the genus Drosophila are generalists for host plant selection, rapid evolution of olfactory receptors, gustatory receptors, odorant-binding proteins, and cytochrome P450s has enabled diverse host specializations of different members of the genus. Here, I review diversification of the chemoreceptor repertoire among members of the genus Drosophila along with co-evolution of detoxification mechanisms that may have enabled occupation of diverse host plant ecological niches.

Evolutionary history and classification of Micropia retroelements in Drosophilidae species

Transposable elements (TEs) have the main role in shaping the evolution of genomes and host species, contributing to the creation of new genes and promoting rearrangements frequently associated with new regulatory networks. Support for these hypotheses frequently results from studies with model species, and Drosophila provides a great model organism to the study of TEs. Micropia belongs to the Ty3/Gypsy group of long terminal repeats (LTR) retroelements and comprises one of the least studied Drosophila transposable elements. In this study, we assessed the evolutionary history of Micropia within Drosophilidae, while trying to assist in the classification of this TE. At first, we performed searches of Micropia presence in the genome of natural populations from several species. Then, based on searches within online genomic databases, we retrieved Micropia-like sequences from the genomes of distinct Drosophilidae species. We expanded the knowledge of Micropia distribution within Drosophila species. The Micropia retroelements we detected consist of an array of divergent sequences, which we subdivided into 20 subfamilies. Even so, a patchy distribution of Micropia sequences within the Drosophilidae phylogeny could be identified, with incongruences between the species phylogeny and the Micropia phylogeny. Comparing the pairwise synonymous distance (dS) values between Micropia and three host nuclear sequences, we found several cases of unexpectedly high levels of similarity between Micropia sequences in divergent species. All these findings provide a hypothesis to the evolution of Micropia within Drosophilidae, which include several events of vertical and horizontal transposon transmission, associated with ancestral polymorphisms and recurrent Micropia sequences diversification.

Analysis of genome-wide SNPs based on 2b-RAD sequencing of pooled samples reveals signature of selection in different populations of Haemonchus contortus

The parasitic nematode Haemonchus contortus is one of the world's most important parasites of small ruminants that causes significant economic losses to the livestock sector. The population structure and selection in its various strains are poorly understood. No study so far compared its different populations using genome-wide data. Here, we focused on different geographic populations of H. contours from China (Tibet, TB; Hubei, HB; Inner Mongolia, IM; Sichuan, SC), UK and Australia (AS), using genome-wide population-genomic approaches, to explore genetic diversity, population structure and selection. We first performed next-generation high-throughput 2b RAD pool sequencing using Illumina technology, and identified single-nucleotide polymorphisms (SNPs) in all the strains. We identified 75,187 SNPs for TB, 82,271 for HB, 82,420 for IM, 79,803 for SC, 83,504 for AS and 78,747 for UK strain. The SNPs revealed low-nucleotide diversity (pi= 0.0092-0.0133) within each strain, and a significant differentiation level (average F_st = 0.34264) among them. Chinese populations TB and SC, along with the UK strain, were more divergent populations. Chinese populations IM and HB showed affinities to the Australian strain. We then analysed signature of selection and detected 44 (UK) and 03 (AS) private selective sweeps containing 49 and 05 genes, respectively. Finally, we performed the functional annotation of selective sweeps and proposed biological significance to signature of selection. Our data suggest that 2b-RAD pool sequencing can be used to assess the signature of selection in H. contortus.

TRF4, the novel TBP-related protein of Drosophila melanogaster, is concentrated at the endoplasmic reticulum and copurifies with proteins participating in the processes associated with endoplasmic reticulum

Understanding the functions of TBP-related factors is essential for studying chromatin assembly and transcription regulation in higher eukaryotes. The novel TBP-related protein-coding gene, trf4, was described in Drosophila melanogaster. trf4 is found only in Drosophila and has likely originated in Drosophila common ancestor. TRF4 protein has a distant homology with TBP and TRF2 in the region of TBP-like domain and is evolutionarily conserved among distinct Drosophila species, which indicates its functional significance. TRF4 is widely expressed in D. melanogaster with high levels of its expression being observed in testes. Interestingly enough, TRF4 has become a cytoplasmic protein having lost nuclear localization signal sequence. TRF4 is concentrated at the endoplasmic reticulum (ER) and copurifies with the proteins participating in the ER-associated processes. We suggest that trf4 gene is an example of homolog neofunctionalization by protein subcellular relocalization pathway, where the subcellular relocalization of gene product of duplicated gene leads to the new functions in ER-associated processes.

Identification and characterization of novel conserved RNA structures in Drosophila

BACKGROUND

Comparative genomics approaches have facilitated the discovery of many novel non-coding and structured RNAs (ncRNAs). The increasing availability of related genomes now makes it possible to systematically search for compensatory base changes - and thus for conserved secondary structures - even in genomic regions that are poorly alignable in the primary sequence. The wealth of available transcriptome data can add valuable insight into expression and possible function for new ncRNA candidates. Earlier work identifying ncRNAs in Drosophila melanogaster made use of sequence-based alignments and employed a sliding window approach, inevitably biasing identification toward RNAs encoded in the more conserved parts of the genome.

RESULTS

To search for conserved RNA structures (CRSs) that may not be highly conserved in sequence and to assess the expression of CRSs, we conducted a genome-wide structural alignment screen of 27 insect genomes including D. melanogaster and integrated this with an extensive set of tiling array data. The structural alignment screen revealed ∼30,000 novel candidate CRSs at an estimated false discovery rate of less than 10%. With more than one quarter of all individual CRS motifs showing sequence identities below 60%, the predicted CRSs largely complement the findings of sliding window approaches applied previously. While a sixth of the CRSs were ubiquitously expressed, we found that most were expressed in specific developmental stages or cell lines. Notably, most statistically significant enrichment of CRSs were observed in pupae, mainly in exons of untranslated regions, promotors, enhancers, and long ncRNAs. Interestingly, cell lines were found to express a different set of CRSs than were found in vivo. Only a small fraction of intergenic CRSs were co-expressed with the adjacent protein coding genes, which suggests that most intergenic CRSs are independent genetic units.

CONCLUSIONS

This study provides a more comprehensive view of the ncRNA transcriptome in fly as well as evidence for differential expression of CRSs during development and in cell lines.

Combined QTL and Genome Scan Analyses With the Help of 2b-RAD Identify Growth-Associated Genetic Markers in a New Fast-Growing Carp Strain

Common carp is one of the oldest and most popular cultured freshwater fish species both globally and in China. In a previous study, we used a carp strain with a long breeding tradition in China, named Huanghe, to create a new fast-growing strain by selection for fast growth for 6 years. The growth performance at 8 months of age has been improved by 20.84%. To achieve this, we combined the best linear unbiased prediction with marker-assisted selection techniques. Recent progress in genome-wide association studies and genomic selection in livestock breeding inspired common carp breeders to consider genome-based breeding approaches. In this study, we developed a 2b-RAD sequence assay as a means of investigating the quantitative trait loci in common carp. A total of 4,953,017,786 clean reads were generated for 250 specimens (average reads/specimen = 19,812,071) with BsaXI Restriction Enzyme. From these, 56,663 SNPs were identified, covering 50 chromosomes and 3,377 scaffolds. Principal component analysis indicated that selection and control groups are relatively clearly distinct. Top 1% of Fst values was selected as the threshold signature of artificial selection. Among the 244 identified loci, genes associated with sex-related factors and nutritional metabolism (especially fat metabolism) were annotated. Eighteen QTL were associated with growth parameters. Body length at 3 months of age and body weight (both at 3 and 8 months) were controlled by polygenic effects, but body size (length, depth, width) at 8 months of age was controlled mainly by several loci with major effects. Importantly, a single shared QTL (IGF2 gene) partially controlled the body length, depth, and width. By merging the above results, we concluded that mainly the genes related to neural pathways, sex and fatty acid metabolism contributed to the improved growth performance of the new Huanghe carp strain. These findings are one of the first investigations into the potential use of genomic selection in the breeding of common carp. Moreover, our results show that combining the Fst, QTL mapping and CRISPR–Cas9 methods can be an effective way to identify important novel candidate molecular markers in economic breeding programs.

Lima A, Nakamura AM, Fernandes F, Sobrinho I, de Brito RA. Evidence for Introgression Among Three Species of the Anastrepha fraterculus Group, a Radiating Species Complex of Fruit Flies

Introgression should no longer be considered as rare a phenomenon as once thought, since several studies have recently documented gene flow between closely related and radiating species. Here, we investigated evolutionary relationships among three closely related species of fruit flies of the Anastrepha fraterculus group (Anastrepha fraterculus, A. obliqua and A. sororcula). We sequenced a set of 20 genes and implemented a combined populational and phylogenetic inference with a model selection approach by an ABC framework in order to elucidate the demographic history of these species. The phylogenetic histories inferred from most genes showed a great deal of discordance and substantial shared polymorphic variation. The analysis of several population and speciation models reveal that this shared variation is better explained by introgression rather than convergence by parallel mutation or incomplete lineage sorting. Our results consistently showed these species evolving under an isolation with migration model experiencing a continuous and asymmetrical pattern of gene flow involving all species pairs, even though still showed a more closely related relationship between A. fraterculus and A. sororcula when compared with A. obliqua. This suggests that these species have been exchanging genes since they split from their common ancestor ∼2.6 MYA ago. We also found strong evidence for recent population expansion that appears to be consequence of anthropic activities affecting host crops of fruit flies. These findings point that the introgression here found may have been driven by genetic drift and not necessary by selection, which has implications for tracking and managing fruit flies.

A genetic linkage map of Pleurotus tuoliensis integrated with physical mapping of the de novo sequenced genome and the mating type loci

Background

Pleurotus tuoliensis (Bailinggu) is a commercially cultivated mushroom species with an increasing popularity in China and other Asian countries. Commercial profits are now low, mainly due to a low yield, long cultivation period and sensitivity to diseases. Breeding efforts are thus required to improve agronomical important traits. Developing saturated genetic linkage and physical maps is a start for applying genetic and molecular approaches to accelerate the precise breeding programs.

Results

Here we present a genetic linkage map for P. tuoliensis constructed by using 115 haploid monokaryons derived from a hybrid strain H6. One thousand one hundred and eighty-two SNP markers developed by 2b–RAD (type IIB restriction-site associated DNA) approach were mapped to 12 linkage groups. The map covers 1073 cM with an average marker spacing of 1.0 cM. The genome of P. tuoliensis was de novo sequenced as 40.8 Mb and consisted of 500 scaffolds (>500 bp), which showed a high level of colinearity to the genome of P. eryngii var. eryngii. A total of 97.4% SNP markers (1151) were physically localized on 78 scaffolds, and the physical length of these anchored scaffolds were 33.9 Mb representing 83.1% of the whole genome. Mating type loci A and B were mapped on separate linkage groups and identified physically on the assembled genomes. Five putative pheromone receptors and two putative pheromone precursors were identified for the mating type B locus.

Conclusions

This study reported a first genetic linkage map integrated with physical mapping of the de novo sequenced genome and the mating type loci of an important cultivated mushroom in China, P. tuoliensis. The de novo sequenced and annotated genome, assembled using a 2b–RAD generated linkage map, provides a basis for marker-assisted breeding of this economic important mushroom species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4421-z) contains supplementary material, which is available to authorized users.

A Horizontally Transferred Autonomous Helitron Became a Full Polydnavirus Segment in Cotesia vestalis

Bracoviruses associate symbiotically with thousands of parasitoid wasp species in the family Braconidae, working as virulence gene vectors, and allowing the development of wasp larvae within hosts. These viruses are composed of multiple DNA circles that are packaged into infective particles, and injected together with wasp's eggs during parasitization. One of the viral segments of Cotesia vestalis bracovirus contains a gene that has been previously described as a helicase of unknown origin. Here, we demonstrate that this gene is a Rep/Helicase from an intact Helitron transposable element that covers the viral segment almost entirely. We also provide evidence that this element underwent at least two horizontal transfers, which appear to have occurred consecutively: first from a Drosophila host ancestor to the genome of the parasitoid wasp C. vestalis and its bracovirus, and then from C. vestalis to a lepidopteran host (Bombyx mori). Our results reinforce the idea of parasitoid wasps as frequent agents of horizontal transfers in eukaryotes. Additionally, this Helitron-bracovirus segment is the first example of a transposable element that effectively became a whole viral circle.

Unveiling the Identity of Wenwan Walnuts and Phylogenetic Relationships of Asian Juglans Species Using Restriction Site-Associated DNA-Sequencing

Juglans species have considerable ecological and economic value worldwide. In China, Wenwan walnuts have been collected by aristocrats and noblemen for more than 2000 years. As a diversity center of Asian Juglans, five species are widely distributed in China. The most famous of these is Mahetao (J. hopeiensis), which is an uncharacterized species that is mostly cultivated. Wild J. hopeiensis individuals are very rare and are endemic to Hebei Province. Because of the minimal variations in previously used molecular markers and the heterogeneity between chloroplast and nuclear genomes, determining the phylogenetic relationships among the Juglans species has been challenging, and has hindered subsequent evolutionary inferences. In this study, we collected enough materials for both cultivated and wild Mahetao to construct well-resolved phylogenetic trees for Asian Juglans species. We used a high-throughput genome-wide restriction site-associated DNA sequencing method. Consequently, the identity of J. hopeiensis has been clearly resolved. Our results indicate that J. hopeiensis is a hybrid of J. regia and J. mandshurica. However, J. hopeiensis, J. regia and J. sigillata should be considered as a single species from section Juglans. Additionally, J. ailantifolia, J. cathayensis, and J. mandshurica likely represent one species from section Cardiocaryon according to morphological and molecular studies. These results are supported by population structure analysis and morphological comparison. We propose that J. hopeiensis trees growing in the wild should be conserved because of the economic value of their nuts. These trees may be of particular importance to impoverished communities. Furthermore, they may serve as a valuable genetic resource relevant for enhancing the production of edible walnuts. The 2b-RAD method is a viable option for future phylogenetic studies of Juglans species as well as other plant species.

Comparative analysis of behavioral and transcriptional variation underlying CO2 sensory neuron function and development in Drosophila

Carbon dioxide is an important environmental cue for many insects, regulating many behaviors including some that have direct human impacts. To further improve our understanding of how this system varies among closely related insect species, we examined both the behavioral response to CO₂ as well as the transcriptional profile of key developmental regulators of CO₂ sensory neurons in the olfactory system across the Drosophila genus. We found that CO₂ generally evokes repulsive behavior across most of the Drosophilids we examined, but this behavior has been lost or reduced in several lineages. Comparisons of transcriptional profiles from the developing and adult antennae for subset these species suggest that behavioral differences in some species may be due to differences in the expression of the CO₂ co-receptor Gr63a. Furthermore, these differences in Gr63a expression are correlated with changes in the expression of a few genes known to be involved in the development of the CO₂ circuit, namely dac, an important regulator of sensilla fate for sensilla that house CO₂ ORNs, and mip120, a member of the MMB/dREAM epigenetic regulatory complex that regulates CO₂ receptor expression. In contrast, most of the other known structural, molecular, and developmental components of the peripheral Drosophila CO₂ olfactory system seem to be well-conserved across all examined lineages. These findings suggest that certain components of CO₂ sensory ORN development may be more evolutionarily labile, and may contribute to differences in CO₂-evoked behavioral responses across species.

Genome-wide SNP analysis using 2b-RAD sequencing identifies the candidate genes putatively associated with resistance to ivermectin in Haemonchus contortus

Background

The excessive and uncontrolled use of anthelmintics, e.g. ivermectin (IVM) for the treatment of livestock parasites has led to widespread resistance in gastrointestinal nematodes, such as Haemonchus contortus. There is an urgent need for better management of drug-use in nematode control and development of novel anthelmintics. Discovery and identification of anthelmintic resistance-associate molecules/markers can provide a basis for rational anthelmintics-use and development of novel drugs. Recent studies have shown that ivermectin resistance in H. contortus is likely to be multi-genic in nature except for several genes coding for IVM target and efflux pump. However, no other IVM resistance-associated genes were characterized by conventional methods or strategies. In the present study we adopted a new strategy, i.e. using genome-wide single nucleotide polymorphism (SNP) analysis based on 2b-RAD sequencing, for discovering SNPs markers across the genomes in both IVM susceptible and resistant isolates of H. contortus and identifying potential IVM resistance-associated genes.

Results

We discovered 2962 and 2667 SNPs within both susceptible and resistant strains of H. contortus, respectively. A relative lower and similar genetic variations were observed within both resistant and susceptible strains (average π values were equal to 0.1883 and 0.1953, respectively); whereas a high genetic variation was found across both strains (average π value was equal to 0.3899). A significant differentiation across 2b-RAD tags nucleotide sites was also observed between the two strains (average F_ST value was equal to 0.3076); the larger differences in average F_ST were observed at SNPs loci between coding and noncoding (including intronic) regions. Comparison between resistant and susceptible strains revealed that 208 SNPs loci exhibited significantly elevated F_ST values, 24 SNPs of those loci were located in the CDS regions of the nine genes and were likely to have signature of IVM directional selection. Seven of the nine candidate genes were predicted to code for some functional proteins such as potential IVM target and/or efflux pump proteins, component proteins of receptor complex in membrane on neuromuscular cells, and transcriptional regulation proteins. Those genes might be involved in resistance to IVM.

Conclusions

Our data suggest that candidate genes putatively associated with resistance to IVM in H. contortus may be identified by genome-wide SNP analysis using 2b-RAD sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1959-6) contains supplementary material, which is available to authorized users.

Gene Ages, Nomenclatures, and Functional Diversification of the Methuselah/Methuselah-Like GPCR Family in Drosophila and Tribolium

Affecting lifespan regulation and oxidative stress resistance, the G-protein coupled receptor (GPCR) gene methuselah (mth) plays important roles in the life history of Drosophila melanogaster. Substantial progress has been made in elucidating the molecular pathways by which mth affects these traits, yet conflicting ideas exist as to how old these genetic interactions are as well as how old the mth gene itself is. Root to these issues is the complex gene family history of the Mth/Mthl GPCR family, which experienced independent expansions in a variety of animal clades, leading to at least six subfamilies in insects. Within insects, drosophilid flies stand out by possessing up to three times more Mth/Mthl receptors due to the expansion of a single subfamily, the mth superclade subfamily, which contains an even younger subfamily introduced here as the melanogaster subgroup subfamily. As a result, most of the 16 Mth/Mthl receptors of D. melanogaster are characterized by n:1 orthology relationships to singleton mth superclade homologs in nondrosophilid species. This challenge is exacerbated by the inconsistent naming of Mth/Mthl orthologs across species. To consolidate this situation, we review established ortholog relationships among insect Mth/Mthl receptors, clarify the gene nomenclatures in two important satellite model species, the fruit fly relative D. virilis and the red flour beetle Tribolium castaneum, and discuss the genetic and functional evolution of the D. melanogaster Mth GPCR.

Serial sequencing of isolength RAD tags for cost-efficient genome-wide profiling of genetic and epigenetic variations

Isolength restriction site-associated DNA (isoRAD) sequencing is a very simple but powerful approach that was originally developed for genome-wide genotyping at minimal labor and cost, and it has recently extended its applicability to allow quantification of DNA methylation levels. The isoRAD method is distinct from other genotyping-by-sequencing (GBS) methods because of its use of special restriction enzymes to produce isolength tags (32-36 bp), and sequencing of these uniform tags can bring many benefits. However, the relatively short tags produced by the original protocol are mostly suited to single-end (SE) sequencing (36-50 bp), and therefore they cannot efficiently match the gradually increased sequencing capacity of next-generation sequencing (NGS) platforms. To address this issue, we describe an advanced protocol that allows the preparation of five concatenated isoRAD tags for Illumina paired-end (PE) sequencing (100-150 bp). The configuration of the five concatenated tags is highly flexible, and can be defined by users to work with a desired combination of samples and/or restriction enzymes to suit specific research purposes. In comparison with the original protocol, the advanced protocol has an additional digestion and ligation step, and library preparation can be completed in ∼8 h.

Genetics on the Fly: A Primer on the Drosophila Model System

Fruit flies of the genus Drosophila have been an attractive and effective genetic model organism since Thomas Hunt Morgan and colleagues made seminal discoveries with them a century ago. Work with Drosophila has enabled dramatic advances in cell and developmental biology, neurobiology and behavior, molecular biology, evolutionary and population genetics, and other fields. With more tissue types and observable behaviors than in other short-generation model organisms, and with vast genome data available for many species within the genus, the fly's tractable complexity will continue to enable exciting opportunities to explore mechanisms of complex developmental programs, behaviors, and broader evolutionary questions. This primer describes the organism's natural history, the features of sequenced genomes within the genus, the wide range of available genetic tools and online resources, the types of biological questions Drosophila can help address, and historical milestones.

How sexual selection can drive the evolution of costly sperm ornamentation

Post-copulatory sexual selection (PSS), fuelled by female promiscuity, is credited with the rapid evolution of sperm quality traits across diverse taxa. Yet, our understanding of the adaptive significance of sperm ornaments and the cryptic female preferences driving their evolution is extremely limited. Here we review the evolutionary allometry of exaggerated sexual traits (for example, antlers, horns, tail feathers, mandibles and dewlaps), show that the giant sperm of some Drosophila species are possibly the most extreme ornaments in all of nature and demonstrate how their existence challenges theories explaining the intensity of sexual selection, mating-system evolution and the fundamental nature of sex differences. We also combine quantitative genetic analyses of interacting sex-specific traits in D. melanogaster with comparative analyses of the condition dependence of male and female reproductive potential across species with varying ornament size to reveal complex dynamics that may underlie sperm-length evolution. Our results suggest that producing few gigantic sperm evolved by (1) Fisherian runaway selection mediated by genetic correlations between sperm length, the female preference for long sperm and female mating frequency, and (2) longer sperm increasing the indirect benefits to females. Our results also suggest that the developmental integration of sperm quality and quantity renders post-copulatory sexual selection on ejaculates unlikely to treat male-male competition and female choice as discrete processes.

Convergent Balancing Selection on an Antimicrobial Peptide in Drosophila

Genes of the immune system often evolve rapidly and adaptively, presumably driven by antagonistic interactions with pathogens [1-4]. Those genes encoding secreted antimicrobial peptides (AMPs), however, have failed to exhibit conventional signatures of strong adaptive evolution, especially in arthropods (e.g., [5, 6]) and often segregate for null alleles and gene deletions [3, 4, 7, 8]. Furthermore, quantitative genetic studies have failed to associate naturally occurring polymorphism in AMP genes with variation in resistance to infection [9-11]. Both the lack of signatures of positive selection in AMPs and lack of association between genotype and immune phenotypes have yielded an interpretation that AMP genes evolve under relaxed evolutionary constraint, with enough functional redundancy that variation in, or even loss of, any particular peptide would have little effect on overall resistance [12, 13]. In stark contrast to the current paradigm, we identified a naturally occurring amino acid polymorphism in the AMP Diptericin that is highly predictive of resistance to bacterial infection in Drosophila melanogaster [13]. The identical amino acid polymorphism arose in parallel in the sister species D. simulans, by independent mutation with equivalent phenotypic effect. Convergent substitutions at the same amino acid residue have evolved at least five times across the Drosophila genus. We hypothesize that the alternative alleles are maintained by balancing selection through context-dependent or fluctuating selection. This pattern of evolution appears to be common in AMPs but is invisible to conventional screens for adaptive evolution that are predicated on elevated rates of amino acid divergence.

Evaluation of the 2b-RAD method for genomic selection in scallop breeding

The recently developed 2b-restriction site-associated DNA (2b-RAD) sequencing method provides a cost-effective and flexible genotyping platform for aquaculture species lacking sufficient genomic resources. Here, we evaluated the performance of this method in the genomic selection (GS) of Yesso scallop (Patinopecten yessoensis) through simulation and real data analyses using six statistical models. Our simulation analysis revealed that the prediction accuracies obtained using the 2b-RAD markers were slightly lower than those obtained using all polymorphic loci in the genome. Furthermore, a small subset of markers obtained from a reduced tag representation (RTR) library presented comparable performance to that obtained using all markers, making RTR be an attractive approach for GS purpose. Six GS models exhibited variable performance in prediction accuracy depending on the scenarios (e.g., heritability, sample size, population structure), but Bayes-alphabet and BLUP-based models generally outperformed other models. Finally, we performed the evaluation using an empirical dataset composed of 349 Yesso scallops that were derived from five families. The prediction accuracy for this empirical dataset could reach 0.4 based on optimal GS models. In summary, the genotyping flexibility and cost-effectiveness make 2b-RAD be an ideal genotyping platform for genomic selection in aquaculture breeding programs.

Extending RAD tag analysis to microbial ecology: a comparison between MultiLocus Sequence Typing and 2b-RAD to investigate Listeria monocytogenes genetic structure

The advent of next-generation sequencing (NGS) has dramatically changed bacterial typing technologies, increasing our ability to differentiate bacterial isolates. Despite it is now possible to sequence a bacterial genome in a few days and at reasonable costs, most genetic analyses do not require whole-genome sequencing, which also remains impractical for large population samples due to the cost of individual library preparation and bioinformatics. More traditional sequencing approaches, however, such as MultiLocus Sequence Typing (mlst) are quite laborious and time-consuming, especially for large-scale analyses. In this study, a genotyping approach based on restriction site-associated (RAD) tag sequencing, 2b-RAD, was applied to characterize Listeria monocytogenes strains. To verify the feasibility of the method, an in silico analysis was performed on 30 available complete genomes. For the same set of strains, in silico mlst analysis was conducted as well. Subsequently, 2b-RAD and mlst analyses were experimentally carried out on 58 isolates collected from food samples or food-processing sites. The obtained results demonstrate that 2b-RAD predicts mlst types and often provides more detailed information on population structure than mlst. Moreover, the majority of variants differentiating identical sequence type isolates mapped against accessory fragments, thus providing additional information to characterize strains. Although mlst still represents a reliable typing method, large-scale studies on molecular epidemiology and public health, as well as bacterial phylogenetics, population genetics and biosafety could benefit of a low cost and fast turnaround time approach such as the 2b-RAD analysis proposed here.

Assessment of rival males through the use of multiple sensory cues in the fruitfly Drosophila pseudoobscura

Environments vary stochastically, and animals need to behave in ways that best fit the conditions in which they find themselves. The social environment is particularly variable, and responding appropriately to it can be vital for an animal's success. However, cues of social environment are not always reliable, and animals may need to balance accuracy against the risk of failing to respond if local conditions or interfering signals prevent them detecting a cue. Recent work has shown that many male Drosophila fruit flies respond to the presence of rival males, and that these responses increase their success in acquiring mates and fathering offspring. In Drosophila melanogaster males detect rivals using auditory, tactile and olfactory cues. However, males fail to respond to rivals if any two of these senses are not functioning: a single cue is not enough to produce a response. Here we examined cue use in the detection of rival males in a distantly related Drosophila species, D. pseudoobscura, where auditory, olfactory, tactile and visual cues were manipulated to assess the importance of each sensory cue singly and in combination. In contrast to D. melanogaster, male D. pseudoobscura require intact olfactory and tactile cues to respond to rivals. Visual cues were not important for detecting rival D. pseudoobscura, while results on auditory cues appeared puzzling. This difference in cue use in two species in the same genus suggests that cue use is evolutionarily labile, and may evolve in response to ecological or life history differences between species.

Genomic repeat abundances contain phylogenetic signal

A large proportion of genomic information, particularly repetitive elements, is usually ignored when researchers are using next-generation sequencing. Here we demonstrate the usefulness of this repetitive fraction in phylogenetic analyses, utilizing comparative graph-based clustering of next-generation sequence reads, which results in abundance estimates of different classes of genomic repeats. Phylogenetic trees are then inferred based on the genome-wide abundance of different repeat types treated as continuously varying characters; such repeats are scattered across chromosomes and in angiosperms can constitute a majority of nuclear genomic DNA. In six diverse examples, five angiosperms and one insect, this method provides generally well-supported relationships at interspecific and intergeneric levels that agree with results from more standard phylogenetic analyses of commonly used markers. We propose that this methodology may prove especially useful in groups where there is little genetic differentiation in standard phylogenetic markers. At the same time as providing data for phylogenetic inference, this method additionally yields a wealth of data for comparative studies of genome evolution.