Scaling-up RADseq methods for large datasets of non-invasive samples: Lessons for library construction and data preprocessing.

Genetic non-invasive sampling (gNIS) is a critical tool for population genetics studies, supporting conservation efforts while imposing minimal impacts on wildlife. However, gNIS often presents variable levels of DNA degradation and non-endogenous contamination, which can incur considerable processing costs. Furthermore, the use of restriction-site-associated DNA sequencing methods (RADseq) for assessing thousands of genetic markers introduces the challenge of obtaining large sets of shared loci with similar coverage across multiple individuals. Here, we present an approach to handling large-scale gNIS-based datasets using data from the spotted hyena population inhabiting the Ngorongoro Crater in Tanzania. We generated 3RADseq data for more than a thousand individuals, mostly from faecal mucus samples collected non-invasively and varying in DNA degradation and contamination level. Using small-scale sequencing, we screened samples for endogenous DNA content, removed highly contaminated samples, confirmed overlap fragment length between libraries, and balanced individual representation in a sequencing pool. We evaluated the impact of (1) DNA degradation and contamination of non-invasive samples, (2) PCR duplicates and (3) different SNP filters on genotype accuracy based on Mendelian error estimated for parent-offspring trio datasets. Our results showed that when balanced for sequencing depth, contaminated samples presented similar genotype error rates to those of non-contaminated samples. We also showed that PCR duplicates and different SNP filters impact genotype accuracy. In summary, we showed the potential of using gNIS for large-scale genetic monitoring based on SNPs and demonstrated how to improve control over library preparation by using a weighted re-pooling strategy that considers the endogenous DNA content.

How genomics can help biodiversity conservation.

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Evidence of backcross inviability and mitochondrial DNA paternal leakage in sea turtle hybrids.

Hybridization is known to be part of many species' evolutionary history. Sea turtles have a fascinating hybridization system in which species separated by as much as 43 million years are still capable of hybridizing. Indeed, the largest nesting populations in Brazil of loggerheads (Caretta caretta) and hawksbills (Eretmochelys imbricata) have a high incidence of hybrids between these two species. A third species, olive ridleys (Lepidochelys olivacea), is also known to hybridize although at a smaller scale. Here, we used restriction site-associated DNA sequencing (RAD-Seq) markers, mitogenomes, and satellite-telemetry to investigate the patterns of hybridization and introgression in the Brazilian sea turtle population and their relationship with the migratory behaviours between feeding and nesting aggregations. We also explicitly test if the mixing of two divergent genomes in sea turtle hybrids causes mitochondrial paternal leakage. We developed a new species-specific PCR-assay capable of detecting mitochondrial DNA (mtDNA) inheritance from both parental species and performed ultra-deep sequencing to estimate the abundance of each mtDNA type. Our results show that all adult hybrids are first generation (F1) and most display a loggerhead migratory behaviour. We detected paternal leakage in F1 hybrids and different proportions of mitochondria from maternal and paternal species. Although previous studies showed no significant fitness decrease in hatchlings, our results support genetically-related hybrid breakdown possibly caused by cytonuclear incompatibility. Further research on hybrids from other populations in addition to Brazil and between different species will show if backcross inviability and mitochondrial paternal leakage is observed across sea turtle species.

Evolution of the ancestral mammalian karyotype and syntenic regions.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

The era of reference genomes in conservation genomics.

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics.

New Genetic Insights About Hybridization and Population Structure of Hawksbill and Loggerhead Turtles From Brazil.

An extremely high incidence of hybridization among sea turtles is found along the Brazilian coast. This atypical phenomenon and its impact on sea turtle conservation can be elucidated through research focused on the evolutionary history of sea turtles. We assessed high-quality multilocus haplotypes of 143 samples of the 5 species of sea turtles that occur along the Brazilian coast to investigate the hybridization process and the population structure of hawksbill (Eretmochelys imbricata) and loggerhead turtles (Caretta caretta). The multilocus data were initially used to characterize interspecific hybrids. Introgression (F2 hybrids) was only confirmed in hatchlings of F1 hybrid females (hawksbill × loggerhead), indicating that introgression was either previously overestimated and F2 hybrids may not survive to adulthood, or the first-generation hybrid females nesting in Brazil were born as recent as few decades ago. Phylogenetic analyses using nuclear markers recovered the mtDNA-based Indo-Pacific and Atlantic lineages for hawksbill turtles, demonstrating a deep genetic divergence dating from the early Pliocene. In addition, loggerhead turtles that share a common feeding area and belong to distinct Indo-Pacific and Atlantic mtDNA clades present no clear genetic differentiation at the nuclear level. Finally, our results indicate that hawksbill and loggerhead rookeries along the Brazilian coast are likely connected by male-mediated gene flow.

