Genomic reconstruction of the successful establishment of a feralized bovine population on the subantarctic island of Amsterdam.

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean zebu in the founder ancestry. Comparative analysis of whole genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population's success in the harsh environment can be attributed to pre-adaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates.

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

Chromosome reciprocal translocations have accompanied subspecies evolution in bananas.

Chromosome rearrangements and the way that they impact genetic differentiation and speciation have long raised questions from evolutionary biologists. They are also a major concern for breeders because of their bearing on chromosome recombination. Banana is a major crop that derives from inter(sub)specific hybridizations between various once geographically isolated Musa species and subspecies. We sequenced 155 accessions, including banana cultivars and representatives of Musa diversity, and genotyped-by-sequencing 1059 individuals from 11 progenies. We precisely characterized six large reciprocal translocations and showed that they emerged in different (sub)species of Musa acuminata, the main contributor to currently cultivated bananas. Most diploid and triploid cultivars analyzed were structurally heterozygous for 1 to 4 M. acuminata translocations, highlighting their complex origin. We showed that all translocations induced a recombination reduction of variable intensity and extent depending on the translocations, involving only the breakpoint regions, a chromosome arm, or an entire chromosome. The translocated chromosomes were found preferentially transmitted in many cases. We explore and discuss the possible mechanisms involved in this preferential transmission and its impact on translocation colonization.

Recombination and Large Structural Variations Shape Interspecific Edible Bananas Genomes.

Admixture and polyploidization are major recognized eukaryotic genome evolutionary processes. Their impacts on genome dynamics vary among systems and are still partially deciphered. Many banana cultivars are triploid (sometimes diploid) interspecific hybrids between Musa acuminata (A genome) and M. balbisiana (B genome). They have no or very low fertility, are vegetatively propagated and have been classified as "AB," "AAB," or "ABB" based on morphological characters. We used NGS sequence data to characterize the A versus B chromosome composition of nine diploid and triploid interspecific cultivars, to compare the chromosome structures of A and B genomes and analyze A/B chromosome segregations in a polyploid context. We showed that interspecific recombination occurred frequently between A and B chromosomes. We identified two large structural variations between A and B genomes, a reciprocal translocation and an inversion that locally affected recombination and led to segregation distortion and aneuploidy in a triploid progeny. Interspecific recombination and large structural variations explained the mosaic genomes observed in edible bananas. The unprecedented resolution in deciphering their genome structure allowed us to start revisiting the origins of banana cultivars and provided new information to gain insight into the impact of interspecificity on genome evolution. It will also facilitate much more effective assessment of breeding strategies.

Two large reciprocal translocations characterized in the disease resistance-rich burmannica genetic group of Musa acuminata.

BACKGROUND AND AIMS: Banana cultivars are derived from hybridizations involving Musa acuminata subspecies. The latter diverged following geographical isolation in distinct South-east Asian continental regions and islands. Observation of chromosome pairing irregularities in meiosis of hybrids between these subspecies suggested the presence of large chromosomal structural variations. The aim of this study was to characterize such rearrangements. METHODS: Marker (single nucleotide polymorphism) segregation in a self-progeny of the 'Calcutta 4' accession and mate-pair sequencing were used to search for chromosomal rearrangements in comparison with the M. acuminata ssp. malaccensis genome reference sequence. Signature segment junctions of the revealed chromosome structures were identified and searched in whole-genome sequencing data from 123 wild and cultivated Musa accessions. KEY RESULTS: Two large reciprocal translocations were characterized in the seedy banana M. acuminata ssp. burmannicoides 'Calcutta 4' accession. One consisted of an exchange of a 240 kb distal region of chromosome 2 with a 7.2 Mb distal region of chromosome 8. The other involved an exchange of a 20.8 Mb distal region of chromosome 1 with a 11.6 Mb distal region of chromosome 9. Both translocations were found only in wild accessions belonging to the burmannicoides/burmannica/siamea subspecies. Only two of the 87 cultivars analysed displayed the 2/8 translocation, while none displayed the 1/9 translocation. CONCLUSION: Two large reciprocal translocations were identified that probably originated in the burmannica genetic group. Accurate characterization of these translocations should enhance the use of this disease resistance-rich burmannica group in breeding programmes.

