Lessons from the deep: mechanisms behind diversification of eukaryotic protein complexes

Genetic variation is the major mechanism behind adaptation and evolutionary change. As most proteins operate through interactions with other proteins, changes in protein complex composition and subunit sequence provide potentially new functions. Comparative genomics can reveal expansions, losses and sequence divergence within protein-coding genes, but in silico analysis cannot detect subunit substitutions or replacements of entire protein complexes. Insights into these fundamental evolutionary processes require broad and extensive comparative analyses, from both in silico and experimental evidence. Here, we combine data from both approaches and consider the gamut of possible protein complex compositional changes that arise during evolution, citing examples of complete conservation to partial and total replacement by functional analogues. We focus in part on complexes in trypanosomes as they represent one of the better studied non-animal/non-fungal lineages, but extend insights across the eukaryotes by extensive comparative genomic analysis. We argue that gene loss plays an important role in diversification of protein complexes and hence enhancement of eukaryotic diversity.

Identification of hidden N4-like viruses and their interactions with hosts

IMPORTANCE

The findings of this study are significant, as N4-like viruses represent a unique viral lineage with a distinct replication mechanism and a conserved core genome. This work has resulted in a comprehensive global map of the entire N4-like viral lineage, including information on their distribution in different biomes, evolutionary divergence, genomic diversity, and the potential for viral-mediated host metabolic reprogramming. As such, this work significantly contributes to our understanding of the ecological function and viral-host interactions of bacteriophages.

Molecular Identification and Evolutionary Divergence of the Sri Lankan Sambar Deer, Rusa unicolor (Kerr 1792)

The Sambar is one of the largest deer species distributed mainly in Asia, and it has been listed as a vulnerable species. Taxonomy based on morphological characterization has been the gold standard method used to identify the Sambar deer species. Yet, morphological identification is challenging and requires expertise. To conduct species identification and taxonomic decisions, we performed the molecular identification of R. unicolor found in Sri Lanka using DNA barcodes, COI, and Cyt b to compare the Sri Lankan R. unicolor with the Indian R. unicolor and other R. unicolor subspecies. We obtained mitochondrial DNA sequences from COI and Cyt b from blood samples collected from the wet zone in Sri Lanka. A phylogenetic tree was constructed based on the Bayesian analyses using MrBayes 3.2.7. Molecular dating was implemented in Bayesian Evolutionary Analysis Sampling Trees (BEAST v1.8.2) on the concatenated sequence using a log-normal relaxed clock and Yule species tree prior, with four categories. The results showed that the Sri Lankan R. unicolor is genetically different from the Indian R. unicolor and other R. unicolor subspecies. The divergence occurred approximately 1.1 MYA (million years ago) in the Pleistocene era. The results are essential for designing new conservation platforms for these Sambar deer species.

Model worms: knowledge gains and risks associated with the use of model species in parasitological research

Model parasite species, whose entire life cycle can be completed in the laboratory and maintained for multiple generations, have played a fundamental role in our understanding of host–parasite interactions. Yet, keeping parasites in laboratory conditions may expose them to unnatural evolutionary pressures, and using laboratory cultures for research is therefore not without limitations. Using 2 widely-used model helminth species, the cestode Hymenolepis diminuta and the nematode Heligmosomoides polygyrus, I illustrate the caution needed when interpreting experimental results on model species. I first review more than 1200 experimental studies published on these species in the past 4 decades, to determine which research areas they have contributed to. This is followed by an examination of the institutional laboratory cultures that have provided the parasites used in these studies. Some of these have persisted for decades and accounted for a substantial proportion of published studies, whereas others have been short-lived. Using information provided by the curators of active cultures, I summarize data on their origins and maintenance conditions. Finally, I discuss how laboratory cultures may have been subject to the influence of evolutionary genetic processes, such as founder effects, genetic drift and inbreeding. I also address the possibility that serial passage through laboratory hosts across multiple generations has exerted artificial selection on several parasite traits, resulting in genetic and phenotypic divergence among laboratory cultures, and between these cultures and natural parasite populations. I conclude with recommendations for the continued usage of laboratory helminth cultures aimed at maximizing their important contribution to parasitological research.

The divergence of mean phenotypes under persistent directional selection

Numerous organismal traits, particularly at the cellular level, are likely to be under persistent directional selection across phylogenetic lineages. Unless all mutations affecting such traits have large enough effects to be efficiently selected in all species, gradients in mean phenotypes are expected to arise as a consequence of differences in the power of random genetic drift, which varies by approximately five orders of magnitude across the Tree of Life. Prior theoretical work examining the conditions under which such gradients can arise focused on the simple situation in which all genomic sites affecting the trait have identical and constant mutational effects. Here, we extend this theory to incorporate the more biologically realistic situation in which mutational effects on a trait differ among nucleotide sites. Pursuit of such modifications leads to the development of semi-analytic expressions for the ways in which selective interference arises via linkage effects in single-effects models, which then extend to more complex scenarios. The theory developed clarifies the conditions under which mutations of different selective effects mutually interfere with each others' fixation and shows how variance in effects among sites can substantially modify and extend the expected scaling relationships between mean phenotypes and effective population sizes.