Long-read DNA metabarcoding of ribosomal RNA in the analysis of fungi from aquatic environments.

DNA metabarcoding is widely used to study prokaryotic and eukaryotic microbial diversity. Technological constraints limit most studies to marker lengths below 600 base pairs (bp). Longer sequencing reads of several thousand bp are now possible with third-generation sequencing. Increased marker lengths provide greater taxonomic resolution and allow for phylogenetic methods of classification, but longer reads may be subject to higher rates of sequencing error and chimera formation. In addition, most bioinformatics tools for DNA metabarcoding were designed for short reads and are therefore unsuitable. Here, we used Pacific Biosciences circular consensus sequencing (CCS) to DNA-metabarcode environmental samples using a ca. 4,500 bp marker that included most of the eukaryote SSU and LSU rRNA genes and the complete ITS region. We developed an analysis pipeline that reduced error rates to levels comparable to short-read platforms. Validation using a mock community indicated that our pipeline detected 98% of chimeras de novo. We recovered 947 OTUs from water and sediment samples from a natural lake, 848 of which could be classified to phylum, 397 to genus and 330 to species. By allowing for the simultaneous use of three databases (Unite, SILVA and RDP LSU), long-read DNA metabarcoding provided better taxonomic resolution than any single marker. We foresee the use of long reads enabling the cross-validation of reference sequences and the synthesis of ribosomal rRNA gene databases. The universal nature of the rRNA operon and our recovery of >100 nonfungal OTUs indicate that long-read DNA metabarcoding holds promise for studies of eukaryotic diversity more broadly.

Demography and selection shape transcriptomic divergence in field crickets.

Gene flow, demography, and selection can result in similar patterns of genomic variation and disentangling their effects is key to understanding speciation. Here, we assess transcriptomic variation to unravel the evolutionary history of Gryllus rubens and Gryllus texensis, cryptic field cricket species with highly divergent mating behavior. We infer their demographic history and screen their transcriptomes for footprints of selection in the context of the inferred demography. We find strong support for a long history of bidirectional gene flow, which ceased during the late Pleistocene, and a bottleneck in G. rubens consistent with a peripatric origin of this species. Importantly, the demographic history has likely strongly shaped patterns of genetic differentiation (empirical FST distribution). Concordantly, FST -based selection detection uncovers a large number of outliers, likely comprising many false positives, echoing recent theoretical insights. Alternative genetic signatures of positive selection, informed by the demographic history of the sibling species, highlighted a smaller set of loci; many of these are candidates for controlling variation in mating behavior. Our results underscore the importance of demography in shaping overall patterns of genetic divergence and highlight that examining both demography and selection facilitates a more complete understanding of genetic divergence during speciation.

Transcriptome profiling of ontogeny in the acridid grasshopper Chorthippus biguttulus.

Acridid grasshoppers (Orthoptera:Acrididae) are widely used model organisms for developmental, evolutionary, and neurobiological research. Although there has been recent influx of orthopteran transcriptomic resources, many use pooled ontogenetic stages obscuring information about changes in gene expression during development. Here we developed a de novo transcriptome spanning 7 stages in the life cycle of the acridid grasshopper Chorthippus biguttulus. Samples from different stages encompassing embryonic development through adults were used for transcriptomic profiling, revealing patterns of differential gene expression that highlight processes in the different life stages. These patterns were validated with semi-quantitative RT-PCR. Embryonic development showed a strongly differentiated expression pattern compared to all of the other stages and genes upregulated in this stage were involved in signaling, cellular differentiation, and organ development. Our study is one of the first to examine gene expression during post-embryonic development in a hemimetabolous insect and we found that only the fourth and fifth instars had clusters of genes upregulated during these stages. These genes are involved in various processes ranging from synthesis of biogenic amines to chitin binding. These observations indicate that post-embryonic ontogeny is not a continuous process and that some instars are differentiated. Finally, genes upregulated in the imago were generally involved in aging and immunity. Our study highlights the importance of looking at ontogeny as a whole and indicates promising directions for future research in orthopteran development.

Canine distemper virus in the Serengeti ecosystem: molecular adaptation to different carnivore species.

