Low-coverage Genome Sequencing Research Articles

How many species are there? Lineage diversification and hidden speciation in Solanaceae from highland grasslands in southern South America.

Species delimitation can be challenging when analysing recently diverged species, especially those taxonomically synonymised due to morphological similarities. We aimed to untangle the relationships between two grassland species, Petunia guarapuavensis and Petunia scheideana, exploring the dynamics of fast divergence and addressing their species delimitation. We used a low-coverage genome sequencing and population genomic approach to distinguish species and populations between P. guarapuavensis and P. scheideana. Our analysis focused on detecting structuration, hybridisation/introgression, and phylogenetic patterns. We employed demographic models to support species delimitation while exploring potential phylogeographic barriers influencing gene flow. Our findings indicated differentiation between the two species and revealed another lineage, which was phylogenetically distinct from the others and had no evidence of gene flow with them. The presence of a river acted as a phylogeographic barrier, limiting gene flow and allowing for structuration between closely related lineages. The optimal species delimitation scenario involved secondary contact between well-established lineages. The rapid divergence observed in these Petunia species explains the lack of significant morphological differences, as floral diagnostic traits in species sharing the pollinators tend to evolve more slowly. This study highlights the complexity of species delimitation in recently diverged groups and emphasises the importance of genomic approaches in understanding evolutionary relationships and speciation dynamics.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States.

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

Skimming genomes for systematics and DNA barcodes of corals.

Numerous genomic methods developed over the past two decades have enabled the discovery and extraction of orthologous loci to help resolve phylogenetic relationships across various taxa and scales. Genome skimming (or low-coverage genome sequencing) is a promising method to not only extract high-copy loci but also 100s to 1000s of phylogenetically informative nuclear loci (e.g., ultraconserved elements [UCEs] and exons) from contemporary and museum samples. The subphylum Anthozoa, including important ecosystem engineers (e.g., stony corals, black corals, anemones, and octocorals) in the marine environment, is in critical need of phylogenetic resolution and thus might benefit from a genome-skimming approach. We conducted genome skimming on 242 anthozoan corals collected from 1886 to 2022. Using existing target-capture baitsets, we bioinformatically obtained UCEs and exons from the genome-skimming data and incorporated them with data from previously published target-capture studies. The mean number of UCE and exon loci extracted from the genome skimming data was 1837 ± 662 SD for octocorals and 1379 ± 476 SD loci for hexacorals. Phylogenetic relationships were well resolved within each class. A mean of 1422 ± 720 loci was obtained from the historical specimens, with 1253 loci recovered from the oldest specimen collected in 1886. We also obtained partial to whole mitogenomes and nuclear rRNA genes from >95% of samples. Bioinformatically pulling UCEs, exons, mitochondrial genomes, and nuclear rRNA genes from genome skimming data is a viable and low-cost option for phylogenetic studies. This approach can be used to review and support taxonomic revisions and reconstruct evolutionary histories, including historical museum and type specimens.

Interspecific cytogenomic comparison reveals a potential chromosomal centromeric marker in Proceratophrys frog species.

Among the repetitive elements, satellite DNA (SatDNA) emerges as extensive arrays of highly similar tandemly repeated units, spanning megabases in length. Given that the satDNA PboSat01-176, previously characterized in P. boiei, prompted our interest for having a high abundance in P. boiei and potential for centromeric satellite, here, we employed various approaches, including low coverage genome sequencing, followed by computational analysis and chromosomal localization techniques in four Proceratophrys species and, investigating the genomic presence and sharing, as well as its potential for chromosomal centromere marker in Proceratophrys frog species. Our findings demonstrate that PboSat01-176 exhibits high abundance across all four Proceratophrys species, displaying distinct characteristics that establish it as the predominant repetitive DNA element in these species. The satellite DNA is prominently clustered in the peri/centromeric region of the chromosomes, particularly in the heterochromatic regions. The widespread presence of PboSat01-176 in closely related Proceratophrys species reinforces the validity of the library hypothesis for repetitive sequences. Thus, this study highlighted the utility of the satDNA family PboSat01-176 as a reliable centromeric marker in Proceratophrys species, with potential to be applied in other species of anuran amphibians. The observed sharing and maintenance of this sequence within the genus suggest possibilities for future research, particularly through expanded sampling to elucidate parameters that underlie the library hypothesis and the evolutionary dynamics of satDNA sequences.

