Population genomics and conservation management of a declining tropical rodent

Conservation management is improved by incorporating information about the spatial distribution of population genetic diversity into planning strategies. Northern Australia is the location of some of the world’s most severe ongoing declines of endemic mammal species, yet we have little genetic information from this regional mammal assemblage to inform a genetic perspective on conservation assessment and planning. We used next-generation sequencing data from remnant populations of the threatened brush-tailed rabbit-rat (Conilurus penicillatus) to compare patterns of genomic diversity and differentiation across the landscape and investigate standardised hierarchical genomic diversity metrics to better understand brush-tailed rabbit-rat population genomic structure. We found strong population structuring, with high levels of differentiation between populations (F_ST = 0.21–0.78). Two distinct genomic lineages between the Tiwi Islands and mainland are also present. Prioritisation analysis showed that one population in both lineages would need to be conserved to retain at least ~80% of alleles for the species. Analysis of standardised genomic diversity metrics showed that approximately half of the total diversity occurs among lineages (δ = 0.091 from grand total γ = 0.184). We suggest that a focus on conserving remnant island populations may not be appropriate for the preservation of species-level genomic diversity and adaptive potential, as these populations represent a small component of the total diversity and a narrow subset of the environmental conditions in which the species occurs. We also highlight the importance of considering both genomic and ecological differentiation between source and receiving populations when considering translocations for conservation purposes.Subject terms: Ecological genetics, Population genetics, Conservation biology, Biogeography     

An environmental gradient dominates ecological and genetic differentiation of marine invertebrates between the North and Baltic Sea

Environmental gradients have emerged as important barriers to structuring populations and species distributions. We set out to test whether the strong salinity gradient from the marine North Sea to the brackish Baltic Sea in northern Europe represents an ecological and genetic break, and to identify life history traits that correlate with the strength of this break. We accumulated mitochondrial cytochrome oxidase subunit 1 sequence data, and data on the distribution, salinity tolerance, and life history for 28 species belonging to the Cnidaria, Crustacea, Echinodermata, Mollusca, Polychaeta, and Gastrotricha. We included seven non‐native species covering a broad range of times since introduction, in order to gain insight into the pace of adaptation and differentiation. We calculated measures of genetic diversity and differentiation across the environmental gradient, coalescent times, and migration rates between North and Baltic Sea populations, and analyzed correlations between genetic and life history data. The majority of investigated species is either genetically differentiated and/or adapted to the lower salinity conditions of the Baltic Sea. Species exhibiting population structure have a range of patterns of genetic diversity in comparison with the North Sea, from lower in the Baltic Sea to higher in the Baltic Sea, or equally diverse in North and Baltic Sea. Two of the non‐native species showed signs of genetic differentiation, their times since introduction to the Baltic Sea being about 80 and >700 years, respectively. Our results indicate that the transition from North Sea to Baltic Sea represents a genetic and ecological break: The diversity of genetic patterns points toward independent trajectories in the Baltic compared with the North Sea, and ecological differences with regard to salinity tolerance are common. The North Sea–Baltic Sea region provides a unique setting to study evolutionary adaptation during colonization processes at different stages by jointly considering native and non‐native species.     

Genomic and phenotypic differentiation among Methanosarcina mazei populations from Columbia River sediment

Methanogenic archaea are genotypically and phenotypically diverse organisms that are integral to carbon cycling in anaerobic environments. Owing to their genetic tractability and ability to be readily cultivated, Methanosarcina spp. have become a powerful model system for understanding methanogen biology at the cellular systems level. However, relatively little is known of how genotypic and phenotypic variation is partitioned in Methanosarcina populations inhabiting natural environments and the possible ecological and evolutionary implications of such variation. Here, we have identified how genomic and phenotypic diversity is partitioned within and between Methanosarcina mazei populations obtained from two different sediment environments in the Columbia River Estuary (Oregon, USA). Population genomic analysis of 56 M. mazei isolates averaging <1% nucleotide divergence revealed two distinct clades, which we refer to as ‘mazei-T'' and ‘mazei-WC''. Genomic analyses showed that these clades differed in gene content and fixation of allelic variants, which point to potential differences in primary metabolism and also interactions with foreign genetic elements. This hypothesis of niche partitioning was supported by laboratory growth experiments that revealed significant differences in trimethylamine utilization. These findings improve our understanding of the ecologically relevant scales of genomic variation in natural systems and demonstrate interactions between genetic and ecological diversity in these easily cultivable and genetically tractable model methanogens.     

