Interlocus Gene Conversion, Natural Selection, and Paralog Homogenization.

Following a duplication, the resulting paralogs tend to diverge. While mutation and natural selection can accelerate this process, they can also slow it. Here, we quantify the paralog homogenization that is caused by point mutations and interlocus gene conversion (IGC). Among 164 duplicated teleost genes, the median percentage of postduplication codon substitutions that arise from IGC rather than point mutation is estimated to be between 7% and 8%. By differentiating between the nonsynonymous codon substitutions that homogenize the protein sequences of paralogs and the nonhomogenizing nonsynonymous substitutions, we estimate the homogenizing nonsynonymous rates to be higher for 163 of the 164 teleost data sets as well as for all 14 data sets of duplicated yeast ribosomal protein-coding genes that we consider. For all 14 yeast data sets, the estimated homogenizing nonsynonymous rates exceed the synonymous rates.

Genetic adaptations in the population history of Arabidopsis thaliana.

A population encounters a variety of environmental stresses, so the full source of its resilience can only be captured by collecting all the signatures of adaptation to the selection of the local environment in its population history. Based on the multiomic data of Arabidopsis thaliana, we constructed a database of phenotypic adaptations (p-adaptations) and gene expression (e-adaptations) adaptations in the population. Through the enrichment analysis of the identified adaptations, we inferred a likely scenario of adaptation that is consistent with the biological evidence from experimental work. We analyzed the dynamics of the allele frequencies at the 23,880 QTLs of 174 traits and 8,618 eQTLs of 1,829 genes with respect to the total SNPs in the genomes and identified 650 p-adaptations and 3,925 e-adaptations [false discovery rate (FDR) = 0.05]. The population underwent large-scale p-adaptations and e-adaptations along 4 lineages. Extremely cold winters and short summers prolonged seed dormancy and expanded the root system architecture. Low temperatures prolonged the growing season, and low light intensity required the increased chloroplast activity. The subtropical and humid environment enhanced phytohormone signaling pathways in response to the biotic and abiotic stresses. Exposure to heavy metals selected alleles for lower heavy metal uptake from soil, lower growth rate, lower resistance to bacteria, and higher expression of photosynthetic genes were selected. The p-adaptations are directly interpretable, while the coadapted gene expressions reflect the physiological requirements for the adaptation. The integration of this information characterizes when and where the population has experienced environmental stress and how the population responded at the molecular level.

Massive Loss of Transcription Factors Promotes the Initial Diversification of Placental Mammals.

As one of the most successful group of organisms, mammals occupy a variety of niches on Earth as a result of macroevolution. Transcription factors (TFs), the fundamental regulators of gene expression, may also have evolved. To examine the relationship between TFs and mammalian macroevolution, we analyzed 140,821 de novo-identified TFs and their birth and death histories from 96 mammalian species. Gene tree vs. species tree reconciliation revealed that placental mammals experienced an upsurge in TF losses around 100 million years ago (Mya) and also near the Cretaceous-Paleogene boundary (K-Pg boundary, 66 Mya). Early Euarchontoglires, Laurasiatheria and marsupials appeared between 100 and 95 Mya and underwent initial diversification. The K-Pg boundary was associated with the massive extinction of dinosaurs, which lead to adaptive radiation of mammals. Surprisingly, TF loss decelerated, rather than accelerated, molecular evolutionary rates of their target genes. As the rate of molecular evolution is affected by the mutation rate, the proportion of neutral mutations and the population size, the decrease in molecular evolution may reflect increased functional constraints to survive target genes.

Molecular Evolutionary Rate Predicts Intraspecific Genetic Polymorphism and Species-Specific Selection.

It is unknown what determines genetic diversity and how genetic diversity is associated with various biological traits. In this work, we provide insight into these issues. By comparing genetic variation of 14,671 mammalian gene trees with thousands of individual human, chimpanzee, gorilla, mouse, and dog/wolf genomes, we found that intraspecific genetic diversity can be predicted by long-term molecular evolutionary rates rather than de novo mutation rates. This relationship was established during the early stage of mammalian evolution. Moreover, we developed a method to detect fluctuations of species-specific selection on genes based on the deviations of intraspecific genetic diversity predicted from long-term rates. We showed that the evolution of epithelial cells, rather than connective tissue, mainly contributed to morphological evolution of different species. For humans, evolution of the immune system and selective sweeps caused by infectious diseases are the most representative examples of adaptive evolution.

