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.

Genome sequence and annotation of Periconia digitata a hopeful biocontrol agent of phytopathogenic oomycetes.

The Periconia fungal genus belongs to the phylum Ascomycota, order Pleosporales, family Periconiaceae. Periconia are found in many habitats, but little is known about their ecology. Several species from this genus produce bioactive molecules. Periconia digitata extracts were shown to be deadly active against the pine wilt nematode. Furthermore, P. digitata was shown to inhibit the plant pathogenic oomycete Phytophthora parasitica. Because P. digitata has great potential as a biocontrol agent and high quality genomic resources are still lacking in the Periconiaceae family, we generated long-read genomic data for P. digitata. Using PacBio Hifi sequencing technology, we obtained a highly-contiguous genome assembled in 13 chromosomes and totaling ca. 39 Mb. In addition, we produced a reference transcriptome, based on 12 different culture conditions, and proteomic data to support the genome annotation. Besides representing a new reference genome within the Periconiaceae, this work will contribute to our better understanding of the Eukaryotic tree of life and opens new possibilities in terms of biotechnological applications.

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.

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.

Comparison of two molecular barcodes for the study of equine strongylid communities with amplicon sequencing.

Basic knowledge on the biology and epidemiology of equine strongylid species still needs to be improved to contribute to the design of better parasite control strategies. Nemabiome metabarcoding is a convenient tool to quantify and identify species in bulk samples that could overcome the hurdle that cyathostomin morphological identification represents. To date, this approach has relied on the internal transcribed spacer 2 (ITS-2) of the ribosomal RNA gene, with a limited investigation of its predictive performance for cyathostomin communities. Using DNA pools of single cyathostomin worms, this study aimed to provide the first elements to compare performances of the ITS-2 and a cytochrome c oxidase subunit I (COI) barcode newly developed in this study. Barcode predictive abilities were compared across various mock community compositions of two, five and 11 individuals from distinct species. The amplification bias of each barcode was estimated. Results were also compared between various types of biological samples, i.e., eggs, infective larvae or adults. Bioinformatic parameters were chosen to yield the closest representation of the cyathostomin community for each barcode, underscoring the need for communities of known composition for metabarcoding purposes. Overall, the proposed COI barcode was suboptimal relative to the ITS-2 rDNA region, because of PCR amplification biases, reduced sensitivity and higher divergence from the expected community composition. Metabarcoding yielded consistent community composition across the three sample types. However, imperfect correlations were found between relative abundances from infective larvae and other life-stages for Cylicostephanus species using the ITS-2 barcode. While the results remain limited by the considered biological material, they suggest that additional improvements are needed for both the ITS-2 and COI barcodes.

A supernumerary "B-sex" chromosome drives male sex determination in the Pachón cavefish, Astyanax mexicanus.

Sex chromosomes are generally derived from a pair of classical type-A chromosomes, and relatively few alternative models have been proposed up to now.1,2 B chromosomes (Bs) are supernumerary and dispensable chromosomes with non-Mendelian inheritance found in many plant and animal species3,4 that have often been considered as selfish genetic elements that behave as genome parasites.5,6 The observation that in some species Bs can be either restricted or predominant in one sex7-14 raised the interesting hypothesis that Bs could play a role in sex determination.15 The characterization of putative B master sex-determining (MSD) genes, however, has not yet been provided to support this hypothesis. Here, in Astyanax mexicanus cavefish originating from Pachón cave, we show that Bs are strongly male predominant. Based on a high-quality genome assembly of a B-carrying male, we characterized the Pachón cavefish B sequence and found that it contains two duplicated loci of the putative MSD gene growth differentiation factor 6b (gdf6b). Supporting its role as an MSD gene, we found that the Pachón cavefish gdf6b gene is expressed specifically in differentiating male gonads, and that its knockout induces male-to-female sex reversal in B-carrying males. This demonstrates that gdf6b is necessary for triggering male sex determination in Pachón cavefish. Altogether these results bring multiple and independent lines of evidence supporting the conclusion that the Pachón cavefish B is a "B-sex" chromosome that contains duplicated copies of the gdf6b gene, which can promote male sex determination in this species.

Correction to: Quantitative trait loci for resistance to Flavobacterium psychrophilum in rainbow trout: effect of the mode of infection and evidence of epistatic interactions.

After publication of this work [1], we noted that there was an error in Table 3 Line 4.