Genome assembly using Nanopore-guided long and error-free DNA reads.

BACKGROUND: Long-read sequencing technologies were launched a few years ago, and in contrast with short-read sequencing technologies, they offered a promise of solving assembly problems for large and complex genomes. Moreover by providing long-range information, it could also solve haplotype phasing. However, existing long-read technologies still have several limitations that complicate their use for most research laboratories, as well as in large and/or complex genome projects. In 2014, Oxford Nanopore released the MinION® device, a small and low-cost single-molecule nanopore sequencer, which offers the possibility of sequencing long DNA fragments. RESULTS: The assembly of long reads generated using the Oxford Nanopore MinION® instrument is challenging as existing assemblers were not implemented to deal with long reads exhibiting close to 30% of errors. Here, we presented a hybrid approach developed to take advantage of data generated using MinION® device. We sequenced a well-known bacterium, Acinetobacter baylyi ADP1 and applied our method to obtain a highly contiguous (one single contig) and accurate genome assembly even in repetitive regions, in contrast to an Illumina-only assembly. Our hybrid strategy was able to generate NaS (Nanopore Synthetic-long) reads up to 60 kb that aligned entirely and with no error to the reference genome and that spanned highly conserved repetitive regions. The average accuracy of NaS reads reached 99.99% without losing the initial size of the input MinION® reads. CONCLUSIONS: We described NaS tool, a hybrid approach allowing the sequencing of microbial genomes using the MinION® device. Our method, based ideally on 20x and 50x of NaS and Illumina reads respectively, provides an efficient and cost-effective way of sequencing microbial or small eukaryotic genomes in a very short time even in small facilities. Moreover, we demonstrated that although the Oxford Nanopore technology is a relatively new sequencing technology, currently with a high error rate, it is already useful in the generation of high-quality genome assemblies.

Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks.

Sexual breeding at the tetraploid level is a promising strategy for rootstock breeding in citrus. Due to the interspecific origin of most of the conventional diploid citrus rootstocks that produced the tetraploid germplasm, the optimization of this strategy requires better knowledge of the meiotic behavior of the tetraploid parents. This work used Genotyping By Sequencing (GBS) data from 103 tetraploid hybrids to study the meiotic behavior and generate a high-density recombination landscape for their tetraploid intergenic Swingle citrumelo and interspecific Volkamer lemon progenitors. A genetic association study was performed with root architecture traits. For citrumelo, high preferential chromosome pairing was revealed and led to an intermediate inheritance with a disomic tendency. Meiosis in Volkamer lemon was more complex than that of citrumelo, with mixed segregation patterns from disomy to tetrasomy. The preferential pairing resulted in low interspecific recombination levels and high interspecific heterozygosity transmission by the diploid gametes. This meiotic behavior affected the efficiency of Quantitative Trait Loci (QTL) detection. Nevertheless, it enabled a high transmission of disease and pest resistance candidate genes from P. trifoliata that are heterozygous in the citrumelo progenitor. The tetrazyg strategy, using doubled diploids of interspecific origin as parents, appears to be efficient in transferring the dominant traits selected at the parental level to the tetraploid progenies.

Genome-wide mapping of individual replication fork velocities using nanopore sequencing.

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

Lifespan prolonging mechanisms and insulin upregulation without fat accumulation in long-lived reproductives of a higher termite.

Kings and queens of eusocial termites can live for decades, while queens sustain a nearly maximal fertility. To investigate the molecular mechanisms underlying their long lifespan, we carried out transcriptomics, lipidomics and metabolomics in Macrotermes natalensis on sterile short-lived workers, long-lived kings and five stages spanning twenty years of adult queen maturation. Reproductives share gene expression differences from workers in agreement with a reduction of several aging-related processes, involving upregulation of DNA damage repair and mitochondrial functions. Anti-oxidant gene expression is downregulated, while peroxidability of membranes in queens decreases. Against expectations, we observed an upregulated gene expression in fat bodies of reproductives of several components of the IIS pathway, including an insulin-like peptide, Ilp9. This pattern does not lead to deleterious fat storage in physogastric queens, while simple sugars dominate in their hemolymph and large amounts of resources are allocated towards oogenesis. Our findings support the notion that all processes causing aging need to be addressed simultaneously in order to prevent it.

