Is speciation an unrelenting march to reproductive isolation?

Speciation is often portrayed as an "incomplete" or "incipient" process if two groups of organisms, technically distinguishable either by morphology or genetics, can exchange genes. The ultimate outcome of diversification, given this perspective, is complete reproductive isolation. But an increasing amount of evidence suggests that speciation is rarely complete and inter-fertility between different taxonomically accepted species is consistently maintained. In this issue of Molecular Ecology, Linan et al. (2021) provide results that bridge evolutionary processes from populations to phylogenies that indicate suites of closely related tree species in the Mascarene Islands actively exchange genes, evolving as a nested set of syngameons with a hierarchical pattern of interfertility. The deep insight into diversification provided by this study is particularly powerful because of the genomic scale of the data and the complete taxonomic sampling of an island clade evolving in situ. The prevalence of syngameon dynamics in a broad range of organisms indicates that we should adopt a fluid and comprehensive approach to defining evolutionary units for conservation and research. We should move beyond focusing on single endangered species in evolutionary and ecological isolation from other species but consider the entire network of potentially interfertile species and the potential for future adaptation and innovation, particularly in a human dominated world.

The oak syngameon: more than the sum of its parts.

One of Anthropocene's most daunting challenges for conservation biology is habitat extinction, caused by rapid global change. Tree diversity has persisted through previous episodes of rapid change, even global extinctions. Given the pace of current change, our management of extant diversity needs to facilitate and even enhance the natural ability of trees to adapt and diversify. Numerous processes contribute to this evolutionary flexibility, including introgression, a widespread yet under-studied process. Reproductive networks, in which species remain distinct despite interspecific gene flow, are called syngameons, a concept largely inspired from work focusing on Quercus. Delineating and analyzing such species groups, empirically and theoretically, will provide insights into the nonadditive effects on evolution of numerous partially interfertile species exchanging genetic material episodically under changing environmental conditions. To conserve tree diversity, crossing experiments designed with an empirical and theoretical understanding of the constituent syngameon should be set up to assist diversification and adaptation in the Anthropocene. Our increasingly detailed knowledge of the oak genome and of oak interspecific and intraspecific phenotypic variation will improve our ability to sustain the diversity of this tree through an unpredictable and unprecedented future.

Fruit development of Lithocarpus (Fagaceae) and the role of heterochrony in their evolution.

Stone oaks, or Lithocarpus species of Fagaceae are ecologically important canopy trees in the tropical and subtropical forests over East Asia, and the fruits of which are important food sources for insects and vertebrates there. The great fruit morphological variation of this genus represents two fruit types, acorn and enclosed receptacle fruit types. However, the evolutionary mechanisms of differentiation into these two fruit types with contrasting morphology remain a puzzle. To reveal the morphogenetic properties of two fruit types, we observed tissue differentiation and development among 20 Lithocarpus species from fruit set to maturity. Unlike in fruits of Quercus, the endocarp differentiation in Lithocarpus fruits occurred later than exocarp and mesocarp. Cupules provided further protection of developing seeds, particularly of acorn-type fruits. Fruits of Lithocarpus and Quercus acorns share similar insect predators. At fruit set, both acorn and enclosed receptacle types were largely identical, with similar tissue morphology and the sequence of differentiation. The distinct difference between two fruit types at maturity came from varied rates and degrees of development between the pericarp and receptacle tissues. We found that heterochrony between two tissues could create substantially divergent ecological strategies for protection and dispersal of their seeds, which is essential for the evolution of two fruit types.

Averting biodiversity collapse in tropical forest protected areas.