Morphology and niche evolution influence hummingbird speciation rates

How traits affect speciation is a long-standing question in evolution. We investigate whether speciation rates are affected by the traits themselves or by the rates of their evolution, in hummingbirds, a clade with great variation in speciation rates, morphology and ecological niches. Further, we test two opposing hypotheses, postulating that speciation rates are promoted by trait conservatism or, alternatively, by trait divergence. To address these questions, we analyse morphological (body mass and bill length) and niche traits (temperature and precipitation position and breadth, and mid-elevation), using a variety of methods to estimate speciation rates and correlate them with traits and their evolutionary rates. When it comes to the traits, we find faster speciation in smaller hummingbirds with shorter bills, living at higher elevations and experiencing greater temperature ranges. As for the trait evolutionary rates, we find that speciation increases with rates of divergence in the niche traits, but not in the morphological traits. Together, these results reveal the interplay of mechanisms through which different traits and their evolutionary rates (conservatism or divergence) influence the origination of hummingbird diversity.

RAD51 and RAD51B Play Diverse Roles in the Repair of DNA Double Strand Breaks in Physcomitrium patens

RAD51 is involved in finding and invading homologous DNA sequences for accurate homologous recombination (HR). Its paralogs have evolved to regulate and promote RAD51 functions. The efficient gene targeting and high HR rates are unique in plants only in the moss Physcomitrium patens (P. patens). In addition to two functionally equivalent RAD51 genes (RAD1-1 and RAD51-2), other RAD51 paralogues were also identified in P. patens. For elucidation of RAD51's involvement during DSB repair, two knockout lines were constructed, one mutated in both RAD51 genes (Pprad51-1-2) and the second with mutated RAD51B gene (Pprad51B). Both lines are equally hypersensitive to bleomycin, in contrast to their very different DSB repair efficiency. Whereas DSB repair in Pprad51-1-2 is even faster than in WT, in Pprad51B, it is slow, particularly during the second phase of repair kinetic. We interpret these results as PpRAD51-1 and -2 being true functional homologs of ancestral RAD51 involved in the homology search during HR. Absence of RAD51 redirects DSB repair to the fast NHEJ pathway and leads to a reduced 5S and 18S rDNA copy number. The exact role of the RAD51B paralog remains unclear, though it is important in damage recognition and orchestrating HR response.

Behavioral plasticity can facilitate evolution in urban environments

Plasticity-led evolution is central to evolutionary theory. Although challenging to study in nature, this process may be particularly apparent in novel environments such as cities. We document abundant evidence of plastic behavioral changes in urban animals, including learning, contextual, developmental, and transgenerational plasticities. Using behavioral drive as a conceptual framework, our analysis of notable case studies suggests that plastic behaviors, such as altered habitat use, migration, diurnal and seasonal activity, and courtship, can have faciliatory and cascading effects on urban evolution via spatial, temporal, and mate-choice mechanisms. Our findings highlight (i) the need to incorporate behavioral plasticity more formally into urban evolutionary research and (ii) the opportunity provided by urban environments to study behavioral mechanisms of plasticity-led processes.

Probing the genomic limits of de-extinction in the Christmas Island rat

Three principal methods are under discussion as possible pathways to “true” de-extinction; i.e., back-breeding, cloning, and genetic engineering.^1,2 Of these, while the latter approach is most likely to apply to the largest number of extinct species, its potential is constrained by the degree to which the extinct species genome can be reconstructed. We explore this question using the extinct Christmas Island rat (Rattus macleari) as a model, an endemic rat species that was driven extinct between 1898 and 1908.3, 4, 5 We first re-sequenced its genome to an average of >60× coverage, then mapped it to the reference genomes of different Rattus species. We then explored how evolutionary divergence from the extant reference genome affected the fraction of the Christmas Island rat genome that could be recovered. Our analyses show that even when the extremely high-quality Norway brown rat (R. norvegicus) is used as a reference, nearly 5% of the genome sequence is unrecoverable, with 1,661 genes recovered at lower than 90% completeness, and 26 completely absent. Furthermore, we find the distribution of regions affected is not random, but for example, if 90% completeness is used as the cutoff, genes related to immune response and olfaction are excessively affected. Ultimately, our approach demonstrates the importance of applying similar analyses to candidates for de-extinction through genome editing in order to provide critical baseline information about how representative the edited form would be of the extinct species.

•

Evolutionary divergence limits the completeness of extinct species genomes

•

The extinct Christmas Island rat was re-sequenced to ca. 60× coverage

•

Nevertheless, 4.85% of the Norway brown rat genome remains absent after mapping

•

Absences are not random; immune response and olfaction are excessively affected

Disentangling the admixed trails of grey wolf evolution

