Reference-free population genomics from next-generation transcriptome data and the vertebrate-invertebrate gap

Development of an infiltration-based RNA preservation method for cryogen-free storage of leaves for gene expression analyses in field-grown plants

BACKGROUND

Gene expression is a fundamental process for plants to express their phenotype, and its analysis is the basis of molecular studies. However, the instability of RNA often poses an obstacle to analyzing plants grown in fields or remote locations where the availability of liquid nitrogen or dry ice is limited. To deepen our understanding of plant phenotypes and tolerance to field-specific stresses, it is crucial to develop methodologies to maintain plant RNA intact and safely transfer it for downstream analyses such as qPCR and RNA-seq.

RESULTS

In this study, the author developed a novel tissue preservation method that involved the infiltration of RNA preservation solution into the leaf apoplast using a syringe and subsequent storage at 4 °C. RNA-seq using samples stored for 5 d and principal component analyses showed that rice leaves treated with the infiltration method maintained the original transcriptome pattern better than those treated with the traditional method when the leaves were simply immersed in the solution. Additionally, it was also found that extracted RNA can be transported with minimum risk of degradation when it is bound to the membrane of RNA extraction kits. The developed infiltration method was applied to rice plants grown in a local farmer's field in northern Madagascar to analyze the expression of nutrient-responsive genes, suggesting nutrient imbalances in some of the fields examined.

CONCLUSIONS

This study showed that the developed infiltration method was effective in preserving the transcriptome status of rice and sorghum leaves when liquid nitrogen or a deep freezer is not available. The developed method was useful for diagnosing plants in the field based on the expression of nutrient-responsive marker genes. Moreover, the method used to protect RNA samples from degradation during transportation offers the possibility to use them for RNA-seq. This novel technique could pave the way for revealing the molecular basis of plant phenotypes by accelerating gene expression analyses using plant samples that are unique in the field.

Mating systems and recombination landscape strongly shape genetic diversity and selection in wheat relatives

How and why genetic diversity varies among species is a long-standing question in evolutionary biology. Life history traits have been shown to explain a large part of observed diversity. Among them, mating systems have one of the strongest impacts on genetic diversity, with selfing species usually exhibiting much lower diversity than outcrossing relatives. Theory predicts that a high rate of selfing amplifies selection at linked sites, reducing genetic diversity genome-wide, but frequent bottlenecks and rapid population turn-over could also explain low genetic diversity in selfers. However, how linked selection varies with mating systems and whether it is sufficient to explain the observed difference between selfers and outcrossers has never been tested. Here, we used the Aegilops/Triticum grass species, a group characterized by contrasted mating systems (from obligate outcrossing to high selfing) and marked recombination rate variation across the genome, to quantify the effects of mating system and linked selection on patterns of neutral and selected polymorphism. By analyzing phenotypic and transcriptomic data of 13 species, we show that selfing strongly affects genetic diversity and the efficacy of selection by amplifying the intensity of linked selection genome-wide. In particular, signatures of adaptation were only found in the highly recombining regions in outcrossing species. These results bear implications for the evolution of mating systems and, more generally, for our understanding of the fundamental drivers of genetic diversity.

Extensive immune receptor repertoire diversity in disease-resistant rice landraces

Plants have powerful defense mechanisms and extensive immune receptor repertoires, yet crop monocultures are prone to epidemic diseases. Rice (Oryza sativa) is susceptible to many diseases, such as rice blast caused by Magnaporthe oryzae. Varietal resistance of rice to blast relies on intracellular nucleotide binding, leucine-rich repeat (NLR) receptors that recognize specific pathogen molecules and trigger immune responses. In the Yuanyang terraces in southwest China, rice landraces rarely show severe losses to disease whereas commercial inbred lines show pronounced field susceptibility. Here, we investigate within-landrace NLR sequence diversity of nine rice landraces and eleven modern varieties using complexity reduction techniques. We find that NLRs display high sequence diversity in landraces, consistent with balancing selection, and that balancing selection at NLRs is more pervasive in landraces than modern varieties. Notably, modern varieties lack many ancient NLR haplotypes that are retained in some landraces. Our study emphasizes the value of standing genetic variation that is maintained in farmer landraces as a resource to make modern crops and agroecosystems less prone to disease. The conservation of landraces is, therefore, crucial for ensuring food security in the face of dynamic biotic and abiotic threats.

De novo transcriptome for Chiloscyllium griseum, a long-tail carpet shark of the Indian waters

Sharks have thrived in the oceans for 400 million years, experienced five extinctions and evolved into today's apex predators. However, enormous genome size, poor karyotyping and limited tissue sampling options are the bottlenecks in shark research. Sharks of the family Orectolobiformes act as model species in transcriptome research with exceptionally high reproductive fecundity, catch prominence and oviparity. The present study illustrates a de novo transcriptome for an adult grey bamboo shark, Chiloscyllium griseum (Chondrichthyes; Hemiscyllidae) using paired-end RNA sequencing. Around 150 million short Illumina reads were obtained from five different tissues and assembled using the Trinity assembler. 70,647 hits on Uniprot by BLASTX was obtained after the transcriptome annotation. The data generated serve as a basis for transcriptome-based population genetic studies and open up new avenues in the field of comparative transcriptomics and conservation biology.

The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Single-cell transcriptomics supports presence of cryptic species and reveals low levels of population genetic diversity in two testate amoebae morphospecies with large population sizes

The enormous population sizes and wide biogeographical distribution of many microbial eukaryotes set the expectation of high levels of intraspecific genetic variation. However, studies investigating protist populations remain scarce, mostly due to limited 'omics data. Instead, most genetics studies of microeukaryotes have thus far relied on single loci, which can be misleading and do not easily allow for detection of recombination, a hallmark of sexual reproduction. Here, we analyze >40 genes from 72 single-cell transcriptomes from two morphospecies-Hyalosphenia papilio and Hyalosphenia elegans-of testate amoebae (Arcellinida, Amoebozoa) to assess genetic diversity in samples collected over four years from New England bogs. We confirm the existence of cryptic species based on our multilocus dataset, which provides evidence of recombination within and high levels of divergence between the cryptic species. At the same time, total levels of genetic diversity within cryptic species are low, suggesting that these abundant organisms have small effective population sizes, perhaps due to extinction and repopulation events coupled with efficient modes of dispersal. This study is one of the first to investigate population genetics in uncultivable heterotrophic protists using transcriptomics data and contributes towards understanding cryptic species of nonmodel microeukaryotes.

