Genome-wide identification and evolution of WNK kinases in Bambusoideae and transcriptional profiling during abiotic stress in Phyllostachys edulis

With-no-lysine (WNK) kinases play vital roles in abiotic stress response, circadian rhythms, and regulation of flowering time in rice, Arabidopsis, and Glycine max. However, there are no previous reports of WNKs in the Bambusoideae, although genome sequences are available for diploid, tetraploid, and hexaploid bamboo species. In the present study, we identified 41 WNK genes in five bamboo species and analysed gene evolution, phylogenetic relationship, physical and chemical properties, cis-elements, and conserved motifs. We predicted the structure of PeWNK proteins of moso bamboo and determined the exposed, buried, structural and functional amino acids. Real-time qPCR analysis revealed that PeWNK5, PeWNK7, PeWNK8, and PeWNK11 genes are involved in circadian rhythms. Analysis of gene expression of different organs at different developmental stages revealed that PeWNK genes are tissue-specific. Analysis of various abiotic stress transcriptome data (drought, salt, SA, and ABA) revealed significant gene expression levels in all PeWNKs except PeWNK11. In particular, PeWNK8 and PeWNK9 were significantly down- and up-regulated, respectively, after abiotic stress treatment. A co-expression network of PeWNK genes also showed that PeWNK2, PeWNK4, PeWNK7, and PeWNK8 were co-expressed with transcriptional regulators related to abiotic stress. In conclusion, our study identified the PeWNKs of moso bamboo involved in circadian rhythms and abiotic stress response. In addition, this study serves as a guide for future functional genomic studies of the WNK genes of the Bambusoideae.

Differential sensory and immune gene evolution in sea turtles with contrasting demographic and life histories

Marine turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 MYA, yet the genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remain largely unknown. Additionally, many populations have declined drastically due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for green (Chelonia mydas) and leatherback (Dermochelys coriacea) turtles, representing the two extant marine turtle families (MRCA ~60 MYA). Generally, these genomes are highly syntenic and homologous. Non-collinearity was associated with higher copy numbers of immune, zinc-finger, or olfactory receptor (OR) genes in green turtles. Gene family analyses suggested that ORs related to waterborne odorants have expanded in green turtles and contracted in leatherbacks, which may underlie immunological and sensory adaptations assisting navigation and occupancy of neritic versus pelagic environments, and diet specialization. Microchromosomes showed reduced collinearity, and greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, demographic history and diversity analyses showed stark contrasts between species, indicating that leatherback turtles have had a low yet stable effective population size, extremely low diversity when compared to other reptiles, and a higher proportion of deleterious variants, reinforcing concern over the persistence of this species under future climate scenarios. These highly contiguous genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Genome-Wide Identification and Comparative Analysis of the ASR Gene Family in the Rosaceae and Expression Analysis of PbrASRs During Fruit Development

The members of the Abscisic Acid (ABA) Stress and Ripening gene family (ASR) encode a class of plant-specific proteins with ABA/WDS domains that play important roles in fruit ripening, abiotic stress tolerance and biotic stress resistance in plants. The ASR gene family has been widely investigated in the monocotyledons and dicotyledons. Although the genome sequence is already available for eight fruit species of the Rosaceae, there is far less information about the evolutionary characteristics and the function of the ASR genes in the Rosaceae than in other plant families. Twenty-seven ASR genes were identified from species in the Rosaceae and divided into four subfamilies (I, II, III, and IV) on the basis of structural characteristics and phylogenetic analysis. Purifying selection was the primary force for ASR family gene evolution in eight Rosaceae species. qPCR experiments showed that the expression pattern of PbrASR genes from Pyrus bretschneideri was organ-specific, being mainly expressed in flower, fruit, leaf, and root. During fruit development, the mRNA abundance levels of different PbrASR genes were either down- or up-regulated, and were also induced by exogenous ABA. Furthermore, subcellular localization results showed that PbrASR proteins were mainly located in the nucleus and cytoplasm. These results provide a theoretical foundation for investigation of the evolution, expression, and functions of the ASR gene family in commercial fruit species of the Rosaceae family.