The predominant phylogenetic patterns within a genome do not always reflect correctly the history of evolutionary divergence and speciation, and the true phylogenetic signal tends to be concentrated within low-recombination regions of the genome. In this issue of Molecular Ecology, Hennelly et al. show that this is also the case for intraspecific relationships that are characterized by considerable gene flow between lineages. The study reconstructs the phylogenetic relationships of Indian and Tibetan wolves with other grey wolf (Canis lupus) populations worldwide, and demonstrates that these two populations represent phylogenetically distinct lineages. This inference was supported by using low-recombination regions of autosomal chromosomes and the X chromosome, which proved to be essential for correct inference of the lineage splitting order. Their study illustrates the power of analytical approaches that implement knowledge of genome evolution patterns to reconstruct complex intraspecific evolutionary relationships. The study also provides a compelling example of the application of modern phylogenomic approaches in the identification of evolutionarily significant units for the purpose of species conservation.

Dual Activity BLEG-1 from Bacillus lehensis G1 Revealed Structural Resemblance to B3 Metallo-β-Lactamase and Glyoxalase II: An Insight into Its Enzyme Promiscuity and Evolutionary Divergence

Metallo-β-lactamases (MBLs) are class B β-lactamases from the metallo-hydrolase-like MBL-fold superfamily which act on a broad range of β-lactam antibiotics. A previous study on BLEG-1 (formerly called Bleg1_2437), a hypothetical protein from Bacillus lehensis G1, revealed sequence similarity and activity to B3 subclass MBLs, despite its evolutionary divergence from these enzymes. Its relatedness to glyoxalase II (GLXII) raises the possibility of its enzymatic promiscuity and unique structural features compared to other MBLs and GLXIIs. This present study highlights that BLEG-1 possessed both MBL and GLXII activities with similar catalytic efficiencies. Its crystal structure revealed highly similar active site configuration to YcbL and GloB GLXIIs from Salmonella enterica, and L1 B3 MBL from Stenotrophomonas maltophilia. However, different from GLXIIs, BLEG-1 has an insertion of an active-site loop, forming a binding cavity similar to B3 MBL at the N-terminal region. We propose that BLEG-1 could possibly have evolved from GLXII and adopted MBL activity through this insertion.

Evolutionary Divergence and Biased Expression of NAC Transcription Factors in Hexaploid Bread Wheat (Triticum aestivum L.)

The NAC genes, a large plant-specific family of transcription factors, regulate a wide range of pathways involved in development and response to biotic and abiotic stress. In this study, the NAC transcription factors were identified in 27 green plants, and the results showed that NAC transcription factors in plants undergo an appearance stage from water to land and a number expansion stage from gymnosperm to angiosperm. Investigating the evolutionary process of the NAC transcription factors from diploid species to hexaploid wheat revealed that tandem replications during the polyploidization process is an important event for increasing the number of NAC transcription factors in wheat. Then, the molecular characteristics, phylogenetic relationships, and expression patterns of 462 NAC transcription factors of hexaploid wheat (TaNACs) were analyzed. The protein structure results showed that TaNAC was relatively conservative at the N-terminal that contains five subdomains. All these TaNACs were divided into Group I and Group II by phylogenetic analysis, and the TaNACs in Group I should undergo strong artificial selection based on single nucleotide polymorphism (SNP) analysis. Through genome synteny and phylogenetic analysis, these TaNACs were classified into 88 groups and 9 clusters. The biased expression results of these TaNACs showed that there are 24 groups and 67 groups of neofunctionalization genes under biotic and abiotic stress, respectively, and 16 groups and 59 groups of subfunctionalization genes. This shows that neofunctionalization plays an important role in coping with different stresses. Our study provides new insights into the evolution of NAC transcription factors in hexaploid wheat.

In silico Evolutionary Divergence Analysis Suggests the Potentiality of Capsid Protein VP2 in Serotype-Independent Foot-and-Mouth Disease Virus Detection

Foot-and-mouth disease (FMD) is an economically devastating disease of the livestock worldwide and caused by the FMD virus (FMDV), which has seven immunologically distinct serotypes (O, A, Asia1, C, and SAT1-SAT3). Studies suggest that VP2 is relatively conserved among three surface-exposed capsid proteins (VP1-VP3) of FMDV, but the level of conservation has not yet been reported. Here we analyzed the comparative evolutionary divergence of VP2 and VP1 to determine the level of conservation in VP2 at different hierarchical levels of three FMDV serotypes (O, A, and Asia1) currently circulating in Asia through an in-depth computational analysis of 14 compiled datasets and designed a consensus VP2 protein that can be used for the development of a serotype-independent FMDV detection tool. The phylogenetic analysis clearly represented a significant level of conservation in VP2 over VP1 at each subgroup level. The protein variability analysis and mutational study showed the presence of 67.4% invariant amino acids in VP2, with the N-terminal end being highly conserved. Nine inter-serotypically conserved fragments located on VP2 have been identified, among which four sites showed promising antigenicity value and surface exposure. The designed 130 amino acid long consensus VP2 protein possessed six surface-exposed B cell epitopes, which suggests the possible potentiality of the protein for the development of a serotype-independent FMDV detection tool in Asia. Conclusively, this is the first study to report the comparative evolutionary divergence between VP2 and VP1, along with proposing the possible potentiality of a designed protein candidate in serotype-independent FMDV detection.