A Population Genomic Assessment of Three Decades of Evolution in a Natural Drosophila Population

Population genetics seeks to illuminate the forces shaping genetic variation, often based on a single snapshot of genomic variation. However, utilizing multiple sampling times to study changes in allele frequencies can help clarify the relative roles of neutral and non-neutral forces on short time scales. This study compares whole-genome sequence variation of recently collected natural population samples of Drosophila melanogaster against a collection made approximately 35 years prior from the same locality—encompassing roughly 500 generations of evolution. The allele frequency changes between these time points would suggest a relatively small local effective population size on the order of 10,000, significantly smaller than the global effective population size of the species. Some loci display stronger allele frequency changes than would be expected anywhere in the genome under neutrality—most notably the tandem paralogs Cyp6a17 and Cyp6a23, which are impacted by structural variation associated with resistance to pyrethroid insecticides. We find a genome-wide excess of outliers for high genetic differentiation between old and new samples, but a larger number of adaptation targets may have affected SNP-level differentiation versus window differentiation. We also find evidence for strengthening latitudinal allele frequency clines: northern-associated alleles have increased in frequency by an average of nearly 2.5% at SNPs previously identified as clinal outliers, but no such pattern is observed at random SNPs. This project underscores the scientific potential of using multiple sampling time points to investigate how evolution operates in natural populations, by quantifying how genetic variation has changed over ecologically relevant timescales.     

Recurrent selection explains parallel evolution of genomic regions of high relative but low absolute differentiation in a ring species

Recent technological developments allow investigation of the repeatability of evolution at the genomic level. Such investigation is particularly powerful when applied to a ring species, in which spatial variation represents changes during the evolution of two species from one. We examined genomic variation among three subspecies of the greenish warbler ring species, using genotypes at 13 013 950 nucleotide sites along a new greenish warbler consensus genome assembly. Genomic regions of low within‐group variation are remarkably consistent between the three populations. These regions show high relative differentiation but low absolute differentiation between populations. Comparisons with outgroup species show the locations of these peaks of relative differentiation are not well explained by phylogenetically conserved variation in recombination rates or selection. These patterns are consistent with a model in which selection in an ancestral form has reduced variation at some parts of the genome, and those same regions experience recurrent selection that subsequently reduces variation within each subspecies. The degree of heterogeneity in nucleotide diversity is greater than explained by models of background selection, but is consistent with selective sweeps. Given the evidence that greenish warblers have had both population differentiation for a long period of time and periods of gene flow between those populations, we propose that some genomic regions underwent selective sweeps over a broad geographic area followed by within‐population selection‐induced reductions in variation. An important implication of this ‘sweep‐before‐differentiation’ model is that genomic regions of high relative differentiation may have moved among populations more recently than other genomic regions.     

Fine-scale comparative phylogeography of a sympatric sister species triplet of subterranean diving beetles from a single calcrete aquifer in Western Australia