A comprehensive framework for detecting copy number variants from single nucleotide polymorphism data: 'rCNV', a versatile r package for paralogue and CNV detection

Recent studies have highlighted the significant role of copy number variants (CNVs) in phenotypic diversity, environmental adaptation and species divergence across eukaryotes. The presence of CNVs also has the potential to introduce genotyping biases, which can pose challenges to accurate population and quantitative genetic analyses. However, detecting CNVs in genomes, particularly in non-model organisms, presents a formidable challenge. To address this issue, we have developed a statistical framework and an accompanying r software package that leverage allelic-read depth from single nucleotide polymorphism (SNP) data for accurate CNV detection. Our framework capitalises on two key principles. First, it exploits the distribution of allelic-read depth ratios in heterozygotes for individual SNPs by comparing it against an expected distribution based on binomial sampling. Second, it identifies SNPs exhibiting an apparent excess of heterozygotes under Hardy-Weinberg equilibrium. By employing multiple statistical tests, our method not only enhances sensitivity to sampling effects but also effectively addresses reference biases, resulting in optimised SNP classification. Our framework is compatible with various NGS technologies (e.g. RADseq, Exome-capture). This versatility enables CNV calling from genomes of diverse complexities. To streamline the analysis process, we have implemented our framework in the user-friendly r package 'rCNV', which automates the entire workflow seamlessly. We trained our models using simulated data and validated their performance on four datasets derived from different sequencing technologies, including RADseq (Chinook salmon-Oncorhynchus tshawytscha), Rapture (American lobster-Homarus americanus), Exome-capture (Norway spruce-Picea abies) and WGS (Malaria mosquito-Anopheles gambiae).

Cosegregation of recombinant chromatids maintains genome-wide heterozygosity in an asexual nematode

In asexual animals, female meiosis is modified to produce diploid oocytes. If meiosis still involves recombination, this is expected to lead to a rapid loss of heterozygosity, with adverse effects on fitness. Many asexuals, however, have a heterozygous genome, the underlying mechanisms being most often unknown. Cytological and population genomic analyses in the nematode Mesorhabditis belari revealed another case of recombining asexual being highly heterozygous genome-wide. We demonstrated that heterozygosity is maintained despite recombination because the recombinant chromatids of each chromosome pair cosegregate during the unique meiotic division. A theoretical model confirmed that this segregation bias is necessary to account for the observed pattern and likely to evolve under a wide range of conditions. Our study uncovers an unexpected type of non-Mendelian genetic inheritance involving cosegregation of recombinant chromatids.

Target enrichment of long open reading frames and ultraconserved elements to link microevolution and macroevolution in non-model organisms

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here, we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced substantial data sets. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites). Variant calling on ORFs and UCEs of Coelaturini from the Malawi Basin produced ~2000 SNPs per population pair. Estimates of nucleotide diversity and population differentiation were similar for ORFs and UCEs. They were low compared to previous estimates in molluscs, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming off-target sequence data from the same enriched libraries of Coelaturini from the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays the gene order inferred for the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for integrative genomic studies of microevolutionary and macroevolutionary dynamics in non-model organisms.

Genetic parallelism between European flat oyster populations at the edge of their natural range

Although all marine ecosystems have experienced global-scale losses, oyster reefs have shown the greatest. Therefore, substantial efforts have been dedicated to restoration of such ecosystems during the last two decades. In Europe, several pilot projects for the restoration of the native European flat oyster, Ostrea edulis, recently begun and recommendations to preserve genetic diversity and to conduct monitoring protocols have been made. In particular, an initial step is to test for genetic differentiation against homogeneity among the oyster populations potentially involved in such programs. Therefore, we conducted a new sampling of wild populations at the European scale and a new genetic analysis with 203 markers to (1) confirm and study in more detail the pattern of genetic differentiation between Atlantic and Mediterranean populations, (2) identify potential translocations that could be due to aquaculture practices and (3) investigate the populations at the fringe of the geographical range, since they seemed related despite their geographic distance. Such information should be useful to enlighten the choice of the animals to be translocated or reproduced in hatcheries for further restocking. After the confirmation of the general geographical pattern of genetic structure and the identification of one potential case of aquaculture transfer at a large scale, we were able to detect genomic islands of differentiation mainly in the form of two groups of linked markers, which could indicate the presence of polymorphic chromosomal rearrangements. Furthermore, we observed a tendency for these two islands and the most differentiated loci to show a parallel pattern of differentiation, grouping the North Sea populations with the Eastern Mediterranean and Black Sea populations, against geography. We discussed the hypothesis that this genetic parallelism could be the sign of a shared evolutionary history of the two groups of populations despite them being at the border of the distribution nowadays.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation

Longevity is a defining, heritable trait that varies dramatically between species. To resolve the genetic regulation of this trait, we have mined genomic variation in rockfishes, which range in longevity from 11 to over 205 years. Multiple shifts in rockfish longevity have occurred independently and in a short evolutionary time frame, thus empowering convergence analyses. Our analyses reveal a common network of genes under convergent evolution, encompassing established aging regulators such as insulin signaling, yet also identify flavonoid (aryl-hydrocarbon) metabolism as a pathway modulating longevity. The selective pressures on these pathways indicate the ancestral state of rockfishes was long lived and that the changes in short-lived lineages are adaptive. These pathways were also used to explore genome-wide association studies of human longevity, identifying the aryl-hydrocarbon metabolism pathway to be significantly associated with human survival to the 99th percentile. This evolutionary intersection defines and cross-validates a previously unappreciated genetic architecture that associates with the evolution of longevity across vertebrates.

Common microgeographical selection patterns revealed in four European conifers

Microgeographical adaptation occurs when the effects of directional selection persist despite gene flow. Traits and genetic loci under selection can then show adaptive divergence, against the backdrop of little differentiation at other traits or loci. How common such events are and how strong the selection is that underlies them remain open questions. Here, we discovered and analysed microgeographical patterns of genomic divergence in four European and Mediterranean conifers with widely differing life-history traits and ecological requirements (Abies alba MIll., Cedrus atlantica [Endl.] Manetti, Pinus halepensis Mill. and Pinus pinaster Aiton) by screening pairs from geographically close forest stands sampled along steep ecological gradients. We inferred patterns of genomic divergence by applying a combination of divergence outlier detection methods, demographic modelling, Approximate Bayesian Computation inferences and genomic annotation to genomic data. Surprisingly for such small geographical scales, we showed that selection is strong in all species but generally affects different loci in each. A clear signature of selection was systematically detected on a fraction of the genome, of the order of 0.1%-1% of the loci depending on the species. The novel modelling method we designed for estimating selection coefficients showed that the microgeographical selection coefficient scaled by population size (Ns) was 2-30. Our results convincingly suggest that selection maintains within-population diversity at microgeographical scales in spatially heterogeneous environments. Such genetic diversity is likely to be a major reservoir of adaptive potential, helping populations to adapt under fluctuating environmental conditions.