Phylogeographic and population genetic analyses of Cucurbita moschata reveal divergence of two mitochondrial lineages linked to an elevational gradient

PREMISE

Domestication usually involves local adaptation to environmental conditions. Cucurbita species are a promising model for studying these processes. Cucurbita moschata is the third major crop in the genus because of its economic value and because it displays high landrace diversity, but research about its genetic diversity, population structure, and phylogeography is limited. We aimed at understanding how geography and elevation shape the distribution of genetic diversity in C. moschata landraces in Mexico.

METHODS

We sampled fruits from 24 localities throughout Mexico. We assessed 11 nuclear microsatellite loci, one mtDNA region, and three cpDNA regions but found no variation in cpDNA. We explored genetic structure with cluster analysis, and phylogeographic relationships with haplotype network analysis.

RESULTS

Mitochondrial genetic diversity was high, and nuclear genetic differentiation among localities was intermediate compared to other domesticated Cucurbita. We found high levels of inbreeding. We recovered two mitochondrial lineages: highland (associated with the Trans-Mexican Volcanic Belt) and lowland. Nuclear microsatellites show that localities from the Yucatan Peninsula constitute a well-differentiated group.

CONCLUSIONS

Mexico is an area of high diversity for C. moschata, and these landraces represent important plant genetic resources. In Mexico this species is characterized by divergence processes linked to an elevational gradient, which could be related to adaptation and may be of value for applications in agriculture. The Isthmus of Tehuantepec may be a partial barrier to gene flow. Morphological variation, agricultural management, and cultural differences may be related to this pattern of genetic structure, but further studies are needed.

Salt-tolerance diversity in diploid and polyploid cotton (Gossypium) species

The development of salt-tolerant genotypes is pivotal for the effective utilization of salinized land and to increase global crop productivity. Several cotton species comprise the most important source of textile fibers globally, and these are increasingly grown on marginal or increasingly saline agroecosystems. The allopolyploid cotton species also provide a model system for polyploid research, of relevance here because polyploidy was suggested to be associated with increased adaptation to stress. To evaluate the genetic variation of salt tolerance among cotton species, 17 diverse accessions of allopolyploid (AD-genome) and diploid (A- and D-genome) Gossypium were evaluated for a total of 29 morphological and physiological traits associated with salt tolerance. For most morphological and physiological traits, cotton accessions showed highly variable responses to 2 weeks of exposure to moderate (50 mm NaCl) and high (100 mm NaCl) hydroponic salinity treatments. Our results showed that the most salt-tolerant species were the allopolyploid Gossypium mustelinum from north-east Brazil, the D-genome diploid Gossypium klotzschianum from the Galapagos Islands, followed by the A-genome diploids of Africa and Asia. Generally, A-genome accessions outperformed D-genome cottons under salinity conditions. Allopolyploid accessions from either diploid genomic group did not show significant differences in salt tolerance, but they were more similar to one of the two progenitor lineages. Our findings demonstrate that allopolyploidy in itself need not be associated with increased salinity stress tolerance and provide information for using the secondary Gossypium gene pool to breed for improved salt tolerance.

Environmental variables associated with Nothophaeocryptopus gaeumannii population structure and Swiss needle cast severity in Western Oregon and Washington

The environment has a strong influence on the abundance and distribution of plant pathogenic organisms and plays a major role in plant disease. Climatological factors may also alter the dynamics of the interactions between plant pathogens and their hosts. Nothophaeocryptopus (=Phaeocryptopus) gaeumannii, the causal agent of Swiss needle cast (SNC) of Douglas-fir, is endemic to western North America where it exists as two sympatric, reproductively isolated lineages. The abundance of this fungus and the severity of SNC are strongly influenced by climate. We used statistical and population genetic analyses to examine relationships between environment, pathogen population structure, and SNC severity. Although N. gaeumannii Lineage 2 in western Oregon and Washington was most abundant where SNC symptoms were most severe, we did not detect a significant relationship between Lineage 2 and disease severity. Warmer winter temperatures were inversely correlated with foliage retention (AFR) and positively correlated with the relative abundance of Lineage 2 (PL2). However when distance inland, which was strongly correlated with both AFR and PL2, was included in the model, there was no significant relationship between Lineage 2 and AFR. Spring/early summer dew point temperatures also were positively associated with total N. gaeumannii abundance (colonization index (CI)) and inversely correlated with AFR. Warmer summer mean temperatures were associated with lower CI and higher AFR. Our results suggest that the two lineages have overlapping environmental optima, but slightly different tolerance ranges. Lineage 2 was absent from more inland sites where winters are colder and summers are warm and dry, while Lineage 1 occurred at most sites across an environmental gradient suggesting broader environmental tolerance. These relationships suggest that climate influences the abundance and distribution of this ecologically important plant pathogen and may have played a role in the evolutionary divergence of these two cryptic fungal lineages.

A Genome-Wide Association Study of Highly Heritable Agronomic Traits in Durum Wheat