The rapid disruption of tropical forests probably imperils global biodiversity more than any other contemporary phenomenon. With deforestation advancing quickly, protected areas are increasingly becoming final refuges for threatened species and natural ecosystem processes. However, many protected areas in the tropics are themselves vulnerable to human encroachment and other environmental stresses. As pressures mount, it is vital to know whether existing reserves can sustain their biodiversity. A critical constraint in addressing this question has been that data describing a broad array of biodiversity groups have been unavailable for a sufficiently large and representative sample of reserves. Here we present a uniquely comprehensive data set on changes over the past 20 to 30 years in 31 functional groups of species and 21 potential drivers of environmental change, for 60 protected areas stratified across the world's major tropical regions. Our analysis reveals great variation in reserve 'health': about half of all reserves have been effective or performed passably, but the rest are experiencing an erosion of biodiversity that is often alarmingly widespread taxonomically and functionally. Habitat disruption, hunting and forest-product exploitation were the strongest predictors of declining reserve health. Crucially, environmental changes immediately outside reserves seemed nearly as important as those inside in determining their ecological fate, with changes inside reserves strongly mirroring those occurring around them. These findings suggest that tropical protected areas are often intimately linked ecologically to their surrounding habitats, and that a failure to stem broad-scale loss and degradation of such habitats could sharply increase the likelihood of serious biodiversity declines.

Inference of Markovian properties of molecular sequences from NGS data and applications to comparative genomics.

MOTIVATION: Next-generation sequencing (NGS) technologies generate large amounts of short read data for many different organisms. The fact that NGS reads are generally short makes it challenging to assemble the reads and reconstruct the original genome sequence. For clustering genomes using such NGS data, word-count based alignment-free sequence comparison is a promising approach, but for this approach, the underlying expected word counts are essential.A plausible model for this underlying distribution of word counts is given through modeling the DNA sequence as a Markov chain (MC). For single long sequences, efficient statistics are available to estimate the order of MCs and the transition probability matrix for the sequences. As NGS data do not provide a single long sequence, inference methods on Markovian properties of sequences based on single long sequences cannot be directly used for NGS short read data. RESULTS: Here we derive a normal approximation for such word counts. We also show that the traditional Chi-square statistic has an approximate gamma distribution ,: using the Lander-Waterman model for physical mapping. We propose several methods to estimate the order of the MC based on NGS reads and evaluate those using simulations. We illustrate the applications of our results by clustering genomic sequences of several vertebrate and tree species based on NGS reads using alignment-free sequence dissimilarity measures. We find that the estimated order of the MC has a considerable effect on the clustering results ,: and that the clustering results that use a N: MC of the estimated order give a plausible clustering of the species. AVAILABILITY AND IMPLEMENTATION: Our implementation of the statistics developed here is available as R package 'NGS.MC' at http://www-rcf.usc.edu/∼fsun/Programs/NGS-MC/NGS-MC.html CONTACT: fsun@usc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Patterns of genomic diversification reflect differences in life history and reproductive biology between figs (Ficus) and the stone oaks (Lithocarpus).

One of the remarkable aspects of the tremendous biodiversity found in tropical forests is the wide range of evolutionary strategies that have produced this diversity, indicating many paths to diversification. We compare two diverse groups of trees with profoundly different biologies to discover whether these differences are reflected in their genomes. Ficus (Moraceae), with its complex co-evolutionary relationship with obligate pollinating wasps, produces copious tiny seeds that are widely dispersed. Lithocarpus (Fagaceae), with generalized insect pollination, produces large seeds that are poorly dispersed. We hypothesize that these different reproductive biologies and life history strategies should have a profound impact on the basic properties of genomic divergence within each genus. Using shallow whole genome sequencing for six species of Ficus, seven species of Lithocarpus, and three outgroups, we examined overall genomic diversity, how it is shared among the species within each genus, and the fraction of this shared diversity that agrees with the major phylogenetic pattern. A substantially larger fraction of the genome is shared among species of Lithocarpus, a considerable amount of this shared diversity was incongruent with the general background history of the genomes, and each fig species possessed a substantially larger fraction of unique diversity than Lithocarpus.

Exploring the potential of gametic reconstruction of parental genotypes by F1 hybrids as a bridge for rapid introgression.