Dioecy and chromosomal sex determination are maintained through allopolyploid speciation in the plant genus Mercurialis

Polyploidization may precipitate dramatic changes to the genome, including chromosome rearrangements, gene loss, and changes in gene expression. In dioecious plants, the sex-determining mechanism may also be disrupted by polyploidization, with the potential evolution of hermaphroditism. However, while dioecy appears to have persisted through a ploidy transition in some species, it is unknown whether the newly formed polyploid maintained its sex-determining system uninterrupted, or whether dioecy re-evolved after a period of hermaphroditism. Here, we develop a bioinformatic pipeline using RNA-sequencing data from natural populations to demonstrate that the allopolyploid plant Mercurialis canariensis directly inherited its sex-determining region from one of its diploid progenitor species, M. annua, and likely remained dioecious through the transition. The sex-determining region of M. canariensis is smaller than that of its diploid progenitor, suggesting that the non-recombining region of M. annua expanded subsequent to the polyploid origin of M. canariensis. Homeologous pairs show partial sexual subfunctionalization. We discuss the possibility that gene duplicates created by polyploidization might contribute to resolving sexual antagonism.

Inter-Specific Genetic Exchange Despite Strong Divergence in Deep-Sea Hydrothermal Vent Gastropods of the Genus Alviniconcha

Deep hydrothermal vents are highly fragmented and unstable habitats at all temporal and spatial scales. Such environmental dynamics likely play a non-negligible role in speciation. Little is, however, known about the evolutionary processes that drive population-level differentiation and vent species isolation and, more specifically, how geography and habitat specialisation interplay in the species history of divergence. In this study, the species range and divergence of Alviniconcha snails that occupy active Western Pacific vent fields was assessed by using sequence variation data of the mitochondrial Cox1 gene, RNAseq, and ddRAD-seq. Combining morphological description and sequence datasets of the three species across five basins, we confirmed that A. kojimai, A. boucheti, and A. strummeri, while partially overlapping over their range, display high levels of divergence in the three genomic compartments analysed that usually encompass values retrieved for reproductively isolated species with divergences rang from 9% to 12.5% (mtDNA) and from 2% to 3.1% (nuDNA). Moreover, the three species can be distinguished on the basis of their external morphology by observing the distribution of bristles and the shape of the columella. According to this sampling, A. boucheti and A. kojimai form an east-to-west species abundance gradient, whereas A. strummeri is restricted to the Futuna Arc/Lau and North Fiji Basins. Surprisingly, population models with both gene flow and population size heterogeneities among genomes indicated that these three species are still able to exchange genes due to secondary contacts at some localities after a long period of isolation.

A population-level statistic for assessing Mendelian behavior of genotyping-by-sequencing data from highly duplicated genomes

Background

Given the economic and environmental importance of allopolyploids and other species with highly duplicated genomes, there is a need for methods to distinguish paralogs, i.e. duplicate sequences within a genome, from Mendelian loci, i.e. single copy sequences that pair at meiosis. The ratio of observed to expected heterozygosity is an effective tool for filtering loci but requires genotyping to be performed first at a high computational cost, whereas counting the number of sequence tags detected per genotype is computationally quick but very ineffective in inbred or polyploid populations. Therefore, new methods are needed for filtering paralogs.

Results

We introduce a novel statistic, H_ind/H_E, that uses the probability that two reads sampled from a genotype will belong to different alleles, instead of observed heterozygosity. The expected value of H_ind/H_E is the same across all loci in a dataset, regardless of read depth or allele frequency. In contrast to methods based on observed heterozygosity, it can be estimated and used for filtering loci prior to genotype calling. In addition to filtering paralogs, it can be used to filter loci with null alleles or high overdispersion, and identify individuals with unexpected ploidy and hybrid status. We demonstrate that the statistic is useful at read depths as low as five to 10, well below the depth needed for accurate genotype calling in polyploid and outcrossing species.

Conclusions

Our methodology for estimating H_ind/H_E across loci and individuals, as well as determining reasonable thresholds for filtering loci, is implemented in polyRAD v1.6, available at https://github.com/lvclark/polyRAD. In large sequencing datasets, we anticipate that the ability to filter markers and identify problematic individuals prior to genotype calling will save researchers considerable computational time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04635-9.

From genotype to phenotype: Genetic redundancy and the maintenance of an adaptive polymorphism in the context of high gene flow

A central question in evolution is how several adaptive phenotypes are maintained within a species. Theory predicts that the genetic determination of a trait, and in particular the amounts of redundancy in the mapping of genotypes to phenotypes, mediates evolutionary outcomes of phenotypic selection. In Mediterranean wild thyme, numerous discrete chemical phenotypes (chemotypes) occur in close geographic proximity. Chemotypes are defined by the predominant monoterpene produced by individual plants in their essential oil. In this study, we analyze the ecological genetics of six chemotypes nested within two well‐established chemical families (hereafter ecotypes). Ecotypes, and chemotypes within ecotypes, are spatially segregated, and their distributions track local differences in the abiotic environment. By combining population genomic, phenotypic, and environmental data from 700 individuals, we show how the genetics of ecotype determination mediates this evolutionary response. Variation in three terpene‐synthase loci explains variation in ecotype identity, with one single locus accounting for as much as 78% of this variation. Phenotypic selection combined with low segregating genotypic redundancy of ecotypes leaves a clear footprint at the genomic level: alleles associated with ecotype identity track environmental variation despite extensive gene flow. Different chemotypes within each ecotype differentially track environmental variation. Their identity is determined by multiple loci and displays a wider range of genotypic redundancy that dilutes phenotypic selection on their characteristic alleles. Our study thus provides a novel illustration of how genetic redundancy of a phenotype modulates the ability of selection to maintain adaptive differentiation. Identifying the precise genetics of the chemical polymorphism in thyme is the next crucial step for our understanding of the origin and maintenance of a polymorphism that is present in many aromatic plants.

A derived ZW chromosome system in Amborella trichopoda, representing the sister lineage to all other extant flowering plants