Exploratory analysis of multi-trait coadaptations in light of population history.

During the process of range expansion, populations encounter a variety of environments. They respond to the local environments by modifying their mutually interacting traits. Common approaches of landscape analysis include first focusing on the genes that undergo diversifying selection or directional selection in response to environmental variation. To understand the whole history of populations, it is ideal to capture the history of their range expansion with reference to the series of surrounding environments and to infer the multitrait coadaptation. To this end, we propose a complementary approach; it is an exploratory analysis using up-to-date methods that integrate population genetic features and features of selection on multiple traits. First, we conduct correspondence analysis of site frequency spectra, traits, and environments with auxiliary information of population-specific fixation index (F ST). This visualizes the structure and the ages of populations and helps infer the history of range expansion, encountered environmental changes, and selection on multiple traits. Next, we further investigate the inferred history using an admixture graph that describes the population split and admixture. Finally, principal component analysis of the selection on edge-by-trait (SET) matrix identifies multitrait coadaptation and the associated edges of the admixture graph. We introduce a newly defined factor loadings of environmental variables in order to identify the environmental factors that caused the coadaptation. A numerical simulation of one-dimensional stepping-stone population expansion showed that the exploratory analysis reconstructed the pattern of the environmental selection that was missed by analysis of individual traits. Analysis of a public dataset of natural populations of black cottonwood in northwestern America identified the first principal component (PC) coadaptation of photosynthesis- vs growth-related traits responding to the geographical clines of temperature and day length. The second PC coadaptation of volume-related traits suggested that soil condition was a limiting factor for aboveground environmental selection.

Population structure of chum salmon and selection on the markers collected for stock identification.

Genetic stock identification (GSI) is a major management tool of Pacific salmon (Oncorhynchus Spp.) that has provided rich genetic baseline data of allozymes, microsatellites, and single-nucleotide polymorphisms (SNPs) across the Pacific Rim. Here, we analyzed published data sets for adult chum salmon (Oncorhynchus keta), namely 10 microsatellites, 53 SNPs, and a mitochondrial DNA locus (mtDNA3, control region, and NADH-3 combined) in samples from 495 locations in the same distribution range (n = 61,813). TreeMix analysis of the microsatellite loci identified the greatest convergence toward Japanese/Korean populations and suggested two admixture events from Japan/Korea to Russia and the Alaskan Peninsula. The SNPs had been purposively collected from rapidly evolving genes to increase the power of GSI. The largest expected heterozygosity was observed in Japanese/Korean populations for microsatellites, whereas it was largest in Western Alaskan populations for SNPs, reflecting the SNP discovery process. A regression of SNP population structures on those of microsatellites indicated the selection of the SNP loci according to deviations from the predicted structures. Specifically, we matched the sampling locations of the SNPs with those of the microsatellites and performed regression analyses of SNP allele frequencies on a 2-dimensional scaling (MDS) of matched locations obtained from microsatellite pairwise F ST values. The MDS first axis indicated a latitudinal cline in American and Russian populations, whereas the second axis showed differentiation of Japanese/Korean populations. The top five outlier SNPs included mtDNA3, U502241 (unknown), GnRH373, ras1362, and TCP178, which were identified by principal component analysis. We summarized the functions of 53 nuclear genes surrounding SNPs and the mtDNA3 locus by referring to a gene database system and propose how they may influence the fitness of chum salmon.

Understanding population structure in an evolutionary context: population-specific FST and pairwise FST.

