Zika virus (ZIKV) is a mosquito-borne flavivirus first isolated in Uganda in 1947. Although entomological and virologic surveillance have reported ZIKV enzootic activity in diverse countries of Africa and Asia, few human cases were reported until 2007, when a Zika fever epidemic took place in Micronesia. In the context of West Africa, the WHO Collaborating Centre for Arboviruses and Hemorrhagic Fever at Institut Pasteur of Dakar (http://www.pasteur.fr/recherche/banques/CRORA/) reports the periodic circulation of ZIKV since 1968. Despite several reports on ZIKV, the genetic relationships among viral strains from West Africa remain poorly understood. To evaluate the viral spread and its molecular epidemiology, we investigated 37 ZIKV isolates collected from 1968 to 2002 in six localities in Senegal and Côte d'Ivoire. In addition, we included strains from six other countries. Our results suggested that these two countries in West Africa experienced at least two independent introductions of ZIKV during the 20th century, and that apparently these viral lineages were not restricted by mosquito vector species. Moreover, we present evidence that ZIKV has possibly undergone recombination in nature and that a loss of the N154 glycosylation site in the envelope protein was a possible adaptive response to the Aedes dalzieli vector.

Within the Haloferax genus, both the surface (S)-layer protein, and the glycans that can decorate it, vary between species, which can potentially result in many different surface types, analogous to bacterial serotypes. This variation may mediate phenotypes, such as sensitivity to different viruses and mating preferences. Here, we describe S-layer glycoproteins found in multiple Haloferax strains and perform comparative genomics analyses of major and alternative glycosylation clusters of isolates from two coastal sites. We analyze the phylogeny of individual glycosylation genes and demonstrate that while the major glycosylation cluster tends to be conserved among closely related strains, the alternative cluster is highly variable. Thus, geographically- and genetically-related strains may exhibit diverse surface structures to such an extent that no two isolates present an identical surface profile.

With the emergence of analytical software for the inference of viral evolution, a number of studies have focused on estimating important parameters such as the substitution rate and the time to the most recent common ancestor (t MRCA) for rapidly evolving viruses. Coupled with an increasing abundance of sequence data sampled under widely different schemes, an effort to keep results consistent and comparable is needed. This study emphasizes commonly disregarded problems in the inference of evolutionary rates in viral sequence data when sampling is unevenly distributed on a temporal scale through a study of the foot-and-mouth (FMD) disease virus serotypes SAT 1 and SAT 2. Our study shows that clustered temporal sampling in phylogenetic analyses of FMD viruses will strongly bias the inferences of substitution rates and t MRCA because the inferred rates in such data sets reflect a rate closer to the mutation rate rather than the substitution rate. Estimating evolutionary parameters from viral sequences should be performed with due consideration of the differences in short-term and longer-term evolutionary processes occurring within sets of temporally sampled viruses, and studies should carefully consider how samples are combined.

With the development of reliable recombination detection tools and an increasing number of available genome sequences, many studies have reported evidence of recombination in a wide range of virus genera. Recombination is apparently a major mechanism in virus evolution, allowing viruses to evolve more quickly by providing immediate direct access to many more areas of a sequence space than are accessible by mutation alone. Recombination has been widely described amongst the insect-transmitted plant viruses in the genus Begomovirus (family Geminiviridae), with potential recombination hot- and cold-spots also having been identified. Nevertheless, because very little is understood about either the biochemical predispositions of different genomic regions to recombine or what makes some recombinants more viable than others, the sources of the evolutionary and biochemical forces shaping distinctive recombination patterns observed in nature remain obscure. Here we present a detailed analysis of unique recombination events detectable in the DNA-A and DNA-A-like genome components of bipartite and monopartite begomoviruses. We demonstrate both that recombination breakpoint hot- and cold-spots are conserved between the two groups of viruses, and that patterns of sequence exchange amongst the genomes are obviously non-random. Using a computational technique designed to predict structural perturbations in chimaeric proteins, we demonstrate that observed recombination events tend to be less disruptive than sets of simulated ones. Purifying selection acting against natural recombinants expressing improperly folded chimaeric proteins is therefore a major determinant of natural recombination patterns in begomoviruses.

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

We have studied the evolution of a type IV secretion system (T4SS), in Bartonella, which is thought to have changed function from conjugation to erythrocyte adherence following a recent horizontal gene transfer event. The system, called Trw, is unique among T4SSs in that genes encoding both exo- and intracellular components are located within the same duplicated fragment. This provides an opportunity to study the influence of selection on proteins involved in host-pathogen interactions. We sequenced the trw locus from several strains of Bartonella henselae and investigated its evolutionary history by comparisons to other Bartonella species. Several instances of recombination and gene conversion events where detected in the 2- to 5-fold duplicated gene fragments encompassing trwJIH, explaining the homogenization of the anchoring protein TrwI and the divergence of the minor pilus protein TrwJ. A phylogenetic analysis of the 7- to 8-fold duplicated gene coding for the major pilus protein TrwL displayed 2 distinct clades, likely representing a subfunctionalization event. The analyses of the B. henselae strains also identified a recent horizontal transfer event of almost the complete trwL region. We suggest that the switch in function of the T4SS was mediated by the duplication of the genes encoding pilus components and their diversification by combinatorial sequence shuffling within and among genomes. We suggest that the pilus proteins have evolved by diversifying selection to match a divergent set of erythrocyte surface structures, consistent with the trench warfare coevolutionary model.

Viruses are the most abundant and the smallest organisms, which are relatively simple to sequence. Genome sequence data of viruses for individual species to populations out‐ number that of other species. Although this offers an opportunity to study viral diversity at varying levels of taxonomic hierarchy, it also poses challenges for systematic and struc‐ tured organization of data and its downstream processing. Extensive computational anal‐ yses using a number of algorithms and programs have opened exciting opportunities for virus discovery and diagnostics, apart from augmenting our understanding of the intri‐ guing world of viruses. Unravelling evolutionary dynamics of viruses permits improved understanding of phenomena such as quasispecies diversity, role of mutations in host switching and drug resistance, which enables the tangible measurements of genotype and phenotype of viruses. Improved understanding of geno-/serotype diversity in corre‐ lation with antigenic diversity will facilitate rational design and development of effica‐ cious vaccines against emerging and re-emerging viruses. Mathematical models developed using the genomic data could be used to predict the spread of viruses due to vector switching and the (re)emergence due to host switching and, thereby, contribute to‐ wards designing public health policies for disease management and control.

Viral particles in stool samples from wild-living chimpanzees were analysed using random PCR amplification and sequencing. Sequences encoding proteins distantly related to the replicase protein of single-stranded circular DNA viruses were identified. Inverse PCR was used to amplify and sequence multiple small circular DNA viral genomes. The viral genomes were related in size and genome organization to vertebrate circoviruses and plant geminiviruses but with a different location for the stem-loop structure involved in rolling circle DNA replication. The replicase genes of these viruses were most closely related to those of the much smaller (approximately 1 kb) plant nanovirus circular DNA chromosomes. Because the viruses have characteristics of both animal and plant viruses, we named them chimpanzee stool-associated circular viruses (ChiSCV). Further metagenomic studies of animal samples will greatly increase our knowledge of viral diversity and evolution.

Viking smallpox diversity Humans have a notable capacity to withstand the ravages of infectious diseases. Smallpox killed millions of people but drove Jenner's invention of vaccination, which eventually led to the annihilation of this virus, declared in 1980. To investigate the history of smallpox, Mühlemann et al. obtained high-throughput shotgun sequencing data from 1867 human remains ranging from >31,000 to 150 years ago (see the Perspective by Alcamí). Thirteen positive samples emerged, 11 of which were northern European Viking Age people (6th to 7th century CE). Although the sequences were patchy and incomplete, four could be used to infer a phylogenetic tree. This showed distinct Viking Age lineages with multiple gene inactivations. The analysis pushes back the date of the earliest variola infection in humans by ∼1000 years and reveals the existence of a previously unknown virus clade. Science, this issue p. eaaw8977; see also p. 376 Variola virus DNA sequences isolated from ancient human remains reveal the presence of smallpox in 7th-century northern Europe. INTRODUCTION Variola virus (VARV), the causative agent of smallpox, is estimated to have killed between 300 million and 500 million people in the 20th century and was responsible for widespread mortality and suffering for at least several preceding centuries. Humans are the only known host of VARV, and smallpox was declared eradicated in 1980. The timeline of the emergence of smallpox in humans is unclear. Based on sequence data up to 360 years old, the most recent common ancestor of VARV has been dated to the 16th or 17th century. This contrasts with written records of possible smallpox infections dating back at least 3000 years and mummified remains suggestive of smallpox dating to 3570 years ago. RATIONALE Ancient virus sequences recovered from archaeological remains provide direct molecular evidence of past infections, give detail of genetic changes that have occurred during the evolution of the virus, and can reveal viable virus sequence diversity not currently present in modern viruses. In the case of VARV, ancient sequences may also reduce the gap between the written historical record of possible early smallpox infections and the dating of the oldest available VARV sequences. We therefore screened high-throughput shotgun sequencing data from skeletal and dental remains of 1867 humans living in Eurasia and the Americas between ~31,630 and ~150 years ago for the presence of sequences matching VARV. RESULTS VARV sequences were recovered from 13 northern European individuals, including 11 dated to ~600–1050 CE, overlapping the Viking Age, and we reconstructed near-complete VARV genomes for four of them. The samples predate the earliest confirmed smallpox cases by ~1000 years. Eleven of the recovered sequences fall into a now-extinct sister clade of the modern VARVs in circulation prior to the eradication of smallpox, while two sequences from the 19th century group with modern VARV. The inferred date of the most recent common ancestor of VARV is ~1700 years ago. The number of functional genes is generally reduced in orthopoxviruses with narrow host ranges. A comparison of the gene content of the Viking Age sequences shows great contrast with that of modern VARV. Three genes that are active in all modern VARV sequences were inactive over 1000 years ago in some or all ancient VARV. Among 10 genes inactive in modern and Viking Age VARV, the mutations causing the inactivations are different and the genes are predicted to be active in the ancestor of both clades, suggesting parallel evolution. Fourteen genes inactivated in modern VARV are active in some or all of the ancient sequences, eight of which encode known virulence factors or immunomodulators. The active gene counts of the four higher-coverage Viking Age viral genomes provide snapshots from an ~350-year period, showing the reduction of gene content during the evolution of VARV. These genomes support suggestions that orthopoxvirus species derive from a common ancestor containing all genes present in orthopoxviruses today, with the reduction in active gene count conjectured to be the result of long-term adaptation within host species. CONCLUSION The Viking Age sequences reported here push the definitive date of the earliest VARV infection in humans back by ~1000 years. These sequences, combined with early written records of VARV epidemics in southern and western Europe, suggest a pan-European presence of smallpox from the late 6th century. The ancient viruses are part of a previously unknown, now-extinct virus clade and were following a genotypic evolutionary path that differs from modern VARV. The reduction in gene content shows that multiple combinations of active genes have led to variola viruses capable of circulating widely within the human population. Sample locations and genome evolution of Viking Age variola virus. (Top left) Sample VK533 with dagger, in situ. (Bottom left) Positive sample locations. Higher- and lower-coverage samples have solid and open circles, respectively. (Right) Diverse gene inactivation patterns in the four higher-coverage 600–1000 CE Viking Age sequences, within clade and as compared to sequences from other human VARVs and the phylogenetically nearest animal poxviruses, camelpox and taterapox. Open circles show intact genes, circle colors (independent for each gene) indicate different inactivations. Gene names at the bottom. Smallpox, one of the most devastating human diseases, killed between 300 million and 500 million people in the 20th century alone. We recovered viral sequences from 13 northern European individuals, including 11 dated to ~600–1050 CE, overlapping the Viking Age, and reconstructed near-complete variola virus genomes for four of them. The samples predate the earliest confirmed smallpox cases by ~1000 years, and the sequences reveal a now-extinct sister clade of the modern variola viruses that were in circulation before the eradication of smallpox. We date the most recent common ancestor of variola virus to ~1700 years ago. Distinct patterns of gene inactivation in the four near-complete sequences show that different evolutionary paths of genotypic host adaptation resulted in variola viruses that circulated widely among humans.

UNLABELLED Tomato yellow leaf curl disease (TYLCD) is one of the most devastating viral diseases affecting tomato crops in tropical, subtropical and temperate regions of the world. Here, we focus on the interactions through recombination between the different begomovirus species causing TYLCD, provide an overview of the interactions with the cellular genes involved in viral replication, and highlight recent progress on the relationships between these viruses and their vector, the whitefly Bemisia tabaci. TAXONOMY The tomato yellow leaf curl virus-like viruses (TYLCVs) are a complex of begomoviruses (family Geminiviridae, genus Begomovirus) including 10 accepted species: Tomato yellow leaf curl Axarquia virus (TYLCAxV), Tomato yellow leaf curl China virus (TYLCCNV), Tomato yellow leaf curl Guangdong virus (TYLCGuV), Tomato yellow leaf curl Indonesia virus (TYLCIDV), Tomato yellow leaf curl Kanchanaburi virus (TYLVKaV), Tomato yellow leaf curl Malaga virus (TYLCMalV), Tomato yellow leaf curl Mali virus (TYLCMLV), Tomato yellow leaf curl Sardinia virus (TYLCSV), Tomato yellow leaf curl Thailand virus (TYLCTHV), Tomato yellow leaf curl Vietnam virus (TYLCVNV) and Tomato yellow leaf curl virus(TYLCV). We follow the species demarcation criteria of the International Committee on Taxonomy of Viruses (ICTV), the most important of which is an 89% nucleotide identity threshold between full-length DNA-A component nucleotide sequences for begomovirus species. Strains of a species are defined by a 93% nucleotide identity threshold. HOST RANGE The primary host of TYLCVs is tomato (Solanum lycopersicum), but they can also naturally infect other crops [common bean (Phaseolus vulgaris), sweet pepper (Capsicum annuum), chilli pepper (C. chinense) and tobacco (Nicotiana tabacum)], a number of ornamentals [petunia (Petuniaxhybrida) and lisianthus (Eustoma grandiflora)], as well as common weeds (Solanum nigrum and Datura stramonium). TYLCVs also infect the experimental host Nicotiana benthamiana. DISEASE SYMPTOMS Infected tomato plants are stunted or dwarfed, with leaflets rolled upwards and inwards; young leaves are slightly chlorotic; in recently infected plants, fruits might not be produced or, if produced, are small and unmarketable. In common bean, some TYLCVs produce the bean leaf crumple disease, with thickening, epinasty, crumpling, blade reduction and upward curling of leaves, as well as abnormal shoot proliferation and internode reduction; the very small leaves result in a bushy appearance.

UNLABELLED RDP2 is a Windows 95/XP program that examines nucleotide sequence alignments and attempts to identify recombinant sequences and recombination breakpoints using 10 published recombination detection methods, including GENECONV, BOOTSCAN, MAXIMUM chi(2), CHIMAERA and SISTER SCANNING. The program enables fast automated analysis of large alignments (up to 300 sequences containing 13 000 sites), and interactive exploration, management and verification of results with different recombination detection and tree drawing methods. AVAILABILITY RDP2 is available free from the RDP2 website (http://darwin.uvigo.es/rdp/rdp.html) CONTACT darren@science.uct.ac.za SUPPLEMENTARY INFORMATION Detailed descriptions of RDP2 and the methods it implements are included in the program manual, which can be downloaded from the RDP2 website.

UNLABELLED Maize streak virus (MSV; Genus Mastrevirus, Family Geminiviridae) occurs throughout Africa, where it causes what is probably the most serious viral crop disease on the continent. It is obligately transmitted by as many as six leafhopper species in the Genus Cicadulina, but mainly by C. mbila Naudé and C. storeyi. In addition to maize, it can infect over 80 other species in the Family Poaceae. Whereas 11 strains of MSV are currently known, only the MSV-A strain is known to cause economically significant streak disease in maize. Severe maize streak disease (MSD) manifests as pronounced, continuous parallel chlorotic streaks on leaves, with severe stunting of the affected plant and, usuallly, a failure to produce complete cobs or seed. Natural resistance to MSV in maize, and/or maize infections caused by non-maize-adapted MSV strains, can result in narrow, interrupted streaks and no obvious yield losses. MSV epidemiology is primarily governed by environmental influences on its vector species, resulting in erratic epidemics every 3-10 years. Even in epidemic years, disease incidences can vary from a few infected plants per field, with little associated yield loss, to 100% infection rates and complete yield loss. TAXONOMY The only virus species known to cause MSD is MSV, the type member of the Genus Mastrevirus in the Family Geminiviridae. In addition to the MSV-A strain, which causes the most severe form of streak disease in maize, 10 other MSV strains (MSV-B to MSV-K) are known to infect barley, wheat, oats, rye, sugarcane, millet and many wild, mostly annual, grass species. Seven other mastrevirus species, many with host and geographical ranges partially overlapping those of MSV, appear to infect primarily perennial grasses. PHYSICAL PROPERTIES MSV and all related grass mastreviruses have single-component, circular, single-stranded DNA genomes of approximately 2700 bases, encapsidated in 22 x 38-nm geminate particles comprising two incomplete T = 1 icosahedra, with 22 pentameric capsomers composed of a single 32-kDa capsid protein. Particles are generally stable in buffers of pH 4-8. DISEASE SYMPTOMS In infected maize plants, streak disease initially manifests as minute, pale, circular spots on the lowest exposed portion of the youngest leaves. The only leaves that develop symptoms are those formed after infection, with older leaves remaining healthy. As the disease progresses, newer leaves emerge containing streaks up to several millimetres in length along the leaf veins, with primary veins being less affected than secondary or tertiary veins. The streaks are often fused laterally, appearing as narrow, broken, chlorotic stripes, which may extend over the entire length of severely affected leaves. Lesion colour generally varies from white to yellow, with some virus strains causing red pigmentation on maize leaves and abnormal shoot and flower bunching in grasses. Reduced photosynthesis and increased respiration usually lead to a reduction in leaf length and plant height; thus, maize plants infected at an early stage become severely stunted, producing undersized, misshapen cobs or giving no yield at all. Yield loss in susceptible maize is directly related to the time of infection: infected seedlings produce no yield or are killed, whereas plants infected at later times are proportionately less affected. DISEASE CONTROL Disease avoidance can be practised by only planting maize during the early season when viral inoculum loads are lowest. Leafhopper vectors can also be controlled with insecticides such as carbofuran. However, the development and use of streak-resistant cultivars is probably the most effective and economically viable means of preventing streak epidemics. Naturally occurring tolerance to MSV (meaning that, although plants become systemically infected, they do not suffer serious yield losses) has been found, which has primarily been attributed to a single gene, msv-1. However, other MSV resistance genes also exist and improved resistance has been achieved by concentrating these within individual maize genotypes. Whereas true MSV immunity (meaning that plants cannot be symptomatically infected by the virus) has been achieved in lines that include multiple small-effect resistance genes together with msv-1, it has proven difficult to transfer this immunity into commercial maize genotypes. An alternative resistance strategy using genetic engineering is currently being investigated in South Africa. USEFUL WEBSITES http://www.mcb.uct.ac.za/MSV/mastrevirus.htm; http://www.danforthcenter.org/iltab/geminiviridae/geminiaccess/mastrevirus/Mastrevirus.htm.

Tuberculosis (TB) and leprosy are ancient human diseases, as demonstrated by an abundant osteoarchaeological record and by molecular analyses of prehistoric and historic human remains affected by the disease. Unfortunately, mycobacterial infections (which include TB and leprosy) are not diseases of the past. Today, mycobacteria are still an enormous burden to humans, with an estimated two billion individuals currently being affected by the disease (in a latent or active stage) worldwide. It has become clear that these successful human pathogens can rapidly adapt to their changing environment and host. Evolutionary theory provides an excellent framework to test hypotheses that aim to explain how these pathogens are able to sustain such changes. Understanding the evolutionary mechanisms of the genus Mycobacterium will improve our understanding, and perhaps our ability to control, the processes that yield to these pathogens' evolution and acquisition of drug resistance. Among these adaptive mechanisms, recombination is without doubt an effective way to acquire foreign genetic material into the mycobacterial genome. The scientific community has not reached a consensus on the occurrence of recombination in the genus Mycobacterium. This dissertation aimed to screen a data set comprising 2354 homologous genome sections from publicly available mycobacterial sequences representing 18 strains from 13 species for the occurrence of recombination. Homologous genomic sets were aligned and concatenated in the chromosomal order of one mycobacterial reference species. The concatenated sequence was screened for recombination using a variety of statistical methods, with each event being evaluated by comparing maximum likelihood phylogenies of each recombinant section with the non-recombining portion of the dataset. Incongruent phylogenies were identified by comparing the site-wise log-likelihoods of each tree using multiple tests. The analyses described here indicate that the hypothesis of recombination in mycobacterial homologues cannot be rejected. While an analysis of non-homologous mycobacterial regions will be required to assess the occurrence of nonhomologous recombination, this study shows that we cannot assume mycobacterial clonality, even in homologous genomic sections.

The recent emergence of Zika virus (ZIKV) has been concentrated in the Caribbean, Southeastern United States, and South- and Central America; resulting in travel-based cases being reported around the globe. As multi-disciplinary collaborations are combatting the ZIKV outbreak, the need to validate the sequence of existing strains has become apparent. Here, we report high-quality sequence data for multiple ZIKV strains made publicly available through the National Institutes of Health- (NIH) funded biorepository, BEI Resources (www.beiresources.org). Next-generation sequencing, 3′ rapid amplification of cDNA ends (RACE), and viral genome annotation pipelines generated GenBank sequence records for 16 BEI Resources strains. Minor variants, consensus mutations, and consensus insertions/deletions were identified within the viral stocks using next-generation sequencing (NGS) and consensus changes were confirmed with Sanger sequencing. Bioinformatics analyses of the sequencing results confirm that the virus stocks available to the scientific research community through BEI Resources adequately represent the viral population diversity of ZIKV.

The genotyping of highly polymorphic multigene families across many individuals used to be a particularly challenging task because of methodological limitations associated with traditional approaches. Next‐generation sequencing (NGS) can overcome most of these limitations, and it is increasingly being applied in population genetic studies of multigene families. Here, we critically review NGS bioinformatic approaches that have been used to genotype the major histocompatibility complex (MHC) immune genes, and we discuss how the significant advances made in this field are applicable to population genetic studies of gene families. Increasingly, approaches are introduced that apply thresholds of sequencing depth and sequence similarity to separate alleles from methodological artefacts. We explain why these approaches are particularly sensitive to methodological biases by violating fundamental genotyping assumptions. An alternative strategy that utilizes ultra‐deep sequencing (hundreds to thousands of sequences per amplicon) to reconstruct genotypes and applies statistical methods on the sequencing depth to separate alleles from artefacts appears to be more robust. Importantly, the ‘degree of change’ (DOC) method avoids using arbitrary cut‐off thresholds by looking for statistical boundaries between the sequencing depth for alleles and artefacts, and hence, it is entirely repeatable across studies. Although the advances made in generating NGS data are still far ahead of our ability to perform reliable processing, analysis and interpretation, the community is developing statistically rigorous protocols that will allow us to address novel questions in evolution, ecology and genetics of multigene families. Future developments in third‐generation single molecule sequencing may potentially help overcome problems that still persist in de novo multigene amplicon genotyping when using current second‐generation sequencing approaches.

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.

The evolution of the mitochondrial (mt) genome is far from being fully understood. Systematic investigations into the modes of inheritance, rates and patterns of recombination, nucleotide substitution, and structural changes in the mt genome are still lacking in many groups of plants. In this study, we sequenced >11kbp mtDNA segments from multiple accessions of 36 pine species to characterize the evolutionary patterns of mtDNA in the genus Pinus. We found extremely low substitution rates and complex repetitive sequences scattered across different genome regions, as well as chimeric structures that were probably generated by multiple intergenomic recombinations. The mtDNA-based phylogeny of the genus differed from that based on chloroplast and nuclear DNA in the placement of several groups of species. Such discordances suggest a series of mtDNA capture events during past range shifts of the pine species and that both vertical and horizontal inheritance are implicated in the evolution of mtDNA in Pinus. MtDNA dating revealed that most extant lineages of the genus originated during Oligocene-Miocene radiation and subgenus Strobus diversified earlier than subgenus Pinus. Our findings illustrate a reticular evolutionary pathway for the mt genome through capture and recombination in the genus Pinus, and provide new insights into the evolution of the genus.

The birth-and-death model of multigene family evolution describes how families can expand by duplication and contract by gene deletion and formation of pseudogenes. The phylogenetic stability of a gene is thought to be related to the degree of functional importance. However, it is unclear how much evolution of a gene in a gene family is driven by adaptive versus neutral processes. The cytochrome P450s are one of the most diverse and well-studied multigene families, involved in both physiological and xenobiotic functions. Bees have a high toxin exposure due to their diet of nectar and pollen, as well as the resin-collecting behavior exhibited by some bees. Here, we describe the P450s of the orchid bee Euglossa dilemma. Orchid bees are a neotropical clade in which males form perfume bouquets used in courtship displays by collecting a diverse set of volatile compounds, resulting in high chemical compound exposure. We conducted phylogenetic and selection analyses across ten bee species encompassing three bee families. We do not find a relationship between the ecology of a bee species and its P450 repertoire. Our analyses reveal that P450 clades can be classified into stable and unstable clades, and that genes involved in xenobiotic metabolism are more likely to belong to unstable clades. Furthermore, we find that unstable clades are under more dynamic evolutionary pressures, with signals of adaptive evolution detected, suggesting that both gene duplication and positive selection driving sequence divergence have played a role in the diversification of bee P450s. Our works highlights the complexity of multigene family evolution which does not always follow generalized predictions. Significance statement Gene family evolution is characterized by deletion and duplication, but it is unclear if this is driven by adaptive or neutral evolutionary processes. To investigate the dynamics of gene family evolution we analyze the P450s of ten bee species, a well-studied gene family involved in both physiological and detoxification roles. We find that genes involved in detoxification are more likely to belong to unstable clades, exhibit dynamic evolutionary pressures, and show signals of adaptive evolution. However, we do not find evidence for a relationship between resin collection or sociality, and the P450 repertoire of a bee species. Our findings highlight the complexity of multigene family evolution which does not always follow generalized predictions.

The annual sunflowers (Helianthus sect. Helianthus) present a formidable challenge for phylogenetic inference because of ancient hybrid speciation, recent introgression, and suspected issues with deep coalescence. Here we analyze sequence data from 11 nuclear DNA (nDNA) genes for multiple genotypes of species within the section to (1) reconstruct the phylogeny of this group, (2) explore the utility of nDNA gene trees for detecting hybrid speciation and introgression; and (3) test an empirical method of hybrid identification based on the phylogenetic congruence of nDNA gene trees from tightly linked genes. We uncovered considerable topological heterogeneity among gene trees with or without three previously identified hybrid species included in the analyses, as well as a general lack of reciprocal monophyly of species. Nonetheless, partitioned Bayesian analyses provided strong support for the reciprocal monophyly of all species except H. annuus (0.89 PP), the most widespread and abundant annual sunflower. Previous hypotheses of relationships among taxa were generally strongly supported (1.0 PP), except among taxa typically associated with H. annuus, apparently due to the paraphyly of the latter in all gene trees. While the individual nDNA gene trees provided a useful means for detecting recent hybridization, identification of ancient hybridization was problematic for all ancient hybrid species, even when linkage was considered. We discuss biological factors that affect the efficacy of phylogenetic methods for hybrid identification.

The ability of pathogens to escape the host's immune response is crucial for the establishment of persistent infections and can influence virulence. Recombination has been observed to contribute to this process by generating novel genetic variants. Although distinctive recombination patterns have been described in many viral pathogens, little is known about the influence of biases in the recombination process itself relative to selective forces acting on newly formed recombinants. Understanding these influences is important for determining how recombination contributes to pathogen genome and proteome evolution. Most previous research on recombination-driven protein evolution has focused on relatively simple proteins, usually in the context of directed evolution experiments. Here, we study recombination in the envelope gene of HIV-1 between primary isolates belonging to subtypes that recombine naturally in the HIV/AIDS pandemic. By characterizing the early steps in the generation of recombinants, we provide novel insights into the evolutionary forces that shape recombination patterns within viral populations. Specifically, we show that the combined effects of mechanistic processes that determine the locations of recombination breakpoints across the HIV-1 envelope gene, and purifying selection acting against dysfunctional recombinants, can explain almost the entire distribution of breakpoints found within this gene in nature. These constraints account for the surprising paucity of recombination breakpoints found in infected individuals within this highly variable gene. Thus, the apparent randomness of HIV evolution via recombination may in fact be relatively more predictable than anticipated. In addition, the dominance of purifying selection in localized areas of the HIV genome defines regions where functional constraints on recombinants appear particularly strong, pointing to vulnerable aspects of HIV biology.