Flavonols improve thermotolerance in tomato pollen during germination and tube elongation by maintaining ROS homeostasis

Elevated temperatures impair pollen performance and reproductive success, resulting in lower crop yields. The Solanum lycopersicum anthocyanin reduced ( are ) mutant has a FLAVANONE 3 HYDROXYLASE ( F3H ) gene mutation resulting in impaired synthesis of flavonol antioxidants. The are mutant has reduced pollen performance and seed set relative to the VF36 parental line, which is accentuated at elevated temperatures. Transformation of are with the wild-type F3H gene, or chemical complementation with flavonols, prevented temperature-dependent ROS accumulation in pollen and reversed are's reduced viability, germination, and tube elongation to VF36 levels. VF36 transformed with an F3H overexpression construct prevented temperature driven ROS increases and impaired pollen performance, revealing thermotolerance results from elevated flavonol synthesis. Although stigmas of are had reduced flavonols and elevated ROS, the growth of are pollen tubes were similarly impaired in both are and VF36 pistils. RNA-Seq was performed at optimal and stress temperatures in are , VF36, and the VF36 F3H overexpression line at multiple timepoints across pollen tube elongation. Differentially expressed gene numbers increased with duration of elevated temperature in all genotypes, with the largest number in are . These findings suggest potential agricultural interventions to combat the negative effects of heat-induced ROS in pollen that leads to reproductive failure.

One sentence summary

Flavonol antioxidants reduce the negative impacts of elevated temperatures on pollen performance by reducing levels of heat induced reactive oxygen species and modulation of heat-induced changes in the pollen transcriptome.

The Catalytic-Dependent and -Independent Roles of Lsd1 and Lsd2 Lysine Demethylases in Heterochromatin Formation in Schizosaccharomyces pombe

In eukaryotes, heterochromatin plays a critical role in organismal development and cell fate acquisition, through regulating gene expression. The evolutionarily conserved lysine-specific demethylases, Lsd1 and Lsd2, remove mono- and dimethylation on histone H3, serving complex roles in gene expression. In the fission yeast Schizosaccharomyces pombe, null mutations of Lsd1 and Lsd2 result in either severe growth defects or inviability, while catalytic inactivation causes minimal defects, indicating that Lsd1 and Lsd2 have essential functions beyond their known demethylase activity. Here, we show that catalytic mutants of Lsd1 or Lsd2 partially assemble functional heterochromatin at centromeres in RNAi-deficient cells, while the C-terminal truncated alleles of Lsd1 or Lsd2 exacerbate heterochromatin formation at all major heterochromatic regions, suggesting that Lsd1 and Lsd2 repress heterochromatic transcripts through mechanisms both dependent on and independent of their catalytic activities. Lsd1 and Lsd2 are also involved in the establishment and maintenance of heterochromatin. At constitutive heterochromatic regions, Lsd1 and Lsd2 regulate one another and cooperate with other histone modifiers, including the class II HDAC Clr3 and the Sirtuin family protein Sir2 for gene silencing, but not with the class I HDAC Clr6. Our findings explore the roles of lysine-specific demethylases in epigenetic gene silencing at heterochromatic regions.

A Deep-sequencing-assisted, Spontaneous Suppressor Screen in the Fission Yeast Schizosaccharomyces pombe