Populations are shaped by their history. It is crucial to interpret population structure in an evolutionary context. Pairwise FST measures population structure, whereas population-specific FST measures deviation from the ancestral population. To understand the current population structure and a population's history of range expansion, we propose a representation method that overlays population-specific FST estimates on a sampling location map, and on an unrooted neighbor-joining tree and a multi-dimensional scaling plot inferred from a pairwise FST distance matrix. We examined the usefulness of our procedure using simulations that mimicked population colonization from an ancestral population and by analyzing published human, Atlantic cod, and wild poplar data. Our results demonstrated that population-specific FST values identify the source population and trace the evolutionary history of its derived populations. Conversely, pairwise FST values represent the current population structure. By integrating the results of both estimators, we obtained a new picture of the population structure that incorporates evolutionary history. The generalized least squares estimate of genome-wide population-specific FST indicated that the wild poplar population expanded its distribution to the north, where daylight hours are long in summer, to coastal areas with abundant rainfall, and to the south where summers are dry. Genomic data highlight the power of the bias-corrected moment estimators of FST, whether global, pairwise, or population-specific, that provide unbiased estimates of FST. All FST moment estimators described in this paper have reasonable processing times and are useful in population genomics studies.

Evolution of Reproductive Life History in Mammals and the Associated Change of Functional Constraints.

Phylogenetic trees based on multiple genomic loci enable us to estimate the evolution of functional constraints that operate on genes based on lineage-specific fluctuation of the evolutionary rate at particular gene loci, "gene-branch interactions". Using this information as predictors, our previous work inferred that the common ancestor of placental mammals was nocturnal, insectivorous, solitary, and bred seasonally. Here, we added seven new continuous traits including lifespan, bodyweight, and five reproduction-related traits and inferred the coevolution network of 14 core life history traits for 89 mammals. In this network, bodyweight and lifespan are not directly connected to each other; instead, their correlation is due to both of them coevolving with gestation period. Diurnal mammals are more likely to be monogamous than nocturnal mammals, while arboreal mammals tend to have a smaller litter size than terrestrial mammals. Coevolution between diet and the seasonal breeding behavior test shows that year-round breeding preceded the dietary change to omnivory, while seasonal breeding preceded the dietary change to carnivory. We also discuss the evolution of reproductive strategy of mammals. Genes selected as predictors were identified as well; for example, genes function as tumor suppressor were selected as predictors of weaning age.

Waves Out of the Korean Peninsula and Inter- and Intra-Species Replacements in Freshwater Fishes in Japan.

The Japanese archipelago is located at the periphery of the continent of Asia. Rivers in the Japanese archipelago, separated from the continent of Asia by about 17 Ma, have experienced an intermittent exchange of freshwater fish taxa through a narrow land bridge generated by lowered sea level. As the Korean Peninsula and Japanese archipelago were not covered by an ice sheet during glacial periods, phylogeographical analyses in this region can trace the history of biota that were, for a long time, beyond the last glacial maximum. In this study, we analyzed the phylogeography of four freshwater fish taxa, Hemibarbus longirostris, dark chub Nipponocypris temminckii, Tanakia ssp. and Carassius ssp., whose distributions include both the Korean Peninsula and Western Japan. We found for each taxon that a small component of diverse Korean clades of freshwater fishes migrated in waves into the Japanese archipelago to form the current phylogeographic structure of biota. The replacements of indigenous populations by succeeding migrants may have also influenced the phylogeography.

Multiple Isolated Transcription Factors Act as Switches and Contribute to Species Uniqueness.

Mammals have variable numbers (1300-2000) of transcription factors (TFs), but the reasons for this large variation are unclear. To investigate general TF patterns, we de novo identified 156,906 TFs from 96 mammalian species. We identified more than 500 human isolated TFs that are rarely reported in human TF-to-TF networks. Mutations in the genes of these TFs were less lethal than those of connected TFs. Consequently, these isolated TFs are more tolerant of changes and have become unique during speciation. They may also serve as a source of variation for TF evolution. Reconciliation of TF-family phylogenetic trees with a mammalian species tree revealed an average of 37.8% TF gains and 15.0% TF losses over 177 million years, which implies that isolated TFs are pervasive in mammals. Compared with non-TF interacting genes, TF-interacting genes have unique TF profiles and have higher expression levels in mice than in humans. Different expression levels of the same TF-interacting gene contribute to species-specific phenotypes. Formation and loss of isolated TFs enabling unique TF profiles may provide variable switches that adjust divergent expression profiles of target genes to generate species-specific phenotypes, thereby making species unique.