Calcrete aquifers in the arid Yilgarn region of central Western Australia are a biodiversity hotspot for stygofauna. A distinct pattern of interspecific size class variation among subterranean dytiscid beetle species has been observed in 29 of these aquifers where either two or three small, medium and/or large sympatric species are found that are in some cases sister species. We used a 3.5 km² grid of bores to sample dytiscids on a fine-scale and employed a comparative phylogeographical and population genetic approach to investigate the origins of a sympatric sister species triplet of diving beetles from a single aquifer. Mitochondrial DNA sequence data from the Cytochrome oxidase c subunit I gene revealed that all three species have high levels of haplotype diversity with ancient (∼1 million years ago) intra-specific coalescence of haplotypes, but low levels of nucleotide diversity. Population analyses provide evidence for multiple expansion events within each species. There was spatial heterogeneity in the distribution of genetic variation and abundance both within and among the three taxa. Population analyses revealed significant fine-scale differentiation with isolation by distance for Paroster macrosturtensis and P. mesosturtensis , but not the smallest species P. microsturtensis . Haplotype network analyses provided limited or no evidence for past population fragmentation within the large and small species, but substantial historical divergence was observed in P. mesosturtensis that was not spatially structured. A patchy population structure with contemporaneous and historical isolation by distance in the three species is likely to have been a significant isolating and diversifying force, preventing us from ruling out a potential role for allopatric divergence during speciation of this beetle sister triplet.     

The population genomic signature of environmental selection in the widespread insect-pollinated tree species Frangula alnus at different geographical scales

The evaluation of the molecular signatures of selection in species lacking an available closely related reference genome remains challenging, yet it may provide valuable fundamental insights into the capacity of populations to respond to environmental cues. We screened 25 native populations of the tree species Frangula alnus subsp. alnus (Rhamnaceae), covering three different geographical scales, for 183 annotated single-nucleotide polymorphisms (SNPs). Standard population genomic outlier screens were combined with individual-based and multivariate landscape genomic approaches to examine the strength of selection relative to neutral processes in shaping genomic variation, and to identify the main environmental agents driving selection. Our results demonstrate a more distinct signature of selection with increasing geographical distance, as indicated by the proportion of SNPs (i) showing exceptional patterns of genetic diversity and differentiation (outliers) and (ii) associated with climate. Both temperature and precipitation have an important role as selective agents in shaping adaptive genomic differentiation in F. alnus subsp. alnus, although their relative importance differed among spatial scales. At the ‘intermediate'' and ‘regional'' scales, where limited genetic clustering and high population diversity were observed, some indications of natural selection may suggest a major role for gene flow in safeguarding adaptability. High genetic diversity at loci under selection in particular, indicated considerable adaptive potential, which may nevertheless be compromised by the combined effects of climate change and habitat fragmentation.     

Genomic Organization of Repetitive DNA Elements and Its Implications for the Chromosomal Evolution of Channid Fishes (Actinopterygii,Perciformes)

Channid fishes, commonly referred to as “snakeheads”, are currently very important in Asian fishery and aquaculture due to the substantial decline in natural populations because of overexploitation. A large degree of chromosomal variation has been found in this family, mainly through the use of conventional cytogenetic investigations. In this study, we analyzed the karyotype structure and the distribution of 7 repetitive DNA sequences in several Channa species from different Thailand river basins. The aim of this study was to investigate the chromosomal differentiation among species and populations to improve upon the knowledge of its biodiversity and evolutionary history. Rearrangements, such as pericentric inversions, fusions and polyploidization, appear to be important events during the karyotypic evolution of this genus, resulting in the chromosomal diversity observed among the distinct species and even among populations of the same species. In addition, such variability is also increased by the genomic dynamism of repetitive elements, particularly by the differential distribution and accumulation of rDNA sequences on chromosomes. This marked diversity is likely linked to the lifestyle of the snakehead fishes and their population fragmentation, as already identified for other fish species. The karyotypic features highlight the biodiversity of the channid fishes and justify a taxonomic revision of the genus Channa, as well as of the Channidae family as a whole, as some nominal species may actually constitute species complexes.     

Environmental metabarcoding reveals heterogeneous drivers of microbial eukaryote diversity in contrasting estuarine ecosystems