The genetic basis and evolution of sex determination in dioecious plants is emerging as an active area of research with exciting advances in genome sequencing and analysis technologies. As the sole species within the sister lineage to all other extant flowering plants, Amborella trichopoda is an important model for understanding the evolution and development of flowers. Plants typically produce only male or female flowers, but sex determination mechanisms are unknown for the species. Sequence data derived from plants of natural origin and an F1 mapping population were used to identify sex-linked genes and the nonrecombining region. Amborella trichopoda has a ZW sex determination system. Analysis of genes in a 4 Mb nonrecombining sex-determination region reveals recent divergence of Z and W gametologs, and few Z- and W-specific genes. The sex chromosomes of A. trichopoda evolved less than 16.5 Myr ago, long after the divergence of the extant angiosperms.

SNP Development in Penaeus vannamei via Next-Generation Sequencing and DNA Pool Sequencing

Next-generation sequencing and pool sequencing have been widely used in SNP (single-nucleotide polymorphism) detection and population genetics research; however, there are few reports on SNPs related to the growth of Penaeus vannamei. The purpose of this study was to call SNPs from rapid-growing (RG) and slow-growing (SG) individuals’ transcriptomes and use DNA pool sequencing to assess the reliability of SNPs. Two parameters were applied to detect SNPs. One parameter was the p-values generated using Fisher’s exact test, which were used to calculate the significance of allele frequency differences between RG and SG. The other one was the AFI (minor allele frequency imbalance), which was defined to highlight the fold changes in MAF (minor allele frequency) values between RG and SG. There were 216,015 hypothetical SNPs, which were obtained based on the transcriptome data. Finally, 104 high-quality SNPs and 96,819 low-quality SNPs were predicted. Then, 18 high-quality SNPs and 17 low-quality SNPs were selected to assess the reliability of the detection process. Here, 72.22% (13/18) accuracy was achieved for high-quality SNPs, while only 52.94% (9/17) accuracy was achieved for low-quality SNPs. These SNPs enrich the data for population genetics studies of P. vannamei and may play a role in the development of SNP markers for future breeding studies.

Plant genera Cannabis and Humulus share the same pair of well-differentiated sex chromosomes

We recently described, in Cannabis sativa, the oldest sex chromosome system documented so far in plants (12-28 Myr old). Based on the estimated age, we predicted that it should be shared by its sister genus Humulus, which is known also to possess XY chromosomes. Here, we used transcriptome sequencing of an F₁ family of H. lupulus to identify and study the sex chromosomes in this species using the probabilistic method SEX-DETector. We identified 265 sex-linked genes in H. lupulus, which preferentially mapped to the C. sativa X chromosome. Using phylogenies of sex-linked genes, we showed that a region of the sex chromosomes had already stopped recombining in an ancestor of both species. Furthermore, as in C. sativa, Y-linked gene expression reduction is correlated to the position on the X chromosome, and highly Y degenerated genes showed dosage compensation. We report, for the first time in Angiosperms, a sex chromosome system that is shared by two different genera. Thus, recombination suppression started at least 21-25 Myr ago, and then (either gradually or step-wise) spread to a large part of the sex chromosomes (c. 70%), leading to a degenerated Y chromosome.

Genomic Approaches for Conservation Management in Australia under Climate Change

Conservation genetics has informed threatened species management for several decades. With the advent of advanced DNA sequencing technologies in recent years, it is now possible to monitor and manage threatened populations with even greater precision. Climate change presents a number of threats and challenges, but new genomics data and analytical approaches provide opportunities to identify critical evolutionary processes of relevance to genetic management under climate change. Here, we discuss the applications of such approaches for threatened species management in Australia in the context of climate change, identifying methods of facilitating viability and resilience in the face of extreme environmental stress. Using genomic approaches, conservation management practices such as translocation, targeted gene flow, and gene-editing can now be performed with the express intention of facilitating adaptation to current and projected climate change scenarios in vulnerable species, thus reducing extinction risk and ensuring the protection of our unique biodiversity for future generations. We discuss the current barriers to implementing conservation genomic projects and the efforts being made to overcome them, including communication between researchers and managers to improve the relevance and applicability of genomic studies. We present novel approaches for facilitating adaptive capacity and accelerating natural selection in species to encourage resilience in the face of climate change.

Detecting sex-linked genes using genotyped individuals sampled in natural populations

We propose a method, SDpop, able to infer sex-linkage caused by recombination suppression typical of sex chromosomes. The method is based on the modeling of the allele and genotype frequencies of individuals of known sex in natural populations. It is implemented in a hierarchical probabilistic framework, accounting for different sources of error. It allows statistical testing for the presence or absence of sex chromosomes, and detection of sex-linked genes based on the posterior probabilities in the model. Furthermore, for gametologous sequences, the haplotype and level of nucleotide polymorphism of each copy can be inferred, as well as the divergence between them. We test the method using simulated data, as well as data from both a relatively recent and an old sex chromosome system (the plant Silene latifolia and humans) and show that, for most cases, robust predictions are obtained with 5 to 10 individuals per sex.

RNA sequencing describes both population structure and plasticity-selection dynamics in a non-model fish

Background

Messenger RNA sequencing is becoming more common in studies of non-model species and is most often used for gene expression-based investigations. However, the method holds potential for numerous other applications as well—including analyses of alternative splicing, population structure, and signatures of selection. To maximize the utility of mRNA data sets, distinct analyses may be combined such as by exploring dynamics between gene expression with signatures of selection in the context of population structure. Here, we compare two published data sets describing two populations of a minnow species endemic to the San Francisco Estuary (Sacramento splittail, Pogonichthys macrolepidotus): a microsatellite data set showing population structure, and an mRNA whole transcriptome data set obtained after the two populations were exposed to a salinity challenge. We compared measures of population structure and genetic variation using single nucleotide polymorphisms (SNPs) called from mRNA from the whole transcriptome sequencing study with those patterns determined from microsatellites. For investigating plasticity and evolution, intra- and inter-population transcriptome plasticity was investigated with differential gene expression, differential exon usage, and gene expression variation. Outlier SNP analysis was also performed on the mRNA data set and signatures of selection and phenotypic plasticity were investigated on an individual-gene basis.

Results

We found that mRNA sequencing revealed patterns of population structure consistent with those found with microsatellites, but with lower magnitudes of genetic variation and population differentiation consistent with widespread purifying selection expected when using mRNA. In addition, within individual genes, phenotypic plasticity or signatures of selection were found in almost mutual exclusion (except heatr6, nfu1, slc22a6, sya, and mmp13).

Conclusions

These results show that an mRNA sequencing data set may have multiple uses, including describing population structure and for investigating the mechanistic interplay of evolution and plasticity in adaptation. MRNA sequencing thus complements traditional sequencing methods used for population genetics, in addition to its utility for describing phenotypic plasticity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07592-4.