Renormalized basal metabolic rate describes the human aging process and longevity.

The question of why we age and finally die has been a central subject in the life, medical, and health sciences. Many aging theories have proposed biomarkers that are related to aging. However, they do not have sufficient power to predict the aging process and longevity. We here propose a new biomarker of human aging based on the mass-specific basal metabolic rate (msBMR). It is well known by the Harris-Benedict equation that the msBMR declines with age but varies among individual persons. We tried to renormalize the msBMR by primarily incorporating the body mass index into this equation. The renormalized msBMR (RmsBMR) which was derived in one cohort of American men (n = 25,425) was identified as one of the best biomarkers of aging, because it could well reproduce the observed respective American, Italian, and Japanese data on the mortality rate and survival curve. A recently observed plateau of the mortality rate in centenarians corresponded to the lowest value (threshold) of the RmsBMR, which stands for the final stage of human life. A universal decline of the RmsBMR with age was associated with the mitochondrial number decay, which was caused by a slight fluctuation of the dynamic fusion/fission system. This decay form was observed by the measurement in mice. Finally, the present approach explained the reason why the BMR in mammals is regulated by the empirical algometric scaling law.

Rigorous monitoring of a large-scale marine stock enhancement program demonstrates the need for comprehensive management of fisheries and nursery habitat.

Release of captively-bred individuals into the wild is one of the most popular tools in fisheries, forestry, and wildlife management, and introgression of hatchery-reared animals into wild populations is of global concern. However, research and monitoring of impacts on wild populations are generally lacking, and the benefit of hatcheries for long-term fisheries and conservation goals is unclear. Using spatio-temporal genetic monitoring and a four-dacade time series of catch data, we quantified the effects on the size and genetic diversity of wild populations of one of the world's largest marine stock enhancement programs-the red sea bream (Pagrus major) in Kagoshima Bay, Japan. Our analyses found that the stock enhancement program reduced genetic diversity of the population, but the genetic effect diminished with increased size of the wild population. Increases to the seaweed communities and reduced fishing efforts were the primary factors associated with the wild population recovery; effects of aquaculture were much smaller. Our results represent crucial evidence that hatcheries for enhancement and conservation of populations cannot be successful over the long term unless sufficient efforts are also made to reduce harvest rates and rehabilitate natural habitats.

Stochastic Variational Inference for Bayesian Phylogenetics: A Case of CAT Model.

The pattern of molecular evolution varies among gene sites and genes in a genome. By taking into account the complex heterogeneity of evolutionary processes among sites in a genome, Bayesian infinite mixture models of genomic evolution enable robust phylogenetic inference. With large modern data sets, however, the computational burden of Markov chain Monte Carlo sampling techniques becomes prohibitive. Here, we have developed a variational Bayesian procedure to speed up the widely used PhyloBayes MPI program, which deals with the heterogeneity of amino acid profiles. Rather than sampling from the posterior distribution, the procedure approximates the (unknown) posterior distribution using a manageable distribution called the variational distribution. The parameters in the variational distribution are estimated by minimizing Kullback-Leibler divergence. To examine performance, we analyzed three empirical data sets consisting of mitochondrial, plastid-encoded, and nuclear proteins. Our variational method accurately approximated the Bayesian inference of phylogenetic tree, mixture proportions, and the amino acid propensity of each component of the mixture while using orders of magnitude less computational time.

A new approach for comprehensively describing heterogametic sex chromosomes.

Notwithstanding the rapid developments in sequencing techniques, Y and W sex chromosomes have still been mostly excluded from whole genome sequencing projects due to their high repetitive DNA content. Therefore, Y and W chromosomes are poorly described in most species despite their biological importance. Several methods were developed for identifying Y or W-linked sequences among unmapped scaffolds. However, it is not enough to discover functional regions from short unmapped scaffolds. Here, we provide a new and simple strategy based on k-mer comparison for comprehensive analysis of the W chromosome in Bombyx mori. Using this novel method, we effectively assembled de novo 1281 W-derived genome contigs (totaling 1.9 Mbp), and identified 156 W-linked transcript RNAs and 345 W-linked small RNAs. This method will help in the elucidation of mechanisms of sexual development and exploration of W chromosome biological functions, and provide insights into the evolution of sex chromosomes. Moreover, we showed this method can be employed in identifying heterogametic sex chromosomes (W and Y chromosomes) in many other species where genomic information is still scarce.