Telomere-to-telomere gapless chromosomes of banana using nanopore sequencing.

Long-read technologies hold the promise to obtain more complete genome assemblies and to make them easier. Coupled with long-range technologies, they can reveal the architecture of complex regions, like centromeres or rDNA clusters. These technologies also make it possible to know the complete organization of chromosomes, which remained complicated before even when using genetic maps. However, generating a gapless and telomere-to-telomere assembly is still not trivial, and requires a combination of several technologies and the choice of suitable software. Here, we report a chromosome-scale assembly of a banana genome (Musa acuminata) generated using Oxford Nanopore long-reads. We generated a genome coverage of 177X from a single PromethION flowcell with near 17X with reads longer than 75 kbp. From the 11 chromosomes, 5 were entirely reconstructed in a single contig from telomere to telomere, revealing for the first time the content of complex regions like centromeres or clusters of paralogous genes.

WFS1 mutations are frequent monogenic causes of juvenile-onset diabetes mellitus in Lebanon.

Most cases of juvenile-onset diabetes (JOD) are diagnosed as type 1 diabetes (T1D), for which genetic studies conducted in outbred Caucasian populations support the concept of multifactorial inheritance. However, this view may be partly challenged in particular population settings. In view of the suggestive evidence for a high prevalence of Wolfram syndrome (WFS) in Lebanon, the phenotypic variability associated with WFS1 mutations, and the high consanguinity rate in Lebanon, we aimed to evaluate the contribution of WFS1 mutations as monogenic determinants to JOD in Lebanon. We performed a family-based genetic study, with linkage analysis followed by systematic mutation screening of WFS1 exons in all JOD probands. The study population consisted of an unbiased recruitment of all juvenile-onset insulin-dependent diabetic patients from a specialized diabetes pediatric clinic in Beirut, Lebanon. Homozygous or compound heterozygous WFS1 mutations were found in 22 of the 399 JOD probands (5.5%), resulting in WFS (17 probands) or in non-syndromic non-autoimmune diabetes mellitus (DM, five probands). These accounted for 12.1% (21/174) of probands in consanguineous families, compared with 0.4% (1/225) in non-consanguineous families. Of the 38 patients identified with homozygous or compound heterozygous WFS1 mutations, 11 (29%) had non-syndromic DM, all of whom carried a particular WFS1 mutation, WFS1(LIB), encoding a protein with an extended C-terminal domain. This mutation resulted in a delayed onset or absence of extrapancreatic features. These results underscore the major impact of population-specific factors, such as population-specific mutations and founder effects, and family structure in the genetic determinism of JOD.

Colonization kinetics and implantation follow-up of the sewage microbiome in an urban wastewater treatment plant.

The Seine-Morée wastewater treatment plant (SM_WWTP), with a capacity of 100,000 population-equivalents, was fed with raw domestic wastewater during all of its start-up phase. Its microbiome resulted from the spontaneous evolution of wastewater-borne microorganisms. This rare opportunity allowed us to analyze the sequential microbiota colonization and implantation follow up during the start-up phase of this WWTP by means of regular sampling carried out over 8 months until the establishment of a stable and functional ecosystem. During the study, biological nitrification-denitrification and dephosphatation occurred 68 days after the start-up of the WWTP, followed by flocs decantation 91 days later. High throughput sequencing of 18S and 16S rRNA genes was performed using Illumina's MiSeq and PGM Ion Torrent platforms respectively, generating 584,647 16S and 521,031 18S high-quality sequence rDNA reads. Analyses of 16S and 18S rDNA datasets show three colonization phases occurring concomitantly with nitrification, dephosphatation and floc development processes. Thus, we could define three microbiota profiles that sequentially colonized the SM_WWTP: the early colonizers, the late colonizers and the continuous spectrum population. Shannon and inverse Simpson diversity indices indicate that the highest microbiota diversity was reached at days 133 and 82 for prokaryotes and eukaryotes respectively; after that, the structure and complexity of the wastewater microbiome reached its functional stability. This study demonstrates that physicochemical parameters and microbial metabolic interactions are the main forces shaping microbial community structure, gradually building up and maintaining a functionally stable microbial ecosystem.