Uncovering the genetic basis of key agronomic traits, and particularly of drought tolerance, addresses an important priority for durum wheat improvement. Here, a genome-wide association study (GWAS) in 493 durum wheat accessions representing a worldwide collection was employed to address the genetic basis of 17 agronomically important traits and a drought wilting score. Using a linear mixed model with 4 inferred subpopulations and a kinship matrix, we identified 90 marker-trait-associations (MTAs) defined by 78 markers. These markers could be merged into 44 genomic loci by linkage disequilibrium (r 2 > 0.2). Based on sequence alignment of the markers to the reference genome of bread wheat, we identified 14 putative candidate genes involved in enzymes, hormone-response, and transcription factors. The GWAS in durum wheat and a previous quantitative trait locus (QTL) analysis in bread wheat identified a consensus QTL locus.4B.1 conferring drought tolerance, which was further scanned for the presence of potential candidate genes. A haplotype analysis of this region revealed that two minor haplotypes were associated with both drought tolerance and reduced plant stature, thought to be the effect of linkage with the semi-dwarfing gene Rht-B1. Haplotype variants in the key chromosome 4B region were informative regarding evolutionary divergence among durum, emmer and bread wheat. Over all, the data are relevant in the context of durum wheat improvement and the isolation of genes underlying variation in some important quantitative traits.

DSH5, a dihydrosphingosine C4 hydroxylase gene family member, shows spatially restricted expression in rice and is lethal when expressed ectopically

Dihydrosphingosine C4 hydroxylase (DSH), a diiron-binding membrane enzyme, catalyzes the hydration of dihydrosphingosine and acyl-sphinganine to produce phytosphingosine and phytoceramide, respectively. Rice has two types of DSH homologs: general DSHs, namely DSH1, DSH2 and DSH4, and others that show spatial expression profiles, namely DSH3 and DSH5. The general DSHs exist in many plant species. These DSHs showed similarity in their functions and complemented the yeast sur2D mutation. In contrast, homologs of DSH3 and DSH5 were found only in monocot plants. Phylogenetic analysis placed these DSHs in different clades that are evolutionarily divergent from those of the general DSHs. DSH3 and DSH5 showed low-level expression. DSH5 expression was specifically in vascular bundle tissues. Ectopic expression of DSH5 induced a dwarf phenotype characterized by severe growth inhibition and an increase in the thickness of the leaf body caused by enlargement of bulliform cells in the leaves. However, no significant difference was observed in the amount of sphingolipid species. DSH5 did not complement the yeast sur2D mutation, implying that DSH5 has little effect on sphingolipid metabolism. These findings suggested that DSH3 and DSH5 originated and diverged in monocot plants.

Comparative Genomic and Transcriptomic Analyses of Family-1 UDP Glycosyltransferase in Prunus Mume

Glycosylation mediated by Family-1 UDP-glycosyltransferases (UGTs) plays crucial roles in plant growth and adaptation to various stress conditions. Prunus mume is an ideal crop for analyzing flowering for its early spring flowering characteristics. Revealing the genomic and transcriptomic portfolio of the UGT family in P. mume, a species in which UGTs have not yet been investigated, is therefore important. In this study, 130 putative UGT genes were identified and phylogenetically clustered into 14 groups. These PmUGTs were distributed unevenly across eight chromosomes and 32 tandem duplication and 8 segmental duplication pairs were revealed. A highly conserved intron insertion event was revealed on the basis of intron/exon patterns within PmUGTs. According to RNA-seq data, these PmUGTs were specifically expressed in different tissues and during the bud dormancy process. In addition, we confirmed the differential expression of some representative genes in response to abscisic acid treatment. Our results will provide important information on the UGT family in P. mume that should aid further characterization of their biological roles in response to environmental stress.

Phylogenetic divergence of cell biological features

Most cellular features have a range of states, but understanding the mechanisms responsible for interspecific divergence is a challenge for evolutionary cell biology. Models are developed for the distribution of mean phenotypes likely to evolve under the joint forces of mutation and genetic drift in the face of constant selection pressures. Mean phenotypes will deviate from optimal states to a degree depending on the effective population size, potentially leading to substantial divergence in the absence of diversifying selection. The steady-state distribution for the mean can even be bimodal, with one domain being largely driven by selection and the other by mutation pressure, leading to the illusion of phenotypic shifts being induced by movement among alternative adaptive domains. These results raise questions as to whether lineage-specific selective pressures are necessary to account for interspecific divergence, providing a possible platform for the establishment of null models for the evolution of cell-biological traits.

Comparative Genomics of Non-TNL Disease Resistance Genes from Six Plant Species

