Adaptation of African swine fever virus to MA-104 cells: Implications of unique genetic variations

African swine fever virus (ASFV) is a large, double-stranded DNA virus that causes a fatal, contagious disease specifically in pigs. However, prevention and control of ASFV outbreaks have been hampered by the lack of an effective vaccine or antiviral treatment for ASFV. Although ASFV has been reported to adapt to a variety of continuous cell lines, the phenotypic and genetic changes associated with ASFV adaptation to MA-104 cells remain poorly understood. Here, we adapted ASFV field isolates to efficiently propagate through serial viral passages in MA-104 cells. The adapted ASFV strain developed a pronounced cytopathic effect and robust infection in MA-104 cells. Interestingly, the adapted variant maintained its tropism in primary porcine kidney macrophages. Whole genome analysis of the adapted virus revealed unique gene deletions in the left and right variable regions of the viral genome compared to other previously reported cell culture-adapted ASFV strains. Notably, gene duplications at the 5' and 3' ends of the viral genome were in reverse complementary alignment with their paralogs. Single point mutations in protein-coding genes and intergenic regions were also observed in the viral genome. Collectively, our results shed light on the significance of these unique genetic changes during adaptation, which facilitate the growth of ASFV in MA-104 cells.

The evolutionary history and spectral tuning of vertebrate visual opsins

Many animals depend on the sense of vision for survival. In eumetazoans, vision requires specialized, light-sensitive cells called photoreceptors. Light reaches the photoreceptors and triggers the excitation of light-detecting proteins called opsins. Here, we describe the story of visual opsin evolution from the ancestral bilaterian to the extant vertebrate lineages. We explain the mechanisms determining color vision of extant vertebrates, focusing on opsin gene losses, duplications, and the expression regulation of vertebrate opsins. We describe the sequence variation both within and between species that has tweaked the sensitivities of opsin proteins towards different wavelengths of light. We provide an extensive resource of wavelength sensitivities and mutations that have diverged light sensitivity in many vertebrate species and predict how these mutations were accumulated in each lineage based on parsimony. We suggest possible natural and sexual selection mechanisms underlying these spectral differences. Understanding how molecular changes allow for functional adaptation of animals to different environments is a major goal in the field, and therefore identifying mutations affecting vision and their relationship to photic selection pressures is imperative. The goal of this review is to provide a comprehensive overview of our current understanding of opsin evolution in vertebrates.

Expansion of viral genomes with viral protein genome linked copies

Virus protein-linked genome (VPg) proteins are required for replication. VPgs are duplicated in a subset of RNA viruses however their roles are not fully understood and the extent of viral genomes containing VPg copies has not been investigated in detail. Here, we generated a novel bioinformatics approach to identify VPg sequences in viral genomes using hidden Markov models (HMM) based on alignments of dicistrovirus VPg sequences. From metagenomic datasets of dicistrovirus genomes, we identified 717 dicistrovirus genomes containing VPgs ranging from a single copy to 8 tandem copies. The VPgs are classified into nine distinct types based on their sequence and length. The VPg types but not VPg numbers per viral genome followed specific virus clades, thus suggesting VPgs co-evolved with viral genomes. We also identified VPg duplications in aquamavirus and mosavirus genomes. This study greatly expands the number of viral genomes that contain VPg copies and indicates that duplicated viral sequences are more widespread than anticipated.

Update of the keratin gene family: evolution, tissue-specific expression patterns, and relevance to clinical disorders

Intermediate filament (IntFil) genes arose during early metazoan evolution, to provide mechanical support for plasma membranes contacting/interacting with other cells and the extracellular matrix. Keratin genes comprise the largest subset of IntFil genes. Whereas the first keratin gene appeared in sponge, and three genes in arthropods, more rapid increases in keratin genes occurred in lungfish and amphibian genomes, concomitant with land animal-sea animal divergence (~ 440 to 410 million years ago). Human, mouse and zebrafish genomes contain 18, 17 and 24 non-keratin IntFil genes, respectively. Human has 27 of 28 type I "acidic" keratin genes clustered at chromosome (Chr) 17q21.2, and all 26 type II "basic" keratin genes clustered at Chr 12q13.13. Mouse has 27 of 28 type I keratin genes clustered on Chr 11, and all 26 type II clustered on Chr 15. Zebrafish has 18 type I keratin genes scattered on five chromosomes, and 3 type II keratin genes on two chromosomes. Types I and II keratin clusters-reflecting evolutionary blooms of keratin genes along one chromosomal segment-are found in all land animal genomes examined, but not fishes; such rapid gene expansions likely reflect sudden requirements for many novel paralogous proteins having divergent functions to enhance species survival following sea-to-land transition. Using data from the Genotype-Tissue Expression (GTEx) project, tissue-specific keratin expression throughout the human body was reconstructed. Clustering of gene expression patterns revealed similarities in tissue-specific expression patterns for previously described "keratin pairs" (i.e., KRT1/KRT10, KRT8/KRT18, KRT5/KRT14, KRT6/KRT16 and KRT6/KRT17 proteins). The ClinVar database currently lists 26 human disease-causing variants within the various domains of keratin proteins.

