Adaptive gene loss in the common bean pan-genome during range expansion and domestication.

The common bean (Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

Soluble CD14 produced by bovine mammary epithelial cells modulates their response to full length LPS.

Bovine mastitis remains a major disease in cattle world-wide. In the mammary gland, mammary epithelial cells (MEC) are sentinels equipped with receptors allowing them to detect and respond to the invasion by bacterial pathogens, in particular Escherichia coli. Lipopolysaccharide (LPS) is the major E. coli motif recognized by MEC through its interaction with the TLR4 receptor and the CD14 co-receptor. Previous studies have highlighted the role of soluble CD14 (sCD14) in the efficient recognition of LPS molecules possessing a full-length O-antigen (LPSS). We demonstrate here that MEC are able to secrete CD14 and are likely to contribute to the presence of sCD14 in milk. We then investigated how sCD14 modulates and is required for the response of MEC to LPSS. This study highlights the key role of sCD14 for the full activation of the Myd88-independent pathway by LPSS. We also identified several lncRNA that are activated in MEC in response to LPS, including one lncRNA showing homologies with the mir-99a-let-7c gene (MIR99AHG). Altogether, our results show that a full response to LPS by mammary epithelial cells requires sCD14 and provide detailed information on how milk sCD14 can contribute to an efficient recognition of LPS from coliform pathogens.

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.

A candidate sex determination locus in amphibians which evolved by structural variation between X- and Y-chromosomes.

Most vertebrates develop distinct females and males, where sex is determined by repeatedly evolved environmental or genetic triggers. Undifferentiated sex chromosomes and large genomes have caused major knowledge gaps in amphibians. Only a single master sex-determining gene, the dmrt1-paralogue (dm-w) of female-heterogametic clawed frogs (Xenopus; ZW♀/ZZ♂), is known across >8740 species of amphibians. In this study, by combining chromosome-scale female and male genomes of a non-model amphibian, the European green toad, Bufo(tes) viridis, with ddRAD- and whole genome pool-sequencing, we reveal a candidate master locus, governing a male-heterogametic system (XX♀/XY♂). Targeted sequencing across multiple taxa uncovered structural X/Y-variation in the 5'-regulatory region of the gene bod1l, where a Y-specific non-coding RNA (ncRNA-Y), only expressed in males, suggests that this locus initiates sex-specific differentiation. Developmental transcriptomes and RNA in-situ hybridization show timely and spatially relevant sex-specific ncRNA-Y and bod1l-gene expression in primordial gonads. This coincided with differential H3K4me-methylation in pre-granulosa/pre-Sertoli cells, pointing to a specific mechanism of amphibian sex determination.

Hyperglycemia and circadian disruption lead to retinal dysfunction in a stabilized colony of the fat sand rat Psammomys obesus.

PURPOSE: The Fat Sand Rat (Psammomys obesus) recapitulates several features of human pre-proliferative diabetic retinopathy, but data are restricted to wild animals, incompatible with stringent biomedical research criteria. To overcome this barrier, we characterized retinal changes in a colony of P. obsesus maintained under strictly controlled housing conditions. METHODS: Animals were maintained on low or high caloric energy diets, and raised under either standard (12 h light/12 h dark) or shortened (5 h light/5 h dark) photoperiods. Visual responses were tested by electroretinography, while structural/molecular changes were assayed by immunochemistry and molecular biology (RNAseq and qPCR). RESULTS: Whereas high calorie diet alone did not induce hyperglycemia, coupled with short photoperiod >80 % animals developed severe hyper-insulinemia by 15 weeks, and 16 % animals further developed hyperglycemia. In these groups, electroretinography showed significant declines in visual responses in both hyper-insulinemic and hyperglycemic animals, especially in photopic (cone) responses. Transcriptomics analysis of hyperglycemic compared to low caloric controls revealed major upregulation in pathways involved in glial activation, extracellular matrix remodeling, inflammation, cytokine production, partial ischemic responses and angiogenesis. Western blotting against rhodopsin and cone opsin also showed decreased levels in both groups, overall decreases being greater for cones than rods in hyperglycemic animals. CONCLUSIONS: P. obesus maintained in rigorously monitored captive conditions, albeit showing attenuated responses to dietary overload compared to wild counterparts, nevertheless do develop some retinal features of diabetic retinopathy-like degeneration. Such a colony with known sanitary status opens their broader use for biomedical research.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues.

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

The genome of the ant Tetramorium bicarinatum reveals a tandem organization of venom peptides genes allowing the prediction of their regulatory and evolutionary profiles.

BACKGROUND: Venoms have evolved independently over a hundred times in the animal kingdom to deter predators and/or subdue prey. Venoms are cocktails of various secreted toxins, whose origin and diversification provide an appealing system for evolutionary researchers. Previous studies of the ant venom of Tetramorium bicarinatum revealed several Myrmicitoxin (MYRTX) peptides that gathered into seven precursor families suggesting different evolutionary origins. Analysis of the T. bicarinatum genome enabling further genomic approaches was necessary to understand the processes underlying the evolution of these myrmicitoxins. RESULTS: Here, we sequenced the genome of Tetramorium bicarinatum and reported the organisation of 44 venom peptide genes (vpg). Of the eleven chromosomes that make up the genome of T. bicarinatum, four carry the vpg which are organized in tandem repeats. This organisation together with the ML evolutionary analysis of vpg sequences, is consistent with evolution by local duplication of ancestral genes for each precursor family. The structure of the vpg into two or three exons is conserved after duplication events while the promoter regions are the least conserved parts of the vpg even for genes with highly identical sequences. This suggests that enhancer sequences were not involved in duplication events, but were recruited from surrounding regions. Expression level analysis revealed that most vpg are highly expressed in venom glands, although one gene or group of genes is much more highly expressed in each family. Finally, the examination of the genomic data revealed that several genes encoding transcription factors (TFs) are highly expressed in the venom glands. The search for binding sites (BS) of these TFs in the vpg promoters revealed hot spots of GATA sites in several vpg families. CONCLUSION: In this pioneering investigation on ant venom genes, we provide a high-quality assembly genome and the annotation of venom peptide genes that we think can fosters further genomic research to understand the evolutionary history of ant venom biochemistry.

OROSOMUCOID PROTEIN 1 regulation of sphingolipid synthesis is required for nodulation in Aeschynomene evenia.

Legumes establish symbiotic interactions with nitrogen-fixing rhizobia that are accommodated in root-derived organs known as nodules. Rhizobial recognition triggers a plant symbiotic signaling pathway that activates 2 coordinated processes: infection and nodule organogenesis. How these processes are orchestrated in legume species utilizing intercellular infection and lateral root base nodulation remains elusive. Here, we show that Aeschynomene evenia OROSOMUCOID PROTEIN 1 (AeORM1), a key regulator of sphingolipid biosynthesis, is required for nodule formation. Using A. evenia orm1 mutants, we demonstrate that alterations in AeORM1 function trigger numerous early aborted nodules, defense-like reactions, and shorter lateral roots. Accordingly, AeORM1 is expressed during lateral root initiation and elongation, including at lateral root bases where nodule primordium form in the presence of symbiotic bradyrhizobia. Sphingolipidomics revealed that mutations in AeORM1 lead to sphingolipid overaccumulation in roots relative to the wild type, particularly for very long-chain fatty acid-containing ceramides. Taken together, our findings reveal that AeORM1-regulated sphingolipid homeostasis is essential for rhizobial infection and nodule organogenesis, as well as for lateral root development in A. evenia.

A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants.

BACKGROUND: The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome. METHODS: We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines. RESULTS: We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference. CONCLUSIONS: We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.

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.

Whole genome sequencing reveals signals of adaptive admixture in Creole cattle.

The Creole cattle from Guadeloupe (GUA) are well adapted to the tropical environment. Its admixed genome likely played an important role in such adaptation. Here, we sought to detect genomic signatures of selection in the GUA genome. For this purpose, we sequenced 23 GUA individuals and combined our data with sequenced genomes of 99 animals representative of European, African and indicine groups. We detect 17,228,983 single nucleotide polymorphisms (SNPs) in the GUA genome, providing the most detailed exploration, to date, of patterns of genetic variation in this breed. We confirm the higher level of African and indicine ancestries, compared to the European ancestry and we highlight the African origin of indicine ancestry in the GUA genome. We identify five strong candidate regions showing an excess of indicine ancestry and consistently supported across the different detection methods. These regions encompass genes with adaptive roles in relation to immunity, thermotolerance and physical activity. We confirmed a previously identified horn-related gene, RXFP2, as a gene under strong selective pressure in the GUA population likely owing to human-driven (socio-cultural) pressure. Findings from this study provide insight into the genetic mechanisms associated with resilience traits in livestock.

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.

Mediator Subunit Med15 Regulates Cell Morphology and Mating in Candida lusitaniae.

Candida lusitaniae is an emerging opportunistic pathogenic yeast capable of shifting from yeast to pseudohyphae form, and it is one of the few Candida species with the ability to reproduce sexually. In this study, we showed that a dpp3Δ mutant, inactivated for a putative pyrophosphatase, is impaired in cell separation, pseudohyphal growth and mating. The defective phenotypes were not restored after the reconstruction of a wild-type DPP3 locus, reinforcing the hypothesis of the presence of an additional mutation that we suspected in our previous study. Genetic crosses and genome sequencing identified an additional mutation in MED15, encoding a subunit of the mediator complex that functions as a general transcriptional co-activator in Eukaryotes. We confirmed that inactivation of MED15 was responsible for the defective phenotypes by rescuing the dpp3Δ mutant with a wild-type copy of MED15 and constructing a med15Δ knockout mutant that mimics the phenotypes of dpp3Δ in vitro. Proteomic analyses revealed the biological processes under the control of Med15 and involved in hyphal growth, cell separation and mating. This is the first description of the functions of MED15 in the regulation of hyphal growth, cell separation and mating, and the pathways involved in C. lusitaniae.

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.

Insights into the molecular mechanisms of pesticide tolerance in the Aporrectodea caliginosa earthworm.

Diffuse pollution of the environment by pesticides has become a major soil threat to non-target organisms, such as earthworms for which declines have been reported. However some endogeic species are still abundant and persist in intensively cultivated fields, suggesting they become tolerant to long-term anthropogenic pressure. We thus considered the working hypothesis that populations of Aporrectodea caliginosa earthworms from conventionally managed fields developed a tolerance to pesticides compared with those from organically managed fields. To investigate this hypothesis, we studied earthworm populations of the same genetic lineage from soils that were either lowly or highly contaminated by pesticides to detect any constitutive expression of differentially expressed molecular pathways between these populations. Earthworm populations were then experimentally exposed to a fungicide-epoxiconazole-in the laboratory to identify different molecular responses when newly exposed to a pesticide. State-of-the-art omics technology (RNA sequencing) and bioinformatics were used to characterize molecular mechanisms of tolerance in a non-targeted way. Additional physiological traits (respirometry, growth, bioaccumulation) were monitored to assess tolerance at higher levels of biological organization. In the present study, we generated the de novo assembly transcriptome of A. caliginosa consisting of 64,556 contigs with N50 = 2862 pb. In total, 43,569 Gene Ontology terms were identified for 21,593 annotated sequences under the three main ontologies (biological processes, cellular components and molecular functions). Overall, we revealed that two same lineage populations of A. caliginosa earthworms, inhabiting similar pedo-climatic environment, have distinct gene expression pathways after they long-lived in differently managed agricultural soils with a contrasted pesticide exposure history for more than 22 years. The main difference was observed regarding metabolism, with upregulated pathways linked to proteolytic activities and the mitochondrial respiratory chain in the highly exposed population. This study improves our understanding of the long-term impact of chronic exposure of soil engineers to pesticide residues.

Venomics survey of six myrmicine ants provides insights into the molecular and structural diversity of their peptide toxins.

Among ants, Myrmicinae represents the most speciose subfamily. The venom composition previously described for these social insects is extremely variable, with alkaloids predominant in some genera while, conversely, proteomics studies have revealed that some myrmicine ant venoms are peptide-rich. Using integrated transcriptomic and proteomic approaches, we characterized the venom peptidomes of six ants belonging to the different tribes of Myrmicinae. We identified a total of 79 myrmicitoxins precursors which can be classified into 38 peptide families according to their mature sequences. Myrmicine ant venom peptidomes showed heterogeneous compositions, with linear and disulfide-bonded monomers as well as dimeric toxins. Several peptide families were exclusive to a single venom whereas some were retrieved in multiple species. A hierarchical clustering analysis of precursor signal sequences led us to divide the myrmicitoxins precursors into eight families, including some that have already been described in other aculeate hymenoptera such as secapin-like peptides and voltage-gated sodium channel (NaV) toxins. Evolutionary and structural analyses of two representatives of these families highlighted variation and conserved patterns that might be crucial to explain myrmicine venom peptide functional adaptations to biological targets.

Mosaic Evolution of Molecular Pathways for Sex Pheromone Communication in a Butterfly.

Unraveling the origin of molecular pathways underlying the evolution of adaptive traits is essential for understanding how new lineages emerge, including the relative contribution of conserved ancestral traits and newly evolved derived traits. Here, we investigated the evolutionary divergence of sex pheromone communication from moths (mostly nocturnal) to butterflies (mostly diurnal) that occurred ~119 million years ago. In moths, it is the females that typically emit pheromones to attract male mates, but in butterflies males emit pheromones that are used by females for mate choice. The molecular bases of sex pheromone communication are well understood in moths, but they have remained relatively unexplored in butterflies. We used a combination of transcriptomics, real time qPCR, and phylogenetics to identify genes involved in the different steps (i.e., production, regulation, and reception) of sex pheromone communication of the butterfly Bicyclus anynana. Our results show that the biosynthesis and reception of sex pheromones relies both on moth-specific gene families (reductases) and on more ancestral insect gene families (desaturases, olfactory receptors, odorant binding proteins). Interestingly, B. anynana appears to use what was believed to be the moth-specific neuropeptide Pheromone Biosynthesis Activating Neuropeptide (PBAN) for regulating sex pheromone production. Altogether, our results suggest that a mosaic pattern best explains how sex pheromone communication evolved in butterflies, with some molecular components derived from moths, and others conserved from more ancient insect ancestors. This is the first large-scale investigation of the genetic pathways underlying sex pheromone communication in a butterfly.

A mutant-based analysis of the establishment of Nod-independent symbiosis in the legume Aeschynomene evenia.

Intensive research on nitrogen-fixing symbiosis in two model legumes has uncovered the molecular mechanisms, whereby rhizobial Nod factors activate a plant symbiotic signaling pathway that controls infection and nodule organogenesis. In contrast, the so-called Nod-independent symbiosis found between Aeschynomene evenia and photosynthetic bradyrhizobia, which does not involve Nod factor recognition nor infection thread formation, is less well known. To gain knowledge on how Nod-independent symbiosis is established, we conducted a phenotypic and molecular characterization of A. evenia lines carrying mutations in different nodulation genes. Besides investigating the effect of the mutations on rhizobial symbiosis, we examined their consequences on mycorrhizal symbiosis and in nonsymbiotic conditions. Analyzing allelic mutant series for AePOLLUX, Ca2+/calmodulin dependent kinase, AeCYCLOPS, nodulation signaling pathway 2 (AeNSP2), and nodule inception demonstrated that these genes intervene at several stages of intercellular infection and during bacterial accommodation. We provide evidence that AeNSP2 has an additional nitrogen-dependent regulatory function in the formation of axillary root hairs at lateral root bases, which are rhizobia-colonized infection sites. Our investigation of the recently discovered symbiotic actor cysteine-rich receptor-like kinase specified that it is not involved in mycorrhization; however, it is essential for both symbiotic signaling and early infection during nodulation. These findings provide important insights on the modus operandi of Nod-independent symbiosis and contribute to the general understanding of how rhizobial-legume symbioses are established by complementing the information acquired in model legumes.

The MtZ Strain: Molecular Characteristics and Outbreak Investigation of the Most Successful Mycobacterium tuberculosis Strain in Aragon Using Whole-Genome Sequencing.

Since 2004, a tuberculosis surveillance protocol has been carried out in Aragon, thereby managing to detect all tuberculosis outbreaks that take place in the community. The largest outbreak was caused by a strain named Mycobacterium tuberculosis Zaragoza (MtZ), causing 242 cases as of 2020. The main objective of this work was to analyze this outbreak and the molecular characteristics of this successful strain that could be related to its greater transmission. To do this, we first applied whole-genome sequencing to 57 of the isolates. This revealed two principal transmission clusters and six subclusters arising from them. The MtZ strain belongs to L4.8 and had eight specific single nucleotide polymorphisms (SNPs) in genes considered to be virulence factors [ptpA, mc3D, mc3F, VapB41, pks15 (two SNPs), virS, and VapC50]. Second, a transcriptomic study was carried out to better understand the multiple IS6110 copies present in its genome. This allowed us to observe three effects of IS6110: the disruption of the gene in which the IS6110 is inserted (desA3), the overexpression of a gene (ppe38), and the absence of transcription of genes (cut1:Rv1765c) due to the recombination of two IS6110 copies. Finally, because of the disruption of ppe38 and ppe71 genes by an IS6110, a study of PE_PGRS secretion was carried out, showing that MtZ secretes these factors in higher amounts than the reference strain, thereby differing from the hypervirulent phenotype described for the Beijing strains. In conclusion, MtZ consists of several SNPs in genes related to virulence, pathogenesis, and survival, as well as other genomic polymorphisms, which may be implicated in its success among our population.