Recurrent gene amplification and soft selective sweeps during evolution of multidrug resistance in malaria parasites

When selection is strong and beneficial alleles have a single origin, local reductions in genetic diversity are expected. However, when beneficial alleles have multiple origins or were segregating in the population prior to a change in selection regime, the impact on genetic diversity may be less clear. We describe an example of such a "soft" selective sweep in the malaria parasite Plasmodium falciparum that involves adaptive genome rearrangements. Amplification in copy number of genome regions containing the pfmdr1 gene on chromosome 5 confer resistance to mefloquine and spread rapidly in the 1990s. Using flanking microsatellite data and real-time polymerase chain reaction determination of copy number, we show that 5-15 independent amplification events have occurred in parasites on the Thailand/Burma border. The amplified genome regions (amplicons) range in size from 14.7 to 49 kb and contain 2-11 genes, with 2-4 copies arranged in tandem. To examine the impact of drug selection on flanking variation, we genotyped 48 microsatellites on chromosome 5 in 326 parasites from a single Thai location. Diversity was reduced in a 170- to 250-kb (10-15 cM) region of chromosomes containing multiple copies of pfmdr1, consistent with hitchhiking resulting from the rapid recent spread of selected chromosomes. However, diversity immediately flanking pfmdr1 is reduced by only 42% on chromosomes bearing multiple amplicons relative to chromosomes carrying a single copy. We highlight 2 features of these results: 1) All amplicon break points occur in monomeric A/T tracts (9-45 bp). Given the abundance of these tracts in P. falciparum, we expect that duplications will occur frequently at multiple genomic locations and have been underestimated as drivers of phenotypic evolution in this pathogen. 2) The signature left by the spread of amplified genome segments is broad, but results in only limited reduction in diversity. If such "soft" sweeps are common in nature, statistical methods based on diversity reduction may be inefficient at detecting evidence for selection in genome-wide marker screens. This may be particularly likely when mutation rate is high, as appears to be the case for gene duplications, and in pathogen populations where effective population sizes are typically very large.     

Drosophila mojoless, a retroposed GSK-3, has functionally diverged to acquire an essential role in male fertility

Retroposition is increasingly recognized as an important mechanism for the acquisition of new genes. We show that a glycogen synthase kinase-3 gene, shaggy (sgg), retroposed at least 50 MYA in the Drosophila genus to generate a new gene, mojoless (mjl). We have extensively analyzed the function of mjl and examined its functional divergence from the parental gene sgg in Drosophila melanogaster. Unlike Sgg, which is expressed in many tissues of both sexes, Mjl is expressed specifically in the male germ line, where it is required for male germ line survival. Our analysis indicates that mjl has acquired a specific function in the maintenance of male germ line viability. However, it has not completely lost its ancestral biochemical function and can partially compensate for loss of the parental gene sgg when ectopically expressed in somatic cells. We postulate that mjl has undergone functional diversification and is now under stabilizing selection in the Drosophila genus.     

Not born equal: increased rate asymmetry in relocated and retrotransposed rodent gene duplicates

Duplicated genes frequently evolve at different rates. This asymmetry is evidence of natural selection's ability to discriminate between the 2 copies, subjecting them to different levels of purifying selection or even permitting adaptive evolution of one or both copies. However, if gene duplication creates pairs of protein-coding sequences that are initially identical, this raises the question of how selection tells the 2 copies apart. Here, we investigated asymmetric sequence divergence of recently duplicated genes in rodents and related this to 2 possible sources of such asymmetry: gene relocation as a consequence of duplication and retrotransposition as a mechanism of gene duplication. We found that most young rodent duplicates that have been relocated were created by retrotransposition. The degree of rate asymmetry in gene pairs where one copy has been relocated (either by retrotransposition or DNA-based duplication) is greater than in pairs formed by local DNA-based duplication events. Furthermore, by considering the direction of transposition for distant duplicates, we found a consistent tendency for retrogenes to undergo accelerated protein evolution relative to their static paralogs, whereas DNA-based transpositions showed no such tendency. Finally, we demonstrate that the faster sequence evolution of retrogenes correlates with the profound alteration of their expression pattern that is precipitated by retrotransposition.     

