Origin of a function by tandem gene duplication limits the evolutionary capability of its sister copy

Chromosomal-level genome of a sheet-web spider provides insight into the composition and evolution of venom

Spiders are the most abundant venomous predators in the world. Previous research related to spider venom has mostly relied on transcriptomes and proteomes, with only a few high-quality genomes available. This is far from consistent with the species diversity of spiders. In this study, we constructed a high-quality chromosome-level genome assembly of Hylyphantes graminicola, which contained 13 chromosomes, with a genome length of 931.68 Mb and scaffold N50 of 77.07 Mb. Integrating genome, transcriptome, and proteome profiling, we identified a total of 59 coding genes among nine toxin gene families. Among them, Group 7 allergen (ALL7) protein was reported in spider venom for the first time. Its coding genes had a predicted signal peptide and maintained high expression levels in the venom, suggesting that ALL7 plays an important role in venom and maybe is a type of newly discovered venom toxin in the spider. By implementing comparative genomics, we found a similar gene number of main toxin gene families in spiders and the scorpion genome with conservative evolutionary rates, indicating that these toxin genes could be an ancient (~400 million years) and a conserved "basic toolkit" for spiders and scorpions to perform primary defense functions. Obtaining high-quality chromosome-level genomes from spiders not only facilitates venom research and toxin resource application, but also can improve comparative genomic analysis in other important traits, like the evolution of silk or behavior.

Evidence for Stabilizing Selection Driving Mutational Turnover of Short Motifs in the Eukaryotic Complementary Sex Determiner (Csd) Protein

Short linear motifs (SLiMs) can play pivotal functional roles in proteins, such as targeting proteins to specific subcellular localizations, modulating the efficiency of translation and tagging proteins for degradation. Until recently we had little knowledge about SLiM evolution. Only a few amino acids in these motifs are functionally important, making them likely to evolve ex nihilo and suggesting that they can play key roles in protein evolution. Several reports now suggest that these motifs can appear and disappear while their function in the protein is preserved, a process sometimes referred to as “turnover”. However, there has been a lack of specific experiments to determine whether independently evolved motifs do indeed have the same function, which would conclusively determine whether the process of turnover actually occurs. In this study, we experimentally detected evidence for such a mutational turnover process for nuclear localization signals (NLS) during the post-duplication divergence of the Complementary sex determiner (Csd) and Feminizer (Fem) proteins in the honeybee (Apis mellifera) lineage. Experiments on the nuclear transport activity of protein segments and those of the most recent common ancestor (MRCA) sequences revealed that three new NLS motifs evolved in the Csd protein during the post-duplication divergence while other NLS motifs were lost that existed before duplication. A screen for essential and newly evolved amino acids revealed that new motifs in the Csd protein evolved by one or two missense mutations coding for lysine. Amino acids that were predating the duplication were also essential in the acquisition of the C1 motif suggesting that the ex nihilo origin was constrained by preexisting amino acids in the physical proximity. Our data support a model in which stabilizing selection maintains the constancy of nuclear transport function but allowed mutational turnover of the encoding NLS motifs.

Uneven distribution of complementary sex determiner (csd) alleles in Apis mellifera population

The complementary sex determiner (csd) gene determines the sex of the western honey bee (Apis mellifera L.). Bees that are heterozygous at the csd locus develop into females; whereas hemizygous bees develop into males. The co-occurrence of two identical csd alleles in a single diploid genome leads to the genetic death of the bee. Thus, the maintenance of csd diversity in the population is favoured. The number and distribution of csd alleles is particularly interesting in light of the recent decline in the honey bee population. In this study, we analysed the distribution of csd alleles in two Polish populations separated by about 100 km. We analysed the maternal alleles of 193 colonies and found 121 different alleles. We also analysed the distribution and frequency of the alleles, and found that they are distributed unevenly. We show that the methods that have been used so far to estimate the total worldwide number of csd alleles have significantly underestimated their diversity. We also show that the uneven distribution of csd alleles is caused by a large number of infrequent alleles, which most likely results from the fact that these alleles are generated very frequently.