Nuclear transport genes recurrently duplicate by means of RNA intermediates in Drosophila but not in other insects

Background The nuclear transport machinery is involved in a well-known male meiotic drive system in Drosophila. Fast gene evolution and gene duplications have been major underlying mechanisms in the evolution of meiotic drive systems, and this might include some nuclear transport genes in Drosophila. So, using a comprehensive, detailed phylogenomic study, we examined 51 insect genomes for the duplication of the same nuclear transport genes. Results We find that most of the nuclear transport duplications in Drosophila are of a few classes of nuclear transport genes, RNA mediated and fast evolving. We also retrieve many pseudogenes for the Ran gene. Some of the duplicates are relatively young and likely contributing to the turnover expected for genes under strong but changing selective pressures. These duplications are potentially revealing what features of nuclear transport are under selection. Unlike in flies, we find only a few duplications when we study the Drosophila duplicated nuclear transport genes in dipteran species outside of Drosophila, and none in other insects. Conclusions These findings strengthen the hypothesis that nuclear transport gene duplicates in Drosophila evolve either as drivers or suppressors of meiotic drive systems or as other male-specific adaptations circumscribed to flies and involving a handful of nuclear transport functions.

Archaic introgression and variation in pharmacogenes and implications for local adaptation

Modern humans carry Neanderthal and Denisovan (archaic) genome elements which may have been a result of environmental adaptation. These effects may be particularly evident in pharmacogenes - genes responsible for the processing of exogenous substances such as food, pollutants, and medications. However, the health implications and contribution of archaic ancestry in pharmacogenes of modern humans remains understudied. We characterize eleven key cytochrome P450 (CYP450) genes involved in drug metabolizing reactions in three Neanderthal and one Denisovan individuals and examine archaic introgression in modern human populations. We infer the metabolizing efficiency of these eleven genes in archaic individuals and show important genetic differences relative to modern human variants. We identify archaic-specific SNVs in each CYP450 gene, including some that are potentially damaging, which may result in altered metabolism in modern human people carrying these variants. We highlight four genes which display interesting patterns of archaic variation: CYP2B6 - we find a large number of unique variants in the Vindija Neanderthal, some of which are shared with a small subset of African modern humans; CYP2C9 - containing multiple variants that are shared between Europeans and Neanderthals; CYP2A6*12 - a variant defined by a hybridization event that was found in humans and Neanderthals, suggesting the recombination event predates both species; and CYP2J2 - in which we hypothesize a Neanderthal variant was re-introduced in non-African populations by archaic admixture. The genetic variation identified in archaic individuals imply environmental pressures that may have driven CYP450 gene evolution.

Origin and evolutionary landscape of Nr2f transcription factors across Metazoa

Background Nuclear Receptor Subfamily 2 Group F (Nr2f) orphan nuclear hormone transcription factors (TFs) are fundamental regulators of many developmental processes in invertebrates and vertebrates. Despite the importance of these TFs throughout metazoan development, previous work has not clearly outlined their evolutionary history. Results We integrated molecular phylogeny with comparisons of intron/exon structure, domain architecture, and syntenic conservation to define critical evolutionary events that distinguish the Nr2f gene family in Metazoa. Our data indicate that a single ancestral eumetazoan Nr2f gene predated six main Bilateria subfamilies, which include single Nr2f homologs, here referred to as Nr2f1/2/5/6, that are present in invertebrate protostomes and deuterostomes, Nr2f1/2 homologs in agnathans, and Nr2f1, Nr2f2, Nr2f5, and Nr2f6 orthologs that are found in gnathostomes. Four cnidarian Nr2f1/2/5/6 and three agnathan Nr2f1/2 members are each due to independent expansions, while the vertebrate Nr2f1/Nr2f2 and Nr2f5/Nr2f6 members each form paralogous groups that arose from the established series of whole-genome duplications (WGDs). Nr2f6 members are the most divergent Nr2f subfamily in gnathostomes. Interestingly, in contrast to the other gnathostome Nr2f subfamilies, Nr2f5 has been independently lost in numerous vertebrate lineages. Furthermore, our analysis shows there are differential expansions and losses of Nr2f genes in teleosts following their additional rounds of WGDs. Conclusion Overall, our analysis of Nr2f gene evolution helps to reveal the origins and previously unrecognized relationships of this ancient TF family, which may allow for greater insights into the conservation of Nr2f functions that shape Metazoan body plans.