Dioecy Is Associated with High Genetic Diversity and Adaptation Rates in the Plant Genus Silene

About 15,000 angiosperm species (∼6%) have separate sexes, a phenomenon known as dioecy. Why dioecious taxa are so rare is still an open question. Early work reported lower species richness in dioecious compared with nondioecious sister clades, raising the hypothesis that dioecy may be an evolutionary dead-end. This hypothesis has been recently challenged by macroevolutionary analyses that detected no or even positive effect of dioecy on diversification. However, the possible genetic consequences of dioecy at the population level, which could drive the long-term fate of dioecious lineages, have not been tested so far. Here, we used a population genomics approach in the Silene genus to look for possible effects of dioecy, especially for potential evidence of evolutionary handicaps of dioecy underlying the dead-end hypothesis. We collected individual-based RNA-seq data from several populations in 13 closely related species with different sexual systems: seven dioecious, three hermaphroditic, and three gynodioecious species. We show that dioecy is associated with increased genetic diversity, as well as higher selection efficacy both against deleterious mutations and for beneficial mutations. The results hold after controlling for phylogenetic inertia, differences in species census population sizes and geographic ranges. We conclude that dioecious Silene species neither show signs of increased mutational load nor genetic evidence for extinction risk. We discuss these observations in the light of the possible demographic differences between dioecious and self-compatible hermaphroditic species and how this could be related to alternatives to the dead-end hypothesis to explain the rarity of dioecy.

De novo Transcriptome Assembly of Myllocerinus aurolineatus Voss in Tea Plants

Myllocerinus aurolineatus Voss is a species of the insecta class in the arthropod. In this study, we first observed and identified M. aurolineatus Voss in tea plants in Guizhou, China, where it caused severe quantity and quality losses in tea plants. Knowledge on M. aurolineatus Voss genome is inadequate, especially for biological or functional research. We performed the first transcriptome sequencing by using the Illumina Hiseq™ technique on M. aurolineatus Voss. Over 55.9 million high-quality paired-end reads were generated and assembled into 69,439 unigenes using the Trinity short read software, resulting in a cluster of 1,207 bp of the N50 length. A total of 69,439 genes were predicted by BLAST to known proteins in the NCBI database and were distributed into Gene Ontology (20,190), eukaryotic complete genomes (12,488), and the Kyoto Encyclopedia of Genes and Genomes (3,170). We also identified 96,790 single-nucleotide polymorphisms and 13,121 simple sequence repeats in these unigenes. Our transcriptome data provide a useful resource for future functional studies of M. aurolineatus Voss for dispersal control in tea plants.

The Road to Sorghum Domestication: Evidence From Nucleotide Diversity and Gene Expression Patterns

Native African cereals (sorghum, millets) ensure food security to millions of low-income people from low fertility and drought-prone regions of Africa and Asia. In spite of their agronomic importance, the genetic bases of their phenotype and adaptations are still not well-understood. Here we focus on Sorghum bicolor, which is the fifth cereal worldwide for grain production and constitutes the staple food for around 500 million people. We leverage transcriptomic resources to address the adaptive consequences of the domestication process. Gene expression and nucleotide variability were analyzed in 11 domesticated and nine wild accessions. We documented a downregulation of expression and a reduction of diversity both in nucleotide polymorphism (30%) and gene expression levels (18%) in domesticated sorghum. These findings at the genome-wide level support the occurrence of a global reduction of diversity during the domestication process, although several genes also showed patterns consistent with the action of selection. Nine hundred and forty-nine genes were significantly differentially expressed between wild and domesticated gene pools. Their functional annotation points to metabolic pathways most likely contributing to the sorghum domestication syndrome, such as photosynthesis and auxin metabolism. Coexpression network analyzes revealed 21 clusters of genes sharing similar expression patterns. Four clusters (totaling 2,449 genes) were significantly enriched in differentially expressed genes between the wild and domesticated pools and two were also enriched in domestication and improvement genes previously identified in sorghum. These findings reinforce the evidence that the combined and intricated effects of the domestication and improvement processes do not only affect the behaviors of a few genes but led to a large rewiring of the transcriptome. Overall, these analyzes pave the way toward the identification of key domestication genes valuable for genetic resources characterization and breeding purposes.

Bedrock radioactivity influences the rate and spectrum of mutation

All organisms on Earth are exposed to low doses of natural radioactivity but some habitats are more radioactive than others. Yet, documenting the influence of natural radioactivity on the evolution of biodiversity is challenging. Here, we addressed whether organisms living in naturally more radioactive habitats accumulate more mutations across generations using 14 species of waterlice living in subterranean habitats with contrasted levels of radioactivity. We found that the mitochondrial and nuclear mutation rates across a waterlouse species’ genome increased on average by 60% and 30%, respectively, when radioactivity increased by a factor of three. We also found a positive correlation between the level of radioactivity and the probability of G to T (and complementary C to A) mutations, a hallmark of oxidative stress. We conclude that even low doses of natural bedrock radioactivity influence the mutation rate possibly through the accumulation of oxidative damage, in particular in the mitochondrial genome.

Genomic inference of complex domestication histories in three Solanaceae species

Domestication is a human-induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species-specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a "protracted period" of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species-specific and supported by the few historical records available.

How Much Does Ne Vary Among Species?

Genetic drift is an important evolutionary force of strength inversely proportional to Ne , the effective population size. The impact of drift on genome diversity and evolution is known to vary among species, but quantifying this effect is a difficult task. Here we assess the magnitude of variation in drift power among species of animals via its effect on the mutation load - which implies also inferring the distribution of fitness effects of deleterious mutations. To this aim, we analyze the nonsynonymous (amino-acid changing) and synonymous (amino-acid conservative) allele frequency spectra in a large sample of metazoan species, with a focus on the primates vs. fruit flies contrast. We show that a Gamma model of the distribution of fitness effects is not suitable due to strong differences in estimated shape parameters among taxa, while adding a class of lethal mutations essentially solves the problem. Using the Gamma + lethal model and assuming that the mean deleterious effects of nonsynonymous mutations is shared among species, we estimate that the power of drift varies by a factor of at least 500 between large-Ne and small-Ne species of animals, i.e., an order of magnitude more than the among-species variation in genetic diversity. Our results are relevant to Lewontin's paradox while further questioning the meaning of the Ne parameter in population genomics.

No evidence for accumulation of deleterious mutations and fitness degradation in clonal fish hybrids: Abandoning sex without regrets

