Search

Obtaining high-quality DNA suitable for long-read sequencing can be difficult for many types of tissues and cells, and it is a key step in current genomic studies. The challenge is even greater when it comes to isolating genomic DNA from mammalian spermatozoa, as DNA is tightly packed into a cell with a robust membrane rich in disulfide bonds. Here we describe a method for isolating high molecular weight DNA from Bovine commercial semen straws. This protocol includes a cleaning step to remove diluents and preservatives used for the long-term storage of the semen, which may affect long read sequencing. It is based on a simple salting-out method and avoid the use of spin columns, strong mixing or intensive centrifugation, in order to limit DNA fragmentation. However, we have adapted this protocol to facilitate the disruption of cell membranes and disulfide bonds with strong chaotropic and reducing agents. The average size of the fragments produced was approximately 49 kb, ranging from 25 to 85 kb, according to the femto pulse profiles.This method was used to isolate DNA from semen straws, more than 80 of them were successfully sequenced using the Continuous Long-Read (CLR) sequencing mode on the PacBio SequelII platform to study genome diversity and notably to detect large structural variations within genomes. [ABSTRACT FROM AUTHOR]

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and their widespread integration as a means of transport remains contentious2–4. Here we assemble a collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged around 2200 bce, through close-kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than approximately 2700 bce, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe around 3000 bce and earlier3,5. Finally, we detect significantly shortened generation times at Botai around 3500 bce, a settlement from central Asia associated with corrals and a subsistence economy centred on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.Analyses of 475 ancient horse genomes show modern horses emerged around 2200 bce, coinciding with sudden expansion across Eurasia, refuting the narrative of large horse herds accompanying earlier migrations of steppe peoples across Europe. [ABSTRACT FROM AUTHOR]

Background: The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses. Results: Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants. Conclusions: Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes. [ABSTRACT FROM AUTHOR]

Background: The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. Results: We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19–27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary. Conclusions: Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species. [ABSTRACT FROM AUTHOR]

Inspired by the production of reference data sets in the Genome in a Bottle project, we sequenced one Charolais heifer with different technologies: Illumina paired-end, Oxford Nanopore, Pacific Biosciences (HiFi and CLR), 10X Genomics linked-reads, and Hi-C. In order to generate haplotypic assemblies, we also sequenced both parents with short reads. From these data, we built two haplotyped trio high quality reference genomes and a consensus assembly, using up-to-date software packages. The assemblies obtained using PacBio HiFi reaches a size of 3.2 Gb, which is significantly larger than the 2.7 Gb ARS-UCD1.2 reference. The BUSCO score of the consensus assembly reaches a completeness of 95.8%, among highly conserved mammal genes. We also identified 35,866 structural variants larger than 50 base pairs. This assembly is a contribution to the bovine pangenome for the "Charolais" breed. These datasets will prove to be useful resources enabling the community to gain additional insight on sequencing technologies for applications such as SNP, indel or structural variant calling, and de novo assembly. [ABSTRACT FROM AUTHOR]

Of American origin, a wide diversity of Xylella fastidiosa strains belonging to different subspecies have been reported in Europe since 2013 and its discovery in Italian olive groves. Strains from the subspecies multiplex (ST6 and ST7) were first identified in France in 2015 in urban and natural areas. To trace back the most probable scenario of introduction in France, the molecular evolution rate of this subspecies was estimated at 3.2165 × 10-7 substitutions per site per year, based on heterochronous genome sequences collected worldwide. This rate allowed the dating of the divergence between French and American strains in 1987 for ST6 and in 1971 for ST7. The development of a new VNTR-13 scheme allowed tracing the spread of the bacterium in France, hypothesizing an American origin. Our results suggest that both sequence types were initially introduced and spread in Provence-Alpes-Côte d'Azur (PACA); then they were introduced in Corsica in two waves from the PACA bridgehead populations. A genomic investigation into strains of the plant pathogen, Xylella fastidiosa, in France indicate an American origin with two waves of invasion. [ABSTRACT FROM AUTHOR]

Fast, accurate sequencing methods are needed to identify new variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) genome. Single‐molecule real‐time (SMRT) Pacific Biosciences (PacBio) provides long, highly accurate sequences by circular consensus reads. This study compares the performance of a target capture SMRT PacBio protocol for whole‐genome sequencing (WGS) of SARS‐CoV‐2 to that of an amplicon PacBio SMRT sequencing protocol. The median genome coverage was higher (p < 0.05) with the target capture protocol (99.3% [interquartile range, IQR: 96.3–99.5]) than with the amplicon protocol (99.3% [IQR: 69.9–99.3]). The clades of 65 samples determined with both protocols were 100% concordant. After adjusting for Ct values, S gene coverage was higher with the target capture protocol than with the amplicon protocol. After stratification on Ct values, higher S gene coverage with the target capture protocol was observed only for samples with Ct > 17 (p < 0.01). PacBio SMRT sequencing protocols appear to be suitable for WGS, genotyping, and detecting mutations of SARS‐CoV‐2. [ABSTRACT FROM AUTHOR]