ABA signaling pathway genes and function during abiotic stress and berry ripening in Vitis vinifera

Abscisic acid (ABA) plays important roles in plant development and tolerance to abiotic stresses. Limited information is available regarding ABA signaling pathway genes in grape. In this study, 9 VvPYR/PYLs, 85 VvPP2Cs, 7 VvABIs, 7 VvSnRK2s, and 8 VvABFs were identified in the grape genome. Duplication analysis indicated that whole genome duplication might contribute to the expansion of these gene families. The comprehensive transcriptome analysis in various organs/tissues implied that most of these genes were tissue-specific, and few were environment-specific genes. Exogenous ABA treatment reduced the grape maturation period. VvPP2C59, VvPP2C60, VvPP2C66, and VvABF8 were all involved in tolerance to cold, heat, and drought stresses, revealing their crucial roles in regulating environmental stress responses. This work provides detailed information of ABA signaling pathway genes and new insights regarding their expression patterns during grape development and abiotic stress treatment.

The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest

BACKGROUND

Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture.

RESULTS

Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world.

CONCLUSIONS

The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

PISTILLATA paralogs in Tarenaya hassleriana have diverged in interaction specificity

BACKGROUND

Floral organs are specified by MADS-domain transcription factors that act in a combinatorial manner, as summarized in the (A)BCE model. However, this evolutionarily conserved model is in contrast to a remarkable amount of morphological diversity in flowers. One of the mechanisms suggested to contribute to this diversity is duplication of floral MADS-domain transcription factors. Although gene duplication is often followed by loss of one of the copies, sometimes both copies are retained. If both copies are retained they will initially be redundant, providing freedom for one of the paralogs to change function. Here, we examine the evolutionary fate and functional consequences of a transposition event at the base of the Brassicales that resulted in the duplication of the floral regulator PISTILLATA (PI), using Tarenaya hassleriana (Cleomaceae) as a model system.

RESULTS

The transposition of a genomic region containing a PI gene led to two paralogs which are located at different positions in the genome. The original PI copy is syntenic in position with most angiosperms, whereas the transposed copy is syntenic with the PI genes in Brassicaceae. The two PI paralogs of T. hassleriana have very similar expression patterns. However, they may have diverged in function, as only one of these PI proteins was able to act heterologously in the first whorl of A. thaliana flowers. We also observed differences in protein complex formation between the two paralogs, and the two paralogs exhibit subtle differences in DNA-binding specificity. Sequence analysis indicates that most of the protein sequence divergence between the two T. hassleriana paralogs emerged in a common ancestor of the Cleomaceae and the Brassicaceae.

CONCLUSIONS

We found that the PI paralogs in T. hassleriana have similar expression patterns, but may have diverged at the level of protein function. Data suggest that most protein sequence divergence occurred rapidly, prior to the origin of the Brassicaceae and Cleomaceae. It is tempting to speculate that the interaction specificities of the Brassicaceae-specific PI proteins are different compared to the PI found in other angiosperms. This could lead to PI regulating partly different genes in the Brassicaceae, and ultimately might result in change floral in morphology.

Correlated duplications and losses in the evolution of palmitoylation writer and eraser families

Background

Protein post-translational modifications (PTMs) change protein properties. Each PTM type is associated with domain families that apply the modification (writers), remove the modification (erasers) and bind to the modified sites (readers) together called toolkit domains. The evolutionary origin and diversification remains largely understudied, except for tyrosine phosphorylation. Protein palmitoylation entails the addition of a palmitoyl fatty acid to a cysteine residue. This PTM functions as a membrane anchor and is involved in a range of cellular processes. One writer family and two erasers families are known for protein palmitoylation.

Results