The Cardamine hirsuta genome offers insight into the evolution of morphological diversity

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

Evolution of Social Insect Polyphenism Facilitated by the Sex Differentiation Cascade

The major transition to eusociality required the evolution of a switch to canalize development into either a reproductive or a helper, the nature of which is currently unknown. Following predictions from the ‘theory of facilitated variation’, we identify sex differentiation pathways as promising candidates because of their pre-adaptation to regulating development of complex phenotypes. We show that conserved core genes, including the juvenile hormone-sensitive master sex differentiation gene doublesex (dsx) and a krüppel homolog 2 (kr-h2) with putative regulatory function, exhibit both sex and morph-specific expression across life stages in the ant Cardiocondyla obscurior. We hypothesize that genes in the sex differentiation cascade evolved perception of alternative input signals for caste differentiation (i.e. environmental or genetic cues), and that their inherent switch-like and epistatic behavior facilitated signal transfer to downstream targets, thus allowing them to control differential development into morphological castes.

Division of labor into reproductive queens and helper workers in the societies of ants, bees and wasps is achieved by phenotypic plasticity, which allows individuals to embark on discrete developmental trajectories in response to variable signals. These signals can be genetic, epigenetic or environmental, thereby resembling the extreme variation in signals for sex determination across multicellular animals. We show that common developmental pathways downstream of these input signals, including the conserved sex differentiation gene doublesex, regulate sex and caste-specific phenotypic differentiation in the ant species Cardiocondyla obscurior. Many different mechanisms of gene regulation have been implicated in controlling caste-specific development in social insects but these all depend on a higher-level genetic switch. We propose that highly conserved hub genes such as dsx, which can translate variable input signals into large transcription differences using intermediate-level regulators, are tightly linked with the repeated evolutionary transition to eusociality and caste polyphenism.

Correlated expression of retrocopies and parental genes in zebrafish

Previous studies of the function and evolution of retrocopies in plants, Drosophila and non-mammalian chordates provided new insights into the origin of novel genes. However, little is known about retrocopies and their parental genes in teleosts, and it remains obscure whether there is any correlation between them. The present study aimed to characterize the spatial and temporal expression profiles of retrogenes and their parental genes based on RNA-Seq data from Danio rerio embryos and tissues from adult. Using a modified pipeline, 306 retrocopies were identified in the zebrafish genome, most of which exhibited ancient retroposition, and 76 of these showed a Ks < 2.0. Expression of a retrocopy is generally expected to present no correlation with its parental gene, as regulatory regions are not part of the retroposition event. Here, this assumption was tested based on RNA-Seq data from eight stages and thirteen tissue types of zebrafish. However, the result suggested that retrocopies displayed correlated expression with their parental genes. The level of correlation was found to decrease during embryogenesis, but to increase slightly within a tissue using Ks as the proxy for the divergence time. Tissue specificity was also observed: retrocopies were found to be expressed at a more specific level compared with their parental genes. Unlike Drosophila, which has sex chromosomes, zebrafish do not show testis-biased expression. Our study elaborated temporal and spatial patterns of expression of retrocopies in zebrafish, examined the correlation between retrocopies and parental genes and analyzed potential source of regulated elements of retrocopies, which lay a foundation for further functional study of retrocopies.

QTL Mapping of Sex Determination Loci Supports an Ancient Pathway in Ants and Honey Bees