Assessing how natural environmental drivers affect biodiversity underpins our understanding of the relationships between complex biotic and ecological factors in natural ecosystems. Of all ecosystems, anthropogenically important estuaries represent a ‘melting pot'' of environmental stressors, typified by extreme salinity variations and associated biological complexity. Although existing models attempt to predict macroorganismal diversity over estuarine salinity gradients, attempts to model microbial biodiversity are limited for eukaryotes. Although diatoms commonly feature as bioindicator species, additional microbial eukaryotes represent a huge resource for assessing ecosystem health. Of these, meiofaunal communities may represent the optimal compromise between functional diversity that can be assessed using morphology and phenotype–environment interactions as compared with smaller life fractions. Here, using 454 Roche sequencing of the 18S nSSU barcode we investigate which of the local natural drivers are most strongly associated with microbial metazoan and sampled protist diversity across the full salinity gradient of the estuarine ecosystem. In order to investigate potential variation at the ecosystem scale, we compare two geographically proximate estuaries (Thames and Mersey, UK) with contrasting histories of anthropogenic stress. The data show that although community turnover is likely to be predictable, taxa are likely to respond to different environmental drivers and, in particular, hydrodynamics, salinity range and granulometry, according to varied life-history characteristics. At the ecosystem level, communities exhibited patterns of estuary-specific similarity within different salinity range habitats, highlighting the environmental sequencing biomonitoring potential of meiofauna, dispersal effects or both.     

Genomic Consequences of Long-Term Population Decline in Brown Eared Pheasant

Population genetic theory and empirical evidence indicate that deleterious alleles can be purged in small populations. However, this viewpoint remains controversial. It is unclear whether natural selection is powerful enough to purge deleterious mutations when wild populations continue to decline. Pheasants are terrestrial birds facing a long-term risk of extinction as a result of anthropogenic perturbations and exploitation. Nevertheless, there are scant genomics resources available for conservation management and planning. Here, we analyzed comparative population genomic data for the three extant isolated populations of Brown eared pheasant (Crossoptilon mantchuricum) in China. We showed that C. mantchuricum has low genome-wide diversity and a contracting effective population size because of persistent declines over the past 100,000 years. We compared genome-wide variation in C. mantchuricum with that of its closely related sister species, the Blue eared pheasant (C. auritum) for which the conservation concern is low. There were detrimental genetic consequences across all C. mantchuricum genomes including extended runs of homozygous sequences, slow rates of linkage disequilibrium decay, excessive loss-of-function mutations, and loss of adaptive genetic diversity at the major histocompatibility complex region. To the best of our knowledge, this study is the first to perform a comprehensive conservation genomic analysis on this threatened pheasant species. Moreover, we demonstrated that natural selection may not suffice to purge deleterious mutations in wild populations undergoing long-term decline. The findings of this study could facilitate conservation planning for threatened species and help recover their population size.     

What does population structure analysis reveal about the Pterostylis longifolia complex (Orchidaceae)?

Morphologically similar groups of species are common and pose significant challenges for taxonomists. Differences in approaches to classifying unique species can result in some species being overlooked, whereas others are wrongly conserved. The genetic diversity and population structure of the Pterostylis longifolia complex (Orchidaceae) in Tasmania was investigated to determine if four species, and potential hybrids, could be distinguished through genomic AFLP and chloroplast restriction‐fragment‐length polymorphism (RFLP) markers. Analysis of molecular variance (AMOVA) results indicated that little genetic variation was present among taxa, whereas PCoA analyses revealed genetic variation at a regional scale irrespective of taxa. Population genetic structure analyses identified three clusters that correspond to regional genetic and single taxon‐specific phenotypic variation. The results from this study suggest that “longifolia” species have persisted throughout the last glacial maximum in Tasmania and that the complex may be best treated as a single taxon with several morphotypes. These results could have serious evolutionary and conservation implications as taxonomic changes could result in the instatement of a single, widespread taxon in which rarer morphotypes are not protected.     

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome‐wide patterns of within‐ and between‐ species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole‐genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long‐term balancing selection have also been crucial components in shaping patterns of genome‐wide variation during the speciation process.     