A High-Quality Grapevine Downy Mildew Genome Assembly Reveals Rapidly Evolving and Lineage-Specific Putative Host Adaptation Genes.

Downy mildews are obligate biotrophic oomycete pathogens that cause devastating plant diseases on economically important crops. Plasmopara viticola is the causal agent of grapevine downy mildew, a major disease in vineyards worldwide. We sequenced the genome of Pl. viticola with PacBio long reads and obtained a new 92.94 Mb assembly with high contiguity (359 scaffolds for a N50 of 706.5 kb) due to a better resolution of repeat regions. This assembly presented a high level of gene completeness, recovering 1,592 genes encoding secreted proteins involved in plant-pathogen interactions. Plasmopara viticola had a two-speed genome architecture, with secreted protein-encoding genes preferentially located in gene-sparse, repeat-rich regions and evolving rapidly, as indicated by pairwise dN/dS values. We also used short reads to assemble the genome of Plasmopara muralis, a closely related species infecting grape ivy (Parthenocissus tricuspidata). The lineage-specific proteins identified by comparative genomics analysis included a large proportion of RxLR cytoplasmic effectors and, more generally, genes with high dN/dS values. We identified 270 candidate genes under positive selection, including several genes encoding transporters and components of the RNA machinery potentially involved in host specialization. Finally, the Pl. viticola genome assembly generated here will allow the development of robust population genomics approaches for investigating the mechanisms involved in adaptation to biotic and abiotic selective pressures in this species.

Quantitative trait loci for resistance to Flavobacterium psychrophilum in rainbow trout: effect of the mode of infection and evidence of epistatic interactions.

BACKGROUND: Bacterial cold-water disease, which is caused by Flavobacterium psychrophilum, is one of the major diseases that affect rainbow trout (Oncorhynchus mykiss) and a primary concern for trout farming. Better knowledge of the genetic basis of resistance to F. psychrophilum would help to implement this trait in selection schemes and to investigate the immune mechanisms associated with resistance. Various studies have revealed that skin and mucus may contribute to response to infection. However, previous quantitative trait loci (QTL) studies were conducted by using injection as the route of infection. Immersion challenge, which is assumed to mimic natural infection by F. psychrophilum more closely, may reveal different defence mechanisms. RESULTS: Two isogenic lines of rainbow trout with contrasting susceptibilities to F. psychrophilum were crossed to produce doubled haploid F2 progeny. Fish were infected with F. psychrophilum either by intramuscular injection (115 individuals) or by immersion (195 individuals), and genotyped for 9654 markers using RAD-sequencing. Fifteen QTL associated with resistance traits were detected and only three QTL were common between the injection and immersion. Using a model that accounted for epistatic interactions between QTL, two main types of interactions were revealed. A "compensation-like" effect was detected between several pairs of QTL for the two modes of infection. An "enhancing-like" interaction effect was detected between four pairs of QTL. Integration of the QTL results with results of a previous transcriptomic analysis of response to F. psychrophilum infection resulted in a list of potential candidate immune genes that belong to four relevant functional categories (bacterial sensors, effectors of antibacterial immunity, inflammatory factors and interferon-stimulated genes). CONCLUSIONS: These results provide new insights into the genetic determinism of rainbow trout resistance to F. psychrophilum and confirm that some QTL with large effects are involved in this trait. For the first time, the role of epistatic interactions between resistance-associated QTL was evidenced. We found that the infection protocol used had an effect on the modulation of defence mechanisms and also identified relevant immune functional candidate genes.

The Characterization of Novel Tissue Microbiota Using an Optimized 16S Metagenomic Sequencing Pipeline.

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.

Genome-wide translational changes induced by the prion [PSI+].

Prions are infectious proteins that can adopt a structural conformation that is then propagated among other molecules of the same protein. [PSI(+)] is an aggregated conformation of the translational release factor eRF3. [PSI(+)] modifies cellular fitness, inducing various phenotypes depending on genetic background. However, the genes displaying [PSI(+)]-controlled expression remain unknown. We used ribosome profiling in isogenic [PSI(+)] and [psi(-)] strains to identify the changes induced by [PSI(+)]. We found 100 genes with stop codon readthrough events and showed that many stress-response genes were repressed in the presence of [PSI(+)]. Surprisingly, [PSI(+)] was also found to affect reading frame selection independently of its effect on translation termination efficiency. These results indicate that [PSI(+)] has a broader impact than initially anticipated, providing explanations for the phenotypic differences between [psi(-)] and [PSI(+)] strains.