Sex determination mechanisms play a central role in life-history characteristics, affecting mating systems, sex ratios, inbreeding tolerance, etc. Downstream components of sex determination pathways are highly conserved, but upstream components evolve rapidly. Evolutionary dynamics of sex determination remain poorly understood, particularly because mechanisms appear so diverse. Here we investigate the origins and evolution of complementary sex determination (CSD) in ants and bees. The honey bee has a well-characterized CSD locus, containing tandemly arranged homologs of the transformer gene [complementary sex determiner (csd) and feminizer (fem)]. Such tandem paralogs appear frequently in aculeate hymenopteran genomes. However, only comparative genomic, but not functional, data support a broader role for csd/fem in sex determination, and whether species other than the honey bee use this pathway remains controversial. Here we used a backcross to test whether csd/fem acts as a CSD locus in an ant (Vollenhovia emeryi). After sequencing and assembling the genome, we computed a linkage map, and conducted a quantitative trait locus (QTL) analysis of diploid male production using 68 diploid males and 171 workers. We found two QTLs on separate linkage groups (CsdQTL1 and CsdQTL2) that jointly explained 98.0% of the phenotypic variance. CsdQTL1 included two tandem transformer homologs. These data support the prediction that the same CSD mechanism has indeed been conserved for over 100 million years. CsdQTL2 had no similarity to CsdQTL1 and included a 236-kb region with no obvious CSD gene candidates, making it impossible to conclusively characterize it using our data. The sequence of this locus was conserved in at least one other ant genome that diverged >75 million years ago. By applying QTL analysis to ants for the first time, we support the hypothesis that elements of hymenopteran CSD are ancient, but also show that more remains to be learned about the diversity of CSD mechanisms.

Whether a developing embryo becomes male or female has significant downstream consequences. Depending on the species, sex can be determined by a wide variety of mechanisms. Sex determination systems can evolve rapidly, but how this occurs, and even how widespread the same mechanism is within a given taxonomic group, remains largely unknown. By experimentally mapping the sex determination architecture in the ant, Vollenhovia emeryi, we found that the well-characterized honey bee sex determination locus originated more than 100 million years ago. However, we also found an additional locus that has no homology to the first. Currently uncharacterized, this locus suggests that different species may use a variety of complementary sex determination mechanisms. Yet, core elements of the complementary sex determination machinery appear to be ancient.

Similar but not the same: insights into the evolutionary history of paralogous sex-determining genes of the dwarf honey bee Apis florea

Studying the fate of duplicated genes provides informative insight into the evolutionary plasticity of biological pathways to which they belong. In the paralogous sex-determining genes complementary sex determiner (csd) and feminizer (fem) of honey bee species (genus Apis), only heterozygous csd initiates female development. Here, the full-length coding sequences of the genes csd and fem of the phylogenetically basal dwarf honey bee Apis florea are characterized. Compared with other Apis species, remarkable evolutionary changes in the formation and localization of a protein-interacting (coiled-coil) motif and in the amino acids coding for the csd characteristic hypervariable region (HVR) are observed. Furthermore, functionally different csd alleles were isolated as genomic fragments from a random population sample. In the predicted potential specifying domain (PSD), a high ratio of πN/πS=1.6 indicated positive selection, whereas signs of balancing selection, commonly found in other Apis species, are missing. Low nucleotide diversity on synonymous and genome-wide, non-coding sites as well as site frequency analyses indicated a strong impact of genetic drift in A. florea, likely linked to its biology. Along the evolutionary trajectory of ~30 million years of csd evolution, episodic diversifying selection seems to have acted differently among distinct Apis branches. Consistently low amino-acid differences within the PSD among pairs of functional heterozygous csd alleles indicate that the HVR is the most important region for determining allele specificity. We propose that in the early history of the lineage-specific fem duplication giving rise to csd in Apis, A. florea csd stands as a remarkable example for the plasticity of initial sex-determining signals.

The evolutionary dynamics of major regulators for sexual development among Hymenoptera species