Anciently duplicated <Emphasis Type="Italic">Broad Complex</Emphasis> exons have distinct temporal functions during tissue morphogenesis

Broad Complex (BRC) is an essential ecdysone-pathway gene required for entry into and progression through metamorphosis in Drosophila melanogaster. Mutations of three BRC complementation groups cause numerous phenotypes, including a common suite of morphogenesis defects involving central nervous system (CNS), adult salivary glands (aSG), and male genitalia. These defects are phenocopied by the juvenile hormone mimic methoprene. Four BRC isoforms are produced by alternative splicing of a protein-binding BTB/POZ-encoding exon (BTB _BRC ) to one of four tandemly duplicated, DNA-binding zinc-finger-encoding exons (Z1 _BRC , Z2 _BRC , Z3 _BRC , Z4 _BRC ). Highly conserved orthologs of BTB _BRC and all four Z _BRC were found among published cDNA sequences or genome databases from Diptera, Lepidoptera, Hymenoptera, and Coleoptera, indicating that BRC arose and underwent internal exon duplication before the split of holometabolous orders. Tramtrack subfamily members, abrupt, tramtrack, fruitless, longitudinals lacking (lola), and CG31666 were characterized throughout Holometabola and used to root phylogenetic analyses of Z _BRC exons, which revealed that the Z _BRC clade includes Z _abrupt . All four Z _BRC domains, including Z4 _BRC , which has no known essential function, are evolving in a manner consistent with selective constraint. We used transgenic rescue to explore how different BRC isoforms contribute to shared tissue-morphogenesis functions. As predicted from earlier studies, the common CNS and aSG phenotypes were rescued by BRC-Z1 in rbp mutants, BRC-Z2 in br mutants, and BRC-Z3 in 2Bc mutants. However, the isoforms are required at two different developmental stages, with BRC-Z2 and -Z3 required earlier than BRC-Z1. The sequential action of BRC isoforms indicates subfunctionalization of duplicated Z _BRC exons even when they contribute to common developmental processes. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Speciation by genome duplication: Repeated origins and genomic composition of the recently formed allopolyploid species Mimulus peregrinus

Whole genome duplication (polyploidization) is a mechanism of “instantaneous” species formation that has played a major role in the evolutionary history of plants. Much of what we know about the early evolution of polyploids is based upon studies of a handful of recently formed species. A new polyploid hybrid (allopolyploid) species Mimulus peregrinus, formed within the last 140 years, was recently discovered on the Scottish mainland and corroborated by chromosome counts. Here, using targeted, high‐depth sequencing of 1200 genic regions, we confirm the parental origins of this new species from M. x robertsii, a sterile triploid hybrid between the two introduced species M. guttatus and M. luteus that are naturalized and widespread in the United Kingdom. We also report a new population of M. peregrinus on the Orkney Islands and demonstrate that populations on the Scottish mainland and Orkney Islands arose independently via genome duplication from local populations of M. x robertsii. Our data raise the possibility that some alleles are already being lost in the evolving M. peregrinus genomes. The recent origins of a new species of the ecological model genus Mimulus via allopolyploidization provide a powerful opportunity to explore the early stages of hybridization and genome duplication in naturally evolved lineages.     

Alternative Transposition Generates New Chimeric Genes and Segmental Duplications at the Maize p1 Locus