Abstract 2411: A deep learning-based multimodal ensemble algorithm for lung cancer early detection with cross-ethnic generalizability

Abstract Background: Lung cancer is the leading cause of cancer-related mortality worldwide, and early detection is crucial for improving patient outcomes. Cell-free DNA (cfDNA) whole genome sequencing (WGS) is a non-invasive technique that has the potential to detect lung cancer early, but it is challenging to develop accurate and generalizable detection algorithms due to the low abundance of tumor-derived cfDNA in early-stage cancer. Methods: We developed a novel deep learning-based multimodal ensemble algorithm that combined two independent models: FEMS (fragment end motif frequency and size) and COV (coverage of genome positions). To assess the performance and ethnic generalizability of our algorithm, we evaluated it on three distinct cohorts: a Korean development cohort, a Korean external cohort, and a Caucasian external cohort. The Korean development cohort for model development comprised 236 lung cancer samples (stage I: 30%, stage II: 12%, stage III: 24%, stage IV: 20%, unknown: 14%) and 2205 normal samples. The Korean external cohort for performance evaluation of the model comprised 288 lung cancer samples (stage I: 31%, stage II: 13%, stage III: 26%, stage IV: 19%, unknown: 12%) and 1463 normal samples. And the Caucasian external cohort for confirm the ethnic generalizability consisted of 126 lung cancer samples (stage I: 17%, stage II: 12%, stage III: 37%, stage IV: 32%, unknown: 3%) and 93 normal samples. All data were generated using Novaseq6000 with a minimum of 40 million reads. Results: Our multimodal ensemble algorithm achieved high sensitivity and specificity in all three cohorts. The Korean development cohort had an overall sensitivity of 94.1% (95% CI: 90.7% to 97%) and specificity of 85%. The Korean external cohort had an overall sensitivity of 90.4% (95% CI: 82.7% to 98.1%) and specificity of 82%. The Caucasian external cohort had an overall sensitivity of 92.1% (95% CI: 87.3% to 96%) and specificity of 82.8%. The Caucasian cohort also had a sensitivity of 80% for early-stage lung cancer. Conclusion: Our deep learning-based multimodal ensemble algorithm is a promising tool for early lung cancer detection using low-coverage cell-free whole genome sequencing (LC-cfWGS). It performed well on three distinct cohorts, including a Caucasian cohort, suggesting that it is robust and generalizable to different populations. This makes it a potential candidate for clinical implementation. Citation Format: Tae-Rim Lee, Jin Mo Ahn, Junnam Lee, Dasom Kim, Byeong-Ho Jeong, Dongryul Oh, Mengchi Wang, Michael Salmans, Andrew Carson, Bryan Leatham, Kristin Fathe, Byung In Lee, Chang-Seok Ki, Young Sik Park, Eun-Hae Cho. A deep learning-based multimodal ensemble algorithm for lung cancer early detection with cross-ethnic generalizability [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 2411.

Genome-wide association study for numbers of vertebrae in Dezhou donkey population reveals new candidate genes

Numbers of vertebrae is an important economic trait associated with body size and meat productivity in animals. However, the genetic basis of vertebrae number in donkey remains to be well understood. The aim of this study was to identify candidate genes affecting the number of thoracic (TVn) and the number of lumbar vertebrae (LVn) in Dezhou donkey. A genome-wide association study was conducted using whole genome sequence data imputed from low-coverage genome sequencing. For TVn, we identified 38 genome-wide significant and 64 suggestive SNPs, which relate to 7 genes (NLGN1, DCC, SLC26A7, TOX, WNT7A, LOC123286078, and LOC123280142). For LVn, we identified 9 genome-wide significant and 38 suggestive SNPs, which relate to 8 genes (GABBR2, FBXO4, LOC123277146, LOC123277359, BMP7, B3GAT1, EML2, and LRP5). The genes involve in the Wnt and TGF-β signaling pathways and may play an important role in embryonic development or bone formation and could be good candidate genes for TVn and LVn.

Chromosomal and genomic data in Neotropical Dorstenia species (Moraceae) and their relationships with African species

Abstract Dorstenia spp. are characterized by a herbaceous habit, the presence of rhizomes and inflorescences being coenanthia. They occur in different forest formations, and it is believed that there are two distribution centres for Dorstenia, one Palaeotropical (Africa) and one Neotropical (Central and South America). The present study aimed to undertake a comparative analysis of the karyotype diversity of Dorstenia, based on chromosome data and DNA C values of 13 Neotropical species, and an analysis of low-coverage genomic sequencing of three African and one Neotropical species. The most frequent chromosome number is 2n = 32, with variations for 2n = 28 and 30, suggesting an important role of polyploidy in the karyotype differentiation of Dorstenia. Karyotypes varied in the number of metacentric, submetacentric and acrocentric chromosomes, with evidence of rearrangements. Genome sizes (C values) varied almost two-fold, from 2C = 2.50 to 5.47 pg. There was no clear relationship between DNA C value variations and karyotype features, such as chromosome numbers and distribution of CMA/DAPI bands or rDNA sites. Proximal heterochromatin was most frequent, but CMA+ and DAPI+ bands predominated in relation to DAPI+ bands. 5S and 35S ribosomal DNA (rDNA) sites were also observed in the proximal regions in almost all species. Based on a chromosome modelling approach, our data suggest that dysploidy events may have been responsible for the chromosome number variations, whereas rearrangements involving repetitive fractions may have been associated with changes in genome size. There was great variability in the amount of rDNA between species, including a remarkable 5S rDNA amplification in D. bahiensis, and a diversity of transposable elements in the four genomes with genomic sequence data. The combination of data from karyotypes, genomes and phylogeny suggests that dysploidy has played a predominant role in the evolution of Neotropical Dorstenia spp.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi.

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

Aegilops crassa Boiss. repeatome characterized using low-coverage NGS as a source of new FISH markers: Application in phylogenetic studies of the Triticeae.

Aegilops crassa Boiss. is polyploid grass species that grows in the eastern part of the Fertile Crescent, Afghanistan, and Middle Asia. It consists of tetraploid (4x) and hexaploid (6x) cytotypes (2n = 4x = 28, D1D (Abdolmalaki et al., 2019) XcrXcr and 2n = 6x = 42, D1D (Abdolmalaki et al., 2019) XcrXcrD2D (Adams and Wendel, 2005), respectively) that are similar morphologically. Although many Aegilops species were used in wheat breeding, the genetic potential of Ae. crassa has not yet been exploited due to its uncertain origin and significant genome modifications. Tetraploid Ae. crassa is thought to be the oldest polyploid Aegilops species, the subgenomes of which still retain some features of its ancient diploid progenitors. The D1 and D2 subgenomes of Ae. crassa were contributed by Aegilops tauschii (2n = 2x = 14, DD), while the Xcr subgenome donor is still unknown. Owing to its ancient origin, Ae. crassa can serve as model for studying genome evolution. Despite this, Ae. crassa is poorly studied genetically and no genome sequences were available for this species. We performed low-coverage genome sequencing of 4x and 6x cytotypes of Ae. crassa, and four Ae. tauschii accessions belonging to different subspecies; diploid wheatgrass Thinopyrum bessarabicum (Jb genome), which is phylogenetically close to D (sub)genome species, was taken as an outgroup. Subsequent data analysis using the pipeline RepeatExplorer2 allowed us to characterize the repeatomes of these species and identify several satellite sequences. Some of these sequences are novel, while others are found to be homologous to already known satellite sequences of Triticeae species. The copy number of satellite repeats in genomes of different species and their subgenome (D1 or Xcr) affinity in Ae. crassa were assessed by means of comparative bioinformatic analysis combined with quantitative PCR (qPCR). Fluorescence in situ hybridization (FISH) was performed to map newly identified satellite repeats on chromosomes of common wheat, Triticum aestivum, 4x and 6x Ae. crassa, Ae. tauschii, and Th. bessarabicum. The new FISH markers can be used in phylogenetic analyses of the Triticeae for chromosome identification and the assessment of their subgenome affinities and for evaluation of genome/chromosome constitution of wide hybrids or polyploid species.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species.

Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

Nyssorhynchus darlingi genome-wide studies related to microgeographic dispersion and blood-seeking behavior

BackgroundIn Brazil, malaria is concentrated in the Amazon Basin, where more than 99% of the annual cases are reported. The main goal of this study was to investigate the population structure and genetic association of the biting behavior of Nyssorhynchus (also known as Anopheles) darlingi, the major malaria vector in the Amazon region of Brazil, using low-coverage genomic sequencing data.MethodsSamples were collected in the municipality of Mâncio Lima, Acre state, Brazil between 2016 and 2017. Different approaches using genotype imputation and no gene imputation for data treatment and low-coverage sequencing genotyping were performed. After the samples were genotyped, population stratification analysis was performed.ResultsWeak but statistically significant stratification signatures were identified between subpopulations separated by distances of approximately 2–3 km. Genome-wide association studies (GWAS) were performed to compare indoor/outdoor biting behavior and blood-seeking at dusk/dawn. A statistically significant association was observed between biting behavior and single nucleotide polymorphism (SNP) markers adjacent to the gene associated with cytochrome P450 (CYP) 4H14, which is associated with insecticide resistance. A statistically significant association between blood-seeking periodicity and SNP markers adjacent to genes associated with the circadian cycle was also observed.ConclusionThe data presented here suggest that low-coverage whole-genome sequencing with adequate processing is a powerful tool to genetically characterize vector populations at a microgeographic scale in malaria transmission areas, as well as for use in GWAS. Female mosquitoes entering houses to take a blood meal may be related to a specific CYP4H14 allele, and female timing of blood-seeking is related to circadian rhythm genes.Graphical

Prevalence and circulation patterns of SARS-CoV-2 variants in European sewage mirror clinical data of 54 European cities

For community-level monitoring, the European Commission under the EU Sewage Sentinel System recommends wastewater-based SARS-CoV-2 surveillance. Tracking SARS-CoV-2 variants in a community is pivotal for appropriate public health response. Genome sequencing of SARS-CoV-2 in wastewater samples for tracking variants is challenging, often resulting in low coverage genome sequences, thereby impeding the detection of the SARS-CoV-2 mutations. Therefore, we aimed at high-coverage SARS-CoV-2 genome sequences from sewage samples which we successfully accomplished. This first pan-European surveillance compared the mutation profiles associated with the variants of concerns: B.1.1.7, P.1, B.1.351 and B.1.617.2 across 20 European countries, including 54 municipalities. The results highlight that SARS-CoV-2 variants detected in the wastewater samples mirror the variants profiles reported in clinical data. This study demonstrated that >98% coverage of SARS-CoV-2 genomic sequences is possible and can be used to track SARS-CoV-2 mutations in wastewater to support identifying variants circulating in a city at the community level.

Oceanic Crustal Fluid Single Cell Genomics Complements Metagenomic and Metatranscriptomic Surveys With Orders of Magnitude Less Sample Volume.

Fluids circulating through oceanic crust play important roles in global biogeochemical cycling mediated by their microbial inhabitants, but studying these sites is challenged by sampling logistics and low biomass. Borehole observatories installed at the North Pond study site on the western flank of the Mid-Atlantic Ridge have enabled investigation of the microbial biosphere in cold, oxygenated basaltic oceanic crust. Here we test a methodology that applies redox-sensitive fluorescent molecules for flow cytometric sorting of cells for single cell genomic sequencing from small volumes of low biomass (approximately 103 cells ml–1) crustal fluid. We compare the resulting genomic data to a recently published paired metagenomic and metatranscriptomic analysis from the same site. Even with low coverage genome sequencing, sorting cells from less than one milliliter of crustal fluid results in similar interpretation of dominant taxa and functional profiles as compared to ‘omics analysis that typically filter orders of magnitude more fluid volume. The diverse community dominated by Gammaproteobacteria, Bacteroidetes, Desulfobacterota, Alphaproteobacteria, and Zetaproteobacteria, had evidence of autotrophy and heterotrophy, a variety of nitrogen and sulfur cycling metabolisms, and motility. Together, results indicate fluorescence activated cell sorting methodology is a powerful addition to the toolbox for the study of low biomass systems or at sites where only small sample volumes are available for analysis.

Comprehensive and accurate genetic variant identification from contaminated and low-coverage Mycobacterium tuberculosis whole genome sequencing data.

Improved understanding of the genomic variants that allow Mycobacterium tuberculosis (Mtb) to acquire drug resistance, or tolerance, and increase its virulence are important factors in controlling the current tuberculosis epidemic. Current approaches to Mtb sequencing, however, cannot reveal Mtb’s full genomic diversity due to the strict requirements of low contamination levels, high Mtb sequence coverage and elimination of complex regions. We have developed the XBS (compleX Bacterial Samples) bioinformatics pipeline, which implements joint calling and machine-learning-based variant filtering tools to specifically improve variant detection in the important Mtb samples that do not meet these criteria, such as those from unbiased sputum samples. Using novel simulated datasets, which permit exact accuracy verification, XBS was compared to the UVP and MTBseq pipelines. Accuracy statistics showed that all three pipelines performed equally well for sequence data that resemble those obtained from culture isolates of high depth of coverage and low-level contamination. In the complex genomic regions, however, XBS accurately identified 9.0 % more SNPs and 8.1 % more single nucleotide insertions and deletions than the WHO-endorsed unified analysis variant pipeline. XBS also had superior accuracy for sequence data that resemble those obtained directly from sputum samples, where depth of coverage is typically very low and contamination levels are high. XBS was the only pipeline not affected by low depth of coverage (5–10×), type of contamination and excessive contamination levels (>50 %). Simulation results were confirmed using whole genome sequencing (WGS) data from clinical samples, confirming the superior performance of XBS with a higher sensitivity (98.8%) when analysing culture isolates and identification of 13.9 % more variable sites in WGS data from sputum samples as compared to MTBseq, without evidence for false positive variants when rRNA regions were excluded. The XBS pipeline facilitates sequencing of less-than-perfect Mtb samples. These advances will benefit future clinical applications of Mtb sequencing, especially WGS directly from clinical specimens, thereby avoiding in vitro biases and making many more samples available for drug resistance and other genomic analyses. The additional genetic resolution and increased sample success rate will improve genome-wide association studies and sequence-based transmission studies.

Genomic analysis of distinct bleaching tolerances among cryptic coral species.

Reef-building coral species are experiencing an unprecedented decline owing to increasing frequency and intensity of marine heatwaves and associated bleaching-induced mortality. Closely related species from the Acropora hyacinthus species complex differ in heat tolerance and in their association with heat-tolerant symbionts. We used low-coverage full genome sequencing of 114 colonies monitored across the 2015 bleaching event in American Samoa to determine the genetic differences among four cryptic species (termed HA, HC, HD and HE) that have diverged in these species traits. Cryptic species differed strongly at thousands of single nucleotide polymorphisms across the genome which are enriched for amino acid changes in the bleaching-resistant species HE. In addition, HE also showed two particularly divergent regions with strong signals of differentiation. One approximately 220 kb locus, HES1, contained the majority of fixed differences in HE. A second locus, HES2, was fixed in HE but polymorphic in the other cryptic species. Surprisingly, non-HE individuals with HE-like haplotypes at HES2 were more likely to bleach. At both loci, HE showed particular sequence similarity to a congener, Acropora millepora. Overall, resilience to bleaching during the third global bleaching event was strongly structured by host cryptic species, buoyed by differences in symbiont associations between these species.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline.

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

Chromosome-specific polymorphic SSR markers in tropical eucalypt species using low coverage whole genome sequences: systematic characterization and validation.

Eucalyptus is one of the major plantation species with wide variety of industrial uses. Polymorphic and informative simple sequence repeats (SSRs) have broad range of applications in genetic analysis. In this study, two individuals of Eucalyptus tereticornis (ET217 and ET86), one individual each from E. camaldulensis (EC17) and E. grandis (EG9) were subjected to whole genome resequencing. Low coverage (10×) genome sequencing was used to find polymorphic SSRs between the individuals. Average number of SSR loci identified was 95,513 and the density of SSRs per Mb was from 157.39 in EG9 to 155.08 in EC17. Among all the SSRs detected, the most abundant repeat motifs were di-nucleotide (59.6%–62.5%), followed by tri- (23.7%–27.2%), tetra- (5.2%–5.6%), penta- (5.0%–5.3%), and hexa-nucleotide (2.7%–2.9%). The predominant SSR motif units were AG/CT and AAG/TTC. Computational genome analysis predicted the SSR length variations between the individuals and identified the gene functions of SSR containing sequences. Selected subset of polymorphic markers was validated in a full-sib family of eucalypts. Additionally, genome-wide characterization of single nucleotide polymorphisms, InDels and transcriptional regulators were carried out. These variations will find their utility in genome-wide association studies as well as understanding of molecular mechanisms involved in key economic traits. The genomic resources generated in this study would provide an impetus to integrate genomics in marker-trait associations and breeding of tropical eucalypts.

CONSULT: accurate contamination removal using locality-sensitive hashing.

A fundamental question appears in many bioinformatics applications: Does a sequencing read belong to a large dataset of genomes from some broad taxonomic group, even when the closest match in the set is evolutionarily divergent from the query? For example, low-coverage genome sequencing (skimming) projects either assemble the organelle genome or compute genomic distances directly from unassembled reads. Using unassembled reads needs contamination detection because samples often include reads from unintended groups of species. Similarly, assembling the organelle genome needs distinguishing organelle and nuclear reads. While k-mer-based methods have shown promise in read-matching, prior studies have shown that existing methods are insufficiently sensitive for contamination detection. Here, we introduce a new read-matching tool called CONSULT that tests whether k-mers from a query fall within a user-specified distance of the reference dataset using locality-sensitive hashing. Taking advantage of large memory machines available nowadays, CONSULT libraries accommodate tens of thousands of microbial species. Our results show that CONSULT has higher true-positive and lower false-positive rates of contamination detection than leading methods such as Kraken-II and improves distance calculation from genome skims. We also demonstrate that CONSULT can distinguish organelle reads from nuclear reads, leading to dramatic improvements in skim-based mitochondrial assemblies.

High diversity, inbreeding and a dynamic Pleistocene demographic history revealed by African buffalo genomes

Genomes retain records of demographic changes and evolutionary forces that shape species and populations. Remnant populations of African buffalo (Syncerus caffer) in South Africa, with varied histories, provide an opportunity to investigate signatures left in their genomes by past events, both recent and ancient. Here, we produce 40 low coverage (7.14×) genome sequences of Cape buffalo (S. c. caffer) from four protected areas in South Africa. Genome-wide heterozygosity was the highest for any mammal for which these data are available, while differences in individual inbreeding coefficients reflected the severity of historical bottlenecks and current census sizes in each population. PSMC analysis revealed multiple changes in Ne between approximately one million and 20 thousand years ago, corresponding to paleoclimatic changes and Cape buffalo colonisation of southern Africa. The results of this study have implications for buffalo management and conservation, particularly in the context of the predicted increase in aridity and temperature in southern Africa over the next century as a result of climate change.

Resolving the Mortierellaceae phylogeny through synthesis of multi-gene phylogenetics and phylogenomics.

Early efforts to classify Mortierellaceae were based on macro- and micromorphology, but sequencing and phylogenetic studies with ribosomal DNA (rDNA) markers have demonstrated conflicting taxonomic groupings and polyphyletic genera. Although some taxonomic confusion in the family has been clarified, rDNA data alone is unable to resolve higher level phylogenetic relationships within Mortierellaceae. In this study, we applied two parallel approaches to resolve the Mortierellaceae phylogeny: low coverage genome (LCG) sequencing and high-throughput, multiplexed targeted amplicon sequencing to generate sequence data for multi-gene phylogenetics. We then combined our datasets to provide a well-supported genome-based phylogeny having broad sampling depth from the amplicon dataset. Resolving the Mortierellaceae phylogeny into monophyletic groups led to the definition of 14 genera, 7 of which are newly proposed. Low-coverage genome sequencing proved to be a relatively cost-effective means of generating a well-resolved phylogeny. The multi-gene phylogenetics approach enabled much greater sampling depth and breadth than the LCG approach, but was unable to resolve higher-level organization of groups. We present this work to resolve some of the taxonomic confusion and provide a genus-level framework to empower future studies on Mortierellaceae diversity, biology, and evolution.