All hymenopteran species, such as bees, wasps and ants, are characterized by the common principle of haplodiploid sex determination in which haploid males arise from unfertilized eggs and females from fertilized eggs. The underlying molecular mechanism has been studied in detail in the western honey bee Apis mellifera, in which the gene complementary sex determiner (csd) acts as primary signal of the sex determining pathway, initiating female development by csd-heterozygotes. Csd arose from gene duplication of the feminizer (fem) gene, a transformer (tra) ortholog, and mediates in conjunction with transformer2 (tra2) sex-specific splicing of fem. Comparative molecular analyses identified fem/tra and its downstream target doublesex (dsx) as conserved unit within the sex determining pathway of holometabolous insects. In this study, we aim to examine evolutionary differences among these key regulators. Our main hypothesis is that sex determining key regulators in Hymenoptera species show signs of coevolution within single phylogenetic lineages. We take advantage of several newly sequenced genomes of bee species to test this hypothesis using bioinformatic approaches. We found evidences that duplications of fem are restricted to certain bee lineages and notable amino acid differences of tra2 between Apis and non-Apis species propose structural changes in Tra2 protein affecting co-regulatory function on target genes. These findings may help to gain deeper insights into the ancestral mode of hymenopteran sex determination and support the common view of the remarkable evolutionary flexibility in this regulatory pathway.

Independent evolutionary origin of fem paralogous genes and complementary sex determination in hymenopteran insects

The primary signal of sex determination in the honeybee, the complementary sex determiner (csd) gene, evolved from a gene duplication event from an ancestral copy of the fem gene. Recently, other paralogs of the fem gene have been identified in several ant and bumblebee genomes. This discovery and the close phylogenetic relationship of the paralogous gene sequences led to the hypothesis of a single ancestry of the csd genetic system of complementary sex determination in the Hymenopteran insects, in which the fem and csd gene copies evolved as a unit in concert with the mutual transfers of sequences (concerted evolution). Here, we show that the paralogous gene copies evolved repeatedly through independent gene duplication events in the honeybee, bumblebee, and ant lineage. We detected no sequence tracts that would indicate a DNA transfer between the fem and the fem1/csd genes between different ant and bee species. Instead, we found tracts of duplication events in other genomic locations, suggesting that gene duplication was a frequent event in the evolution of these genes. These and other evidences suggest that the fem1/csd gene originated repeatedly through gene duplications in the bumblebee, honeybee, and ant lineages in the last 100 million years. Signatures of concerted evolution were not detectable, implicating that the gene tree based on neutral synonymous sites represents the phylogenetic relationships and origins of the fem and fem1/csd genes. Our results further imply that the fem1 and csd gene in bumblebees, honeybees, and ants are not orthologs, because they originated independently from the fem gene. Hence, the widely shared and conserved complementary sex determination mechanism in Hymenopteran insects is controlled by different genes and molecular processes. These findings highlight the limits of comparative genomics and emphasize the requirement to study gene functions in different species and major hymenopteran lineages.

Phylogenetic distribution and evolutionary dynamics of the sex determination genes doublesex and transformer in insects

Sex determination in insects is characterized by a gene cascade that is conserved at the bottom but contains diverse primary signals at the top. The bottom master switch gene doublesex is found in all insects. Its upstream regulator transformer is present in the orders Hymenoptera, Coleoptera and Diptera, but has thus far not been found in Lepidoptera and in the basal lineages of Diptera. transformer is presumed to be ancestral to the holometabolous insects based on its shared domains and conserved features of autoregulation and sex-specific splicing. We interpret that its absence in basal lineages of Diptera and its order-specific conserved domains indicate multiple independent losses or recruitments into the sex determination cascade. Duplications of transformer are found in derived families within the Hymenoptera, characterized by their complementary sex determination mechanism. As duplications are not found in any other insect order, they appear linked to the haplodiploid reproduction of the Hymenoptera. Further phylogenetic analyses combined with functional studies are needed to understand the evolutionary history of the transformer gene among insects.

Sex determination in insects: variations on a common theme

Recent studies in a representative selection of holometabolous insects suggest that, despite diversity at the instructive level, the signal-relaying part of the sex-determining pathway is remarkably well conserved. In principle, it is composed of the transformer gene (tra), which acts as a common binary switch that transduces the selected sexual fate, female when ON, male when OFF, to the downstream effector doublesex(dsx) that controls overt sexual differentiation. An interesting recurrent feature is that tra is switched ON in the early zygote by maternally provisioned tra activity. Different male-determining signals evolved, which prevent maternal activation of zygotic tra to allow for male development. In some species, where lack of maternal activation leaves tra in the OFF state, novel female-determining signals were deployed to activate zygotic tra. It appears that both the instructive end of the pathway upstream of tra as well as the executive end downstream of dsx are primary targets of evolutionary divergence, while the transduction part seems less prone to changes. We propose that this is a feature shared with many other signaling cascades that regulate developmental fates.