A genetic screen for mutant alleles that suppress phenotypic defects caused by a mutation is a powerful approach to identify genes that belong to closely related biochemical pathways. Previous methods such as the Synthetic Genetic Array (SGA) analysis, and random mutagenesis techniques using ultraviolet (UV) or chemicals like ethyl methanesulfonate (EMS) or N-ethyl-N- nitrosourea (ENU), have been widely used but are often costly and laborious. Also, these mutagen-based screening methods are frequently associated with severe side effects on the organism, inducing multiple mutations that add to the complexity of isolating the suppressors. Here, we present a simple and effective protocol to identify suppressor mutations in mutants which confer a growth defect in Schizosaccharomyces pombe. The fitness of cells with a growth deficiency in standard rich liquid media or synthetic liquid media can be monitored for recovery using an automated 96-well plate reader over an extended period. Once a cell acquires a suppressor mutation in the culture, its descendants outcompete those of the parental cells. The recovered cells that have a competitive growth advantage over the parental cells can then be isolated and backcrossed with the parental cells. The suppressor mutations are then identified using whole-genome sequencing. Using this approach, we have successfully isolated multiple suppressors that alleviate the severe growth defects caused by loss of Elf1, an AAA+ family ATPase that is important in nuclear mRNA transport and maintenance of genomic stability. There are currently over 400 genes in S. pombe with mutants conferring a growth defect. As many of these genes are uncharacterized, we propose that our method will hasten the identification of novel functional interactions with this user-friendly, high-throughput approach.

Encoding Data Using Biological Principles: The Multisample Variant Format for Phylogenomics and Population Genomics

Rapid progress in the fields of phylogenomics and population genomics has driven increases in both the size of multi-genomic datasets and the number and complexity of genome-wide analyses. We present the Multisample Variant Format, specifically designed to store multiple sequence alignments for phylogenomics and population genomic analysis. The signature feature of MVF is a distinctive encoding of aligned sites with specific biological information content (e.g., invariant, low-coverage). This biological pattern-based encoding of sequence data allows for rapid filtering and quality control of data and speeds up computation for many analyses. Similar to other modern formats, MVF has a simple data structure and flexible header structure to accommodate project metadata, allowing to also serve as an effective data publication and sharing format. We also propose several variants of the MVF format to accommodate protein and codon alignments, quality scores, and a mix of de novo and reference-aligned data. Using the MVFtools package, MVF files can be converted from other common sequence formats. MVFtools completes tasks ranging from simple transformation and filtering operations to complex genome-wide visualizations in only a few minutes, even on large datasets. In addition to presentation of MVF and MVFtools, we also discuss the application both in MVF and other existing data formats of the broader concept of using biological principles and patterns to inform sequence data encoding.

Quartet Sampling distinguishes lack of support from conflicting support in the green plant tree of life

PREMISE OF THE STUDY

Phylogenetic support has been difficult to evaluate within the green plant tree of life partly due to a lack of specificity between conflicted versus poorly informed branches. As data sets continue to expand in both breadth and depth, new support measures are needed that are more efficient and informative.

METHODS

We describe the Quartet Sampling (QS) method, a quartet-based evaluation system that synthesizes several phylogenetic and genomic analytical approaches. QS characterizes discordance in large-sparse and genome-wide data sets, overcoming issues of alignment sparsity and distinguishing strong conflict from weak support. We tested QS with simulations and recent plant phylogenies inferred from variously sized data sets.

KEY RESULTS

QS scores demonstrated convergence with increasing replicates and were not strongly affected by branch depth. Patterns of QS support from different phylogenies led to a coherent understanding of ancestral branches defining key disagreements, including the relationships of Ginkgo to cycads, magnoliids to monocots and eudicots, and mosses to liverworts. The relationships of ANA-grade angiosperms (Amborella, Nymphaeales, Austrobaileyales), major monocot groups, bryophytes, and fern families are likely highly discordant in their evolutionary histories, rather than poorly informed. QS can also detect discordance due to introgression in phylogenomic data.

CONCLUSIONS

Quartet Sampling is an efficient synthesis of phylogenetic tests that offers more comprehensive and specific information on branch support than conventional measures. The QS method corroborates growing evidence that phylogenomic investigations that incorporate discordance testing are warranted when reconstructing complex evolutionary histories, in particular those surrounding ANA-grade, monocots, and nonvascular plants.

Quartet Sampling distinguishes lack of support from conflicting support in the green plant tree of life