In this work we unravel the evolutionary history of these writer and eraser families. We constructed a high-quality profile hidden Markov model (HMM) of each family, searched for protein family members in fully sequenced genomes and subsequently constructed phylogenetic distributions of the families. We constructed Maximum Likelihood phylogenetic trees and using gene tree rearrangement and tree reconciliation inferred their evolutionary histories in terms of duplication and loss events. We identified lineages where the families expanded or contracted and found that the evolutionary histories of the families are correlated. The results show that the erasers were invented first, before the origin of the eukaryotes. The writers first arose in the eukaryotic ancestor. The writers and erasers show co-expansions in several eukaryotic ancestral lineages. These expansions often seem to be followed by contractions in some or all of the lineages further in evolution.

Conclusions

A general pattern of correlated evolution appears between writer and eraser domains. These co-evolution patterns could be used in new methods for interaction prediction based on phylogenies.

Electronic supplementary material

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

An insight into the evolutionary history of human MHC paralogon

The vertebrate genome contains several closely spaced sets of paralogous genes from distinct gene families on typically two, three or four different chromosomes (paralogons). These four fold paralogy regions have been considered as historical remnants of whole genome duplication events (WGDs/2R hypothesis). To examine the 2R hypothesis, a robust phylogenetic analysis of 40 multigene families with triplicated or quadruplicated distribution on human MHC bearing chromosomes (1/6/9/19) was conducted. Topology comparison approach categorized the members of 40 families into six distinct co-duplicated groups. Genes belonging to a particular co-duplicated group are duplicated concurrently, whereas genes of two different co-duplicated groups do not share their evolutionary history and have not duplicated in harmony. Our results based on this large scale phylogenetic data set contradict the polyploidization model and are indicative of small-scale duplications and rearrangement events that cover the entire span of animal history.

Living apart together: crosstalk between the core and supernumerary genomes in a fungal plant pathogen

BACKGROUND

Eukaryotes display remarkable genome plasticity, which can include supernumerary chromosomes that differ markedly from the core chromosomes. Despite the widespread occurrence of supernumerary chromosomes in fungi, their origin, relation to the core genome and the reason for their divergent characteristics are still largely unknown. The complexity of genome assembly due to the presence of repetitive DNA partially accounts for this.

RESULTS

Here we use single-molecule real-time (SMRT) sequencing to assemble the genome of a prominent fungal wheat pathogen, Fusarium poae, including at least one supernumerary chromosome. The core genome contains limited transposable elements (TEs) and no gene duplications, while the supernumerary genome holds up to 25 % TEs and multiple gene duplications. The core genome shows all hallmarks of repeat-induced point mutation (RIP), a defense mechanism against TEs, specific for fungi. The absence of RIP on the supernumerary genome accounts for the differences between the two (sub)genomes, and results in a functional crosstalk between them. The supernumerary genome is a reservoir for TEs that migrate to the core genome, and even large blocks of supernumerary sequence (>200 kb) have recently translocated to the core. Vice versa, the supernumerary genome acts as a refuge for genes that are duplicated from the core genome.

CONCLUSIONS

For the first time, a mechanism was determined that explains the differences that exist between the core and supernumerary genome in fungi. Different biology rather than origin was shown to be responsible. A "living apart together" crosstalk exists between the core and supernumerary genome, accelerating chromosomal and organismal evolution.

Practical interpretation of CYP2D6 haplotypes: Comparison and integration of automated and expert calling

We describe a population genetic approach to compare samples interpreted with expert calling (EC) versus automated calling (AC) for CYP2D6 haplotyping. The analysis represents 4812 haplotype calls based on signal data generated by the Luminex xMap analyzers from 2406 patients referred to a high-complexity molecular diagnostics laboratory for CYP450 testing. DNA was extracted from buccal swabs. We compared the results of expert calls (EC) and automated calls (AC) with regard to haplotype number and frequency. The ratio of EC to AC was 1:3. Haplotype frequencies from EC and AC samples were convergent across haplotypes, and their distribution was not statistically different between the groups. Most duplications required EC, as only expansions with homozygous or hemizygous haplotypes could be automatedly called. High-complexity laboratories can offer equivalent interpretation to automated calling for non-expanded CYP2D6 loci, and superior interpretation for duplications. We have validated scientific expert calling specified by scoring rules as standard operating procedure integrated with an automated calling algorithm. The integration of EC with AC is a practical strategy for CYP2D6 clinical haplotyping.

The phenotypic impact of the male-specific region of chromosome-Y in inbred mating: the role of genetic variants and gene duplications in multiple inbred rat strains

Backgound

The male-specific region of chromosome-Y (MSY) contributes to phenotypes outside of testis development and has a high rate of evolution between mammalian species. With a lack of genomic crossover, MSY is one of the few genomic areas under similar variation and evolutionary selection in inbred and outbred animal populations, allowing for an assessment of evolutionary mechanisms to translate between the populations.