Despite its inherent costs, sexual reproduction is ubiquitous in nature, and the mechanisms to protect it from a competitive displacement by asexuality remain unclear. Popular mutation-based explanations, like the Muller's ratchet and the Kondrashov's hatchet, assume that purifying selection may not halt the accumulation of deleterious mutations in the nonrecombining genomes, ultimately leading to their degeneration. However, empirical evidence is scarce and it remains particularly unclear whether mutational degradation proceeds fast enough to ensure the decay of clonal organisms and to prevent them from outcompeting their sexual counterparts. To test this hypothesis, we jointly analysed the exome sequences and the fitness-related phenotypic traits of the sexually reproducing fish species and their clonal hybrids, whose evolutionary ages ranged from F1 generations to 300 ky. As expected, mutations tended to accumulate in the clonal genomes in a time-dependent manner. However, contrary to the predictions, we found no trend towards increased nonsynonymity of mutations acquired by clones, nor higher radicality of their amino acid substitutions. Moreover, there was no evidence for fitness degeneration in the old clones compared with that in the younger ones. In summary, although an efficacy of purifying selection may still be reduced in the asexual genomes, our data indicate that its efficiency is not drastically decreased. Even the oldest investigated clone was found to be too young to suffer fitness consequences from a mutation accumulation. This suggests that mechanisms other than mutation accumulation may be needed to explain the competitive advantage of sex in the short term.

Evidence for Dosage Compensation in Coccinia grandis, a Plant with a Highly Heteromorphic XY System

About 15,000 angiosperms are dioecious, but the mechanisms of sex determination in plants remain poorly understood. In particular, how Y chromosomes evolve and degenerate, and whether dosage compensation evolves as a response, are matters of debate. Here, we focus on Coccinia grandis, a dioecious cucurbit with the highest level of X/Y heteromorphy recorded so far. We identified sex-linked genes using RNA sequences from a cross and a model-based method termed SEX-DETector. Parents and F1 individuals were genotyped, and the transmission patterns of SNPs were then analyzed. In the >1300 sex-linked genes studied, maximum X-Y divergence was 0.13-0.17, and substantial Y degeneration is implied by an average Y/X expression ratio of 0.63 and an inferred gene loss on the Y of ~40%. We also found reduced Y gene expression being compensated by elevated expression of corresponding genes on the X and an excess of sex-biased genes on the sex chromosomes. Molecular evolution of sex-linked genes in C. grandis is thus comparable to that in Silene latifolia, another dioecious plant with a strongly heteromorphic XY system, and cucurbits are the fourth plant family in which dosage compensation is described, suggesting it might be common in plants.

Transcriptome of Sphaerospora molnari (Cnidaria, Myxosporea) blood stages provides proteolytic arsenal as potential therapeutic targets against sphaerosporosis in common carp

Background

Parasites employ proteases to evade host immune systems, feed and replicate and are often the target of anti-parasite strategies to disrupt these interactions. Myxozoans are obligate cnidarian parasites, alternating between invertebrate and fish hosts. Their genes are highly divergent from other metazoans, and available genomic and transcriptomic datasets are limited. Some myxozoans are important aquaculture pathogens such as Sphaerospora molnari replicating in the blood of farmed carp before reaching the gills for sporogenesis and transmission. Proliferative stages cause a massive systemic lymphocyte response and the disruption of the gill epithelia by spore-forming stages leads to respiratory problems and mortalities. In the absence of a S. molnari genome, we utilized a de novo approach to assemble the first transcriptome of proliferative myxozoan stages to identify S. molnari proteases that are upregulated during the first stages of infection when the parasite multiplies massively, rather than in late spore-forming plasmodia. Furthermore, a subset of orthologs was used to characterize 3D structures and putative druggable targets.

Results

An assembled and host filtered transcriptome containing 9436 proteins, mapping to 29,560 contigs was mined for protease virulence factors and revealed that cysteine proteases were most common (38%), at a higher percentage than other myxozoans or cnidarians (25–30%). Two cathepsin Ls that were found upregulated in spore-forming stages with a presenilin like aspartic protease and a dipeptidyl peptidase. We also identified downregulated proteases in the spore-forming development when compared with proliferative stages including an astacin metallopeptidase and lipases (qPCR). In total, 235 transcripts were identified as putative proteases using a MEROPS database. In silico analysis of highly transcribed cathepsins revealed potential drug targets within this data set that should be prioritised for development.

Conclusions

In silico surveys for proteins are essential in drug discovery and understanding host-parasite interactions in non-model systems. The present study of S. molnari’s protease arsenal reveals previously unknown proteases potentially used for host exploitation and immune evasion. The pioneering dataset serves as a model for myxozoan virulence research, which is of particular importance as myxozoan diseases have recently been shown to emerge and expand geographically, due to climate change.

Is adaptation limited by mutation? A timescale-dependent effect of genetic diversity on the adaptive substitution rate in animals

Whether adaptation is limited by the beneficial mutation supply is a long-standing question of evolutionary genetics, which is more generally related to the determination of the adaptive substitution rate and its relationship with species effective population size (Ne) and genetic diversity. Empirical evidence reported so far is equivocal, with some but not all studies supporting a higher adaptive substitution rate in large-Ne than in small-Ne species. We gathered coding sequence polymorphism data and estimated the adaptive amino-acid substitution rate ωa, in 50 species from ten distant groups of animals with markedly different population mutation rate θ. We reveal the existence of a complex, timescale dependent relationship between species adaptive substitution rate and genetic diversity. We find a positive relationship between ωa and θ among closely related species, indicating that adaptation is indeed limited by the mutation supply, but this was only true in relatively low-θ taxa. In contrast, we uncover no significant correlation between ωa and θ at a larger taxonomic scale, suggesting that the proportion of beneficial mutations scales negatively with species' long-term Ne.

An efficient RNA-seq-based segregation analysis identifies the sex chromosomes of Cannabis sativa

Cannabissativa–derived tetrahydrocannabinol (THC) production is increasing very fast worldwide. C. sativa is a dioecious plant with XY Chromosomes, and only females (XX) are useful for THC production. Identifying the sex chromosome sequence would improve early sexing and better management of this crop; however, the C. sativa genome projects have failed to do so. Moreover, as dioecy in the Cannabaceae family is ancestral, C. sativa sex chromosomes are potentially old and thus very interesting to study, as little is known about old plant sex chromosomes. Here, we RNA-sequenced a C. sativa family (two parents and 10 male and female offspring, 576 million reads) and performed a segregation analysis for all C. sativa genes using the probabilistic method SEX-DETector. We identified >500 sex-linked genes. Mapping of these sex-linked genes to a C. sativa genome assembly identified the largest chromosome pair being the sex chromosomes. We found that the X-specific region (not recombining between X and Y) is large compared to other plant systems. Further analysis of the sex-linked genes revealed that C. sativa has a strongly degenerated Y Chromosome and may represent the oldest plant sex chromosome system documented so far. Our study revealed that old plant sex chromosomes can have large, highly divergent nonrecombining regions, yet still be roughly homomorphic.