Chloroplast diversity in Vouacapoua americana (Caesalpiniaceae), a neotropical forest tree

The chloroplast genome has been widely used to describe genetic diversity in plant species. Its maternal inheritance in numerous angiosperm species and low mutation rate are suitable characters when inferring historical events such as possible recolonization routes. Here we have studied chloroplast DNA variation using PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) with seven pairs of primers and four restriction enzymes in 14 populations of Vouacapoua americana (Caesalpiniaceae) a neotropical tree sampled throughout French Guiana. Population diversity (Hs), total gene diversity (Ht) and differentiation among populations (GST) were estimated using Nei's method as 0.09, 0.87 and 0.89, respectively. This is consistent with the limited gene flow associated with synzoochory in this species. The genetic structure observed in the north of French Guiana suggests that historical events such as contractions and recent recolonizations have had a large impact on the distribution of genetic diversity in this species.     

A meta-analysis revealed insights into the sources,conservation and impact of microRNA 5′-isoforms in four model species

MicroRNA (miRNA) 5′-isoforms, or 5′-isomiRs, are small-RNA species that originate from the same genomic loci as the major miRNAs with their 5′ ends shifted from the 5′ ends of the miRNAs by a few nucleotides. Although 5′-isomiRs have been reported, their origins, properties and potential functions remain to be examined. We systematically studied 5′-isomiRs in human, mouse, fruitfly and worm by analysing a large collection of small non-coding RNA and mRNA profiling data. The results revealed a broad existence of 5′-isomiRs in the four species, many of which were conserved and could arise from genomic loci of canonical and non-canonical miRNAs. The well-conserved 5′-isomiRs have several features, including a preference of the 3p over the 5p arms of hairpins of conserved mammalian miRNAs, altered 5′-isomiRs across species and across tissues, and association with structural variations of miRNA hairpins. Importantly, 5′-isomiRs and their major miRNAs may have different mRNA targets and thus potentially play distinct roles of gene regulation, as shown by an integrative analysis combining miRNA and mRNA profiling data from psoriatic and normal human skin and from murine miRNA knockout assays. Indeed, 18 5′-isomiRs had aberrant expression in psoriatic human skin, suggesting their potential function in psoriasis pathogenesis. The results of the current study deepened our understanding of the diversity and conservation of miRNAs, their plasticity in gene regulation and potential broad function in complex diseases.     

Data-driven identification of major axes of functional variation in bacteria

The discovery of major axes of correlated functional variation among species and habitats has revealed the fundamental trade-offs structuring both functional and taxonomic diversity in eukaryotes such as plants. Whether such functional axes exist in the bacterial realm and whether they could explain bacterial taxonomic turnover among ecosystems remains unknown. Here, we use a data-driven approach to leverage global genomic and metagenomic datasets to reveal the existence of major axes of functional variation explaining both evolutionary differentiation within Bacteria and their ecological sorting across diverse habitats. We show that metagenomic variation among bacterial communities from various ecosystems is structured along a few axes of correlated functional pathways. Similar clusters of traits explained phylogenetic trait variation among >16,000 bacterial genomes, suggesting that functional turnover among bacterial communities from distinct habitats does not only result from the differential filtering of similar functions among communities, but also from phylogenetic correlations among these functions. Concordantly, functional pathways associated with trait clusters that were most important for defining functional turnover among bacterial communities were also those that had the highest phylogenetic signal in the bacterial genomic phylogeny. This study overall underlines the important role of evolutionary history in shaping contemporary distributions of bacteria across ecosystems.     

A comparative hierarchical analysis of bacterioplankton and biofilm metacommunity structure in an interconnected pond system