The Pacific halibut (Hippoglossus stenolepis) is a key species in the North Pacific Ocean and Bering Sea ecosystems, where it also supports important fisheries. However, the lack of genomic resources limits our understanding of evolutionary, environmental and anthropogenic forces affecting key life history characteristics of Pacific halibut and prevents the application of genomic tools in fisheries management and conservation efforts. In the present study, we report on the first generation of a high‐quality chromosome‐level assembly of the Pacific halibut genome, with an estimated size of 602 Mb, 24 chromosome‐length scaffolds that contain 99.8% of the assembly and a N50 scaffold length of 27.3 Mb. In the first application of this important resource, we conducted genome‐wide analyses of sex‐specific genetic variation by pool sequencing and characterized a potential sex‐determining region in chromosome 9 with a high density of female‐specific SNPs. Within this region, we identified the bmpr1ba gene as a potential candidate for master sex‐determining (MSD) gene. bmpr1ba is a member of the TGF‐β family that in teleosts has provided the largest number of MSD genes, including a paralogue of this gene in Atlantic herring. The genome assembly constitutes an essential resource for future studies on Pacific halibut population structure and dynamics, evolutionary history and responses to environmental and anthropogenic influences. Furthermore, the genomic location of the sex‐determining region in Pacific halibut has been identified and a putative candidate MSD gene has been proposed, providing further support for the rapid evolution of sex‐determining mechanisms in teleost fish. [ABSTRACT FROM AUTHOR]

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex‐determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome‐scale genome assembly, we identified an anti‐Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex‐linked in males and expressed only in testes. These results point to a Y chromosome male‐specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N‐terminal domain, lacking the cysteine‐rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference‐guided assembly of 11 additional Pangasiids, along with sex‐linkage studies, revealed that this truncated amhr2by duplication is a male‐specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex‐determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems. [ABSTRACT FROM AUTHOR]

Objective: Causal mutations for major genes that underlie a broad range of morphological traits are often located within exons of genes that then affect protein functions. Non-model organism genetic studies are not easy to perform due to the lack of genome-wide molecular tools such as SNP genotyping array. Genotyping-By-Sequencing (GBS) methods offer an alternative. Consequently, we used this approach that is focused on the exome to target and identify major genes in rabbit populations. Data description We used a heterologous enrichment method before sequencing, allowing us to capture the rabbit exome using the marketed human panel since mammal protein coding genes are well conserved across the phylogenic tree of species. This targeted strategy was performed on 52 French rabbits from 5 different French strains (Californian, New-Zealand, Castor, Chinchilla and Laghmere). We generated 3.4 billion sequencing reads and approximately 29–140 million of reads per DNA sample. The expected exome coverage per sample ranged between 118 and 566X. The present dataset could be useful for the scientific community working on rabbit species in order to (i) improve the annotation of the rabbit reference genome Oryctolagus cuniculus (OryCun2.0), (ii) enrich the characterization of polymorphisms segregating in rabbits and (iii) evaluate the genetic biodiversity in different rabbit strains. Raw sequences were deposited in the European Nucleotide Archive (ENA) at the European Molecular Biology Laboratory- European Bioinformatics Institute (EMBL-EBI) data portal under bioproject accession number PRJEB37917. [ABSTRACT FROM AUTHOR]

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long‐reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi‐C, we generated a high‐continuity chromosome‐scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome‐size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male‐specific duplicate of the anti‐Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex‐specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by+) from XX genetic females (amhr2by−). Our high‐quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex‐determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes. [ABSTRACT FROM AUTHOR]

In third generation sequencing, the production of quality data requires the selection of molecules longer than ∼20 kbp, but the size selection threshold of most purification technologies is smaller than this target. Here, we describe a technology operated in a capillary with a tunable selection threshold in the range of 3 to 40 kbp controlled by an electric field. We demonstrate that the selection cut-off is sharp, the purification yield is high, and the purification throughput is scalable. We also provide an analytical model that the actuation settings of the filter. The selection of high molecular weight genomic DNA from the melon Cucumis melo L., a diploid organism of ∼0.45 Gbp, is then reported. Linked-read sequencing data show that the N50 phase block size, which scores the correct representation of two chromosomes, is enhanced by a factor of 2 after size selection, establishing the relevance and versatility of our technology. [ABSTRACT FROM AUTHOR]