Rates of Molecular Evolution Suggest Natural History of Life History Traits and a Post-K-Pg Nocturnal Bottleneck of Placentals.

Life history and behavioral traits are often difficult to discern from the fossil record, but evolutionary rates of genes and their changes over time can be inferred from extant genomic data. Under the neutral theory, molecular evolutionary rate is a product of mutation rate and the proportion of neutral mutations [1, 2]. Mutation rates may be shared across the genome, whereas proportions of neutral mutations vary among genes because functional constraints vary. By analyzing evolutionary rates of 1,185 genes in a phylogeny of 89 mammals, we extracted historical profiles of functional constraints on these rates in the form of gene-branch interactions. By applying a novel statistical approach to these profiles, we reconstructed the history of ten discrete traits related to activity, diet, and social behaviors. Our results indicate that the ancestor of placental mammals was solitary, seasonally breeding, insectivorous, and likely nocturnal. The results suggest placental diversification began 10-20 million years before the K-Pg boundary (66 million years ago), with some ancestors of extant placental mammals becoming diurnal and adapted to different diets. However, from the Paleocene to the Eocene-Oligocene transition (EOT, 33.9 mya), we detect a post-K-Pg nocturnal bottleneck where all ancestral lineages of extant placentals were nocturnal. Although diurnal placentals may have existed during the elevated global temperatures of the Paleocene-Eocene thermal maximum [3], we hypothesize that diurnal placentals were selectively extirpated during or after the global cooling of the EOT, whereas some nocturnal lineages survived due to preadaptations to cold environments [4].

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest.

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect's natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India-South China-Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

The empirical Bayes estimators of fine-scale population structure in high gene flow species.

An empirical Bayes (EB) pairwise FST estimator was previously introduced and evaluated for its performance by numerical simulation. In this study, we conducted coalescent simulations and generated genetic population structure mechanistically, and compared the performance of the EBFST with Nei's GST , Nei and Chesser's bias-corrected GST (GST_NC ), Weir and Cockerham's θ (θWC ) and θ with finite sample correction (θWC_F ). We also introduced EB estimators for Hedrick' G'ST and Jost' D. We applied these estimators to publicly available SNP genotypes of Atlantic herring. We also examined the power to detect the environmental factors causing the population structure. Our coalescent simulations revealed that the finite sample correction of θWC is necessary to assess population structure using pairwise FST values. For microsatellite markers, EBFST performed the best among the present estimators regarding both bias and precision under high gene flow scenarios (FST≤0.032). For 300 SNPs, EBFST had the highest precision in all cases, but the bias was negative and greater than those for GST_NC and θWC_F in all cases. GST_NC and θWC_F performed very similarly at all levels of FST . As the number of loci increased up to 10 000, the precision of GST_NC and θWC_F became slightly better than for EBFST for cases with FST≥0.004, even though the size of the bias remained constant. The EB estimators described the fine-scale population structure of the herring and revealed that ~56% of the genetic differentiation was caused by sea surface temperature and salinity. The R package finepop for implementing all estimators used here is available on CRAN.

Expression map of a complete set of gustatory receptor genes in chemosensory organs of Bombyx mori.