Gradual molecular evolution of a sex determination switch through incomplete penetrance of femaleness

Some genes regulate phenotypes that are either present or absent. They are often important regulators of developmental switches and are involved in morphological evolution. We have little understanding of the molecular mechanisms by which these absence/presence gene functions have evolved, because the phenotype and fitness of molecular intermediate forms are unknown. Here, we studied the sex-determining switch of 14 natural sequence variants of the csd gene among 76 genotypes of the honeybee (Apis mellifera). Heterozygous genotypes (different specificities) of the csd gene determine femaleness, while hemizygous genotypes (single specificity) determine maleness. Homozygous genotypes of the csd gene (same specificity) are lethal. We found that at least five amino acid differences and length variation between Csd specificities in the specifying domain (PSD) were sufficient to regularly induce femaleness. We estimated that, on average, six pairwise amino acid differences evolved under positive selection. We also identified a natural evolutionary intermediate that showed only three amino acid length differences in the PSD relative to its parental allele. This genotype showed an intermediate fitness because it implemented lethality regularly and induced femaleness infrequently (i.e., incomplete penetrance). We suggest incomplete penetrance as a mechanism through which new molecular switches can gradually and adaptively evolve.

Nucleotide variability at its limit? Insights into the number and evolutionary dynamics of the sex-determining specificities of the honey bee Apis mellifera

Deciphering the evolutionary processes driving nucleotide variation in multiallelic genes is limited by the number of genetic systems in which such genes occur. The complementary sex determiner (csd) gene in the honey bee Apis mellifera is an informative example for studying allelic diversity and the underlying evolutionary forces in a well-described model of balancing selection. Acting as the primary signal of sex determination, diploid individuals heterozygous for csd develop into females, whereas csd homozygotes are diploid males that have zero fitness. Examining 77 of the functional heterozygous csd allele pairs, we established a combinatorical criteria that provide insights into the minimum number of amino acid differences among those pairs. Given a data set of 244 csd sequences, we show that the total number of csd alleles found in A. mellifera ranges from 53 (locally) to 87 (worldwide), which is much higher than was previously reported (20). Using a coupon-collector model, we extrapolate the presence of in total 116–145 csd alleles worldwide. The hypervariable region (HVR) is of particular importance in determining csd allele specificity, and we provide for this region evidence of high evolutionary rate for length differences exceeding those of microsatellites. The proportion of amino acids driven by positive selection and the rate of nonsynonymous substitutions in the HVR-flanking regions reach values close to 1 but differ with respect to the HVR length. Using a model of csd coalescence, we identified the high originating rate of csd specificities as a major evolutionary force, leading to an origin of a novel csd allele every 400,000 years. The csd polymorphism frequencies in natural populations indicate an excess of new mutations, whereas signs of ancestral transspecies polymorphism can still be detected. This study provides a comprehensive view of the enormous diversity and the evolutionary forces shaping a multiallelic gene.

Evolution of recombination and genome structure in eusocial insects

Eusocial Hymenoptera, such as the European honey bee, Apis mellifera, have the highest recombination rates of multicellular animals.(1) Recently, we showed(2) that a side-effect of recombination in the honey bee, GC biased gene conversion (bGC), helps maintain the unusual bimodal GC-content distribution of the bee genome by increasing GC-content in high recombination areas while low recombination areas are losing GC-content because of biased AT mutations and low rates of bGC. Although the very high recombination rate of A. mellifera makes GC-content evolution easier to study, the pattern is consistent with results found in many other species including mammals and yeast.(3) Also consistent across phyla is the association of higher genetic diversity and divergence with high GC and high recombination areas.(4) (,) (5) Finally, we showed that genes overexpressed in the brains of workers cluster in GC-rich genomic areas with the highest rates of recombination and molecular evolution.(2) In this Addendum we present a conceptual model of how eusociality and high recombination rates may co-evolve.