FORK-seq: replication landscape of the Saccharomyces cerevisiae genome by nanopore sequencing.

Genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200-nucleotide resolution. By quantifying BrdU incorporation along pulse-chased replication intermediates from Saccharomyces cerevisiae, we orient 58,651 replication tracks reproducing population-based replication directionality profiles and map 4964 and 4485 individual initiation and termination events, respectively. Although most events cluster at known origins and fork merging zones, 9% and 18% of initiation and termination events, respectively, occur at many locations previously missed. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

Long-read assembly of the Brassica napus reference genome Darmor-bzh.

BACKGROUND: The combination of long reads and long-range information to produce genome assemblies is now accepted as a common standard. This strategy not only allows access to the gene catalogue of a given species but also reveals the architecture and organization of chromosomes, including complex regions such as telomeres and centromeres. The Brassica genus is not exempt, and many assemblies based on long reads are now available. The reference genome for Brassica napus, Darmor-bzh, which was published in 2014, was produced using short reads and its contiguity was extremely low compared with current assemblies of the Brassica genus. FINDINGS: Herein, we report the new long-read assembly of Darmor-bzh genome (Brassica napus) generated by combining long-read sequencing data and optical and genetic maps. Using the PromethION device and 6 flowcells, we generated ∼16 million long reads representing 93× coverage and, more importantly, 6× with reads longer than 100 kb. This ultralong-read dataset allows us to generate one of the most contiguous and complete assemblies of a Brassica genome to date (contig N50 > 10 Mb). In addition, we exploited all the advantages of the nanopore technology to detect modified bases and sequence transcriptomic data using direct RNA to annotate the genome and focus on resistance genes. CONCLUSION: Using these cutting-edge technologies, and in particular by relying on all the advantages of the nanopore technology, we provide the most contiguous Brassica napus assembly, a resource that will be valuable to the Brassica community for crop improvement and will facilitate the rapid selection of agronomically important traits.

A novel immunodeficiency associated with hypomorphic RAG1 mutations and CMV infection.

Amorphic mutations in the recombination activating genes RAG1 and RAG2 have been reported to cause T- B- SCID, whereas hypomorphic mutations led to the expansion of a few autoimmune T cell clones responsible for the Omenn syndrome phenotype. We report here a novel clinical and immunological phenotype associated with recessive RAG1 hypomorphic mutations in 4 patients from 4 different families. The immunological phenotype consists of the oligoclonal expansion of TCR gammadelta T cells combined with TCR alphabeta T cell lymphopenia. The clinical phenotype consists of severe, disseminated CMV infection and autoimmune blood cell manifestations. Repertoire studies suggest that CMV infection, in the setting of this particular T cell immunodeficiency, may have driven the TCR gammadelta T cell clonal expansion. This observation extends the range of clinical and immunological phenotypes associated with RAG mutations, emphasizing the role of the genetic background and microbial environment in determining disease phenotype.

Evidence for linkage of a new region (11p14) to eczema and allergic diseases.

Asthma, allergic rhinitis (AR) and atopic dermatitis also called eczema are allergic co-morbidites, which are likely to depend on pleiotropic genetic effects as well as on specific genetic factors. After a previous genome-wide linkage screen conducted for asthma and AR in a sample of 295 French EGEA families ascertained through asthmatic subjects, the aim here was to search for genetic factors involved in eczema and more particularly the ones shared by the three allergic diseases using the same EGEA data. In this sake, eczema and phenotypes of "allergic disease" accounting for the joint information on the presence/absence of the three diseases were examined by linkage analyses using the maximum likelihood binomial method. A fine mapping was carried out in regions detected for potential linkage, followed by association studies using the family-based association test (FBAT). Evidence for linkage to 11p14 region was shown for "allergic disease" and eczema. Linkage was also indicated between eczema and 5q13 and between "allergic disease" and both 5p15 and 17q21 regions. Fine mapping supported the evidence of linkage to 11p14 and FBAT analyses showed the association between "allergic disease" and a marker located at the linkage peak on 11p14. Further investigations in this region will allow identifying genetic factor(s) which could have pleiotropic effect in the three allergic diseases.

