scholarly journals Chromosome-level assembly of the Atlantic silverside genome reveals extreme levels of sequence diversity and structural genetic variation

2020 ◽  
Author(s):  
Anna Tigano ◽  
Arne Jacobs ◽  
Aryn P. Wilder ◽  
Ankita Nand ◽  
Ye Zhan ◽  
...  

AbstractThe levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10X Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia) - an established ecological model for studying the phenotypic effects of natural and artificial selection - and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is the highest reported for a fish, or any vertebrate, to date (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ~23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

2021 ◽  
Author(s):  
Lydia R. Heasley ◽  
Juan Lucas Argueso

The budding yeast Saccharomyces cerevisiae has been extensively characterized for many decades and is a critical resource for the study of numerous facets of eukaryotic biology. Recently, the analysis of whole genome sequencing data from over 1000 natural isolates of S. cerevisiae has provided critical insights into the evolutionary landscape of this species by revealing a population structure comprised of numerous genomically diverse lineages. These survey-level analyses have been largely devoid of structural genomic information, mainly because short read sequencing is not suitable for detailed characterization of genomic architecture. Consequently, we still lack a complete perspective of the genomic variation the exists within the species. Single molecule long read sequencing technologies, such as Oxford Nanopore and PacBio, provide sequencing-based approaches with which to rigorously define the structure of a genome, and have empowered yeast geneticists to explore this poorly described realm of eukaryotic genomics. Here, we present the comprehensive genomic structural analysis of a pathogenic isolate of S. cerevisiae, YJM311. We used long read sequence analysis to construct a haplotype-phased, telomere-to-telomere length assembly of the YJM311 diploid genome and characterized the structural variations (SVs) therein. We discovered that the genome of YJM311 contains significant intragenomic structural variation, some of which imparts notable consequences to the genomic stability and developmental biology of the strain. Collectively, we outline a new methodology for creating accurate haplotype-phased genome assemblies and highlight how such genomic analyses can define the structural architectures of S. cerevisiae isolates. It is our hope that through continued structural characterization of S. cerevisiae genomes, such as we have reported here for YJM311, we will comprehensively advance our understanding of eukaryotic genome structure-function relationships, structural diversity, and evolution.


Genetics ◽  
2021 ◽  
Author(s):  
Lydia R Heasley ◽  
Juan Lucas Argueso

Abstract The budding yeast Saccharomyces cerevisiae has been extensively characterized for many decades and is a critical resource for the study of numerous facets of eukaryotic biology. Recently, whole genome sequence analysis of over 1000 natural isolates of S. cerevisiae has provided critical insights into the evolutionary landscape of this species by revealing a population structure comprised of numerous genomically diverse lineages. These survey-level analyses have been largely devoid of structural genomic information, mainly because short read sequencing is not suitable for detailed characterization of genomic architecture. Consequently, we still lack a complete perspective of the genomic variation the exists within the species. Single molecule long read sequencing technologies, such as Oxford Nanopore and PacBio, provide sequencing-based approaches with which to rigorously define the structure of a genome, and have empowered yeast geneticists to explore this poorly described realm of eukaryotic genomics. Here, we present the comprehensive genomic structural analysis of a wild diploid isolate of S. cerevisiae, YJM311. We used long read sequence analysis to construct a haplotype-phased, telomere-to-telomere length assembly of the YJM311 genome and characterized the structural variations (SVs) therein. We discovered that the genome of YJM311 contains significant intragenomic structural variation, some of which imparts notable consequences to the genomic stability and developmental biology of the strain. Collectively, we outline a new methodology for creating accurate haplotype-phased genome assemblies and highlight how such genomic analyses can define the structural architectures of S. cerevisiae isolates. It is our hope that continued structural characterization of S. cerevisiae genomes, such as we have reported here for YJM311, will comprehensively advance our understanding of eukaryotic genome structure-function relationships, structural genomic diversity, and evolution.


2018 ◽  
Author(s):  
Sojung Han ◽  
Aida M. Andrés ◽  
Tomas Marques-Bonet ◽  
Martin Kuhlwilm

AbstractChimpanzees (Pan troglodytes) and bonobos (Pan paniscus) are the closest living relatives of humans, but they show distinct behavioral and physiological differences, particularly regarding female reproduction. Despite their recent rapid decline, the demographic histories of the two species have been different during the past one to two million years, likely having an impact on their genomic diversity. Here, we analyze the inferred functional consequences of genetic variation across 69 individuals, making use of the most complete dataset of genomic variation in the Pan clade to date. We test to which extent the demographic history influences the efficacy of purifying selection in these species. We find that small historical effective population sizes (Ne) correlate not only with small genetic diversity, but also with more homozygous deleterious alleles, and an increased proportion of deleterious changes at low frequencies. Furthermore, we exploit the catalog of deleterious protein-coding changes on each lineage to investigate the putative genetic basis for phenotypic differences between chimpanzees and bonobos. We show that bonobo-specific non-synonymous changes are enriched in genes related to age at menarche in humans, suggesting that the prominent physiological differences in the female reproductive system between chimpanzees and bonobos might be explained, in part, by putatively adaptive changes on the bonobo lineage.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Kelly B. Klingler ◽  
Joshua P. Jahner ◽  
Thomas L. Parchman ◽  
Chris Ray ◽  
Mary M. Peacock

Abstract Background Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada. Results Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θW = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θW = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θW = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide. Conclusions Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.


2020 ◽  
Author(s):  
Kelly Brie Klingler ◽  
Joshua P Jahner ◽  
Thomas L Parchman ◽  
Chris Ray ◽  
Mary Peacock

Abstract Background: Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada. Results: Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (𝜋=0.0006–0.0009; 𝜃W=0.0005–0.0007) relative to populations in California (𝜋=0.0014–0.0019; 𝜃W=0.0011–0.0017) and the Rocky Mountains (𝜋=0.0025–0.0027; 𝜃W=0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D=0.240-0.811), consistent with recent contraction in population sizes range-wide. Conclusions: Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.


2021 ◽  
Vol 8 (3) ◽  
Author(s):  
M. V. Westbury ◽  
K. F. Thompson ◽  
M. Louis ◽  
A. A. Cabrera ◽  
M. Skovrind ◽  
...  

The deep oceans of the Southern Hemisphere are home to several elusive and poorly studied marine megafauna. In the absence of robust observational data for these species, genetic data can aid inferences on population connectivity, demography and ecology. A previous investigation of genetic diversity and population structure in Gray's beaked whale ( Mesoplodon grayi ) from Western Australia and New Zealand found high levels of mtDNA diversity, no geographic structure and stable demographic history. To further investigate phylogeographic and demographic patterns across their range, we generated complete mitochondrial and partial nuclear genomes of 16 of the individuals previously analysed and included additional samples from South Africa ( n = 2) and South Australia ( n = 4), greatly expanding the spatial range of genomic data for the species. Gray's beaked whales are highly elusive and rarely observed, and our data represents a unique and geographically broad dataset. We find relatively high levels of diversity in the mitochondrial genome, despite an absence of population structure at the mitochondrial and nuclear level. Demographic analyses suggest these whales existed at stable levels over at least the past 1.1 million years, with an approximately twofold increase in female effective population size approximately 250 thousand years ago, coinciding with a period of increased Southern Ocean productivity, sea surface temperature and a potential expansion of suitable habitat. Our results suggest that Gray's beaked whales are likely to be resilient to near-future ecosystem changes, facilitating their conservation. Our study demonstrates the utility of low-effort shotgun sequencing in providing ecological information on highly elusive species.


Author(s):  
Anthony M. Gacita ◽  
Dominic E. Fullenkamp ◽  
Joyce Ohiri ◽  
Tess Pottinger ◽  
Megan J. Puckelwartz ◽  
...  

Background: Inherited cardiomyopathy associates with a range of phenotype, mediated by genetic and non-genetic factors. Non-inherited cardiomyopathy also displays varying progression and outcomes. Expression of cardiomyopathy genes is under the regulatory control of promoters and enhancers, and human genetic variation in promoters and enhancers may contribute to this variability. Methods: We superimposed epigenomic profiling from hearts and cardiomyocytes, including promoter-capture chromatin conformation information, to identify enhancers for two cardiomyopathy genes, MYH7 and LMNA . Enhancer function was validated in human cardiomyocytes derived from induced pluripotent stem cells. We also conducted a genome-wide search to ascertain genomic variation in enhancers positioned to alter cardiac expression and correlated one of these variants to cardiomyopathy progression using biobank data. Results: Multiple enhancers were identified and validated for LMNA and MYH7 , including a key enhancer that regulates the switch from MYH6 expression to MYH7 expression. Deletion of this enhancer resulted in a dose-dependent increase in MYH6 and faster contractile rate in engineered heart tissues. We searched for genomic variation in enhancer sequences across the genome, with focus on nucleotide changes that create or interrupt transcription factor binding sites. rs875908 disrupts a TBX5 binding motif and maps to an enhancer region 2KB from the transcriptional start site of MYH7 . Gene editing to remove the enhancer harboring this variant markedly reduced MYH7 expression in human cardiomyocytes. Using biobank-derived data, rs875908 associated with longitudinal echocardiographic features with cardiomyopathy. Conclusions: Enhancers regulate cardiomyopathy gene expression, and genomic variation within these enhancer regions associates with cardiomyopathic progression over time. This integrated approach identified noncoding modifiers of cardiomyopathy and is applicable to other cardiac genes.


2020 ◽  
Author(s):  
Tlou Caswell Chokoe ◽  
Khanyi Hadebe ◽  
Farai Muchadeyi ◽  
Khathutshelo Nephawe ◽  
Edgar Dzomba ◽  
...  

Abstract Background: Indigenous goats forms the majority of populations in smallholder; low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support communal farming. Effective population size (Ne), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Results: Across populations, the current Ne of Gauteng was the lowest at 371 animals, while the historical Ne across populations suggests that the ancestor Ne has decreased by 53.86%, 44.58%, 42.16% and 41.16% in Free State (FS), North West (NW), Limpopo (LP) and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (FROH >5Mb) was highest in FS (0.08±0.09) and lowest for Eastern Cape (EC) (0.02±0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), NW (n = 2) populations only. The mean length of ROH across populations was 7.76Mb and ranged from 1.61Mb KwaZulu-Natal (KZN) to 98.05Mb (GP and NW). Distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were >5Mb in length than the other categories. Assuming two hypothetical ancestral populations, the population from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP shared an origin but have substantial admixture from the EC and NW populations.Conclusions: These findings indicated a greater negative impact of inbreeding in recent times which is important for planning conservation strategies. It was revealed that the occurrence of high Ne and autozygosity varied largely across breeds in communal indigenous goat populations at different recent and ancient events when a genome-wide SNP marker was used.


Genes ◽  
2021 ◽  
Vol 12 (9) ◽  
pp. 1342
Author(s):  
Francesca Bernini ◽  
Alessandro Bagnato ◽  
Stefano Paolo Marelli ◽  
Luisa Zaniboni ◽  
Silvia Cerolini ◽  
...  

Italian autochthonous turkey breeds are an important reservoir of genetic biodiversity that should be maintained with an in vivo approach. The aim of this study, part of the TuBAvI national project on biodiversity, was to use run of homozygosity (ROH), together with others statistical approaches (e.g., Wright’s F-statistics, principal component analysis, ADMIXTURE analysis), to investigate the genomic diversity in several heritage turkey breeds. We performed a genome-wide characterization of ROH-rich regions in seven autochthonous turkey breeds, i.e., Brianzolo (Brzl), Bronzato Comune Italiano (BrCI), Bronzato dei Colli Euganei (CoEu), Parma e Piacenza (PrPc), Nero d’Italia (NeIt), Ermellinato di Rovigo (ErRo) and Romagnolo (Roma). ROHs were detected based on a 650K SNP genotyping. ROH_islands were identified as homozygous ROH regions shared by at least 75% of birds (within breed). Annotation of genes was performed with DAVID. The admixture analyses revealed that six breeds are unique populations while the Roma breed consists in an admixture of founder populations. Effective population size estimated on genomic data shows a numeric contraction. ROH_islands harbour genes that may be interesting for target selection in commercial populations also. Among them the PTGS2 and PLA2G4A genes on chr10 were related to reproduction efficiency. This is the first study mapping genetic variation in autochthonous turkey populations. Breeds were genetically different among them, with the Roma breed proving to be a mixture of the other breeds. The ROH_islands identified harboured genes peculiar to the selection that occurred in heritage breeds. Finally, this study releases previously undisclosed information on existing genetic variation in the turkey species.


Sign in / Sign up

Export Citation Format

Share Document