It is unknown whether bacterioplankton and biofilm communities are structured by the same ecological processes, and whether they influence each other through continuous dispersal (known as mass effects). Using a hierarchical sampling approach we compared the relative importance of ecological processes structuring the dominant fraction (relative abundance ≥0.1%) of bacterioplankton and biofilm communities from three microhabitats (open water, Nuphar and Phragmites sites) at within‐ and among‐pond scale in a set of 14 interconnected shallow ponds. Our results demonstrate that while bacterioplankton and biofilm communities are highly distinct, a similar hierarchy of ecological processes is acting on them. For both community types, most variation in community composition was determined by pond identity and environmental variables, with no effect of space. The highest β‐diversity within each community type was observed among ponds, while microhabitat type (Nuphar, Phragmites, open water) significantly influenced biofilm communities but not bacterioplankton. Mass effects among bacterioplankton and biofilm communities were not detected, as suggested by the absence of within‐site covariation of biofilm and bacterioplankton communities. Both biofilm and plankton communities were thus highly structured by environmental factors (i.e., species sorting), with among‐lake variation being more important than within‐lake variation, whereas dispersal limitation and mass effects were not observed.     

Evolutionary Conservation Genomics Reveals Recent Speciation and Local Adaptation in Threatened Takins

Incorrect species delimitation will lead to inappropriate conservation decisions, especially for threatened species. The takin (Budorcas taxicolor) is a large artiodactyl endemic to the Himalayan–Hengduan–Qinling Mountains and is well known for its threatened status and peculiar appearance. However, the speciation, intraspecies taxonomy, evolutionary history, and adaptive evolution of this species still remain unclear, which greatly hampers its scientific conservation. Here, we de novo assembled a high-quality chromosome-level genome of takin and resequenced the genomes of 75 wild takins. Phylogenomics revealed that takin was positioned at the root of Caprinae. Population genomics based on the autosome, X chromosome, and Y chromosome SNPs and mitochondrial genomes consistently revealed the existence of two phylogenetic species and recent speciation in takins: the Himalayan takin (B. taxicolor) and the Chinese takin (B. tibetana), with the support of morphological evidence. Two genetically divergent subspecies were identified in both takin species, rejecting three previously proposed taxonomical viewpoints. Furthermore, their distribution boundaries were determined, suggesting that large rivers play important roles in shaping the genetic partition. Compared with the other subspecies, the Qinling subspecies presented the lowest genomic diversity, higher linkage disequilibrium, inbreeding, and genetic load, thus is in urgent need of genetic management and protection. Moreover, coat color gene (PMEL) variation may be responsible for the adaptive coat color difference between the two species following Gloger’s rule. Our findings provide novel insights into the recent speciation, local adaptation, scientific conservation of takins, and biogeography of the Himalaya–Hengduan biodiversity hotspot.     

Subspecies in the global human gut microbiome

Population genomics of prokaryotes has been studied in depth in only a small number of primarily pathogenic bacteria, as genome sequences of isolates of diverse origin are lacking for most species. Here, we conducted a large‐scale survey of population structure in prevalent human gut microbial species, sampled from their natural environment, with a culture‐independent metagenomic approach. We examined the variation landscape of 71 species in 2,144 human fecal metagenomes and found that in 44 of these, accounting for 72% of the total assigned microbial abundance, single‐nucleotide variation clearly indicates the existence of sub‐populations (here termed subspecies). A single subspecies (per species) usually dominates within each host, as expected from ecological theory. At the global scale, geographic distributions of subspecies differ between phyla, with Firmicutes subspecies being significantly more geographically restricted. To investigate the functional significance of the delineated subspecies, we identified genes that consistently distinguish them in a manner that is independent of reference genomes. We further associated these subspecies‐specific genes with properties of the microbial community and the host. For example, two of the three Eubacterium rectale subspecies consistently harbor an accessory pro‐inflammatory flagellum operon that is associated with lower gut community diversity, higher host BMI, and higher blood fasting insulin levels. Using an additional 676 human oral samples, we further demonstrate the existence of niche specialized subspecies in the different parts of the oral cavity. Taken together, we provide evidence for subspecies in the majority of abundant gut prokaryotes, leading to a better functional and ecological understanding of the human gut microbiome in conjunction with its host.