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.

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.

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.

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.

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.

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.

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.

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.

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.

Whole genome sequencing and methylome analysis of the wild guinea pig.

BACKGROUND: DNA methylation is a heritable mechanism that acts in response to environmental changes, lifestyle and diseases by influencing gene expression in eukaryotes. Epigenetic studies of wild organisms are mandatory to understand their role in e.g. adaptational processes in the great variety of ecological niches. However, strategies to address those questions on a methylome scale are widely missing. In this study we present such a strategy and describe a whole genome sequence and methylome analysis of the wild guinea pig. RESULTS: We generated a full Wild guinea pig (Cavia aperea) genome sequence with enhanced coverage of methylated regions, benefiting from the available sequence of the domesticated relative Cavia porcellus. This new genome sequence was then used as reference to map the sequence reads of bisulfite treated Wild guinea pig sequencing libraries to investigate DNA-methylation patterns at nucleotide-specific level, by using our here described method, named 'DNA-enrichment-bisulfite-sequencing' (MEBS). The results achieved using MEBS matched those of standard methods in other mammalian model species. The technique is cost efficient, and incorporates both methylation enrichment results and a nucleotide-specific resolution even without a whole genome sequence available. Thus MEBS can be easily applied to extend methylation enrichment studies to a nucleotide-specific level. CONCLUSIONS: The approach is suited to study methylomes of not yet sequenced mammals at single nucleotide resolution. The strategy is transferable to other mammalian species by applying the nuclear genome sequence of a close relative. It is therefore of interest for studies on a variety of wild species trying to answer evolutionary, adaptational, ecological or medical questions by epigenetic mechanisms.

A haplotype-resolved reference genome of a long-distance migratory bat, Pipistrellus nathusii (Keyserling & Blasius, 1839).

We present a complete, chromosome-scale reference genome for the long-distance migratory bat Pipistrellus nathusii. The genome encompasses both haplotypic sets of autosomes and the separation of both sex chromosomes by utilizing highly accurate long-reads and preserving long-range phasing information through the use of three-dimensional chromatin conformation capture sequencing (Hi-C). This genome, accompanied by a comprehensive protein-coding sequence annotation, provides a valuable genomic resource for future investigations into the genomic bases of long-distance migratory flight in bats as well as uncovering the genetic architecture, population structure and evolutionary history of Pipistrellus nathusii. The reference-quality genome presented here gives a fundamental resource to further our understanding of bat genetics and evolution, adding to the growing number of high-quality genetic resources in this field. Here, we demonstrate its use in the phylogenetic reconstruction of the order Chiroptera, and in particular, we present the resources to allow detailed investigations into the genetic drivers and adaptations related to long-distance migration.

MHC genotyping of non-model organisms using next-generation sequencing: a new methodology to deal with artefacts and allelic dropout.

BACKGROUND: The Major Histocompatibility Complex (MHC) is the most important genetic marker to study patterns of adaptive genetic variation determining pathogen resistance and associated life history decisions. It is used in many different research fields ranging from human medical, molecular evolutionary to functional biodiversity studies. Correct assessment of the individual allelic diversity pattern and the underlying structural sequence variation is the basic requirement to address the functional importance of MHC variability. Next-generation sequencing (NGS) technologies are likely to replace traditional genotyping methods to a great extent in the near future but first empirical studies strongly indicate the need for a rigorous quality control pipeline. Strict approaches for data validation and allele calling to distinguish true alleles from artefacts are required. RESULTS: We developed the analytical methodology and validated a data processing procedure which can be applied to any organism. It allows the separation of true alleles from artefacts and the evaluation of genotyping reliability, which in addition to artefacts considers for the first time the possibility of allelic dropout due to unbalanced amplification efficiencies across alleles. Finally, we developed a method to assess the confidence level per genotype a-posteriori, which helps to decide which alleles and individuals should be included in any further downstream analyses. The latter method could also be used for optimizing experiment designs in the future. CONCLUSIONS: Combining our workflow with the study of amplification efficiency offers the chance for researchers to evaluate enormous amounts of NGS-generated data in great detail, improving confidence over the downstream analyses and subsequent applications.

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.

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.

Evolution and population structure of Salmonella enterica serovar Newport.