Duplication and concerted evolution in a master sex determiner under balancing selection

The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.

Evolutionary mechanisms driving the evolution of a large polydnavirus gene family coding for protein tyrosine phosphatases

Background

Gene duplications have been proposed to be the main mechanism involved in genome evolution and in acquisition of new functions. Polydnaviruses (PDVs), symbiotic viruses associated with parasitoid wasps, are ideal model systems to study mechanisms of gene duplications given that PDV genomes consist of virulence genes organized into multigene families. In these systems the viral genome is integrated in a wasp chromosome as a provirus and virus particles containing circular double-stranded DNA are injected into the parasitoids’ hosts and are essential for parasitism success. The viral virulence factors, organized in gene families, are required collectively to induce host immune suppression and developmental arrest. The gene family which encodes protein tyrosine phosphatases (PTPs) has undergone spectacular expansion in several PDV genomes with up to 42 genes.

Results

Here, we present strong indications that PTP gene family expansion occurred via classical mechanisms: by duplication of large segments of the chromosomally integrated form of the virus sequences (segmental duplication), by tandem duplications within this form and by dispersed duplications. We also propose a novel duplication mechanism specific to PDVs that involves viral circle reintegration into the wasp genome. The PTP copies produced were shown to undergo conservative evolution along with episodes of adaptive evolution. In particular recently produced copies have undergone positive selection in sites most likely involved in defining substrate selectivity.

Conclusion

The results provide evidence about the dynamic nature of polydnavirus proviral genomes. Classical and PDV-specific duplication mechanisms have been involved in the production of new gene copies. Selection pressures associated with antagonistic interactions with parasitized hosts have shaped these genes used to manipulate lepidopteran physiology with evidence for positive selection involved in adaptation to host targets.

Complex patterns of divergence among green-sensitive (RH2a) African cichlid opsins revealed by Clade model analyses

BACKGROUND

Gene duplications play an important role in the evolution of functional protein diversity. Some models of duplicate gene evolution predict complex forms of paralog divergence; orthologous proteins may diverge as well, further complicating patterns of divergence among and within gene families. Consequently, studying the link between protein sequence evolution and duplication requires the use of flexible substitution models that can accommodate multiple shifts in selection across a phylogeny. Here, we employed a variety of codon substitution models, primarily Clade models, to explore how selective constraint evolved following the duplication of a green-sensitive (RH2a) visual pigment protein (opsin) in African cichlids. Past studies have linked opsin divergence to ecological and sexual divergence within the African cichlid adaptive radiation. Furthermore, biochemical and regulatory differences between the RH2aα and RH2aβ paralogs have been documented. It thus seems likely that selection varies in complex ways throughout this gene family.

RESULTS

Clade model analysis of African cichlid RH2a opsins revealed a large increase in the nonsynonymous-to-synonymous substitution rate ratio (ω) following the duplication, as well as an even larger increase, one consistent with positive selection, for Lake Tanganyikan cichlid RH2aβ opsins. Analysis using the popular Branch-site models, by contrast, revealed no such alteration of constraint. Several amino acid sites known to influence spectral and non-spectral aspects of opsin biochemistry were found to be evolving divergently, suggesting that orthologous RH2a opsins may vary in terms of spectral sensitivity and response kinetics. Divergence appears to be occurring despite intronic gene conversion among the tandemly-arranged duplicates.

CONCLUSIONS

Our findings indicate that variation in selective constraint is associated with both gene duplication and divergence among orthologs in African cichlid RH2a opsins. At least some of this variation may reflect an adaptive response to differences in light environment. Interestingly, these patterns only became apparent through the use of Clade models, not through the use of the more widely employed Branch-site models; we suggest that this difference stems from the increased flexibility associated with Clade models. Our results thus bear both on studies of cichlid visual system evolution and on studies of gene family evolution in general.

Extensive pyrosequencing reveals frequent intra-genomic variations of internal transcribed spacer regions of nuclear ribosomal DNA