PREMISE OF THE STUDY

Phylogenetic support has been difficult to evaluate within the green plant tree of life partly due to a lack of specificity between conflicted versus poorly informed branches. As data sets continue to expand in both breadth and depth, new support measures are needed that are more efficient and informative.

METHODS

We describe the Quartet Sampling (QS) method, a quartet-based evaluation system that synthesizes several phylogenetic and genomic analytical approaches. QS characterizes discordance in large-sparse and genome-wide data sets, overcoming issues of alignment sparsity and distinguishing strong conflict from weak support. We tested QS with simulations and recent plant phylogenies inferred from variously sized data sets.

KEY RESULTS

QS scores demonstrated convergence with increasing replicates and were not strongly affected by branch depth. Patterns of QS support from different phylogenies led to a coherent understanding of ancestral branches defining key disagreements, including the relationships of Ginkgo to cycads, magnoliids to monocots and eudicots, and mosses to liverworts. The relationships of ANA-grade angiosperms (Amborella, Nymphaeales, Austrobaileyales), major monocot groups, bryophytes, and fern families are likely highly discordant in their evolutionary histories, rather than poorly informed. QS can also detect discordance due to introgression in phylogenomic data.

CONCLUSIONS

Quartet Sampling is an efficient synthesis of phylogenetic tests that offers more comprehensive and specific information on branch support than conventional measures. The QS method corroborates growing evidence that phylogenomic investigations that incorporate discordance testing are warranted when reconstructing complex evolutionary histories, in particular those surrounding ANA-grade, monocots, and nonvascular plants.

Heterogeneous molecular processes among the causes of how sequence similarity scores can fail to recapitulate phylogeny

Sequence similarity tools like Basic Local Alignment Search Tool (BLAST) are essential components of many functional genetic, genomic, phylogenetic and bioinformatic studies. Many modern analysis pipelines use significant sequence similarity scores (p- or E-values) and the ranked order of BLAST matches to test a wide range of hypotheses concerning homology, orthology, the timing of de novo gene birth/death and gene family expansion/contraction. Despite significant contrary findings, many of these tests still implicitly assume that stronger or higher-ranked E-value scores imply closer phylogenetic relationships between sequences. Here, we demonstrate that even though a general relationship does exist between the phylogenetic distance of two sequences and their E-value, significant and misleading errors occur in both the completeness and the order of results under realistic evolutionary scenarios. These results provide additional details to past evidence showing that studies should avoid drawing direct inferences of evolutionary relatedness from measures of sequence similarity alone, and should instead, where possible, use more rigorous phylogeny-based methods.

Powerful methods for detecting introgressed regions from population genomic data

Understanding the types and functions of genes that are able to cross species boundaries-and those that are not-is an important step in understanding the forces maintaining species as largely independent lineages across the remainder of the genome. With large next-generation sequencing data sets we are now able to ask whether introgression has occurred across the genome, and multiple methods have been proposed to detect the signature of such events. Here, we introduce a new summary statistic that can be used to test for introgression, RNDmin , that makes use of the minimum pairwise sequence distance between two population samples relative to divergence to an outgroup. We find that our method offers a modest increase in power over other, related tests, but that all such tests have high power to detect introgressed loci when migration is recent and strong. RNDmin is robust to variation in the mutation rate, and remains reliable even when estimates of the divergence time between sister species are inaccurate. We apply RNDmin to population genomic data from the African mosquitoes Anopheles quadriannulatus and A. arabiensis, identifying three novel candidate regions for introgression. Interestingly, one of the introgressed loci is on the X chromosome, but outside of an inversion separating these two species. Our results suggest that significant, but rare, sharing of alleles is occurring between species that diverged more than 1 million years ago, and that application of these methods to additional systems are likely to reveal similar results.

Phylogenomics Reveals Three Sources of Adaptive Variation during a Rapid Radiation