Interspecific hybridization and genetic introgression are commonly observed in natural populations of many species, especially trees. Among oaks, gene flow between closely related species has been well documented. And yet, hybridization does not lead to a "melting pot", i.e., the homogenization of phenotypic traits. Here, we explore how the combination of several common reproductive and genomic traits could create an avenue for interspecific gene flow that partially explains this apparent paradox. During meiosis, F1 hybrids will produce approximately (½)n "reconstructed" parental gametes, where n equals the number of chromosomes. Crossing over would introduce a small amount of introgressive material. The resulting parental-type gametophytes would probably possess a similar fertilization advantage as conspecific pollen. The resulting "backcross" would actually be the genetic equivalent of a conspecific out-cross, with a small amount of heterospecific DNA captured through crossing over. Even with detailed genomic analysis, the resulting offspring would not appear to be a backcross. This avenue for rapid introgression between species through the F1 hybrid will be viable for organisms that meet certain conditions: low base chromosome number, conserved genomic structure and size, production of billions of gametes/gametophytes during each reproductive event, and conspecific fertilization advantage.

Historical distribution of Sundaland's Dipterocarp rainforests at Quaternary glacial maxima.

The extent of Dipterocarp rainforests on the emergent Sundaland landmass in Southeast Asia during Quaternary glaciations remains a key question. A better understanding of the biogeographic history of Sundaland could help explain current patterns of biodiversity and support the development of effective forest conservation strategies. Dipterocarpaceae trees dominate the rainforests of Sundaland, and their distributions serve as a proxy for rainforest extent. We used species distribution models (SDMs) of 317 Dipterocarp species to estimate the geographic extent of appropriate climatic conditions for rainforest on Sundaland at the last glacial maximum (LGM). The SDMs suggest that the climate of central Sundaland at the LGM was suitable to sustain Dipterocarp rainforest, and that the presence of a previously suggested transequatorial savannah corridor at that time is unlikely. Our findings are supported by palynologic evidence, dynamic vegetation models, extant mammal and termite communities, vascular plant fatty acid stable isotopic compositions, and stable carbon isotopic compositions of cave guano profiles. Although Dipterocarp species richness was generally lower at the LGM, areas of high species richness were mostly found off the current islands and on the emergent Sunda Shelf, indicating substantial species migration and mixing during the transitions between the Quaternary glacial maxima and warm periods such as the present.

Pollinator sharing and gene flow among closely related sympatric dioecious fig taxa.

Hybridization and insect pollination are widely believed to increase rates of plant diversification. The extreme diversity of figs (Ficus) and their obligate pollinators, fig wasps (Agaonidae), provides an opportunity to examine the possible role of pollinator-mediated hybridization in plant diversification. Increasing evidence suggests that pollinator sharing and hybridization occurs among fig taxa, despite relatively strict coevolution with the pollinating wasp. Using five sympatric dioecious fig taxa and their pollinators, we examine the degree of pollinator sharing and inter-taxa gene flow. We experimentally test pollinator preference for floral volatiles, the main host recognition signal, from different figs. All five fig taxa shared pollinators with other taxa, and gene flow occurred between fig taxa within and between sections. Floral volatiles of each taxon attracted more than one pollinator species. Floral volatiles were more similar between closely related figs, which experienced higher levels of pollinator sharing and inter-taxa gene flow. This study demonstrates that pollinator sharing and inter-taxa gene flow occurs among closely related sympatric dioecious fig taxa and that pollinators choose the floral volatiles of multiple fig taxa. The implications of pollinator sharing and inter-taxa gene flow on diversification, occurring even in this highly specialized obligate pollination system, require further study.

An assembly and alignment-free method of phylogeny reconstruction from next-generation sequencing data.