BACKGROUND

Internal transcribed spacer of nuclear ribosomal DNA (nrDNA) is already one of the most popular phylogenetic and DNA barcoding markers. However, the existence of its multiple copies has complicated such usage and a detailed characterization of intra-genomic variations is critical to address such concerns.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we used sequence-tagged pyrosequencing and genome-wide analyses to characterize intra-genomic variations of internal transcribed spacer 2 (ITS2) regions from 178 plant species. We discovered that mutation of ITS2 is frequent, with a mean of 35 variants per species. And on average, three of the most abundant variants make up 91% of all ITS2 copies. Moreover, we found different congeneric species share identical variants in 13 genera. Interestingly, different species across different genera also share identical variants. In particular, one minor variant of ITS2 in Eleutherococcus giraldii was found identical to the ITS2 major variant of Panax ginseng, both from Araliaceae family. In addition, DNA barcoding gap analysis showed that the intra-genomic distances were markedly smaller than those of the intra-specific or inter-specific variants. When each of 5543 variants were examined for its species discrimination efficiency, a 97% success rate was obtained at the species level.

CONCLUSIONS

Identification of identical ITS2 variants across intra-generic or inter-generic species revealed complex species evolutionary history, possibly, horizontal gene transfer and ancestral hybridization. Although intra-genomic multiple variants are frequently found within each genome, the usage of the major variants alone is sufficient for phylogeny construction and species determination in most cases. Furthermore, the inclusion of minor variants further improves the resolution of species identification.

Molecular evolution of Drosophila Sex-lethal and related sex determining genes

Background

Sex determining mechanisms are evolutionarily labile and related species often use different primary signals and gene regulatory networks. This is well illustrated by the sex determining cascade of Drosophila fruitflies, which have recruited Sex-lethal as the master switch and cellular memory of sexual identity, a role performed in other insects by the gene transformer. Here we investigate the evolutionary change in the coding sequences of sex determining genes associated with the recruitment of Sex-lethal. We analyze sequences of Sex-lethal itself, its Drosophila paralogue sister-or-Sex-lethal and downstream targets transformer and doublesex.

Results

We find that the recruitment of sister-or-Sex-lethal was associated with a number of adaptive amino acid substitutions, followed by a tightening of purifying selection within the Drosophila clade. Sequences of the paralogue sister-or-Sex-lethal, in contrast, show a signature of rampant positive selection and relaxation of purifying selection. The recruitment of Sex-lethal as top regulator and memory gene is associated with a significant release from purifying selection in transformer throughout the Drosophila clade. In addition, doublesex shows a signature of positive selection and relaxation of purifying selection in the Drosophila clade. A similar pattern is seen in sequences from the sister Tephritidae clade.

Conclusions

The pattern of molecular evolution we observe for Sex-lethal and its paralogue sister-or-Sex-lethal is not characteristic of a duplication followed by neo-functionalization. Rather, evidence suggests a sub-functionalization scenario achieved through the evolution of sophisticated splicing. As expected, we find that transformer evolves under relaxed purifying selection after the recruitment of Sex-lethal in Drosophila. Finally, the observation of doublesex adaptation in both Drosophila and Tephritidae suggests that these changes are due to ongoing adaptation of downstream sex-specific regulation, rather than being associated the recruitment of Sex-lethal and the resulting change in the topology of the sex determining cascade.

Function and evolution of sex determination mechanisms, genes and pathways in insects

Animals have evolved a bewildering diversity of mechanisms to determine the two sexes. Studies of sex determination genes – their history and function – in non-model insects and Drosophila have allowed us to begin to understand the generation of sex determination diversity. One common theme from these studies is that evolved mechanisms produce activities in either males or females to control a shared gene switch that regulates sexual development. Only a few small-scale changes in existing and duplicated genes are sufficient to generate large differences in sex determination systems. This review summarises recent findings in insects, surveys evidence of how and why sex determination mechanisms can change rapidly and suggests fruitful areas of future research.