Widespread horse-based mobility arose around 2,200 BCE in Eurasia.

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and widespread integration as a means of transportation remains contentious2-4. Here we assemble a large 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 ~2,200 BCE (Before Common Era), through close kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than ~2,700 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 ~3,000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai ~3,500 BCE, a settlement from Central Asia associated with corrals and a subsistence economy centered on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

Large-scale detection and characterization of interchromosomal rearrangements in normozoospermic bulls using massive genotype and phenotype data sets.

In this paper, we developed a highly sensitive approach to detect interchromosomal rearrangements in cattle by searching for abnormal linkage disequilibrium patterns between markers located on different chromosomes in large paternal half-sib families genotyped as part of routine genomic evaluations. We screened 5571 families of artificial insemination sires from 15 breeds and revealed 13 putative interchromosomal rearrangements, 12 of which were validated by cytogenetic analysis and long-read sequencing. These consisted of one Robertsonian fusion, 10 reciprocal translocations, and the first case of insertional translocation reported in cattle. Taking advantage of the wealth of data available in cattle, we performed a series of complementary analyses to define the exact nature of these rearrangements, investigate their origins, and search for factors that may have favored their occurrence. We also evaluated the risks to the livestock industry and showed significant negative effects on several traits in the sires and in their balanced or aneuploid progeny compared with wild-type controls. Thus, we present the most comprehensive and thorough screen for interchromosomal rearrangements compatible with normal spermatogenesis in livestock species. This approach is readily applicable to any population that benefits from large genotype data sets, and will have direct applications in animal breeding. Finally, it also offers interesting prospects for basic research by allowing the detection of smaller and rarer types of chromosomal rearrangements than GTG banding, which are interesting models for studying gene regulation and the organization of genome structure.

The genomics of linkage drag in inbred lines of sunflower.

Crop wild relatives represent valuable sources of alleles for crop improvement, including adaptation to climate change and emerging diseases. However, introgressions from wild relatives might have deleterious effects on desirable traits, including yield, due to linkage drag. Here, we analyzed the genomic and phenotypic impacts of wild introgressions in inbred lines of cultivated sunflower to estimate the impacts of linkage drag. First, we generated reference sequences for seven cultivated and one wild sunflower genotype, as well as improved assemblies for two additional cultivars. Next, relying on previously generated sequences from wild donor species, we identified introgressions in the cultivated reference sequences, as well as the sequence and structural variants they contain. We then used a ridge-regression best linear unbiased prediction (BLUP) model to test the effects of the introgressions on phenotypic traits in the cultivated sunflower association mapping population. We found that introgression has introduced substantial sequence and structural variation into the cultivated sunflower gene pool, including >3,000 new genes. While introgressions reduced genetic load at protein-coding sequences, they mostly had negative impacts on yield and quality traits. Introgressions found at high frequency in the cultivated gene pool had larger effects than low-frequency introgressions, suggesting that the former likely were targeted by artificial selection. Also, introgressions from more distantly related species were more likely to be maladaptive than those from the wild progenitor of cultivated sunflower. Thus, breeding efforts should focus, as far as possible, on closely related and fully compatible wild relatives.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae.

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.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes.

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.

Size Fractionation of Milliliter DNA Samples in Minutes Controlled by an Electric Field of ∼10 V.

DNA size fractionation is an essential tool in molecular biology and is used to isolate targets in a mixture characterized by a broad molecular-weight distribution. Microfluidics was thought to provide the opportunity to create devices capable of enhancing and speeding up the classical fractionation processes. However, this conjecture met limited success due to the low mass or volume throughput of these technologies. We describe the µLAF (µ-laboratory for DNA fractionation) technology for DNA size selection based on the stacking of molecules on films of ∼100 µm in thickness with 105 cm-2 pores ∼2 µm in diameter. Size selection is achieved by controlling the regime of electrohydrodynamic migration through the temporal modulation of an electric field. This technology allows the processing of milliliter-scale samples containing a DNA mass of several hundreds of ng within ∼10 min and the selection of DNA in virtually any size window spanning 200 to 1000 bp. We demonstrate that one operation suffices to fractionate sheared genomic DNA in up to six fractions with collection efficiencies of ∼20-40% and enrichment factors of ∼1.5-3-fold. These performances compare favorably in terms of speed and versatility to those of the current standards.

Whole-genome sequencing of SARS-CoV-2: Comparison of target capture and amplicon single molecule real-time sequencing protocols.

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.

The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution.

The domesticated sunflower, Helianthus annuus L., is a global oil crop that has promise for climate change adaptation, because it can maintain stable yields across a wide variety of environmental conditions, including drought. Even greater resilience is achievable through the mining of resistance alleles from compatible wild sunflower relatives, including numerous extremophile species. Here we report a high-quality reference for the sunflower genome (3.6 gigabases), together with extensive transcriptomic data from vegetative and floral organs. The genome mostly consists of highly similar, related sequences and required single-molecule real-time sequencing technologies for successful assembly. Genome analyses enabled the reconstruction of the evolutionary history of the Asterids, further establishing the existence of a whole-genome triplication at the base of the Asterids II clade and a sunflower-specific whole-genome duplication around 29 million years ago. An integrative approach combining quantitative genetics, expression and diversity data permitted development of comprehensive gene networks for two major breeding traits, flowering time and oil metabolism, and revealed new candidate genes in these networks. We found that the genomic architecture of flowering time has been shaped by the most recent whole-genome duplication, which suggests that ancient paralogues can remain in the same regulatory networks for dozens of millions of years. This genome represents a cornerstone for future research programs aiming to exploit genetic diversity to improve biotic and abiotic stress resistance and oil production, while also considering agricultural constraints and human nutritional needs.