Disease resistance genes (R genes), as part of the plant defense system, have coevolved with corresponding pathogen molecules. The main objectives of this project were to identify non-Toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (nTNL) genes and elucidate their evolutionary divergence across six plant genomes. Using reference sequences from Arabidopsis, we investigated nTNL orthologs in the genomes of common bean, Medicago, soybean, poplar, and rice. We used Hidden Markov Models for sequence identification, performed model-based phylogenetic analyses, visualized chromosomal positioning, inferred gene clustering, and assessed gene expression profiles. We analyzed 908 nTNL R genes in the genomes of the six plant species, and classified them into 12 subgroups based on the presence of coiled-coil (CC), nucleotide binding site (NBS), leucine rich repeat (LRR), resistance to Powdery mildew 8 (RPW8), and BED type zinc finger domains. Traditionally classified CC-NBS-LRR (CNL) genes were nested into four clades (CNL A-D) often with abundant, well-supported homogeneous subclades of Type-II R genes. CNL-D members were absent in rice, indicating a unique R gene retention pattern in the rice genome. Genomes from Arabidopsis, common bean, poplar and soybean had one chromosome without any CNL R genes. Medicago and Arabidopsis had the highest and lowest number of gene clusters, respectively. Gene expression analyses suggested unique patterns of expression for each of the CNL clades. Differential gene expression patterns of the nTNL genes were often found to correlate with number of introns and GC content, suggesting structural and functional divergence.

Genome-Wide Identification of the PHD-Finger Family Genes and Their Responses to Environmental Stresses in Oryza sativa L

The PHD-finger family has been demonstrated to be involved in regulating plant growth and development. However, little information is given for its role in environmental stress responses. Here, we identified a total of 59 PHD family genes in the rice genome. These OsPHDs genes were located on eleven chromosomes and synteny analysis only revealed nine duplicated pairs within the rice PHD family. Phylogenetic analysis of all OsPHDs and PHDs from other species revealed that they could be grouped into two major clusters. Furthermore, OsPHDs were clustered into eight groups and members from different groups displayed a great divergence in terms of gene structure, functional domains and conserved motifs. We also found that with the exception of OsPHD6, all OsPHDs were expressed in at least one of the ten tested tissues and OsPHDs from certain groups were expressed in specific tissues. Moreover, our results also uncovered differential responses of OsPHDs expression to environmental stresses, including ABA (abscisic acid), water deficit, cold and high Cd. By using quantitative real-time PCR, we further confirmed the differential expression of OsPHDs under these stresses. OsPHD1/7/8/13/33 were differentially expressed under water deficit and Cd stresses, while OsPHD5/17 showed altered expression under water deficit and cold stresses. Moreover, OsPHD3/44/28 displayed differential expression under ABA and Cd stresses. In conclusion, our results provide valuable information on the rice PHD family in plant responses to environmental stress, which will be helpful for further characterizing their biological roles in responding to environmental stresses.

Molecular and pathobiological characterization of 61 Potato mop-top virus full-length cDNAs reveals great variability of the virus in the centre of potato domestication, novel genotypes and evidence for recombination

The evolutionary divergence of Potato mop-top virus (PMTV), a tri-partite, single-stranded RNA virus, is exceptionally low, based on the analysis of sequences obtained from isolates from Europe, Asia and North America. In general, RNA viruses exist as dynamic populations of closely related and recombinant genomes that are subjected to continuous genetic variation. The reason behind the low genetic variation of PMTV remains unclear. The question remains as to whether the low variability is a shared property of all PMTV isolates or is a result of the limited number of isolates characterized so far. We hypothesized that higher divergence of the virus might exist in the Andean regions of South America, the centre of potato domestication. Here, we report high variability of PMTV isolates collected from 12 fields in three locations in the Andean region of Peru. To evaluate PMTV genetic variation in Peru, we generated full-length cDNA clones, which allowed reliable comparative molecular and pathobiological characterization of individual isolates. We found significant divergence of the CP-RT and 8K sequences. The 8K cistron, which encodes a viral suppressor of RNA silencing, was found to be under diversifying selection. Phylogenetic analysis determined that, based on the CP-RT sequence, all PMTV isolates could be categorized into three separate lineages (clades). Moreover, we found evidence for recombination between two clades. Using infectious cDNA clones of the representatives of these two clades, as well as reassortants for the RNA-CP genomic component, we determined the pathobiological differences between the lineages, which we coined as S (for severe) and M (for mild) types. Interestingly, all isolates characterized previously (from Europe, Asia and North America) fall into the S-type clade, whereas most of the Peruvian isolates belong to the M-type. Taken together, our results support the notion of the single introduction of PMTV from the centre of potato origin to Europe, and subsequent spread of the S-type into Asia and USA. This is also supported by the suggested novel classification of isolates based on genetic constellations.

Drosophila maleless gene counteracts X global aneuploid effects in males