Most lepidopteran species are herbivores, and interaction with host plants affects their gene expression and behavior as well as their genome evolution. Gustatory receptors (Grs) are expected to mediate host plant selection, feeding, oviposition and courtship behavior. However, due to their high diversity, sequence divergence and extremely low level of expression it has been difficult to identify precisely a complete set of Grs in Lepidoptera. By manual annotation and BAC sequencing, we improved annotation of 43 gene sequences compared with previously reported Grs in the most studied lepidopteran model, the silkworm, Bombyx mori, and identified 7 new tandem copies of BmGr30 on chromosome 7, bringing the total number of BmGrs to 76. Among these, we mapped 68 genes to chromosomes in a newly constructed chromosome distribution map and 8 genes to scaffolds; we also found new evidence for large clusters of BmGrs, especially from the bitter receptor family. RNA-seq analysis of diverse BmGr expression patterns in chemosensory organs of larvae and adults enabled us to draw a precise organ specific map of BmGr expression. Interestingly, most of the clustered genes were expressed in the same tissues and more than half of the genes were expressed in larval maxillae, larval thoracic legs and adult legs. For example, BmGr63 showed high expression levels in all organs in both larval and adult stages. By contrast, some genes showed expression limited to specific developmental stages or organs and tissues. BmGr19 was highly expressed in larval chemosensory organs (especially antennae and thoracic legs), the single exon genes BmGr53 and BmGr67 were expressed exclusively in larval tissues, the BmGr27-BmGr31 gene cluster on chr7 displayed a high expression level limited to adult legs and the candidate CO2 receptor BmGr2 was highly expressed in adult antennae, where few other Grs were expressed. Transcriptional analysis of the Grs in B. mori provides a valuable new reference for finding genes involved in plant-insect interactions in Lepidoptera and establishing correlations between these genes and vital insect behaviors like host plant selection and courtship for mating.

Phylogenomics and Morphology of Extinct Paleognaths Reveal the Origin and Evolution of the Ratites.

The Palaeognathae comprise the flightless ratites and the volant tinamous, and together with the Neognathae constitute the extant members of class Aves. It is commonly believed that Palaeognathae originated in Gondwana since most of the living species are found in the Southern Hemisphere [1-3]. However, this hypothesis has been questioned because the fossil paleognaths are mostly from the Northern Hemisphere in their earliest time (Paleocene) and possessed many putative ancestral characters [4]. Uncertainties regarding the origin and evolution of Palaeognathae stem from the difficulty in estimating their divergence times [1, 2] and their remarkable morphological convergence. Here, we recovered nuclear genome fragments from extinct elephant birds, which enabled us to reconstruct a reliable phylogenomic time tree for the Palaeognathae. Based on the tree, we identified homoplasies in morphological traits of paleognaths and reconstructed their morphology-based phylogeny including fossil species without molecular data. In contrast to the prevailing theories, the fossil paleognaths from the Northern Hemisphere were placed as the basal lineages. Combined with our stable divergence time estimates that enabled a valid argument regarding the correlation with geological events, we propose a new evolutionary scenario that contradicts the traditional view. The ancestral Palaeognathae were volant, as estimated from their molecular evolutionary rates, and originated during the Late Cretaceous in the Northern Hemisphere. They migrated to the Southern Hemisphere and speciated explosively around the Cretaceous-Paleogene boundary. They then extended their distribution to the Gondwana-derived landmasses, such as New Zealand and Madagascar, by overseas dispersal. Gigantism subsequently occurred independently on each landmass.

Cretaceous origin of giant rhinoceros beetles (Dynastini; Coleoptera) and correlation of their evolution with the Pangean breakup.

The giant rhinoceros beetles (Dynastini, Scarabaeidae, Coleoptera) are distributed in tropical and temperate regions in Asia, America and Africa. Recent molecular phylogenetic studies have revealed that the giant rhinoceros beetles can be divided into three clades representing Asia, America and Africa. Although a correlation between their evolution and the continental drift during the Pangean breakup was suggested, there is no accurate divergence time estimation among the three clades based on molecular data. Moreover, there is a long chronological gap between the timing of the Pangean breakup (Cretaceous: 110-148 Ma) and the emergence of the oldest fossil record (Oligocene: 33 Ma). In this study, we estimated their divergence times based on molecular data, using several combinations of fossil calibration sets, and obtained robust estimates. The inter-continental divergence events among the clades were estimated to have occurred about 99 Ma (Asian clade and others) and 78 Ma (American clade and African clade), both of which are after the Pangean breakup. These estimates suggest their inter-continental divergences occurred by overseas sweepstakes dispersal, rather than by vicariances of the population caused by the Pangean breakup.