Analysis of Polycerate Mutants Reveals the Evolutionary Co-option of HOXD1 for Horn Patterning in Bovidae.

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as "headgear," which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

A Chromosome-Scale Genome Assembly Resource for Myriosclerotinia sulcatula Infecting Sedge Grass (Carex sp.).

The fungus Myriosclerotinia sulcatula is a close relative of the notorious polyphagous plant pathogens Botrytis cinerea and Sclerotinia sclerotiorum but exhibits a host range restricted to plants from the Carex genus (Cyperaceae family). To date, there are no genomic resources available for fungi in the Myriosclerotinia genus. Here, we present a chromosome-scale reference genome assembly for M. sulcatula. The assembly contains 24 contigs with a total length of 43.53 Mbp, with scaffold N50 of 2,649.7 kbp and N90 of 1,133.1 kbp. BRAKER-predicted gene models were manually curated using WebApollo, resulting in 11,275 protein-coding genes that we functionally annotated. We provide a high-quality reference genome assembly and annotation for M. sulcatula as a resource for studying evolution and pathogenicity in fungi from the Sclerotiniaceae family.

Novel avian coronavirus and fulminating disease in guinea fowl, France.

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.

Detection of DNA methylation signatures through the lens of genomic imprinting.

Genomic imprinting represents an original model of epigenetic regulation resulting in a parent-of-origin expression. Despite the critical role of imprinted genes in mammalian growth, metabolism and neuronal function, there is no molecular tool specifically targeting them for a systematic evaluation. We show here that enzymatic methyl-seq consistently outperforms the bisulfite-based standard in capturing 165 candidate regions for genomic imprinting in the pig. This highlights the potential for a turnkey, fully customizable and reliable capture tool of genomic regions regulated by cytosine methylation in any population of interest. For the field of genomic imprinting, it opens up the possibility of detecting multilocus imprinting variations across the genome, with implications for basic research, agrigenomics and clinical practice.

Imputed genomes of historical horses provide insights into modern breeding.

Historical genomes can provide important insights into recent genomic changes in horses, especially the development of modern breeds. In this study, we characterized 8.7 million genomic variants from a panel of 430 horses from 73 breeds, including newly sequenced genomes from 20 Clydesdales and 10 Shire horses. We used this modern genomic variation to impute the genomes of four historically important horses, consisting of publicly available genomes from 2 Przewalski's horses, 1 Thoroughbred, and a newly sequenced Clydesdale. Using these historical genomes, we identified modern horses with higher genetic similarity to those in the past and unveiled increased inbreeding in recent times. We genotyped variants associated with appearance and behavior to uncover previously unknown characteristics of these important historical horses. Overall, we provide insights into the history of Thoroughbred and Clydesdale breeds and highlight genomic changes in the endangered Przewalski's horse following a century of captive breeding.

Suspicions of two bridgehead invasions of Xylella fastidiosa subsp. multiplex in France.

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 Bos taurus sequencing methods benchmark for assembly, haplotyping, and variant calling.

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.

Genome structures resolve the early diversification of teleost fishes.

Accurate species phylogenies are a prerequisite for all evolutionary research. Teleosts are the largest and most diversified group of extant vertebrates, but relationships among their three oldest extant lineages remain unresolved. On the basis of seven high-quality new genome assemblies in Elopomorpha (tarpons, eels), we revisited the topology of the deepest branches of the teleost phylogeny using independent gene sequence and chromosomal rearrangement phylogenomic approaches. These analyses converged to a single scenario that unambiguously places the Elopomorpha and Osteoglossomorpha (arapaima, elephantnose fish) in a monophyletic sister group to all other teleosts, i.e., the Clupeocephala lineage (zebrafish, medaka). This finding resolves more than 50 years of controversy on the evolutionary relationships of these lineages and highlights the power of combining different levels of genome-wide information to solve complex phylogenies.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae.

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. 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 duplications (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 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 variants (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 testis than ovary. 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.

Rabbit targeted genomic sequences after heterologous hybridization using human exome.

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.

Generation of a chromosome-level genome assembly for Pacific halibut (Hippoglossus stenolepis) and characterization of its sex-determining genomic region.

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.

Transcriptional response to host chemical cues underpins the expansion of host range in a fungal plant pathogen lineage.

The host range of parasites is an important factor in assessing the dynamics of disease epidemics. The evolution of pathogens to accommodate new hosts may lead to host range expansion, a process the molecular bases of which are largely enigmatic. The fungus Sclerotinia sclerotiorum has been reported to parasitize more than 400 plant species from diverse eudicot families while its close relative, S. trifoliorum, is restricted to plants from the Fabaceae family. We analyzed S. sclerotiorum global transcriptome reprogramming on hosts from six botanical families and reveal a flexible, host-specific transcriptional program. We generated a chromosome-level genome assembly for S. trifoliorum and found near-complete gene space conservation in two representative strains of broad and narrow host range Sclerotinia species. However, S. trifoliorum showed increased sensitivity to the Brassicaceae defense compound camalexin. Comparative analyses revealed a lack of transcriptional response to camalexin in the S. trifoliorum strain and suggest that regulatory variation in detoxification and effector genes at the population level may associate with the genetic accommodation of Brassicaceae in the Sclerotinia host range. Our work proposes transcriptional plasticity and the co-existence of signatures for generalist and polyspecialist adaptive strategies in the genome of a plant pathogen.