The loss of the entire X chromosome in Drosophila males generates a genome-wide aneuploid effect. We have employed a systems biology approach (microarray) to investigate the global aneuploid effect of the maleless (mle) mutation that disrupts the binding of male specific lethal (MSL) proteins that function in dosage compensation. A large number of the genes (144) that encode a broad spectrum of cellular transport proteins and transcription factors are located mainly in the autosomes of Drosophila melanogaster. We found several such targets to be sensitive to the aneuploid effect and conserved with the X chromosome in primitive Drosophila species and Anopheles gambiae. During evolution, they shifted gradually from their X-chromosomal positions to the autosomes in the species D. melanogaster, suggesting that the counteraction of trans-acting regulatory modifiers and their targets in two separate chromosomes is evolutionarily advantageous. These findings suggest a remarkable and previously unexpected level of complexity favoring natural compensation for the aneuploid effect. We propose that the MSL complex functions in dosage compensation in two separate steps. Initially, X-linked genes including X-linked trans-acting modifiers are hyperactivated in dosage compensated males. Later, these compensated modifiers eventually mitigate inverse aneuploid effects for the equality in expression of their autosomal targets in males and females.

DATABASE

Microarray data are available at the Gene Expression Omnibus (GEO) web deposit of National Centre for Biotechnology Information (NCBI) with the accession number GSE78227.

Complete Genome Sequence of Elephant Endotheliotropic Herpesvirus 4, the First Example of a GC-Rich Branch Proboscivirus

A novel group of mammalian DNA viruses called elephant endotheliotropic herpesviruses (EEHVs) belonging to the Proboscivirus genus has been associated with nearly 100 cases of highly lethal acute hemorrhagic disease in young Asian elephants worldwide. The complete 180-kb genomes of prototype strains from three AT-rich branch viruses, EEHV1A, EEHV1B, and EEHV5, have been published. However, less than 6 kb of DNA sequence each from EEHV3, EEHV4, and EEHV7 showed them to be a hugely diverged second major branch with GC-rich characteristics. Here, we determined the complete 206-kb genome of EEHV4(Baylor) directly from trunk wash DNA by next-generation sequencing and de novo assembly procedures. Among a total of 119 genes with an overall colinear organization similar to those of the AT-rich EEHVs, major features of EEHV4 include a family of 26 paralogous 7xTM and vGPCR-like genes plus 25 novel or missing genes. The genome also contains an unusual distribution of tracts of 5 to 11 successive A or T nucleotides in intergenic domains between the mostly much higher GC content protein coding regions. Furthermore, an extremely high GC-rich bias in the third wobble position of codons clearly delineates the coding regions for many but not all proteins. There are also two novel captured cellular genes, including a C-type lectin (vECTL) and an O-linked acetylglucosamine transferase (vOGT), as well as an unusually large and complex Ori-Lyt dyad symmetry domain. Finally, 30 kb from a second strain proved to include three small chimeric domains, indicating the existence of distinct EEHV4A and EEHV4B subtypes. IMPORTANCE Multiple species of herpesviruses from three different lineages of the Proboscivirus genus (EEHV1/6, EEHV2/5, and EEHV3/4/7) infect both Asian and African elephants, but lethal hemorrhagic disease is largely confined to Asian elephant calves and is predominantly associated with EEHV1. Milder disease caused by EEHV5 or EEHV4 is being increasingly recognized as well, but little is known about the latter, which is estimated to have diverged at least 35 million years ago from the others within a distinctive GC-rich branch of the Proboscivirus genus. Here, we have determined the complete genomic DNA sequence of a strain of EEHV4 obtained from a trunk wash sample collected from a surviving Asian elephant calf undergoing asymptomatic shedding during convalescence after an acute hemorrhagic disease episode. This represents the first example from among the three known GC-rich branch Proboscivirus species to be assembled and fully annotated. Several distinctive features of EEHV4 compared to AT-rich branch genomes are described.

Comparison of the Gene Coding Contents and Other Unusual Features of the GC-Rich and AT-Rich Branch Probosciviruses