An efficient RNA-seq-based segregation analysis identifies the sex chromosomes of Cannabis sativa

Cannabis sativa-derived tetrahydrocannabinol (THC) production is increasing very fast worldwide. C. sativa is a dioecious plant with XY Chromosomes, and only females (XX) are useful for THC production. Identifying the sex chromosome sequence would improve early sexing and better management of this crop; however, the C. sativa genome projects have failed to do so. Moreover, as dioecy in the Cannabaceae family is ancestral, C. sativa sex chromosomes are potentially old and thus very interesting to study, as little is known about old plant sex chromosomes. Here, we RNA-sequenced a C. sativa family (two parents and 10 male and female offspring, 576 million reads) and performed a segregation analysis for all C. sativa genes using the probabilistic method SEX-DETector. We identified >500 sex-linked genes. Mapping of these sex-linked genes to a C. sativa genome assembly identified the largest chromosome pair being the sex chromosomes. We found that the X-specific region (not recombining between X and Y) is large compared to other plant systems. Further analysis of the sex-linked genes revealed that C. sativa has a strongly degenerated Y Chromosome and may represent the oldest plant sex chromosome system documented so far. Our study revealed that old plant sex chromosomes can have large, highly divergent nonrecombining regions, yet still be roughly homomorphic.

The de novo transcriptome of workers head of the higher group termite Globitermes sulphureus Haviland (Blattodea: Termitidae)

The subterranean termite Globitermus sulphureus is an important Southeast Asian pest with limited genomic resources that causes damages to agriculture crops and building structures. Therefore, the main goal of this study was to survey the G. sulphureus transcriptome composition. Here, we performed de novo transcriptome for G. sulphureus workers’ heads using Illumina HiSeq paired-end sequencing technology. A total of 88, 639, 408 clean reads were collected and assembled into 243, 057 transcripts and 193, 344 putative genes. The transcripts were annotated with the Trinotate pipeline. In total, 27, 061 transcripts were successfully annotated using BLASTX against the SwissProt database and 17, 816 genes were assigned to 47, 598 GO terms. We classified 14, 223 transcripts into COG classification, resulting in 25 groups of functional annotations. Next, a total of 12, 194 genes were matched in the KEGG pathway and 392 metabolic pathways were predicted based on the annotation. Moreover, we detected two endogenous cellulases in the sequences. The RT-qPCR analysis showed that there were significant differences in the expression levels of two genes β-glucosidase and endo-β-1,4-glucanase between worker and soldier heads of G. sulphureus. This is the first study to characterize the complete head transcriptome of a higher termite G. sulphureus using a high-throughput sequencing. Our study may provide an overview and comprehensive molecular resource for comparative studies of the transcriptomics and genomics of termites.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Evolutionary transcriptomics reveals the origins of olives and the genomic changes associated with their domestication

The olive (Olea europaea L. subsp. europaea) is one of the oldest and most socio-economically important cultivated perennial crop in the Mediterranean region. Yet, its origins are still under debate and the genetic bases of the phenotypic changes associated with its domestication are unknown. We generated RNA-sequencing data for 68 wild and cultivated olive trees to study the genetic diversity and structure both at the transcription and sequence levels. To localize putative genes or expression pathways targeted by artificial selection during domestication, we employed a two-step approach in which we identified differentially expressed genes and screened the transcriptome for signatures of selection. Our analyses support a major domestication event in the eastern part of the Mediterranean basin followed by dispersion towards the West and subsequent admixture with western wild olives. While we found large changes in gene expression when comparing cultivated and wild olives, we found no major signature of selection on coding variants and weak signals primarily affected transcription factors. Our results indicated that the domestication of olives resulted in only moderate genomic consequences and that the domestication syndrome is mainly related to changes in gene expression, consistent with its evolutionary history and life history traits.

A transcriptome-based phylogenetic study of hard ticks (Ixodidae)

Hard ticks are widely distributed across temperate regions, show strong variation in host associations, and are potential vectors of a diversity of medically important zoonoses, such as Lyme disease. To address unresolved issues with respect to the evolutionary relationships among certain species or genera, we produced novel RNA-Seq data sets for nine different Ixodes species. We combined this new data with 18 data sets obtained from public databases, both for Ixodes and non-Ixodes hard tick species, using soft ticks as an outgroup. We assembled transcriptomes (for 27 species in total), predicted coding sequences and identified single copy orthologues (SCO). Using Maximum-likelihood and Bayesian frameworks, we reconstructed a hard tick phylogeny for the nuclear genome. We also obtained a mitochondrial DNA-based phylogeny using published genome sequences and mitochondrial sequences derived from the new transcriptomes. Our results confirm previous studies showing that the Ixodes genus is monophyletic and clarify the relationships among Ixodes sub-genera. This work provides a baseline for studying the evolutionary history of ticks: we indeed found an unexpected acceleration of substitutions for mitochondrial sequences of Prostriata, and for nuclear and mitochondrial genes of two species of Rhipicephalus, which we relate with patterns of genome architecture and changes of life-cycle, respectively.

Codon Usage Bias in Animals: Disentangling the Effects of Natural Selection, Effective Population Size, and GC-Biased Gene Conversion

Selection on codon usage bias is well documented in a number of microorganisms. Whether codon usage is also generally shaped by natural selection in large organisms, despite their relatively small effective population size (Ne), is unclear. In animals, the population genetics of codon usage bias has only been studied in a handful of model organisms so far, and can be affected by confounding, nonadaptive processes such as GC-biased gene conversion and experimental artefacts. Using population transcriptomics data, we analyzed the relationship between codon usage, gene expression, allele frequency distribution, and recombination rate in 30 nonmodel species of animals, each from a different family, covering a wide range of effective population sizes. We disentangled the effects of translational selection and GC-biased gene conversion on codon usage by separately analyzing GC-conservative and GC-changing mutations. We report evidence for effective translational selection on codon usage in large-Ne species of animals, but not in small-Ne ones, in agreement with the nearly neutral theory of molecular evolution. C- and T-ending codons tend to be preferred over synonymous G- and A-ending ones, for reasons that remain to be determined. In contrast, we uncovered a conspicuous effect of GC-biased gene conversion, which is widespread in animals and the main force determining the fate of AT↔GC mutations. Intriguingly, the strength of its effect was uncorrelated with Ne.