The maize Ac/Ds transposon family was the first transposable element system identified and characterized by Barbara McClintock. Ac/Ds transposons belong to the hAT family of class II DNA transposons. We and others have shown that Ac/Ds elements can undergo a process of alternative transposition in which the Ac/Ds transposase acts on the termini of two separate, nearby transposons. Because these termini are present in different elements, alternative transposition can generate a variety of genome alterations such as inversions, duplications, deletions, and translocations. Moreover, Ac/Ds elements transpose preferentially into genic regions, suggesting that structural changes arising from alternative transposition may potentially generate chimeric genes at the rearrangement breakpoints. Here we identified and characterized 11 independent cases of gene fusion induced by Ac alternative transposition. In each case, a functional chimeric gene was created by fusion of two linked, paralogous genes; moreover, each event was associated with duplication of the ∼70-kb segment located between the two paralogs. An extant gene in the maize B73 genome that contains an internal duplication apparently generated by an alternative transposition event was also identified. Our study demonstrates that alternative transposition-induced duplications may be a source for spontaneous creation of diverse genome structures and novel genes in maize.     

Molecular evolution and expression of the CRAL_TRIO protein family in insects

CRAL_TRIO domain proteins are known to bind small lipophilic molecules such as retinal, inositol and Vitamin E and include such gene family members as PINTA, α-tocopherol transfer (ATT) proteins, retinoid binding proteins, and clavesins. In insects, very little is known about either the molecular evolution of this family of proteins or their ligand specificity. Here we characterize insect CRAL_TRIO domain proteins and present the first insect CRAL_TRIO protein phylogeny constructed by performing reciprocal BLAST searches of the reference genomes of Drosophila melanogaster, Anopheles gambiae, Apis mellifera, Tribolium castaneum, Bombyx mori, Manduca sexta and Danaus plexippus. We find several highly conserved amino acid residues in the CRAL_TRIO domain-containing genes across insects and a gene expansion resulting in more than twice as many gene family members in lepidopterans than in other surveyed insect species, but no lepidopteran homolog of the PINTA gene in Drosophila. In addition, we examined the expression pattern of CRAL_TRIO domain genes in Manduca sexta heads using RNA-Seq data. Of the 42 gene family members found in the M. sexta reference genome, we found 30 expressed in the head tissue with similar expression profiles between males and females. Our results suggest this gene family underwent a large expansion in Lepidoptera, making the lepidopteran CRAL_TRIO domain family distinct from other holometabolous insect lineages.     

Species-specific duplications driving the recent expansion of NBS-LRR genes in five Rosaceae species

Background

Disease resistance (R) genes from different Rosaceae species have been identified by map-based cloning for resistance breeding. However, there are few reports describing the pattern of R-gene evolution in Rosaceae species because several Rosaceae genome sequences have only recently become available.

Results

Since most disease resistance genes encode NBS-LRR proteins, we performed a systematic genome-wide survey of NBS-LRR genes between five Rosaceae species, namely Fragaria vesca (strawberry), Malus × domestica (apple), Pyrus bretschneideri (pear), Prunus persica (peach) and Prunus mume (mei) which contained 144, 748, 469, 354 and 352 NBS-LRR genes, respectively. A high proportion of multi-genes and similar Ks peaks (Ks = 0.1- 0.2) of gene families in the four woody genomes were detected. A total of 385 species-specific duplicate clades were observed in the phylogenetic tree constructed using all 2067 NBS-LRR genes. High percentages of NBS-LRR genes derived from species-specific duplication were found among the five genomes (61.81% in strawberry, 66.04% in apple, 48.61% in pear, 37.01% in peach and 40.05% in mei). Furthermore, the Ks and Ka/Ks values of TIR-NBS-LRR genes (TNLs) were significantly greater than those of non-TIR-NBS-LRR genes (non-TNLs), and most of the NBS-LRRs had Ka/Ks ratios less than 1, suggesting that they were evolving under a subfunctionalization model driven by purifying selection.

Conclusions

Our results indicate that recent duplications played an important role in the evolution of NBS-LRR genes in the four woody perennial Rosaceae species. Based on the phylogenetic tree produced, it could be inferred that species-specific duplication has mainly contributed to the expansion of NBS-LRR genes in the five Rosaceae species. In addition, the Ks and Ka/Ks ratios suggest that the rapidly evolved TNLs have different evolutionary patterns to adapt to different pathogens compared with non-TNL resistant genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1291-0) contains supplementary material, which is available to authorized users.