Mitochondrial genetic markers have been extensively used to study the phylogenetics and phylogeography of many birds, including seabirds of the order Procellariiformes. Evidence suggests that part of the mitochondrial genome of Procellariiformes, especially albatrosses, is duplicated, but no DNA fragment covering the entire duplication has been sequenced. We sequenced the complete mitochondrial genome of a non-albatross species of Procellariiformes, Puffinus lherminieri (Audubon's shearwater) using the long-read MinION (ONT) technology. Two mitogenomes were assembled from the same individual, differing by 52 SNPs and in length. The shorter was 19 kb long while the longer was 21 kb, due to the presence of two identical copies of nad6, three tRNA, and two dissimilar copies of the control region (CR). Contrary to albatrosses, cob was not duplicated. We further detected a complex repeated region of undetermined length between the CR and 12S. Long-read sequencing suggests heteroplasmy and a novel arrangement within the duplicated region, indicating a complex evolution of the mitogenome in Procellariiformes. [ABSTRACT FROM AUTHOR]

Background: Substantial progress in high-throughput metagenomic sequencing methodologies has enabled the characterisation of bacteria from various origins (for example gut and skin). However, the recently-discovered bacterial microbiota present within animal internal tissues has remained unexplored due to technical difficulties associated with these challenging samples. Results: We have optimized a specific 16S rDNA-targeted metagenomics sequencing (16S metabarcoding) pipeline based on the Illumina MiSeq technology for the analysis of bacterial DNA in human and animal tissues. This was successfully achieved in various mouse tissues despite the high abundance of eukaryotic DNA and PCR inhibitors in these samples. We extensively tested this pipeline on mock communities, negative controls, positive controls and tissues and demonstrated the presence of novel tissue specific bacterial DNA profiles in a variety of organs (including brain, muscle, adipose tissue, liver and heart). Conclusion: The high throughput and excellent reproducibility of the method ensured exhaustive and precise coverage of the 16S rDNA bacterial variants present in mouse tissues. This optimized 16S metagenomic sequencing pipeline will allow the scientific community to catalogue the bacterial DNA profiles of different tissues and will provide a database to analyse host/bacterial interactions in relation to homeostasis and disease. [ABSTRACT FROM AUTHOR]

For decades, French guinea fowl have been affected by fulminating enteritis of unclear origin. By using metagenomics, we identified a novel avian gammacoronavirus associated with this disease that is distantly related to turkey coronaviruses. Fatal respiratory diseases in humans have recently been caused by coronaviruses of animal origin. [ABSTRACT FROM AUTHOR]

Background: Single Nucleotide Polymorphisms (SNPs) can be used as genetic markers for applications such as genetic diversity studies or genetic mapping. New technologies now allow genotyping hundreds to thousands of SNPs in a single reaction. In order to evaluate the potential of these technologies in pea, we selected a custom 384-SNP set using SNPs discovered in Pisum through the resequencing of gene fragments in different genotypes and by compiling genomic sequence data present in databases. We then designed an Illumina GoldenGate assay to genotype both a Pisum germplasm collection and a genetic mapping population with the SNP set. Results: We obtained clear allelic data for more than 92% of the SNPs (356 out of 384). Interestingly, the technique was successful for all the genotypes present in the germplasm collection, including those from species or subspecies different from the P. sativum ssp sativum used to generate sequences. By genotyping the mapping population with the SNP set, we obtained a genetic map and map positions for 37 new gene markers. Conclusion: Our results show that the Illumina GoldenGate assay can be used successfully for high-throughput SNP genotyping of diverse germplasm in pea. This genotyping approach will simplify genotyping procedures for association mapping or diversity studies purposes and open new perspectives in legume genomics. [ABSTRACT FROM AUTHOR]

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is the causal agent of the COVID-19 pandemic that emerged in late 2019. The outbreak of variants with mutations in the region encoding the spike protein S1 sub-unit that can make them more resistant to neutralizing or monoclonal antibodies is the main point of the current monitoring. This study examines the feasibility of predicting the variant lineage and monitoring the appearance of reported mutations by sequencing only the region encoding the S1 domain by Pacific Bioscience Single Molecule Real-Time sequencing (PacBio SMRT). Using the PacBio SMRT system, we successfully sequenced 186 of the 200 samples previously sequenced with the Illumina COVIDSeq (whole genome) system. PacBio SMRT detected mutations in the S1 domain that were missed by the COVIDseq system in 27/186 samples (14.5%), due to amplification failure. These missing positions included mutations that are decisive for lineage assignation, such as G142D (n = 11), N501Y (n = 6), or E484K (n = 2). The lineage of 172/186 (92.5%) samples was accurately determined by analyzing the region encoding the S1 domain with a pipeline that uses key positions in S1. Thus, the PacBio SMRT protocol is appropriate for determining virus lineages and detecting key mutations. [ABSTRACT FROM AUTHOR]

The spread of SARS-CoV-2 and the resulting disease COVID-19 has killed over 2.6 million people as of 18 March 2021. We have used a modified susceptible, infected, recovered (SIR) epidemiological model to predict how the spread of the virus in regions of France will vary depending on the proportions of variants and on the public health strategies adopted, including anti-COVID-19 vaccination. The proportion of SARS-CoV-2 variant B.1.1.7, which was not detected in early January, increased to become 60% of the forms of SARS-CoV-2 circulating in the Toulouse urban area at the beginning of February 2021, but there was no increase in positive nucleic acid tests. Our prediction model indicates that maintaining public health measures and accelerating vaccination are efficient strategies for the sustained control of SARS-CoV-2. [ABSTRACT FROM AUTHOR]