Parallel pattern of differentiation at a genomic island shared between clinal and mosaic hybrid zones in a complex of cryptic seahorse lineages

Diverging semi-isolated lineages either meet in narrow clinal hybrid zones, or have a mosaic distribution associated with environmental variation. Intrinsic reproductive isolation is often emphasized in the former and local adaptation in the latter, although both reduce gene flow between groups. Rarely are these two patterns of spatial distribution reported in the same study system. Here, we report that the long-snouted seahorse Hippocampus guttulatus is subdivided into discrete panmictic entities by both types of hybrid zones. Along the European Atlantic coasts, a northern and a southern lineage meet in the southwest of France where they coexist in sympatry-i.e., in the same geographical zone-with little hybridization. In the Mediterranean Sea, two lineages have a mosaic distribution, associated with lagoon-like and marine habitats. A fifth lineage was identified in the Black Sea. Genetic homogeneity over large spatial scales contrasts with isolation maintained in sympatry or close parapatry at a fine scale. A high variation in locus-specific introgression rates provides additional evidence that partial reproductive isolation must be maintaining the divergence. We find that fixed differences between lagoon and marine populations in the Mediterranean Sea belong to the most differentiated SNPs between the two Atlantic lineages, against the genome-wide pattern of structure that mostly follow geography. These parallel outlier SNPs cluster on a single chromosome-wide island of differentiation. Since Atlantic lineages do not map to lagoon-sea habitat variation, genetic parallelism at the genomic island suggests a shared genetic barrier contributes to reproductive isolation in contrasting contexts-i.e., spatial versus ecological. We discuss how a genomic hotspot of parallel differentiation could have evolved and become associated both with space and with a patchy environment in a single study system.

Evolution of Young Sex Chromosomes in Two Dioecious Sister Plant Species with Distinct Sex Determination Systems

In the last decade, progress has been made in methods to identify the sex determination system in plants. This gives the opportunity to study sex chromosomes that arose independently at different phylogenetic scales, and thus allows the discovery and the understanding of early stages of sex chromosome evolution. In the genus Silene, sex chromosomes have evolved independently in at least two clades from a nondioecious ancestor, the Melandrium and Otites sections. In the latter, sex chromosomes could be younger than in the section Melandrium, based on phylogenetic studies and as no heteromorphic sex chromosomes have been detected. This section might also exhibit lability in sex determination, because male heterogamy and female heterogamy have been suggested to occur.

In this study, we investigated the sex determination system of two dioecious species in the section Otites (Silene otites and its close relative Silene pseudotites). Applying the new probabilistic method SEX-DETector on RNA-seq data from cross-controlled progenies, we inferred their most likely sex determination system and a list of putative autosomal and sex-linked contigs. We showed that the two phylogenetically close species differed in their sex determination system (XY versus ZW) with sex chromosomes that derived from two different pairs of autosomes. We built a genetic map of the sex chromosomes and showed that both pairs exhibited a large region with lack of recombination. However, the sex-limited chromosomes exhibited no strong degeneration. Finally, using the “ancestral” autosomal expression of sex-linked orthologs of nondioecious S. nutans, we found a slight signature of dosage compensation in the heterogametic females of S. otites.

Genomic Selection in Aquaculture: Application, Limitations and Opportunities With Special Reference to Marine Shrimp and Pearl Oysters

Within aquaculture industries, selection based on genomic information (genomic selection) has the profound potential to change genetic improvement programs and production systems. Genomic selection exploits the use of realized genomic relationships among individuals and information from genome-wide markers in close linkage disequilibrium with genes of biological and economic importance. We discuss the technical advances, practical requirements, and commercial applications that have made genomic selection feasible in a range of aquaculture industries, with a particular focus on molluscs (pearl oysters, Pinctada maxima) and marine shrimp (Litopenaeus vannamei and Penaeus monodon). The use of low-cost genome sequencing has enabled cost-effective genotyping on a large scale and is of particular value for species without a reference genome or access to commercial genotyping arrays. We highlight the pitfalls and offer the solutions to the genotyping by sequencing approach and the building of appropriate genetic resources to undertake genomic selection from first-hand experience. We describe the potential to capture large-scale commercial phenotypes based on image analysis and artificial intelligence through machine learning, as inputs for calculation of genomic breeding values. The application of genomic selection over traditional aquatic breeding programs offers significant advantages through being able to accurately predict complex polygenic traits including disease resistance; increasing rates of genetic gain; minimizing inbreeding; and negating potential limiting effects of genotype by environment interactions. Further practical advantages of genomic selection through the use of large-scale communal mating and rearing systems are highlighted, as well as presenting rate-limiting steps that impact on attaining maximum benefits from adopting genomic selection. Genomic selection is now at the tipping point where commercial applications can be readily adopted and offer significant short- and long-term solutions to sustainable and profitable aquaculture industries.

Nonrandom RNAseq gene expression associated with RNAlater and flash freezing storage methods

RNA sequencing is a popular next-generation sequencing technique for assaying genome-wide gene expression profiles. Nonetheless, it is susceptible to biases that are introduced by sample handling prior gene expression measurements. Two of the most common methods for preserving samples in both field-based and laboratory conditions are submersion in RNAlater and flash freezing in liquid nitrogen. Flash freezing in liquid nitrogen can be impractical, particularly for field collections. RNAlater is a solution for stabilizing tissue for longer-term storage as it rapidly permeates tissue to protect cellular RNA. In this study, we assessed genome-wide expression patterns in 30-day-old fry collected from the same brood at the same time point that were flash-frozen in liquid nitrogen and stored at -80°C or submerged and stored in RNAlater at room temperature, simulating conditions of fieldwork. We show that sample storage is a significant factor influencing observed differential gene expression. In particular, genes with elevated GC content exhibit higher observed expression levels in liquid nitrogen flash-freezing relative to RNAlater storage. Further, genes with higher expression in RNAlater relative to liquid nitrogen experience disproportionate enrichment for functional categories, many of which are involved in RNA processing. This suggests that RNAlater may elicit a physiological response that has the potential to bias biological interpretations of expression studies. The biases introduced to observed gene expression arising from mimicking many field-based studies are substantial and should not be ignored.

Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species

Population genetic theory predicts that selection should be more effective when the effective population size (N_e) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger N_e. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.