A novel X-linked recessive form of Mendelian susceptibility to mycobaterial disease.

BACKGROUND: Mendelian susceptibility to mycobacterial disease (MSMD) is associated with infection caused by weakly virulent mycobacteria in otherwise healthy people. Causal germline mutations in five autosomal genes (IFNGR1, IFNGR2, STAT1, IL12RB1, IL12B) and one X-linked (NEMO) gene have been described. The gene products are physiologically related, as they are involved in interleukin 12/23-dependent, interferon gamma-mediated immunity. However, no genetic aetiology has yet been identified for about half the patients with MSMD. METHODS: A large kindred was studied, including four male maternal relatives with recurrent mycobacterial disease, suggesting X-linked recessive inheritance. Three patients had recurrent disease caused by the bacille Calmette-Guérin vaccine, and the fourth had recurrent tuberculosis. The infections showed tropism for the peripheral lymph nodes. RESULTS: Known autosomal and X-linked genetic aetiologies of MSMD were excluded through genetic and immunological investigations. Genetic linkage analysis of the X-chromosome identified two candidate regions, on Xp11.4-Xp21.2 and Xq25-Xq26.3, with a maximum LOD score of 2. CONCLUSION: A new X-linked recessive form of MSMD is reported, paving the way for the identification of a new MSMD-causing gene.

De novo assembly and annotation of three Leptosphaeria genomes using Oxford Nanopore MinION sequencing.

Leptosphaeria maculans and Leptosphaeria biglobosa are ascomycete phytopathogens of Brassica napus (oilseed rape, canola). Here we report the complete sequence of three Leptosphaeria genomes (L. maculans JN3, L. maculans Nz-T4 and L. biglobosa G12-14). Nz-T4 and G12-14 genome assemblies were generated de novo and the reference JN3 genome assembly was improved using Oxford Nanopore MinION reads. The new assembly of L. biglobosa showed the existence of AT rich regions and pointed to a genome compartmentalization previously unsuspected following Illumina sequencing. Moreover nanopore sequencing allowed us to generate a chromosome-level assembly for the L. maculans reference isolate, JN3. The genome annotation was supported by integrating conserved proteins and RNA sequencing from Leptosphaeria-infected samples. The newly produced high-quality assemblies and annotations of those three Leptosphaeria genomes will allow further studies, notably focused on the tripartite interaction between L. maculans, L. biglobosa and oilseed rape. The discovery of as yet unknown effectors will notably allow progress in B. napus breeding towards L. maculans resistance.

Chromosome-scale assemblies of plant genomes using nanopore long reads and optical maps.

Plant genomes are often characterized by a high level of repetitiveness and polyploid nature. Consequently, creating genome assemblies for plant genomes is challenging. The introduction of short-read technologies 10 years ago substantially increased the number of available plant genomes. Generally, these assemblies are incomplete and fragmented, and only a few are at the chromosome scale. Recently, Pacific Biosciences and Oxford Nanopore sequencing technologies were commercialized that can sequence long DNA fragments (kilobases to megabase) and, using efficient algorithms, provide high-quality assemblies in terms of contiguity and completeness of repetitive regions1-4. However, even though genome assemblies based on long reads exhibit high contig N50s (>1 Mb), these methods are still insufficient to decipher genome organization at the chromosome level. Here, we describe a strategy based on long reads (MinION or PromethION sequencers) and optical maps (Saphyr system) that can produce chromosome-level assemblies and demonstrate applicability by generating high-quality genome sequences for two new dicotyledon morphotypes, Brassica rapa Z1 (yellow sarson) and Brassica oleracea HDEM (broccoli), and one new monocotyledon, Musa schizocarpa (banana). All three assemblies show contig N50s of >5 Mb and contain scaffolds that represent entire chromosomes or chromosome arms.

Oak genome reveals facets of long lifespan.

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes1 but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times2. With 450 species spread throughout Asia, Europe and America3, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical4 and modelling5 approaches have shown that intra-organismal genetic heterogeneity can be selected for6 and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes7. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.