Multiple species of herpesviruses from three different lineages of the Proboscivirus genus (EEHV1/6, EEHV2/5, and EEHV3/4/7) infect either Asian or African elephants, but the highly lethal hemorrhagic disease is largely confined to Asian elephant calves and is predominantly associated with EEHV1. In the accompanying paper [P. D. Ling et al., mSphere 1(3):e00081-15, http://dx.doi.org/10.1128/mSphere.00081-15], we report the complete 206-kb genome of EEHV4, the third different species causing disease in Asian elephants and the first example of a GC-rich branch proboscivirus. To gain insights into the nature and differential properties of these two very anciently diverged lineages of elephant herpesviruses, we describe here several additional unusual features found in the complete GC-rich genome of EEHV4 with particular emphasis on patterns of divergence as well as common unique features that are distinct from those of all other herpesviruses, such as the enlarged AT-rich intergenic domains and gene families, including the large number of vGPCR-like proteins.

Nearly 100 cases of lethal acute hemorrhagic disease in young Asian elephants have been reported worldwide. All tested cases contained high levels of elephant endotheliotropic herpesvirus (EEHV) DNA in pathological blood or tissue samples. Seven known major types of EEHVs have been partially characterized and shown to all belong to the novel Proboscivirus genus. However, the recently determined 206-kb EEHV4 genome proved to represent the prototype of a GC-rich branch virus that is very distinct from the previously published 180-kb EEHV1A, EEHV1B, and EEHV5A genomes, which all fall within an alternative AT-rich branch. Although EEHV4 retains the large family of 7xTM and vGPCR-like genes, six are unique to either just one or the other branch. While both branches display a highly enriched distribution of A and T tracts in intergenic domains, they are generally much larger within the GC-rich branch. Both branches retain the vGCNT1 acetylglucosamine transferase and at least one vOX-2 gene, but the two branches differ by 25 genes overall, with the AT-rich branch encoding a fucosyl transferase (vFUT9) plus two or three more vOX2 proteins and an immunoglobulin-like gene family that are all absent from the GC-rich branch. Several envelope glycoproteins retain only 15 to 20% protein identity or less across the two branches. Finally, the two plausible predicted transcriptional regulatory proteins display no homology at all to those in the alpha-, beta-, or gammaherpesvirus subfamilies. These results reinforce our previous proposal that the probosciviruses should be designated a new subfamily of mammalian herpesviruses.

IMPORTANCE Multiple species of herpesviruses from three different lineages of the Proboscivirus genus (EEHV1/6, EEHV2/5, and EEHV3/4/7) infect either Asian or African elephants, but the highly lethal hemorrhagic disease is largely confined to Asian elephant calves and is predominantly associated with EEHV1. In the accompanying paper [P. D. Ling et al., mSphere 1(3):e00081-15, 10.1128/mSphere.00081-15], we report the complete 206-kb genome of EEHV4, the third different species causing disease in Asian elephants and the first example of a GC-rich branch proboscivirus. To gain insights into the nature and differential properties of these two very anciently diverged lineages of elephant herpesviruses, we describe here several additional unusual features found in the complete GC-rich genome of EEHV4 with particular emphasis on patterns of divergence as well as common unique features that are distinct from those of all other herpesviruses, such as the enlarged AT-rich intergenic domains and gene families, including the large number of vGPCR-like proteins.

Sperm competition and the coevolution of pre- and postcopulatory traits: Weapons evolve faster than testes among onthophagine dung beetles

Reproductive competition generates episodes of both pre- and postcopulatory sexual selection. Theoretical models of sperm competition predict that as the fitness gains from expenditure on the weapons of male combat increase, males should increase their expenditure on weapons and decrease their expenditure on traits that contribute to competitive fertilization success. Although traits subject to sexual selection are known to have accelerated evolutionary rates of phenotypic divergence, it is not known whether the competing demands of investment into pre- and postcopulatory traits affect their relative rates of evolutionary divergence. We use a comparative approach to estimate the rates of divergence in pre- and postcopulatory traits among onthophagine dung beetles. Weapons evolved faster than body size while testes mass and sperm length evolved more slowly than body size, suggesting that precopulatory competition is the stronger episode of sexual selection acting on these beetles. Although horns evolved faster than testes, evolutionary increases in horn length were not associated with evolutionary reductions in testes mass. Our data for onthophagines support the notion that in taxa where males are unable to monopolize paternity, expenditure on both weapons and testes should both be favored.

CNL Disease Resistance Genes in Soybean and Their Evolutionary Divergence

Disease resistance genes (R-genes) encode proteins involved in detecting pathogen attack and activating downstream defense molecules. Recent availability of soybean genome sequences makes it possible to examine the diversity of gene families including disease-resistant genes. The objectives of this study were to identify coiled-coil NBS-LRR (= CNL) R-genes in soybean, infer their evolutionary relationships, and assess structural as well as functional divergence of the R-genes. Profile hidden Markov models were used for sequence identification and model-based maximum likelihood was used for phylogenetic analysis, and variation in chromosomal positioning, gene clustering, and functional divergence were assessed. We identified 188 soybean CNL genes nested into four clades consistent to their orthologs in Arabidopsis. Gene clustering analysis revealed the presence of 41 gene clusters located on 13 different chromosomes. Analyses of the K s-values and chromosomal positioning suggest duplication events occurring at varying timescales, and an extrapericentromeric positioning may have facilitated their rapid evolution. Each of the four CNL clades exhibited distinct patterns of gene expression. Phylogenetic analysis further supported the extrapericentromeric positioning effect on the divergence and retention of the CNL genes. The results are important for understanding the diversity and divergence of CNL genes in soybean, which would have implication in soybean crop improvement in future.

Composition and expression of conserved microRNA genes in diploid cotton (Gossypium) species

MicroRNAs are ubiquitous in plant genomes but vary greatly in their abundance within and conservation among plant lineages. To gain insight into the evolutionary birth/death dynamics of microRNA families, we sequenced small RNA and 5′-end PARE libraries generated from two closely related species of Gossypium. Here, we demonstrate that 33 microRNA families, with similar copy numbers and average evolutionary rates, are conserved in the two congeneric cottons. Analysis of the presence/absence of these microRNA families in other land plants sheds light on their depth of phylogenetic origin and lineage-specific loss/gain. Conserved microRNA families in Gossypium exhibit a striking interspecific asymmetry in expression, potentially connected to relative proximity to neighboring transposable elements. A complex correlated expression pattern of microRNA target genes with their controlling microRNAs indicates that possible functional divergence of conserved microRNA families can also exist even within a single plant genus.