Was the 1993/1994 fatal canine distemper virus (CDV) epidemic in lions and spotted hyaenas in the Serengeti ecosystem caused by the recent spillover of a virulent domestic dog strain or one well adapted to these noncanids? We examine this question using sequence data from 13 'Serengeti' strains including five complete genomes obtained between 1993 and 2011. Phylogenetic and haplotype network analyses reveal that strains from noncanids during the epidemic were more closely related to each other than to those from domestic or wild canids. All noncanid 'Serengeti' strains during the epidemic encoded: (1) one novel substitution G134S in the CDV-V protein; and (2) the rare amino acid combination 519I/549H at two sites under positive selection in the region of the CDV-H protein that binds to SLAM (CD 150) host cell receptors. Worldwide, only a few noncanid strains in the America II lineage encode CDV-H 519I/549H. All canid 'Serengeti' strains during the epidemic coded CDV-V 134G, and CDV-H 519R/549Y, or 519R/549H. A functional assay of cell entry revealed the highest performance by CDV-H proteins encoding 519I/549H in cells expressing lion SLAM receptors, and the highest performance by proteins encoding 519R/549Y, typical of dog strains worldwide, in cells expressing dog SLAM receptors. Our findings are consistent with an epidemic in lions and hyaenas caused by CDV variants better adapted to noncanids than canids, but not with the recent spillover of a dog strain. Our study reveals a greater complexity of CDV molecular epidemiology in multihost environments than previously thought.

PCR Strategies for Complete Allele Calling in Multigene Families Using High-Throughput Sequencing Approaches.

The characterization of multigene families with high copy number variation is often approached through PCR amplification with highly degenerate primers to account for all expected variants flanking the region of interest. Such an approach often introduces PCR biases that result in an unbalanced representation of targets in high-throughput sequencing libraries that eventually results in incomplete detection of the targeted alleles. Here we confirm this result and propose two different amplification strategies to alleviate this problem. The first strategy (called pooled-PCRs) targets different subsets of alleles in multiple independent PCRs using different moderately degenerate primer pairs, whereas the second approach (called pooled-primers) uses a custom-made pool of non-degenerate primers in a single PCR. We compare their performance to the common use of a single PCR with highly degenerate primers using the MHC class I of the Iberian lynx as a model. We found both novel approaches to work similarly well and better than the conventional approach. They significantly scored more alleles per individual (11.33 ± 1.38 and 11.72 ± 0.89 vs 7.94 ± 1.95), yielded more complete allelic profiles (96.28 ± 8.46 and 99.50 ± 2.12 vs 63.76 ± 15.43), and revealed more alleles at a population level (13 vs 12). Finally, we could link each allele's amplification efficiency with the primer-mismatches in its flanking sequences and show that ultra-deep coverage offered by high-throughput technologies does not fully compensate for such biases, especially as real alleles may reach lower coverage than artefacts. Adopting either of the proposed amplification methods provides the opportunity to attain more complete allelic profiles at lower coverages, improving confidence over the downstream analyses and subsequent applications.

Evidence for selection maintaining MHC diversity in a rodent species despite strong density fluctuations.

Strong spatiotemporal variation in population size often leads to reduced genetic diversity limiting the adaptive potential of individual populations. Key genes of adaptive variation are encoded by the immune genes of the major histocompatibility complex (MHC) playing an essential role in parasite resistance. How MHC variation persists in rodent populations that regularly experience population bottlenecks remains an important topic in evolutionary genetics. We analysed the consequences of strong population fluctuations on MHC class II DRB exon 2 diversity in two distant common vole (Microtus arvalis) populations in three consecutive years using a high-throughput sequencing approach. In 143 individuals, we detected 25 nucleotide alleles translating into 14 unique amino acid MHC alleles belonging to at least three loci. Thus, the overall allelic diversity and amino acid distance among the remaining MHC alleles, used as a surrogate for the range of pathogenic antigens that can be presented to T-cells, are still remarkably high. Both study populations did not show significant population differentiation between years, but significant differences were found between sites. We concluded that selection processes seem to be strong enough to maintain moderate levels of MHC diversity in our study populations outcompeting genetic drift, as the same MHC alleles were conserved between years. Differences in allele frequencies between populations might be the outcome of different local parasite pressures and/or genetic drift. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in cyclic populations.

A genes eye view of ontogeny: de novo assembly and profiling of the Gryllus rubens transcriptome.

Crickets (Orthoptera:Gryllidae) are widely used model organisms for developmental, evolutionary, neurobiological and behavioural research. Here, we developed a de novo transcriptome from pooled RNA-seq Illumina data spanning seven stages in the life cycle of Gryllus rubens. Approximately 705 Mbp of data was assembled and filtered to form 27 312 transcripts. We were able to annotate 52% of our transcripts using BLAST and assign at least one gene ontology term to 41%. Pooled samples from three different ontogenetic stages were used for transcriptomic profiling revealing patterns of differential gene expression that highlight processes in the different life stages. Embryonic and early instar development was enriched for ecdysteroid metabolism, cytochrome P450s and glutathione production. Late instar development was enriched for regulation of gene expression and many of the genes highly expressed during this stage were involved in conserved developmental signalling pathways suggesting that these developmental pathways are active beyond embryonic development. Adults were enriched for fat transport (mostly relating to egg production) and production of octopamine, an important neurohormone. We also identified genes involved in conserved developmental pathways (Hedgehog, Hippo, Wnt, JAK/STAT, TGF-beta, Notch, and MEK/ERK). This is the first transcriptome spanning ontogeny in Gryllus rubens and a valuable resource for future work on development and evolution in Orthoptera.