Speciation events often occur in rapid bursts of diversification, but the ecological and genetic factors that promote these radiations are still much debated. Using whole transcriptomes from all 13 species in the ecologically and reproductively diverse wild tomato clade (Solanum sect. Lycopersicon), we infer the species phylogeny and patterns of genetic diversity in this group. Despite widespread phylogenetic discordance due to the sorting of ancestral variation, we date the origin of this radiation to approximately 2.5 million years ago and find evidence for at least three sources of adaptive genetic variation that fuel diversification. First, we detect introgression both historically between early-branching lineages and recently between individual populations, at specific loci whose functions indicate likely adaptive benefits. Second, we find evidence of lineage-specific de novo evolution for many genes, including loci involved in the production of red fruit color. Finally, using a "PhyloGWAS" approach, we detect environment-specific sorting of ancestral variation among populations that come from different species but share common environmental conditions. Estimated across the whole clade, small but substantial and approximately equal fractions of the euchromatic portion of the genome are inferred to contribute to each of these three sources of adaptive genetic variation. These results indicate that multiple genetic sources can promote rapid diversification and speciation in response to new ecological opportunity, in agreement with our emerging phylogenomic understanding of the complexity of both ancient and recent species radiations.

Detection and Polarization of Introgression in a Five-Taxon Phylogeny

When multiple speciation events occur rapidly in succession, discordant genealogies due to incomplete lineage sorting (ILS) can complicate the detection of introgression. A variety of methods, including the [Formula: see text]-statistic (a.k.a. the "ABBA-BABA test"), have been proposed to infer introgression in the presence of ILS for a four-taxon clade. However, no integrated method exists to detect introgression using allelic patterns for more complex phylogenies. Here we explore the issues associated with previous systems of applying [Formula: see text]-statistics to a larger tree topology, and propose new [Formula: see text] tests as an integrated framework to infer both the taxa involved in and the direction of introgression for a symmetric five-taxon phylogeny. Using theory and simulations, we show that the [Formula: see text] statistics correctly identify the introgression donor and recipient lineages, even at low rates of introgression. [Formula: see text] is also shown to have extremely low false-positive rates. The [Formula: see text] tests are computationally inexpensive to calculate and can easily be applied to phylogenomic data sets, both genome-wide and in windows of the genome. In addition, we explore both the principles and problems of introgression detection in even more complex phylogenies.

Mosquito genomics. Extensive introgression in a malaria vector species complex revealed by phylogenomics

Introgressive hybridization is now recognized as a widespread phenomenon, but its role in evolution remains contested. Here, we use newly available reference genome assemblies to investigate phylogenetic relationships and introgression in a medically important group of Afrotropical mosquito sibling species. We have identified the correct species branching order to resolve a contentious phylogeny and show that lineages leading to the principal vectors of human malaria were among the first to split. Pervasive autosomal introgression between these malaria vectors means that only a small fraction of the genome, mainly on the X chromosome, has not crossed species boundaries. Our results suggest that traits enhancing vectorial capacity may be gained through interspecific gene flow, including between nonsister species.

Quantitative genetic analysis indicates natural selection on leaf phenotypes across wild tomato species (Solanum sect. Lycopersicon; Solanaceae)

Adaptive evolution requires both raw genetic material and an accessible path of high fitness from one fitness peak to another. In this study, we used an introgression line (IL) population to map quantitative trait loci (QTL) for leaf traits thought to be associated with adaptation to precipitation in wild tomatoes (Solanum sect. Lycopersicon; Solanaceae). A QTL sign test showed that several traits likely evolved under directional natural selection. Leaf traits correlated across species do not share a common genetic basis, consistent with a scenario in which selection maintains trait covariation unconstrained by pleiotropy or linkage disequilibrium. Two large effect QTL for stomatal distribution colocalized with key genes in the stomatal development pathway, suggesting promising candidates for the molecular bases of adaptation in these species. Furthermore, macroevolutionary transitions between vastly different stomatal distributions may not be constrained when such large-effect mutations are available. Finally, genetic correlations between stomatal traits measured in this study and data on carbon isotope discrimination from the same ILs support a functional hypothesis that the distribution of stomata affects the resistance to CO2 diffusion inside the leaf, a trait implicated in climatic adaptation in wild tomatoes. Along with evidence from previous comparative and experimental studies, this analysis indicates that leaf traits are an important component of climatic niche adaptation in wild tomatoes and demonstrates that some trait transitions between species could have involved few, large-effect genetic changes, allowing rapid responses to new environmental conditions.