Salmonellosis caused by Salmonella enterica serovar Newport is a major global public health concern, particularly because S. Newport isolates that are resistant to multiple drugs (MDR), including third-generation cephalosporins (MDR-AmpC phenotype), have been commonly isolated from food animals. We analyzed 384 S. Newport isolates from various sources by a multilocus sequence typing (MLST) scheme to study the evolution and population structure of the serovar. These were compared to the population structure of S. enterica serovars Enteritidis, Kentucky, Paratyphi B, and Typhimurium. Our S. Newport collection fell into three lineages, Newport-I, Newport-II, and Newport-III, each of which contained multiple sequence types (STs). Newport-I has only a few STs, unlike Newport-II or Newport-III, and has possibly emerged recently. Newport-I is more prevalent among humans in Europe than in North America, whereas Newport-II is preferentially associated with animals. Two STs of Newport-II encompassed all MDR-AmpC isolates, suggesting recent global spread after the acquisition of the bla(CMY-2) gene. In contrast, most Newport-III isolates were from humans in North America and were pansusceptible to antibiotics. Newport was intermediate in population structure to the other serovars, which varied from a single monophyletic lineage in S. Enteritidis or S. Typhimurium to four discrete lineages within S. Paratyphi B. Both mutation and homologous recombination are responsible for diversification within each of these lineages, but the relative frequencies differed with the lineage. We conclude that serovars of S. enterica provide a variety of different population structures.

Estimation of divergence time between two sibling species of the Anopheles (Kerteszia) cruzii complex using a multilocus approach.

BACKGROUND: Anopheles cruzii is the primary human Plasmodium vector in southern and southeastern Brazil. The distribution of this mosquito follows the coast of the Brazilian Atlantic Forest. Previous studies indicated that An. cruzii is a complex of cryptic species. RESULTS: A multilocus approach using six loci, three circadian clock genes and three encoding ribosomal proteins, was implemented to investigate in more detail the genetic differentiation between the An. cruzii populations from Santa Catarina (southern Brazil) and Bahia States (northeastern Brazil) that represent two sibling species. The analysis revealed very high FST values and fixed differences between the two An. cruzii sibling species in all loci, irrespective of their function. An Isolation with Migration model was fit to the data using the IM program. The results reveal no migration in either direction and allowed a rough estimate of the divergence time between the two sibling species. CONCLUSIONS: Population genetics analysis of An. cruzii samples from two Brazilian localities using a multilocus approach confirmed that they represent two different sibling species in this complex. The results suggest that the two species have not exchanged migrants since their separation and that they possibly diverged between 1.1 and 3.6 million years ago, a period of intense climatic changes.

A multiplex single nucleotide polymorphism typing assay for detecting mutations that result in decreased fluoroquinolone susceptibility in Salmonella enterica serovars Typhi and Paratyphi A.

OBJECTIVES: Decreased susceptibility to fluoroquinolones has become a major problem for the successful therapy of human infections caused by Salmonella enterica, especially the life-threatening typhoid and paratyphoid fevers. METHODS: By using Luminex xTAG beads, we developed a rapid, reliable and cost-effective multiplexed genotyping assay for simultaneously detecting 11 mutations in gyrA, gyrB and parE of S. enterica serovars Typhi and Paratyphi A that result in nalidixic acid resistance (Nal(R)) and/or decreased susceptibility to fluoroquinolones. RESULTS: This assay yielded unambiguous single nucleotide polymorphism calls on extracted DNA from 292 isolates of Salmonella Typhi (Nal(R) = 223 and Nal(S) = 69) and 106 isolates of Salmonella Paratyphi A (Nal(R) = 24 and Nal(S) = 82). All of the 247 Nal(R) Salmonella Typhi and Salmonella Paratyphi A isolates were found to harbour at least one of the target mutations, with GyrA Phe-83 as the most common one (143/223 for Salmonella Typhi and 18/24 for Salmonella Paratyphi A). We also identified three GyrB mutations in eight Nal(S) Salmonella Typhi isolates (six for GyrB Phe-464, one for GyrB Leu-465 and one for GyrB Asp-466), and mutations GyrB Phe-464 and GyrB Asp-466 seem to be related to the decreased ciprofloxacin susceptibility phenotype in Salmonella Typhi. This assay can also be used directly on boiled single colonies. CONCLUSIONS: The assay presented here would be useful for clinical and reference laboratories to rapidly screen quinolone-resistant isolates of Salmonella Typhi and Salmonella Paratyphi A, and decipher the underlying genetic changes for epidemiological purposes.