Sequence diversity between class I MHC loci of African native and introduced Bos taurus cattle in Theileria parva endemic regions: in silico peptide binding prediction identifies distinct functional clusters.

There is strong evidence that the immunity induced by live vaccination for control of the protozoan parasite Theileria parva is mediated by class I MHC-restricted CD8(+) T cells directed against the schizont stage of the parasite that infects bovine lymphocytes. The functional competency of class I MHC genes is dependent on the presence of codons specifying certain critical amino acid residues that line the peptide binding groove. Compared with European Bos taurus in which class I MHC allelic polymorphisms have been examined extensively, published data on class I MHC transcripts in African taurines in T. parva endemic areas is very limited. We utilized the multiplexing capabilities of 454 pyrosequencing to make an initial assessment of class I MHC allelic diversity in a population of Ankole cattle. We also typed a population of exotic Holstein cattle from an African ranch for class I MHC and investigated the extent, if any, that their peptide-binding motifs overlapped with those of Ankole cattle. We report the identification of 18 novel allelic sequences in Ankole cattle and provide evidence of positive selection for sequence diversity, including in residues that predominantly interact with peptides. In silico functional analysis resulted in peptide binding specificities that were largely distinct between the two breeds. We also demonstrate that CD8(+) T cells derived from Ankole cattle that are seropositive for T. parva do not recognize vaccine candidate antigens originally identified in Holstein and Boran (Bos indicus) cattle breeds.

Spatial pattern of adaptive and neutral genetic diversity across different biomes in the lesser anteater (Tamandua tetradactyla).

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and activation of the adaptive immune response and are thought to be regulated by natural selection, especially due to pathogen-driven selective pressure. In this study, we investigated the spatial distribution of MHC class II DRB exon 2 gene diversity of the lesser anteater (Tamandua tetradactyla) across five Brazilian biomes using next-generation sequencing and compared the MHC pattern with that of neutral markers (microsatellites). We found a noticeable high level of diversity in DRB (60 amino acid alleles in 65 individuals) and clear signatures of historical positive selection acting on this gene. Higher allelic richness and proportion of private alleles were found in rain forest biomes, especially Amazon forest, a megadiverse biome, possibly harboring greater pathogen richness as well. Neutral markers, however, showed a similar pattern to DRB, demonstrating the strength of demography as an additional force to pathogen-driven selection in shaping MHC diversity and structure. This is the first characterization and description of diversity of a MHC gene for any member of the magna-order Xenarthra, one of the basal lineages of placental mammals.

Endogenous murine leukemia retroviral variation across wild European and inbred strains of house mouse.

BACKGROUND: Endogenous murine leukemia retroviruses (MLVs) are high copy number proviral elements difficult to comprehensively characterize using standard low throughput sequencing approaches. However, high throughput approaches generate data that is challenging to process, interpret and present. RESULTS: Next generation sequencing (NGS) data was generated for MLVs from two wild caught Mus musculus domesticus (from mainland France and Corsica) and for inbred laboratory mouse strains C3H, LP/J and SJL. Sequence reads were grouped using a novel sequence clustering approach as applied to retroviral sequences. A Markov cluster algorithm was employed, and the sequence reads were queried for matches to specific xenotropic (Xmv), polytropic (Pmv) and modified polytropic (Mpmv) viral reference sequences. CONCLUSIONS: Various MLV subtypes were more widespread than expected among the mice, which may be due to the higher coverage of NGS, or to the presence of similar sequence across many different proviral loci. The results did not correlate with variation in the major MLV receptor Xpr1, which can restrict exogenous MLVs, suggesting that endogenous MLV distribution may reflect gene flow more than past resistance to infection.

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence.

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA-seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using FST -based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.

A novel endogenous betaretrovirus in the common vampire bat (Desmodus rotundus) suggests multiple independent infection and cross-species transmission events.

The Desmodus rotundus endogenous betaretrovirus (DrERV) is fixed in the vampire bat D. rotundus population and in other phyllostomid bats but is not present in all species from this family. DrERV is not phylogenetically related to Old World bat betaretroviruses but to betaretroviruses from rodents and New World primates, suggesting recent cross-species transmission. A recent integration age estimation of the provirus in some taxa indicates that an exogenous counterpart might have been in recent circulation.