More accurate phylogenies inferred from low-recombination regions in the presence of incomplete lineage sorting

When speciation events occur in rapid succession, incomplete lineage sorting (ILS) can cause disagreement among individual gene trees. The probability that ILS affects a given locus is directly related to its effective population size (Ne ), which in turn is proportional to the recombination rate if there is strong selection across the genome. Based on these expectations, we hypothesized that low-recombination regions of the genome, as well as sex chromosomes and nonrecombining chromosomes, should exhibit lower levels of ILS. We tested this hypothesis in phylogenomic datasets from primates, the Drosophila melanogaster clade, and the Drosophila simulans clade. In all three cases, regions of the genome with low or no recombination showed significantly stronger support for the putative species tree, although results from the X chromosome differed among clades. Our results suggest that recurrent selection is acting in these low-recombination regions, such that current levels of diversity also reflect past decreases in the effective population size at these same loci. The results also demonstrate how considering the genomic context of a gene tree can assist in more accurate determination of the true species phylogeny, especially in cases where a whole-genome phylogeny appears to be an unresolvable polytomy.

No excess gene movement is detected off the avian or lepidopteran Z chromosome

Most of our knowledge of sex-chromosome evolution comes from male heterogametic (XX/XY) taxa. With the genome sequencing of multiple female heterogametic (ZZ/ZW) taxa, we can now ask whether there are patterns of evolution common to both sex chromosome systems. In all XX/XY systems examined to date, there is an excess of testis-biased retrogenes moving from the X chromosome to the autosomes, which is hypothesized to result from either sexually antagonistic selection or escape from meiotic sex chromosome inactivation (MSCI). We examined RNA-mediated (retrotransposed) and DNA-mediated gene movement in two independently evolved ZZ/ZW systems, birds (chicken and zebra finch) and lepidopterans (silkworm). Even with sexually antagonistic selection likely operating in both taxa and MSCI having been identified in the chicken, we find no evidence for an excess of genes moving from the Z chromosome to the autosomes in either lineage. We detected no excess for either RNA- or DNA-mediated duplicates, across a range of approaches and methods. We offer some potential explanations for this difference between XX/XY and ZZ/ZW sex chromosome systems, but further work is needed to distinguish among these hypotheses. Regardless of the root causes, we have identified an additional, potentially inherent, difference between XX/XY and ZZ/ZW systems.

On the possible role of robustness in the evolution of infectious diseases

Robustness describes the capacity for a biological system to remain canalized despite perturbation. Genetic robustness affords maintenance of phenotype despite mutational input, necessarily involving the role of epistasis. Environmental robustness is phenotypic constancy in the face of environmental variation, where epistasis may be uninvolved. Here we discuss genetic and environmental robustness, from the standpoint of infectious disease evolution, and suggest that robustness may be a unifying principle for understanding how different disease agents evolve. We focus especially on viruses with RNA genomes due to their importance in the evolution of emerging diseases and as model systems to test robustness theory. We present new data on adaptive constraints for a model RNA virus challenged to evolve in response to UV radiation. We also draw attention to other infectious disease systems where robustness theory may prove useful for bridging evolutionary biology and biomedicine, especially the evolution of antibiotic resistance in bacteria, immune evasion by influenza, and malaria parasite infections.