Methods

Using next-generation sequencing, MSY consomic strains, molecular characterization, and large-scale phenotyping, we present here regions of MSY that contribute to inbred strain phenotypes.

Results

We have shown that (1) MSY of rat has nine autosomal gene transposition events with strain-specific selection; (2) sequence variants in MSY occur with a 1.98-fold higher number of variants than other chromosomes in seven sequenced rat strains; (3) Sry, the most studied MSY gene, has undergone extensive gene duplications, driving ubiquitous expression not seen in human or mouse; (4) the expression profile of Sry in the rat is driven by the insertion of the Sry2 copy into an intron of the ubiquitously expressed Kdm5d gene in antisense orientation, but due to several loss of function mutations in the Sry2 protein, nuclear localization and transcriptional control are decreased; (5) expression of Sry copies other than Sry2 in the rat overlaps with the expression profile for human SRY; (6) gene duplications and sequence variants (P76T) of Sry can be selected for phenotypes such as high blood pressure and androgen receptor signaling within inbred mating; and most importantly, (7) per chromosome size, MSY contributes to higher strain-specific phenotypic variation relative to all other chromosomes, with 53 phenotypes showing both a male to female and consomic cross significance.

Conclusion

The data presented supports a high probability of MSY genetic variation altering a broad range of inbred rat phenotypes.

Electronic supplementary material

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

Phylogenomic analysis reveals ancient segmental duplications in the human genome

Evolution of organismal complexity and origin of novelties during vertebrate history has been widely explored in context of both regulation of gene expression and gene duplication events. Ohno (1970) for the first time put forward the idea of two rounds whole genome duplication events as the most plausible explanation for evolutionarizing the vertebrate lineage (2R hypothesis). To test the validity of 2R hypothesis, a robust phylogenomic analysis of multigene families with triplicated or quadruplicated representation on human FGFR bearing chromosomes (4/5/8/10) was performed. Topology comparison approach categorized members of 80 families into five distinct co-duplicated groups. Genes belonging to one co-duplicated group are duplicated concurrently, whereas genes of two different co-duplicated groups do not share their duplication history and have not duplicated in congruency. Our findings contradict the 2R model and are indicative of small-scale duplications and rearrangements that cover the entire span of animal's history.

Discovery of novel plastid phenylalanine (trnF) pseudogenes defines a distinctive clade in Solanaceae

BACKGROUND

The plastome of embryophytes is known for its high degree of conservation in size, structure, gene content and linear order of genes. The duplication of entire tRNA genes or their arrangement in a tandem array composed by multiple pseudogene copies is extremely rare in the plastome. Pseudogene repeats of the trnF gene have rarely been described from the chloroplast genome of angiosperms.

FINDINGS

We report the discovery of duplicated copies of the original phenylalanine (trnFGAA) gene in Solanaceae that are specific to a larger clade within the Solanoideae subfamily. The pseudogene copies are composed of several highly structured motifs that are partial residues or entire parts of the anticodon, T- and D-domains of the original trnF gene.

CONCLUSIONS

The Pseudosolanoid clade consists of 29 genera and includes many economically important plants such as potato, tomato, eggplant and pepper.

Plastid trnF pseudogenes are present in Jaltomata, the sister genus of Solanum (Solanaceae): molecular evolution of tandemly repeated structural mutations

Extensive gene duplication arranged in a tandem array is rare in the plastome of embryophytes. Interestingly, we found pseudogene copies of the trnF gene in the genus Jaltomata, the sister genus of Solanum where such gene duplication has been previously reported. In each Jaltomata sequence available we found two pseudogene copies in close 5'-proximity to the original functional gene. The size of each pseudogene copy ranged between 17 and 48 bp and the anticodon domain was identified as the most conserved element. A common ATT(G)n motif is particularly interesting and its modifications were found to border the 3' of the duplicated regions. Other motifs were partial residues, or entire parts of the T- and D-domains, and both domains proved to be variable in length among the pseudogenes identified. The residues of the 3' and 5' acceptor stem were not found among the copies. We further compared the newly discovered copies of Jaltomata with those ones previously described from Solanum and inferred phylogenetic relationships of the copies aligned. The evolution of Solanum copies, in contrast to Jaltomata, is hard to explain as resulting only in parsimonious changes since reticulate evolutionary patterns were detected among the copies. The dynamic evolutionary patterns of Solanum might be explained by possible inter- or intrachromosomal recombination.