Protein innovation through template switching in the Saccharomyces cerevisiae lineage

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

BOPAL 2.0 and a study of tRNA and rRNA gene evolution in Clostridium

We present BOPAL 2.0, an improved version of the BOPAL algorithm for the evolutionary history inference of tRNA and rRNA genes in bacterial genomes. Our approach can infer complete evolutionary scenarios and ancestral gene orders on a phylogeny and considers a wide range of events such as duplications, deletions, substitutions, inversions and transpositions. It is based on the fact that tRNA and rRNA genes are often organized in operons/clusters in bacteria, and this information is used to help identify orthologous genes for each genome comparison. BOPAL 2.0 introduces new features, such as a triple-wise alignment step, context-aware singleton matching and a second pass of the algorithm. Evaluation on simulated datasets shows that BOPAL 2.0 outperforms the original BOPAL in terms of the accuracy of inferred events and ancestral genomes. We also present a study of the tRNA/rRNA gene evolution in the Clostridium genus, in which the organization of these genes is very divergent. Our results indicate that tRNA and rRNA genes in Clostridium have evolved through numerous duplications, losses, transpositions and substitutions, but very few inversions were inferred.

Gcorn fungi: A Web Tool for Detecting Biases between Gene Evolution and Speciation in Fungi

(1) Background: Fungi contain several millions of species, and the diversification of fungal genes has been achieved by speciation, gene duplication, and horizontal gene transfer. Although several databases provide information on orthologous and paralogous events, these databases show no information on biases between gene mutation and speciation. Here, we designed the Gcorn fungi database to better understand such biases. (2) Methods: Amino acid sequences of fungal genes in 249 species, which contain 2,345,743 sequences, were used for this database. Homologous genes were grouped at various thresholds of the homology index, which was based on the percentages of gene mutations. By grouping genes that showed highly similar homology indices to each other, we showed functional and evolutionary traits in the phylogenetic tree depicted for the gene of interest. (3) Results: Gcorn fungi provides well-summarized information on the evolution of a gene lineage and on the biases between gene evolution and speciation, which are quantitatively identified by the Robinson–Foulds metric. The database helps users visualize these traits using various depictions. (4) Conclusions: Gcorn fungi is an open access database that provides a variety of information with which to understand gene function and evolution.

The Evolution of Hemocyanin Genes in Caenogastropoda: Gene Duplications and Intron Accumulation in Highly Diverse Gastropods

Hemocyanin is the oxygen transport protein of most molluscs and represents an important physiological factor that has to be well-adapted to their environments because of the strong influences of abiotic factors on its oxygen affinity. Multiple independent gene duplications and intron gains have been reported for hemocyanin genes of Tectipleura (Heterobranchia) and the caenogastropod species Pomacea canaliculata, which contrast with the uniform gene architectures of hemocyanins in Vetigastropoda. The goal of this study was to analyze hemocyanin gene evolution within the diverse group of Caenogastropoda in more detail. Our findings reveal multiple gene duplications and intron gains and imply that these represent general features of Apogastropoda hemocyanins. Whereas hemocyanin exon–intron structures are identical within different Tectipleura lineages, they differ strongly within Caenogastropoda among phylogenetic groups as well as between paralogous hemocyanin genes of the same species. Thus, intron accumulation took place more gradually within Caenogastropoda but finally led to a similar consequence, namely, a multitude of introns. Since both phenomena occurred independently within Heterobranchia and Caenogastropoda, the results support the hypothesis that introns may contribute to adaptive radiation by offering new opportunities for genetic variability (multiple paralogs that may evolve differently) and regulation (multiple introns). Our study indicates that adaptation of hemocyanin genes may be one of several factors that contributed to the evolution of the large diversity of Apogastropoda. While questions remain, this hypothesis is presented as a starting point for the further study of hemocyanin genes and possible correlations between hemocyanin diversity and adaptive radiation.