Maintaining two mating types: structure of the mating type locus and its role in heterokaryosis in Podospora anserina.

Pseudo-homothallism is a reproductive strategy elected by some fungi producing heterokaryotic sexual spores containing genetically different but sexually compatible nuclei. This lifestyle appears as a compromise between true homothallism (self-fertility with predominant inbreeding) and complete heterothallism (with exclusive outcrossing). However, pseudohomothallic species face the problem of maintaining heterokaryotic mycelia to fully benefit from this lifestyle, as homokaryons are self-sterile. Here, we report on the structure of chromosome 1 in mat+ and mat- isolates of strain S of the pseudohomothallic fungus Podospora anserina. Chromosome 1 contains either one of the mat+ and mat- mating types of P. anserina, which is mostly found in nature as a mat+/mat- heterokaryotic mycelium harboring sexually compatible nuclei. We identified a "mat" region ∼0.8 Mb long, devoid of meiotic recombination and containing the mating-type idiomorphs, which is a candidate to be involved in the maintenance of the heterokaryotic state, since the S mat+ and S mat- strains have different physiology that may enable hybrid-vigor-like phenomena in the heterokaryons. The mat region contains 229 coding sequences. A total of 687 polymorphisms were detected between the S mat+ and S mat- chromosomes. Importantly, the mat region is colinear between both chromosomes, which calls for an original mechanism of recombination inhibition. Microarray analyses revealed that 10% of the P. anserina genes have different transcriptional profiles in S mat+ and S mat-, in line with their different phenotypes. Finally, we show that the heterokaryotic state is faithfully maintained during mycelium growth of P. anserina, yet mat+/mat+ and mat-/mat- heterokaryons are as stable as mat+/mat- ones, evidencing a maintenance of heterokaryosis that does not rely on fitness-enhancing complementation between the S mat+ and S mat- strains.

Detection of non-coding RNA in bacteria and archaea using the DETR'PROK Galaxy pipeline.

RNA-seq experiments are now routinely used for the large scale sequencing of transcripts. In bacteria or archaea, such deep sequencing experiments typically produce 10-50 million fragments that cover most of the genome, including intergenic regions. In this context, the precise delineation of the non-coding elements is challenging. Non-coding elements include untranslated regions (UTRs) of mRNAs, independent small RNA genes (sRNAs) and transcripts produced from the antisense strand of genes (asRNA). Here we present a computational pipeline (DETR'PROK: detection of ncRNAs in prokaryotes) based on the Galaxy framework that takes as input a mapping of deep sequencing reads and performs successive steps of clustering, comparison with existing annotation and identification of transcribed non-coding fragments classified into putative 5' UTRs, sRNAs and asRNAs. We provide a step-by-step description of the protocol using real-life example data sets from Vibrio splendidus and Escherichia coli.

Transcriptomic profiling of the oyster pathogen Vibrio splendidus opens a window on the evolutionary dynamics of the small RNA repertoire in the Vibrio genus.

Work in recent years has led to the recognition of the importance of small regulatory RNAs (sRNAs) in bacterial regulation networks. New high-throughput sequencing technologies are paving the way to the exploration of an expanding sRNA world in nonmodel bacteria. In the Vibrio genus, compared to the enterobacteriaceae, still a limited number of sRNAs have been characterized, mostly in Vibrio cholerae, where they have been shown to be important for virulence, as well as in Vibrio harveyi. In addition, genome-wide approaches in V. cholerae have led to the discovery of hundreds of potential new sRNAs. Vibrio splendidus is an oyster pathogen that has been recently associated with massive mortality episodes in the French oyster growing industry. Here, we report the first RNA-seq study in a Vibrio outside of the V. cholerae species. We have uncovered hundreds of candidate regulatory RNAs, be it cis-regulatory elements, antisense RNAs, and trans-encoded sRNAs. Conservation studies showed the majority of them to be specific to V. splendidus. However, several novel sRNAs, previously unidentified, are also present in V. cholerae. Finally, we identified 28 trans sRNAs that are conserved in all the Vibrio genus species for which a complete genome sequence is available, possibly forming a Vibrio "sRNA core."