BACKGROUND: Next-generation sequencing technologies are rapidly generating whole-genome datasets for an increasing number of organisms. However, phylogenetic reconstruction of genomic data remains difficult because de novo assembly for non-model genomes and multi-genome alignment are challenging. RESULTS: To greatly simplify the analysis, we present an Assembly and Alignment-Free (AAF) method ( https://sourceforge.net/projects/aaf-phylogeny ) that constructs phylogenies directly from unassembled genome sequence data, bypassing both genome assembly and alignment. Using mathematical calculations, models of sequence evolution, and simulated sequencing of published genomes, we address both evolutionary and sampling issues caused by direct reconstruction, including homoplasy, sequencing errors, and incomplete sequencing coverage. From these results, we calculate the statistical properties of the pairwise distances between genomes, allowing us to optimize parameter selection and perform bootstrapping. As a test case with real data, we successfully reconstructed the phylogeny of 12 mammals using raw sequencing reads. We also applied AAF to 21 tropical tree genome datasets with low coverage to demonstrate its effectiveness on non-model organisms. CONCLUSION: Our AAF method opens up phylogenomics for species without an appropriate reference genome or high sequence coverage, and rapidly creates a phylogenetic framework for further analysis of genome structure and diversity among non-model organisms.

Borneo and Indochina are major evolutionary hotspots for Southeast Asian biodiversity.

Tropical Southeast (SE) Asia harbors extraordinary species richness and in its entirety comprises four of the Earth's 34 biodiversity hotspots. Here, we examine the assembly of the SE Asian biota through time and space. We conduct meta-analyses of geological, climatic, and biological (including 61 phylogenetic) data sets to test which areas have been the sources of long-term biological diversity in SE Asia, particularly in the pre-Miocene, Miocene, and Plio-Pleistocene, and whether the respective biota have been dominated by in situ diversification, immigration and/or emigration, or equilibrium dynamics. We identify Borneo and Indochina, in particular, as major "evolutionary hotspots" for a diverse range of fauna and flora. Although most of the region's biodiversity is a result of both the accumulation of immigrants and in situ diversification, within-area diversification and subsequent emigration have been the predominant signals characterizing Indochina and Borneo's biota since at least the early Miocene. In contrast, colonization events are comparatively rare from younger volcanically active emergent islands such as Java, which show increased levels of immigration events. Few dispersal events were observed across the major biogeographic barrier of Wallace's Line. Accelerated efforts to conserve Borneo's flora and fauna in particular, currently housing the highest levels of SE Asian plant and mammal species richness, are critically required.

Deforestation and fragmentation of natural forests in the upper Changhua watershed, Hainan, China: implications for biodiversity conservation.

Hainan, the largest tropical island in China, belongs to the Indo-Burma biodiversity hotspot. The Changhua watershed is a center of endemism for plants and birds and the cradle of Hainan's main rivers. However, this area has experienced recent and ongoing deforestation and habitat fragmentation. To quantify habitat loss and fragmentation of natural forests, as well as the land-cover changes in the Changhua watershed, we analyzed Landsat images obtained in 1988, 1995, and 2005. Land-cover dynamics analysis showed that natural forests increased in area (97,909 to 104,023 ha) from 1988 to 1995 but decreased rapidly to 76,306 ha over the next decade. Rubber plantations increased steadily throughout the study period while pulp plantations rapidly expanded after 1995. Similar patterns of land cover change were observed in protected areas, indicating a lack of enforcement. Natural forests conversion to rubber and pulp plantations has a general negative effect on biodiversity, primarily through habitat fragmentation. The fragmentation analysis showed that natural forests area was reduced and patch number increased, while patch size and connectivity decreased. These land-cover changes threatened local biodiversity, especially island endemic species. Both natural forests losses and fragmentation should be stopped by strict enforcement to prevent further damage. Preserving the remaining natural forests and enforcing the status of protected areas should be a management priority to maximize the watershed's biodiversity conservation value.

Soils on exposed Sunda shelf shaped biogeographic patterns in the equatorial forests of Southeast Asia.

The marked biogeographic difference between western (Malay Peninsula and Sumatra) and eastern (Borneo) Sundaland is surprising given the long time that these areas have formed a single landmass. A dispersal barrier in the form of a dry savanna corridor during glacial maxima has been proposed to explain this disparity. However, the short duration of these dry savanna conditions make it an unlikely sole cause for the biogeographic pattern. An additional explanation might be related to the coarse sandy soils of central Sundaland. To test these two nonexclusive hypotheses, we performed a floristic cluster analysis based on 111 tree inventories from Peninsular Malaysia, Sumatra, and Borneo. We then identified the indicator genera for clusters that crossed the central Sundaland biogeographic boundary and those that did not cross and tested whether drought and coarse-soil tolerance of the indicator genera differed between them. We found 11 terminal floristic clusters, 10 occurring in Borneo, 5 in Sumatra, and 3 in Peninsular Malaysia. Indicator taxa of clusters that occurred across Sundaland had significantly higher coarse-soil tolerance than did those from clusters that occurred east or west of central Sundaland. For drought tolerance, no such pattern was detected. These results strongly suggest that exposed sandy sea-bed soils acted as a dispersal barrier in central Sundaland. However, we could not confirm the presence of a savanna corridor. This finding makes it clear that proposed biogeographic explanations for plant and animal distributions within Sundaland, including possible migration routes for early humans, need to be reevaluated.

Pseudo-Sanger sequencing: massively parallel production of long and near error-free reads using NGS technology.

BACKGROUND: Usually, next generation sequencing (NGS) technology has the property of ultra-high throughput but the read length is remarkably short compared to conventional Sanger sequencing. Paired-end NGS could computationally extend the read length but with a lot of practical inconvenience because of the inherent gaps. Now that Illumina paired-end sequencing has the ability of read both ends from 600 bp or even 800 bp DNA fragments, how to fill in the gaps between paired ends to produce accurate long reads is intriguing but challenging. RESULTS: We have developed a new technology, referred to as pseudo-Sanger (PS) sequencing. It tries to fill in the gaps between paired ends and could generate near error-free sequences equivalent to the conventional Sanger reads in length but with the high throughput of the Next Generation Sequencing. The major novelty of PS method lies on that the gap filling is based on local assembly of paired-end reads which have overlaps with at either end. Thus, we are able to fill in the gaps in repetitive genomic region correctly. The PS sequencing starts with short reads from NGS platforms, using a series of paired-end libraries of stepwise decreasing insert sizes. A computational method is introduced to transform these special paired-end reads into long and near error-free PS sequences, which correspond in length to those with the largest insert sizes. The PS construction has 3 advantages over untransformed reads: gap filling, error correction and heterozygote tolerance. Among the many applications of the PS construction is de novo genome assembly, which we tested in this study. Assembly of PS reads from a non-isogenic strain of Drosophila melanogaster yields an N50 contig of 190 kb, a 5 fold improvement over the existing de novo assembly methods and a 3 fold advantage over the assembly of long reads from 454 sequencing. CONCLUSIONS: Our method generated near error-free long reads from NGS paired-end sequencing. We demonstrated that de novo assembly could benefit a lot from these Sanger-like reads. Besides, the characteristic of the long reads could be applied to such applications as structural variations detection and metagenomics.

Virtually the Same? Evaluating the Effectiveness of Remote Undergraduate Research Experiences.

In-person undergraduate research experiences (UREs) promote students' integration into careers in life science research. In 2020, the COVID-19 pandemic prompted institutions hosting summer URE programs to offer them remotely, raising questions about whether undergraduates who participate in remote research can experience scientific integration and whether they might perceive doing research less favorably (i.e., not beneficial or too costly). To address these questions, we examined indicators of scientific integration and perceptions of the benefits and costs of doing research among students who participated in remote life science URE programs in Summer 2020. We found that students experienced gains in scientific self-efficacy pre- to post-URE, similar to results reported for in-person UREs. We also found that students experienced gains in scientific identity, graduate and career intentions, and perceptions of the benefits of doing research only if they started their remote UREs at lower levels on these variables. Collectively, students did not change in their perceptions of the costs of doing research despite the challenges of working remotely. Yet students who started with low cost perceptions increased in these perceptions. These findings indicate that remote UREs can support students' self-efficacy development, but may otherwise be limited in their potential to promote scientific integration.

Exploiting sparseness in de novo genome assembly.

BACKGROUND: The very large memory requirements for the construction of assembly graphs for de novo genome assembly limit current algorithms to super-computing environments. METHODS: In this paper, we demonstrate that constructing a sparse assembly graph which stores only a small fraction of the observed k-mers as nodes and the links between these nodes allows the de novo assembly of even moderately-sized genomes (~500 M) on a typical laptop computer. RESULTS: We implement this sparse graph concept in a proof-of-principle software package, SparseAssembler, utilizing a new sparse k-mer graph structure evolved from the de Bruijn graph. We test our SparseAssembler with both simulated and real data, achieving ~90% memory savings and retaining high assembly accuracy, without sacrificing speed in comparison to existing de novo assemblers.

Development of high-throughput SNP-based genotyping in Acacia auriculiformis x A. mangium hybrids using short-read transcriptome data.

BACKGROUND: Next Generation Sequencing has provided comprehensive, affordable and high-throughput DNA sequences for Single Nucleotide Polymorphism (SNP) discovery in Acacia auriculiformis and Acacia mangium. Like other non-model species, SNP detection and genotyping in Acacia are challenging due to lack of genome sequences. The main objective of this study is to develop the first high-throughput SNP genotyping assay for linkage map construction of A. auriculiformis x A. mangium hybrids. RESULTS: We identified a total of 37,786 putative SNPs by aligning short read transcriptome data from four parents of two Acacia hybrid mapping populations using Bowtie against 7,839 de novo transcriptome contigs. Given a set of 10 validated SNPs from two lignin genes, our in silico SNP detection approach is highly accurate (100%) compared to the traditional in vitro approach (44%). Further validation of 96 SNPs using Illumina GoldenGate Assay gave an overall assay success rate of 89.6% and conversion rate of 37.5%. We explored possible factors lowering assay success rate by predicting exon-intron boundaries and paralogous genes of Acacia contigs using Medicago truncatula genome as reference. This assessment revealed that presence of exon-intron boundary is the main cause (50%) of assay failure. Subsequent SNPs filtering and improved assay design resulted in assay success and conversion rate of 92.4% and 57.4%, respectively based on 768 SNPs genotyping. Analysis of clustering patterns revealed that 27.6% of the assays were not reproducible and flanking sequence might play a role in determining cluster compression. In addition, we identified a total of 258 and 319 polymorphic SNPs in A. auriculiformis and A. mangium natural germplasms, respectively. CONCLUSION: We have successfully discovered a large number of SNP markers in A. auriculiformis x A. mangium hybrids using next generation transcriptome sequencing. By using a reference genome from the most closely related species, we converted most SNPs to successful assays. We also demonstrated that Illumina GoldenGate genotyping together with manual clustering can provide high quality genotypes for a non-model species like Acacia. These SNPs markers are not only important for linkage map construction, but will be very useful for hybrid discrimination and genetic diversity assessment of natural germplasms in the future.

The current refugial rainforests of Sundaland are unrepresentative of their biogeographic past and highly vulnerable to disturbance.

Understanding the historical dynamics of forest communities is a critical element for accurate prediction of their response to future change. Here, we examine evergreen rainforest distribution in the Sunda Shelf region at the last glacial maximum (LGM), using a spatially explicit model incorporating geographic, paleoclimatic, and geologic evidence. Results indicate that at the LGM, Sundaland rainforests covered a substantially larger area than currently present. Extrapolation of the model over the past million years demonstrates that the current "island archipelago" setting in Sundaland is extremely unusual given the majority of its history and the dramatic biogeographic transitions caused by global deglaciation were rapid and brief. Compared with dominant glacial conditions, lowland forests were probably reduced from approximately 1.3 to 0.8 x 10(6) km(2) while upland forests were probably reduced by half, from approximately 2.0 to 1.0 x 10(5) km(2). Coastal mangrove and swamp forests experienced the most dramatic change during deglaciations, going through a complete and major biogeographic relocation. The Sundaland forest dynamics of fragmentation and contraction and subsequent expansion, driven by glacial cycles, occur in the opposite phase as those in the northern hemisphere and equatorial Africa, indicating that Sundaland evergreen rainforest communities are currently in a refugial stage. Widespread human-mediated reduction and conversion of these forests in their refugial stage, when most species are passing through significant population bottlenecks, strongly emphasizes the urgency of conservation and management efforts. Further research into the natural process of fragmentation and contraction during deglaciation is necessary to understand the long-term effect of human activity on forest species.