A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

Genes controlling hydrolysate toxin tolerance identified by QTL analysis of the natural Saccharomyces cerevisiae BCC39850

Lignocellulosic material can be converted to valorized products such as fuels. Pretreatment is an essential step in conversion, which is needed to increase the digestibility of the raw material for microbial fermentation. However, pretreatment generates by-products (hydrolysate toxins) that are detrimental to microbial growth. In this study, natural Saccharomyces strains isolated from habitats in Thailand were screened for their tolerance to synthetic hydrolysate toxins (synHTs). The Saccharomyces cerevisiae natural strain BCC39850 (toxin-tolerant) was crossed with the laboratory strain CEN.PK2-1C (toxin-sensitive), and quantitative trait locus (QTL) analysis was performed on the segregants using phenotypic scores of growth (OD600) and glucose consumption. VMS1, DET1, KCS1, MRH1, YOS9, SYO1, and YDR042C were identified from QTLs as candidate genes associated with the tolerance trait. CEN.PK2-1C knockouts of the VMS1, YOS9, KCS1, and MRH1 genes exhibited significantly greater hydrolysate toxin sensitivity to growth, whereas CEN.PK2-1C knock-ins with replacement of VMS1 and MRH1 genes from the BCC39850 alleles showed significant increased ethanol production titers compared with the CEN.PK2-1C parental strain in the presence of synHTs. The discovery of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin tolerance in S. cerevisiae indicates the roles of the endoplasmic-reticulum-associated protein degradation pathway, plasma membrane protein association, and the phosphatidylinositol signaling system in this trait. KEY POINTS: • QTL analysis was conducted using a hydrolysate toxin-tolerant S. cerevisiae natural strain • Deletion of VMS1, YOS9, MRH1, and KCS1 genes associated with hydrolysate toxin-sensitivity • Replacement of VMS1 and MRH1 with natural strain alleles increased ethanol production titers in the presence of hydrolysate toxins.

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

OBJECTIVES

We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China.

METHODS

We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis.

RESULT

The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China.

CONCLUSION

MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

Resistance mechanism of Escherichia coli strains with different ampicillin resistance levels

Antibiotic resistance is an important problem that threatens medical treatment. Differences in the resistance levels of microorganisms cause great difficulties in understanding the mechanisms of antibiotic resistance. Therefore, the molecular reasons underlying the differences in the level of antibiotic resistance need to be clarified. For this purpose, genomic and transcriptomic analyses were performed on three Escherichia coli strains with varying degrees of adaptive resistance to ampicillin. Whole-genome sequencing of strains with different levels of resistance detected five mutations in strains with 10-fold resistance and two additional mutations in strains with 95-fold resistance. Overall, three of the seven mutations occurred as a single base change, while the other four occurred as insertions or deletions. While it was thought that 10-fold resistance was achieved by the effect of mutations in the ftsI, marAR, and rpoC genes, it was found that 95-fold resistance was achieved by the synergistic effect of five mutations and the ampC mutation. In addition, when the general transcriptomic profiles were examined, it was found that similar transcriptomic responses were elicited in strains with different levels of resistance. This study will improve our view of resistance mechanisms in bacteria with different levels of resistance and provide the basis for our understanding of the molecular mechanism of antibiotic resistance in ampicillin-resistant E. coli strains. KEY POINTS: •The mutation of the ampC promoter may act synergistically with other mutations and lead to higher resistance. •Similar transcriptomic responses to ampicillin are induced in strains with different levels of resistance. •Low antibiotic concentrations are the steps that allow rapid achievement of high antibiotic resistance.

Conserved A-to-I RNA editing with non-conserved recoding expands the candidates of functional editing sites

Adenosine-to-inosine (A-to-I) RNA editing recodes the genome and confers flexibility for the organisms to adapt to the environment. It is believed that RNA recoding sites are well suited for facilitating adaptive evolution by increasing the proteomic diversity in a temporal-spatial manner. The function and essentiality of a few conserved recoding sites are recognized. However, the experimentally discovered functional sites only make up a small corner of the total sites, and there is still the need to expand the repertoire of such functional sites with bioinformatic approaches. In this study, we define a new category of RNA editing sites termed 'conserved editing with non-conserved recoding' and systematically identify such sites in Drosophila editomes, figuring out their selection pressure and signals of adaptation at inter-species and intra-species levels. Surprisingly, conserved editing sites with non-conserved recoding are not suppressed and are even slightly overrepresented in Drosophila. DNA mutations leading to such cases are also favoured during evolution, suggesting that the function of those recoding events in different species might be diverged, specialized, and maintained. Finally, structural prediction suggests that such recoding in potassium channel Shab might increase ion permeability and compensate the effect of low temperature. In conclusion, conserved editing with non-conserved recoding might be functional as well. Our study provides novel aspects in considering the adaptive evolution of RNA editing sites and meanwhile expands the candidates of functional recoding sites for future validation.

Genome assembly and variant analysis of two Saccharomyces cerevisiae strains isolated from stingless bee pollen

Products from stingless bees are rich reservoirs of microbial diversity, including yeasts with fermentative potential. Previously, two Saccharomyces cerevisiae strains, JP14 and IP9, were isolated from Jataí (Tetragonisca angustula) and Iraí (Nannotrigona testaceicornis) bees, respectively, aiming at mead production. Both strains presented great osmotic and sulfite tolerance, and ethanol production, although they have a high free amino nitrogen demand. Herein, their genomes were sequenced, assembled, and annotated, and the variants were compared to the S. cerevisiae S288c reference strain. The final assembly of IP9 and JP14 presented high N50 and BUSCO scores, and more than 6430 protein-coding genes. Additionally, nQuire predicted the ploidy of IP9 as diploid, but the results were not enough to determine the ploidy of JP14. The mitochondrial genomes of IP9 and JP14 presented the same gene content as S288c but the genes were rearranged and fragmented in different patterns. Meanwhile, the genes with mutations of high impact (e.g., indels, gain of stop codon) for both yeasts were enriched for transmembrane transport, electron transfer, oxidoreductase, heme binding, fructose, mannose, and glucose transport, activities related to the respiratory chain and sugar metabolism. The IP9 strain presented copy number gains in genes related to sugar transport and cell morphogenesis; in JP14, genes were enriched for disaccharide metabolism and transport, response to reactive oxygen species, and polyamine transport. On the other hand, IP9 presented copy number losses related to disaccharide, thiamine, and aldehyde metabolism, while JP14 presented depletions related to disaccharide, oligosaccharide, asparagine, and aspartate metabolism. Notably, both strains presented a killer toxin gene, annotated from the assembling of unmapped reads, representing a potential mechanism for the control of other microorganisms population in the environment. Therefore, the annotated genomes of JP14 and IP9 presented a high selective pressure for sugar and nitrogen metabolism and stress response, consistent with their isolation source and fermentative properties.

Whole genome resequencing revealed genomic variants and functional pathways related to adaptation in Indian yak populations

The present study aims to identify genomic variants through a whole genome sequencing (WGS) approach and uncover biological pathways associated with adaptation and fitness in Indian yak populations. A total of 30 samples (10 from each population) were included from Arunachali, Himachali and Ladakhi yak populations. WGS analysis revealed a total of 32171644, 27260825, and 32632460 SNPs and 4865254, 4429941, and 4847513 Indels in the Arunachali, Himachali, and Ladakhi yaks, respectively. Genes such as RYR2, SYNE2, BOLA, HF1, and the novel transcript ENSBGRG00000011079 were found to have the maximum number of high impact variants in all three yak populations, and might play a major role in local adaptation. Functional enrichment analysis of genes harboring high impact SNPs revealed overrepresented pathways related to response to stress, immune system regulation, and high-altitude adaptation. This study provides comprehensive information about genomic variants and their annotation in Indian yak populations, thus would serve as a data resource for researchers working on the yaks. Furthermore, it could be well exploited for better yak conservation strategies by estimating population genetics parameters viz., effective population size, inbreeding, and observed and expected heterozygosity.

Evolution of the regulation of developmental gene expression in blind Mexican cavefish

Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.

Differential adaptation of the yeast Candida anglica to fermented food

A single strain of Candida anglica, isolated from cider, is available in international yeast collections. We present here seven new strains isolated from French PDO cheeses. For one of the cheese strains, we achieved a high-quality genome assembly of 13.7 Mb with eight near-complete telomere-to-telomere chromosomes. The genomes of two additional cheese strains and of the cider strain were also assembled and annotated, resulting in a core genome of 5966 coding sequences. Phylogenetic analysis showed that the seven cheese strains clustered together, away from the cider strain. Mating-type locus analysis revealed the presence of a MATa locus in the cider strain but a MATalpha locus in all cheese strains. The presence of LINE retrotransposons at identical genome position in the cheese strains, and two different karyotypic profiles resulting from chromosomal rearrangements were observed. Together, these findings are consistent with clonal propagation of the cheese strains. Phenotypic trait variations were observed within the cheese population under stress conditions whereas the cider strain was found to have a much greater capacity for growth in all conditions tested.

Convergent evolution of allele-specific gene expression that leads to non-small cell lung cancer in different human populations

Phenotypical innovations during evolution are caused by novel mutations, which are usually heterozygous at the beginning. The gene expressions on two alleles of these mutation sites are not necessarily identical, leading to flexible allele-specific regulation in cell systems. We retrieve the transcriptome data of normal and non-small cell lung cancer (NSCLC) tissues from 47 African Americans (AA) and 50 European Americans (EA). We analyze the differentially expressed genes (DEGs) in NSCLC as well as the tumor-specific mutations. Expression and mutation profiles show convergent evolution in AA and EA populations. The tumor-specific mutations are poorly overlapped, but many of them are located in the same genes, mainly oncogenes and tumor suppressor genes. The DEGs in tumors are majorly caused by the mutated alleles rather than normal alleles. The relative expressions of mutated alleles are highly correlated between AA and EA. The differential expression in NSCLC is predominantly mediated by the mutated alleles on heterozygous sites. This molecular mechanism underlying NSCLC oncogenesis is conserved across different human populations, exhibiting convergent evolution. We present this novel angle that differential expression analysis should be performed separately for different alleles. Our ideas should greatly benefit the cancer community.

C. elegans LIN-66 mediates EIF-3/eIF3-dependent protein translation via a cold-shock domain

Protein translation initiation is a conserved process involving many proteins acting in concert. The 13 subunit eukaryotic initiation factor 3 (eIF3) complex is essential for assembly of the pre-initiation complex that scans mRNA and positions ribosome at the initiation codon. We previously reported that a gain-of-function (gf) mutation affecting the G subunit of the Caenorhabditis elegans eIF3 complex, eif-3.g(gf), selectively modulates protein translation in the ventral cord cholinergic motor neurons. Here, through unbiased genetic suppressor screening, we identified that the gene lin-66 mediates eif-3.g(gf)-dependent protein translation in motor neurons. LIN-66 is composed largely of low-complexity amino acid sequences with unknown functional domains. We combined bioinformatics analysis with in vivo functional dissection and identified a cold-shock domain in LIN-66 critical for its function. In cholinergic motor neurons, LIN-66 shows a close association with EIF-3.G in the cytoplasm. The low-complexity amino acid sequences of LIN-66 modulate its subcellular pattern. As cold-shock domains function broadly in RNA regulation, we propose that LIN-66 mediates stimulus-dependent protein translation by facilitating the interaction of mRNAs with EIF-3.G.

Qualitative and quantitative analysis of MED12 c.887G>A causing both missense and splicing variants in X-linked Ohdo syndrome

The phenotypes associated with MED12 pathogenic variants are diverse. Male patients usually have missense variants, but the effects of base substitutions on mRNA splicing have not been investigated. Here, we report a Japanese brother with intellectual disability, characteristic facial appearance with blepharophimosis, cleft palate, Fallot tetralogy, vesicoureteral reflux, and deafness. A known missense pathogenic variant was detected in MED12, NM_005120.3:c.887G>A p.(Arg296Gln), and X-linked Ohdo syndrome was diagnosed in combination with their phenotype. mRNA splicing of MED12 was evaluated qualitatively and quantitatively using long-range PCR-based targeted RNA sequencing (reverse transcribed long amplicon sequencing), and it was shown that this missense variant simultaneously causes aberrant splicing of the 42-bp in-frame deletion in exon 7, r.847_888del, which accounts for approximately 30% of the mRNAs in both siblings. The X chromosome inactivation study showed that the X chromosome carrying the mutant allele was 100% inactivated in the carrier mothers. mRNA level analysis is essential for the accurate interpretation of the effects of variants. In this case, the MED12 protein function may be reduced by more than just an amino acid substitution, resulting in the patients with the most severe phenotype of MED12-related syndrome in males.

Century-long timelines of herbarium genomes predict plant stomatal response to climate change

Dissecting plant responses to the environment is key to understanding whether and how plants adapt to anthropogenic climate change. Stomata, plants' pores for gas exchange, are expected to decrease in density following increased CO2 concentrations, a trend already observed in multiple plant species. However, it is unclear whether such responses are based on genetic changes and evolutionary adaptation. Here we make use of extensive knowledge of 43 genes in the stomatal development pathway and newly generated genome information of 191 Arabidopsis thaliana historical herbarium specimens collected over 193 years to directly link genetic variation with climate change. While we find that the essential transcription factors SPCH, MUTE and FAMA, central to stomatal development, are under strong evolutionary constraints, several regulators of stomatal development show signs of local adaptation in contemporary samples from different geographic regions. We then develop a functional score based on known effects of gene knock-out on stomatal development that recovers a classic pattern of stomatal density decrease over the past centuries, suggesting a genetic component contributing to this change. This approach combining historical genomics with functional experimental knowledge could allow further investigations of how different, even in historical samples unmeasurable, cellular plant phenotypes may have already responded to climate change through adaptive evolution.

Human genetic structure in Northwest France provides new insights into West European historical demography

The demographical history of France remains largely understudied despite its central role toward understanding modern population structure across Western Europe. Here, by exploring publicly available Europe-wide genotype datasets together with the genomes of 3234 present-day and six newly sequenced medieval individuals from Northern France, we found extensive fine-scale population structure across Brittany and the downstream Loire basin and increased population differentiation between the northern and southern sides of the river Loire, associated with higher proportions of steppe vs. Neolithic-related ancestry. We also found increased allele sharing between individuals from Western Brittany and those associated with the Bell Beaker complex. Our results emphasise the need for investigating local populations to better understand the distribution of rare (putatively deleterious) variants across space and the importance of common genetic legacy in understanding the sharing of disease-related alleles between Brittany and people from western Britain and Ireland.

Whole-genome sequencing analysis of soybean diversity across different countries and selection signature of Korean soybean accession

Soybean is an important agricultural crop known for its high protein and oil content, contributing to essential nutritional and health benefits for humans. Domesticated in China over 5,000 years ago, soybean has since adapted to diverse environments and spread worldwide. This study aimed to investigate the genomic characteristics and population structures of 2,317 publicly available soybean whole-genome sequences from diverse geographical regions, including China, Korea, Japan, Europe, North America, and South America. We used large-scale whole-genome sequencing data to perform high-resolution analyses to reveal the genetic characteristics of soybean accessions. Soybean accessions from China and Korea exhibited landrace characteristics, indicating higher genetic diversity and adaptation to local environments. On the other hand, soybean accessions from Japan, the European Union, and South America were found to have low genetic diversity due to artificial selection and breeding for agronomic traits. We also identified key variants and genes associated with the ability to adapt to different environments. In Korean soybean accessions, we observed strong selection signals for isoflavone synthesis, an adaptive trait critical for improving soybean adaptability, survival, and reproductive success by mitigating environmental stress. Identifying specific genomic regions showing unique patterns of selective sweeps for genes such as HIDH, CYP73A11, IFS1, and CYP81E11 associated with isoflavone synthesis provided valuable insights into potential adaptation mechanisms. Our research has significantly improved our understanding of soybean diversity at the genetic level. We have identified key genetic variants and genes influencing adaptability, laying the foundation for future advances in genomics-based breeding programs and crop improvement efforts.

Genome assembly of Hibiscus sabdariffa L. provides insights into metabolisms of medicinal natural products

Hibiscus sabdariffa L. is a widely cultivated herbaceous plant with diverse applications in food, tea, fiber, and medicine. In this study, we present a high-quality genome assembly of H. sabdariffa using more than 33 Gb of high-fidelity (HiFi) long-read sequencing data, corresponding to ∼20× depth of the genome. We obtained 3 genome assemblies of H. sabdariffa: 1 primary and 2 partially haplotype-resolved genome assemblies. These genome assemblies exhibit N50 contig lengths of 26.25, 11.96, and 14.50 Mb, with genome coverage of 141.3, 86.0, and 88.6%, respectively. We also utilized 26 Gb of total RNA sequencing data to predict 154k, 79k, and 87k genes in the respective assemblies. The completeness of the primary genome assembly and its predicted genes was confirmed by the benchmarking universal single-copy ortholog analysis with a completeness rate of 99.3%. Based on our high-quality genomic resources, we constructed genetic networks for phenylpropanoid and flavonoid metabolism and identified candidate biosynthetic genes, which are responsible for producing key intermediates of roselle-specific medicinal natural products. Our comprehensive genomic and functional analysis opens avenues for further exploration and application of valuable natural products in H. sabdariffa.

Genome-wide association analyses using multilocus models on bananas (Musa spp.) reveal candidate genes related to morphology, fruit quality, and yield

Bananas (Musa spp.) are an essential fruit worldwide and rank as the fourth most significant food crop for addressing malnutrition due to their rich nutrients and starch content. The potential of their genetic diversity remains untapped due to limited molecular breeding tools. Our study examined a phenotypically diverse group of 124 accessions from the Colombian Musaceae Collection conserved in AGROSAVIA. We assessed 12 traits categorized into morphology, fruit quality, and yield, alongside sequence data. Our sequencing efforts provided valuable insights, with an average depth of about 7× per accession, resulting in 187,133 single-nucleotide polymorphisms (SNPs) against Musa acuminata (A genome) and 220,451 against Musa balbisiana (B genome). Population structure analysis grouped samples into four and five clusters based on the reference genome. By using different association models, we identified marker-trait associations (MTAs). The mixed linear model revealed four MTAs, while the Bayesian-information and linkage-disequilibrium iteratively nested keyway and fixed and random model for circulating probability unification models identified 82 and 70 MTAs, respectively. We identified 38 and 40 candidate genes in linkage proximity to significant MTAs for the A genome and B genome, respectively. Our findings provide insights into the genetic underpinnings of morphology, fruit quality, and yield. Once validated, the SNP markers and candidate genes can potentially drive advancements in genomic-guided breeding strategies to enhance banana crop improvement.

Insights into the molecular basis of reduced vancomycin susceptibility among three prominent Staphylococcus aureus clonal complexes

Staphylococcus aureus is a leading cause of healthcare-associated infections globally. Vancomycin-resistant S. aureus (VRSA), those with high-level resistance [minimum inhibitory concentration (MIC) of 16-32 µg/mL vancomycin], are uncommon, whereas vancomycin-intermediate S. aureus (VISA; MIC of 4-8 µg/mL), are isolated more frequently and develop during long-term and/or repeated use of the antibiotic. VISA can be difficult to eradicate and infections may persist. Our knowledge of mechanisms that underlie the development of VISA is incomplete. We used a genomics approach to investigate the VISA phenotype in three prominent S. aureus lineages. All VISA clinical isolates tested had increased cell wall thickness compared with vancomycin-susceptible S. aureus strains. Growth rates of clonal complex (CC) 5, CC8, and CC45 clinical isolates were reduced in 2 µg/mL vancomycin compared to media alone. Culture in 2 and 4 µg/mL vancomycin sequentially for two weeks reduced susceptibility to daptomycin, televancin, tigecycline, and vancomycin in a majority of CC5, CC8, and CC45 isolates tested. We identified alleles reported previously to contribute to the VISA phenotype, but unexpectedly, these alleles were unique to each CC. A subtherapeutic concentration of vancomycin elicited changes in the VISA transcriptome-common and unique-among the three CCs tested. Multiple genes, including those encoding a glycerate kinase, an M50 family metallopeptidase, and an uncharacterized membrane protein, were upregulated among all three lineages and not reported previously as associated with VISA. Although there are lineage-specific changes in DNA sequence, our findings suggest changes in the VISA transcriptome constitute a general response to stress that confers reduced susceptibility to multiple antibiotics.

IMPORTANCE

Our understanding of the mechanisms that underlie the development of vancomycin-intermediate Staphylococcus aureus (VISA) is incomplete. To provide a more comprehensive view of this process, we compared genome sequences of clonal complex (CC) 5, CC8, and CC45 VISA clinical isolates and measured changes in the transcriptomes of these isolates during culture with a subtherapeutic concentration of vancomycin. Notably, we identified differentially expressed genes that were lineage-specific or common to the lineages tested, including genes that have not been previously reported to contribute to a VISA phenotype. Changes in gene expression were accompanied by reduced growth rate, increased cell wall thickness, and reduced susceptibility to daptomycin, televancin, tigecycline, and vancomycin. Our results provide support to the idea that changes in gene expression contribute to the development of VISA among three CCs that are a prominent cause of human infections.

CRISPR-Cas9-mediated barcode insertion into Bacillus thuringiensis for surrogate tracking

The use of surrogate organisms can enable researchers to safely conduct research on pathogens and in a broader set of conditions. Being able to differentiate between the surrogates used in the experiments and background contamination as well as between different experiments will further improve research efforts. One effective approach is to introduce unique genetic barcodes into the surrogate genome and track their presence using the quantitative polymerase chain reaction (qPCR). In this report, we utilized the CRISPR-Cas9 methodology, which employs a single plasmid and a transformation step to insert five distinct barcodes into Bacillus thuringiensis, a well-established surrogate for Bacillus anthracis when Risk Group 1 organisms are needed. We subsequently developed qPCR assays for barcode detection and successfully demonstrated the stability of the barcodes within the genome through five cycles of sporulation and germination. Additionally, we conducted whole-genome sequencing on these modified strains and analyzed 187 potential Cas9 off-target sites. We found no correlation between the mutations observed in the engineered strains and the predicted off-target sites, suggesting this genome engineering strategy did not directly result in off-target mutations in the genome. This simple approach has the potential to streamline the creation of barcoded B. thuringiensis strains for use in future studies on surrogate genomes.

IMPORTANCE

The use of Bacillus anthracis as a biothreat agent poses significant challenges for public health and national security. Bacillus anthracis surrogates, like Bacillus thuringiensis, are invaluable tools for safely understanding Bacillus anthracis properties without the safety concerns that would arise from using a virulent strain of Bacillus anthracis. We report a simple method for barcode insertion into Bacillus thuringiensis using the CRISPR-Cas9 methodology and subsequent tracking by quantitative polymerase chain reaction (qPCR). Moreover, whole-genome sequencing data and CRISPR-Cas9 off-target analyses in Bacillus thuringiensis suggest that this gene-editing method did not directly cause unwanted mutations in the genome. This study should assist in the facile development of barcoded Bacillus thuringiensis surrogate strains, among other biotechnological applications in Bacillus species.

Comparative genomic sequencing to characterize Mycoplasma pneumoniae genome, typing, and drug resistance

To analyze the characteristics of Mycoplasma pneumoniae as well as macrolide antibiotic resistance through whole-genome sequencing and comparative genomics. Thirteen clinical strains isolated from 2003 to 2019 were selected, 10 of which were resistant to erythromycin (MIC >64 µg/mL), including 8 P1-type I and 2 P1-type II. Three were sensitive (<1 µg/mL) and P1-type II. One resistant strain had an A→G point mutation at position 2064 in region V of the 23S rRNA, the others had it at position 2063, while the three sensitive strains had no mutation here. Genome assembly and comparative genome analysis revealed a high level of genome consistency within the P1 type, and the primary differences in genome sequences concentrated in the region encoding the P1 protein. In P1-type II strains, three specific gene mutations were identified: C162A and A430G in L4 gene and T1112G mutation in the CARDS gene. Clinical information showed seven cases were diagnosed with severe pneumonia, all of which were infected with drug-resistant strains. Notably, BS610A4 and CYM219A1 exhibited a gene multi-copy phenomenon and shared a conserved functional domain with the DUF31 protein family. Clinically, the patients had severe refractory pneumonia, with pleural effusion, necessitating treatment with glucocorticoids and bronchoalveolar lavage. The primary variations between strains occur among different P1-types, while there is a high level of genomic consistency within P1-types. Three mutation loci associated with specific types were identified, and no specific genetic alterations directly related to clinical presentation were observed.IMPORTANCEMycoplasma pneumoniae is an important pathogen of community-acquired pneumonia, and macrolide resistance brings difficulties to clinical treatment. We analyzed the characteristics of M. pneumoniae as well as macrolide antibiotic resistance through whole-genome sequencing and comparative genomics. The work addressed primary variations between strains that occur among different P1-types, while there is a high level of genomic consistency within P1-types. In P1-type II strains, three specific gene mutations were identified: C162A and A430G in L4 gene and T1112G mutation in the CARDS gene. All the strains isolated from severe pneumonia cases were drug-resistant, two of which exhibited a gene multi-copy phenomenon, sharing a conserved functional domain with the DUF31 protein family. Three mutation loci associated with specific types were identified, and no specific genetic alterations directly related to clinical presentation were observed.

Whole-exome and whole-genome sequencing of 1064 individuals with type 1 diabetes reveals novel genes for diabetic kidney disease

AIMS/HYPOTHESIS

Diabetic kidney disease (DKD) is a severe diabetic complication that affects one third of individuals with type 1 diabetes. Although several genes and common variants have been shown to be associated with DKD, much of the predicted inheritance remains unexplained. Here, we performed next-generation sequencing to assess whether low-frequency variants, extending to a minor allele frequency (MAF) ≤10% (single or aggregated) contribute to the missing heritability in DKD.

METHODS

We performed whole-exome sequencing (WES) of 498 individuals and whole-genome sequencing (WGS) of 599 individuals with type 1 diabetes. After quality control, next-generation sequencing data were available for a total of 1064 individuals, of whom 541 had developed either severe albuminuria or end-stage kidney disease, and 523 had retained normal albumin excretion despite a long duration of type 1 diabetes. Single-variant and gene-aggregate tests for protein-altering variants (PAV) and protein-truncating variants (PTV) were performed separately for WES and WGS data and combined in a meta-analysis. We also performed genome-wide aggregate analyses on genomic windows (sliding window), promoters and enhancers using the WGS dataset.

RESULTS

In the single-variant meta-analysis, no variant reached genome-wide significance, but a suggestively associated common THAP7 rs369250 variant (p=1.50 × 10-5, MAF=49%) was replicated in the FinnGen general population genome-wide association study (GWAS) data for chronic kidney disease and DKD phenotypes. The gene-aggregate meta-analysis provided suggestive evidence (p<4.0 × 10-4) at four genes for DKD, of which NAT16 (MAFPAV≤10%) and LTA (also known as TNFβ, MAFPAV≤5%) are replicated in the FinnGen general population GWAS data. The LTA rs2229092 C allele was associated with significantly lower TNFR1, TNFR2 and TNFR3 serum levels in a subset of FinnDiane participants. Of the intergenic regions suggestively associated with DKD, the enhancer on chromosome 18q12.3 (p=3.94 × 10-5, MAFvariants≤5%) showed interaction with the METTL4 gene; the lead variant was replicated, and predicted to alter binding of the MafB transcription factor.

CONCLUSIONS/INTERPRETATION

Our sequencing-based meta-analysis revealed multiple genes, variants and regulatory regions that were suggestively associated with DKD. However, as no variant or gene reached genome-wide significance, further studies are needed to validate the findings.

Exome sequencing reveals neurodevelopmental genes in simplex consanguineous Iranian families with syndromic autism

BACKGROUND AND OBJECTIVE

Autosomal recessive genetic disorders pose significant health challenges in regions where consanguineous marriages are prevalent. The utilization of exome sequencing as a frequently employed methodology has enabled a clear delineation of diagnostic efficacy and mode of inheritance within multiplex consanguineous families. However, these aspects remain less elucidated within simplex families.

METHODS

In this study involving 12 unrelated simplex Iranian families presenting syndromic autism, we conducted singleton exome sequencing. The identified genetic variants were validated using Sanger sequencing, and for the missense variants in FOXG1 and DMD, 3D protein structure modeling was carried out to substantiate their pathogenicity. To examine the expression patterns of the candidate genes in the fetal brain, adult brain, and muscle, RT-qPCR was employed.

RESULTS

In four families, we detected an autosomal dominant gene (FOXG1), an autosomal recessive gene (CHKB), and two X-linked autism genes (IQSEC2 and DMD), indicating diverse inheritance patterns. In the remaining eight families, we were unable to identify any disease-associated genes. As a result, our variant detection rate stood at 33.3% (4/12), surpassing rates reported in similar studies of smaller cohorts. Among the four newly identified coding variants, three are de novo (heterozygous variant p.Trp546Ter in IQSEC2, heterozygous variant p.Ala188Glu in FOXG1, and hemizygous variant p.Leu211Met in DMD), while the homozygous variant p.Glu128Ter in CHKB was inherited from both healthy heterozygous parents. 3D protein structure modeling was carried out for the missense variants in FOXG1 and DMD, which predicted steric hindrance and spatial inhibition, respectively, supporting the pathogenicity of these human mutants. Additionally, the nonsense variant in CHKB is anticipated to influence its dimerization - crucial for choline kinase function - and the nonsense variant in IQSEC2 is predicted to eliminate three functional domains. Consequently, these distinct variants found in four unrelated individuals with autism are likely indicative of loss-of-function mutations.

CONCLUSIONS

In our two syndromic autism families, we discovered variants in two muscular dystrophy genes, DMD and CHKB. Given that DMD and CHKB are recognized for their participation in the non-cognitive manifestations of muscular dystrophy, it indicates that some genes transcend the boundary of apparently unrelated clinical categories, thereby establishing a novel connection between ASD and muscular dystrophy. Our findings also shed light on the complex inheritance patterns observed in Iranian consanguineous simplex families and emphasize the connection between autism spectrum disorder and muscular dystrophy. This underscores a likely genetic convergence between neurodevelopmental and neuromuscular disorders.

Evolution of Phytophthora infestans on its potato host since the Irish potato famine

Phytophthora infestans is a major oomycete plant pathogen, responsible for potato late blight, which led to the Irish Potato Famine from 1845-1852. Since then, potatoes resistant to this disease have been bred and deployed worldwide. Their resistance (R) genes recognize pathogen effectors responsible for virulence and then induce a plant response stopping disease progression. However, most deployed R genes are quickly overcome by the pathogen. We use targeted sequencing of effector and R genes on herbarium specimens to examine the joint evolution in both P. infestans and potato from 1845-1954. Currently relevant effectors are historically present in P. infestans, but with alternative alleles compared to modern reference genomes. The historic FAM-1 lineage has the virulent Avr1 allele and the ability to break the R1 resistance gene before breeders deployed it in potato. The FAM-1 lineage is diploid, but later, triploid US-1 lineages appear. We show that pathogen virulence genes and host resistance genes have undergone significant changes since the Famine, from both natural and artificial selection.

Evolution of the biosynthetic pathways of terpene scent compounds in roses

It is unknown why roses are terpene-rich, what the terpene biosynthetic pathways in roses are, and why only a few rose species produce the major components of rose essential oil. Here, we assembled two high-quality chromosome-level genomes for Rosa rugosa and Rosa multiflora. We also re-sequenced 132 individuals from the F1 progeny of Rosa chinensis and Rosa wichuraiana and 36 of their related species. Comparative genomics revealed that expansions of the 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) and terpene synthases (TPSs) gene families led to the enrichment of terpenes in rose scent components. We constructed a terpene biosynthesis network and discovered a TPS-independent citronellol biosynthetic pathway in roses through gene functional identification, genome-wide association studies (GWASs), and multi-omic analysis. Heterologous co-expression of rose citronellol biosynthetic genes in Nicotiana benthamiana led to citronellol production. Our genomic and metabolomic analyses suggested that the copy number of NUDX1-1a determines the citronellol content in different rose species. Our findings not only provide additional genome and gene resources and reveal the evolution of the terpene biosynthetic pathways but also present a nearly complete scenario for terpenoid metabolism that will facilitate the breeding of fragrant roses and the production of rose oil.

Allelic variability in the Rpp1 locus conferring resistance to Asian soybean rust revealed by genome-wide association

Soybean is a crucial crop for the Brazilian economy, but it faces challenges from the biotrophic fungus Phakopsora pachyrhizi, which causes Asian Soybean Rust (ASR). In this study, we aimed to identify SNPs associated with resistance within the Rpp1 locus, which is effective against Brazilian ASR populations. We employed GWAS and re-sequencing analyzes to pinpoint SNP markers capable of differentiating between soybean accessions harboring the Rpp1, Rpp1-b and other alternative alleles in the Rpp1 locus and from susceptible soybean cultivars. Seven SNP markers were found to be associated with ASR resistance through GWAS, with three of them defining haplotypes that efficiently distinguished the accessions based on their ASR resistance and source of the Rpp gene. These haplotypes were subsequently validated using a bi-parental population and a diverse set of Rpp sources, demonstrating that the GWAS markers co-segregate with ASR resistance. We then examined the presence of these haplotypes in a diverse set of soybean genomes worldwide, finding a few new potential sources of Rpp1/Rpp1-b. Further genomic sequence analysis revealed nucleotide differences within the genes present in the Rpp1 locus, including the ULP1-NBS-LRR genes, which are potential R gene candidates. These results provide valuable insights into ASR resistance in soybean, thus helping the development of resistant soybean varieties through genetic breeding programs.

European farmhouse brewing yeasts form a distinct genetic group

The brewing industry is constantly evolving, driven by the quest for novel flavours and fermentation characteristics that cater to evolving consumer preferences. This study explores the genetic and phenotypic diversity of European farmhouse yeasts, traditionally used in rural brewing practices and maintained outside of pure culture industrial yeast selection. We isolated landrace brewing yeast strains from diverse geographical locations across Europe, including Norway, Lithuania, Latvia, and Russia, and also included African farmhouse brewing strains from Ghana. Our genomic analysis using long-read and short-read whole genome sequencing uncovered a genetically distinct group that diverges from industrial brewing yeasts. This group, which is closely related to ale brewing strains, is preliminarily named the 'European Farmhouse' group and shows greater predicted admixture from Asian fermentation strains. Through genomic and phenotypic analyses, including flavour metabolite analysis via headspace gas chromatography-mass spectrometry, sugar metabolite analysis via high-performance liquid chromatography, and wort fermentation analysis, we found a broad spectrum of fermentation capabilities, from rapid and efficient fermentation to unique aroma and flavour compound profiles, potentially offering novel traits for brewing applications. This study highlights the importance of preservation of brewing cultural heritage knowledge and resources including yeast cultures. KEY POINTS: • A large set of geographically diverse farmhouse brewing strains were characterized • Norwegian and Baltic farmhouse brewing strains form a distinct genetic group • Farmhouse strains show considerable diversity in fermentation and flavour formation.

Antifungal heteroresistance causes prophylaxis failure and facilitates breakthrough Candida parapsilosis infections

Breakthrough fungal infections in patients on antimicrobial prophylaxis during allogeneic hematopoietic cell transplantation (allo-HCT) represent a major and often unexplained cause of morbidity and mortality. Candida parapsilosis is a common cause of invasive candidiasis and has been classified as a high-priority fungal pathogen by the World Health Organization. In high-risk allo-HCT recipients on micafungin prophylaxis, we show that heteroresistance (the presence of a phenotypically unstable, low-frequency subpopulation of resistant cells (~1 in 10,000)) underlies breakthrough bloodstream infections by C. parapsilosis. By analyzing 219 clinical isolates from North America, Europe and Asia, we demonstrate widespread micafungin heteroresistance in C. parapsilosis. Standard antimicrobial susceptibility tests, such as broth microdilution or gradient diffusion assays, which guide drug selection for invasive infections, fail to detect micafungin heteroresistance in C. parapsilosis. To facilitate rapid detection of micafungin heteroresistance in C. parapsilosis, we constructed a predictive machine learning framework that classifies isolates as heteroresistant or susceptible using a maximum of ten genomic features. These results connect heteroresistance to unexplained antifungal prophylaxis failure in allo-HCT recipients and demonstrate a proof-of-principle diagnostic approach with the potential to guide clinical decisions and improve patient care.

The re-emergence of sexually transmissible multidrug resistant Shigella flexneri 3a, England, United Kingdom

Shigellosis is an enteric infection that transmits through the faecal-oral route, which can occur during sex between men who have sex with men (MSM). Between 2009 and 2014, an epidemic of sexually transmissible Shigella flexneri 3a occurred in England that subsequently declined. However, from 2019 to 2021, despite SARS-CoV-2 restrictions, S. flexneri 3a continued to re-emerge. We explored possible drivers of re-emergence by comparing host demography and pathogen genomics. Cases were primarily among 35-64 year old men in London. Genomic analyses of 502 bacterial isolates showed that the majority (58%) of re-emerging MSM strains were a clonal replacement of the original, with reduced antimicrobial resistance, conservation of plasmid col156_1, and two SNPs with 19 predicted effects. The absence of major changes in the pathogen or host demographics suggest that other factors may have driven the re-emergence of S. flexneri 3a and highlight the need for further work in the area.

Scots pine - panmixia and the elusive signal of genetic adaptation

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Orchidinae-205: A new genome-wide custom bait set for studying the evolution, systematics, and trade of terrestrial orchids

Terrestrial orchids are a group of genetically understudied, yet culturally and economically important plants. The Orchidinae tribe contains many species that produce edible tubers that are used for the production of traditional delicacies collectively called 'salep'. Overexploitation of wild orchids in the Eastern Mediterranean and Western Asia threatens to drive many of these species to extinction, but cost-effective tools for monitoring their trade are currently lacking. Here we present a custom bait kit for target enrichment and sequencing of 205 novel genetic markers that are tailored to phylogenomic applications in Orchidinae s.l. A subset of 31 markers capture genes putatively involved in the production of glucomannan, a water-soluble polysaccharide that gives salep its distinctive properties. We tested the kit on 73 taxa native to the area, demonstrating universally high locus recovery irrespective of species identity, that exceeds the total sequence length obtained with alternative kits currently available. Phylogenetic inference with concatenation and coalescent approaches was robust and showed high levels of support for most clades, including some which were previously unresolved. Resolution for hybridizing and recently radiated lineages remains difficult, but could be further improved by analysing multiple haplotypes and the non-exonic sequences captured by our kit, with the promise to shed new light on the evolution of enigmatic taxa with a complex speciation history. Offering a step-up from traditional barcoding and universal markers, the genome-wide custom loci targeted by Orchidinae-205 are a valuable new resource to study the evolution, systematics and trade of terrestrial orchids.

Genotypic and phenotypic characterisation of a nosocomial outbreak of Candida auris in Spain during 5 years

OBJECTIVES

The investigation of Candida auris outbreaks is needed to provide insights into its population structure and transmission dynamics. We genotypically and phenotypically characterised a C. auris nosocomial outbreak occurred in Consorcio Hospital General Universitario de Valencia (CHGUV), Spain.

METHODS

Data and isolates were collected from CHGUV from September 2017 (first case) until September 2021. Thirty-five isolates, including one from an environmental source, were randomly selected for whole genome sequencing (WGS), and the genomes were analysed along with a database with 335 publicly available genomes, assigning them to one of the five major clades. In order to identify polymorphisms associated with drug resistance, we used the fully susceptible GCA_003014415.1 strain as reference sequence. Known mutations in genes ERG11 and FKS1 conferring resistance to fluconazole and echinocandins, respectively, were investigated. Isolates were classified into aggregating or non-aggregating.

RESULTS

All isolates belonged to clade III and were from an outbreak with a single origin. They clustered close to three publicly available genomes from a hospital from where the first patient was transferred, being the probable origin. The mutation VF125AL in the ERG11 gene, conferring resistance to fluconazole, was present in all the isolates and one isolate also carried the mutation S639Y in the FKS1 gene. All the isolates had a non-aggregating phenotype (potentially more virulent).

CONCLUSIONS

Isolates are genotypically related and phenotypically identical but one with resistance to echinocandins, which seems to indicate that they all belong to an outbreak originated from a single isolate, remaining largely invariable over the years. This result stresses the importance of implementing infection control practices as soon as the first case is detected or when a patient is transferred from a setting with known cases.

Forest genomics in the Caucasus through the lens of its dominant tree species - Fagus orientalis

The last glacial period is known to have greatly influenced the demographic history of temperate forest trees, with important range contractions and post-glacial expansions that led to the formation of multiple genetic lineages and secondary contact zones in the Northern Hemisphere. These dynamics have been extensively studied for European and North American species but are still poorly understood in other temperate regions of rich biodiversity such as the Caucasus. Our study helps filling that gap by deciphering the genomic landscapes of F. orientalis across the South Caucasus. The use of genome-wide data confirmed a past demographic history strongly influenced by the Last Glacial Maximum, revealing two disjunct glacial refugia in the Colchis and Hyrcanian regions. The resulting patterns of genetic diversity, load and differentiation are not always concordant across the region, with genetic load pinpointing the location of the glacial refugia more efficiently than genetic diversity alone. The Hyrcanian forests show depleted genetic diversity and substantial isolation, even if long-distance gene flow is still present with the main centre of diversity in the Greater Caucasus. Finally, we characterize a strong heterogeneity of genetic diversity and differentiation along the species chromosomes, with noticeably a first chromosome showing low diversity and weak differentiation.

Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits

The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveals that those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.

DDX41 haploinsufficiency causes inefficient hematopoiesis under stress and cooperates with p53 mutations to cause hematologic malignancy

Germline heterozygous mutations in DDX41 predispose individuals to hematologic malignancies in adulthood. Most of these DDX41 mutations result in a truncated protein, leading to loss of protein function. To investigate the impact of these mutations on hematopoiesis, we generated mice with hematopoietic-specific knockout of one Ddx41 allele. Under normal steady-state conditions, there was minimal effect on lifelong hematopoiesis, resulting in a mild yet persistent reduction in red blood cell counts. However, stress induced by transplantation of the Ddx41+/- BM resulted in hematopoietic stem/progenitor cell (HSPC) defects and onset of hematopoietic failure upon aging. Transcriptomic analysis of HSPC subsets from the transplanted BM revealed activation of cellular stress responses, including upregulation of p53 target genes in erythroid progenitors. To understand how the loss of p53 affects the phenotype of Ddx41+/- HSPCs, we generated mice with combined Ddx41 and Trp53 heterozygous deletions. The reduction in p53 expression rescued the fitness defects in HSPC caused by Ddx41 heterozygosity. However, the combined Ddx41 and Trp53 mutant mice were prone to developing hematologic malignancies that resemble human myelodysplastic syndrome and acute myeloid leukemia. In conclusion, DDX41 heterozygosity causes dysregulation of the response to hematopoietic stress, which increases the risk of transformation with a p53 mutation.

Historical museum samples reveal signals of selection and drift in response to changing insecticide use in an agricultural pest moth

In response to environmental and human-imposed selective pressures, agroecosystem pests frequently undergo rapid evolution, with some species having a remarkable capacity to rapidly develop pesticide resistance. Temporal sampling of genomic data can comprehensively capture such adaptive changes over time, for example, by elucidating allele frequency shifts in pesticide resistance loci in response to different pesticides. Here, we leveraged museum specimens spanning over a century of collections to generate temporal contrasts between pre- and post-insecticide populations of an agricultural pest moth, Helicoverpa armigera. We used targeted exon sequencing of 254 samples collected across Australia from the pre-1950s (prior to insecticide introduction) to the 1990s, encompassing decades of changing insecticide use. Our sequencing approach focused on genes that are known to be involved in insecticide resistance, environmental sensation, and stress tolerance. We found an overall lack of spatial and temporal population structure change across Australia. In some decades (e.g., 1960s and 1970s), we found a moderate reduction of genetic diversity, implying stochasticity in evolutionary trajectories due to genetic drift. Temporal genome scans showed extensive evidence of selection following insecticide use, although the majority of selected variants were low impact. Finally, alternating trajectories of allele frequency change were suggestive of potential antagonistic pleiotropy. Our results provide new insights into recent evolutionary responses in an agricultural pest and show how temporal contrasts using museum specimens can improve mechanistic understanding of rapid evolution.

Association mapping with a diverse population of Puccinia graminis f. sp. tritici identified avirulence loci interacting with the barley Rpg1 stem rust resistance gene

BACKGROUND

Wheat stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is an important disease of barley and wheat. A diverse sexual Pgt population from the Pacific Northwest (PNW) region of the US contains a high proportion of individuals with virulence on the barley stem rust resistance (R) gene, Rpg1. However, the evolutionary mechanisms of this virulence on Rpg1 are mysterious considering that Rpg1 had not been deployed in the region and the gene had remained remarkably durable in the Midwestern US and prairie provinces of Canada.

METHODS AND RESULTS

To identify AvrRpg1 effectors, genome wide association studies (GWAS) were performed using 113 Pgt isolates collected from the PNW (n = 89 isolates) and Midwest (n = 24 isolates) regions of the US. Disease phenotype data were generated on two barley lines Morex and the Golden Promise transgenic (H228.2c) that carry the Rpg1 gene. Genotype data was generated by whole genome sequencing (WGS) of 96 isolates (PNW = 89 isolates and Midwest = 7 isolates) and RNA sequencing (RNAseq) data from 17 Midwestern isolates. Utilizing ~1.2 million SNPs generated from WGS and phenotype data (n = 96 isolates) on the transgenic line H228.2c, 53 marker trait associations (MTAs) were identified. Utilizing ~140 K common SNPs generated from combined analysis of WGS and RNAseq data, two significant MTAs were identified using the cv Morex phenotyping data. The 55 MTAs defined two distinct avirulence loci, on supercontig 2.30 and supercontig 2.11 of the Pgt reference genome of Pgt isolate CRL 75-36-700-3. The major avirulence locus designated AvrRpg1A was identified with the GWAS using both barley lines and was delimited to a 35 kb interval on supercontig 2.30 containing four candidate genes (PGTG_10878, PGTG_10884, PGTG_10885, and PGTG_10886). The minor avirulence locus designated AvrRpg1B identified with cv Morex contained a single candidate gene (PGTG_05433). AvrRpg1A haplotype analysis provided strong evidence that a dominant avirulence gene underlies the locus.

CONCLUSIONS

The association analysis identified strong candidate AvrRpg1 genes. Further analysis to validate the AvrRpg1 genes will fill knowledge gaps in our understanding of rust effector biology and the evolution and mechanism/s of Pgt virulence on Rpg1.

FOXM1 c.1205 C > A mutation is associated with unilateral Moyamoya disease and inhibits angiogenesis in human brain endothelial cells

Unilateral moyamoya disease (MMD) represents a distinct subtype characterised by occlusive changes in the circle of Willis and abnormal vascular network formation. However, the aetiology and pathogenesis of unilateral MMD remain unclear. In this study, genetic screening of a family with unilateral MMD using whole-genome sequencing helped identify the c.1205 C > A variant of FOXM1, which encodes the transcription factor FOXM1 and plays a crucial role in angiogenesis and cell proliferation, as a susceptibility gene mutation. We demonstrated that this mutation significantly attenuated the proangiogenic effects of FOXM1 in human brain endothelial cells, leading to reduced proliferation, migration, and tube formation. Furthermore, FOXM1 c.1205 C > A results in increased apoptosis of human brain endothelial cells, mediated by the downregulation of the transcription of the apoptosis-inhibiting protein BCL2. These results suggest a potential role for the FOXM1 c.1205 C > A mutation in the pathogenesis of unilateral MMD and may contribute to the understanding and treatment of this condition.

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Isoform-resolved genome annotation enables mapping of tissue-specific betalain regulation in amaranth

Betalains are coloring pigments produced in some families of the order Caryophyllales, where they replace anthocyanins as coloring pigments. While the betalain pathway itself is well studied, the tissue-specific regulation of the pathway remains mostly unknown. We enhance the high-quality Amaranthus hypochondriacus reference genome and produce a substantially more complete genome annotation, incorporating isoform details. We annotate betalain and anthocyanin pathway genes along with their regulators in amaranth and map the genetic control and tissue-specific regulation of the betalain pathway. Our improved genome annotation allowed us to identify causal mutations that lead to a knock-out of red betacyanins in natural accessions of amaranth. We reveal the tissue-specific regulation of flower color via a previously uncharacterized MYB transcription factor, AhMYB2. Downregulation of AhMYB2 in the flower leads to reduced expression of key betalain enzyme genes and loss of red flower color. Our improved amaranth reference genome represents the most complete genome of amaranth to date and is a valuable resource for betalain and amaranth research. High similarity of the flower betalain regulator AhMYB2 to anthocyanin regulators and a partially conserved interaction motif support the co-option of anthocyanin regulators for the betalain pathway as a possible reason for the mutual exclusiveness of the two pigments.

Computational tools for plant genomics and breeding

Plant genomics and crop breeding are at the intersection of biotechnology and information technology. Driven by a combination of high-throughput sequencing, molecular biology and data science, great advances have been made in omics technologies at every step along the central dogma, especially in genome assembling, genome annotation, epigenomic profiling, and transcriptome profiling. These advances further revolutionized three directions of development. One is genetic dissection of complex traits in crops, along with genomic prediction and selection. The second is comparative genomics and evolution, which open up new opportunities to depict the evolutionary constraints of biological sequences for deleterious variant discovery. The third direction is the development of deep learning approaches for the rational design of biological sequences, especially proteins, for synthetic biology. All three directions of development serve as the foundation for a new era of crop breeding where agronomic traits are enhanced by genome design.

The potential influence of genome-wide adaptive divergence on conservation translocation outcome in an isolated greater sage-grouse population

Conservation translocations are an important conservation tool commonly employed to augment declining or reestablish extirpated populations. One goal of augmentation is to increase genetic diversity and reduce the risk of inbreeding depression (i.e., genetic rescue). However, introducing individuals from significantly diverged populations risks disrupting coadapted traits and reducing local fitness (i.e., outbreeding depression). Genetic data are increasingly more accessible for wildlife species and can provide unique insight regarding the presence and retention of introduced genetic variation from augmentation as an indicator of effectiveness and adaptive similarity as an indicator of source and recipient population suitability. We used 2 genetic data sets to evaluate augmentation of isolated populations of greater sage-grouse (Centrocercus urophasianus) in the northwestern region of the species range (Washington, USA) and to retrospectively evaluate adaptive divergence among source and recipient populations. We developed 2 statistical models for microsatellite data to evaluate augmentation outcomes. We used one model to predict genetic diversity after augmentation and compared these predictions with observations of genetic change. We used the second model to quantify the amount of observed reproduction attributed to transplants (proof of population integration). We also characterized genome-wide adaptive divergence among source and recipient populations. Observed genetic diversity (HO = 0.65) was higher in the recipient population than predicted had no augmentation occurred (HO = 0.58) but less than what was predicted by our model (HO = 0.75). The amount of shared genetic variation between the 2 geographically isolated resident populations increased, which is evidence of periodic gene flow previously assumed to be rare. Among candidate adaptive genes associated with elevated fixation index (FST) (143 genes) or local environmental variables (97 and 157 genes for each genotype-environment association method, respectively), we found clusters of genes with related functions that may influence the ability of transplants to use local resources and navigate unfamiliar environments and their reproductive potential, all possible reasons for low genetic retention from augmentation.

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

Genome resequencing reveals population divergence and local adaptation of blacklegged ticks in the United States

Tick vectors and tick-borne disease are increasingly impacting human populations globally. An important challenge is to understand tick movement patterns, as this information can be used to improve management and predictive modelling of tick population dynamics. Evolutionary analysis of genetic divergence, gene flow and local adaptation provides insight on movement patterns at large spatiotemporal scales. We develop low coverage, whole genome resequencing data for 92 blacklegged ticks, Ixodes scapularis, representing range-wide variation across the United States. Through analysis of population genomic data, we find that tick populations are structured geographically, with gradual isolation by distance separating three population clusters in the northern United States, southeastern United States and a unique cluster represented by a sample from Tennessee. Populations in the northern United States underwent population contractions during the last glacial period and diverged from southern populations at least 50 thousand years ago. Genome scans of selection provide strong evidence of local adaptation at genes responding to host defences, blood-feeding and environmental variation. In addition, we explore the potential of low coverage genome sequencing of whole-tick samples for documenting the diversity of microbial pathogens and recover important tick-borne pathogens such as Borrelia burgdorferi. The combination of isolation by distance and local adaptation in blacklegged ticks demonstrates that gene flow, including recent expansion, is limited to geographical scales of a few hundred kilometres.

SDHAF1 confers metabolic resilience to aging hematopoietic stem cells by promoting mitochondrial ATP production

Aging generally predisposes stem cells to functional decline, impairing tissue homeostasis. Here, we report that hematopoietic stem cells (HSCs) acquire metabolic resilience that promotes cell survival. High-resolution real-time ATP analysis with glucose tracing and metabolic flux analysis revealed that old HSCs reprogram their metabolism to activate the pentose phosphate pathway (PPP), becoming more resistant to oxidative stress and less dependent on glycolytic ATP production at steady state. As a result, old HSCs can survive without glycolysis, adapting to the physiological cytokine environment in bone marrow. Mechanistically, old HSCs enhance mitochondrial complex II metabolism during stress to promote ATP production. Furthermore, increased succinate dehydrogenase assembly factor 1 (SDHAF1) in old HSCs, induced by physiological low-concentration thrombopoietin (TPO) exposure, enables rapid mitochondrial ATP production upon metabolic stress, thereby improving survival. This study provides insight into the acquisition of resilience through metabolic reprogramming in old HSCs and its molecular basis to ameliorate age-related hematopoietic abnormalities.

Comparative genome-wide association study on body weight in Chinese native ducks using four models

The Jinling White duck represents a newly developed breed characterized by a rapid growth rate and a superior meat quality, offering significant economic value and research potential; however, the genetic basis underlying their body weight traits remains less understood. Here, we performed whole-genome resequencing for 201 diverse Jinling White male ducks and conducted population genomic analyses, suggesting a rich genetic diversity within the Jinling White duck population. Equipped with our genomic resources, we applied genome-wide association analysis for body weight on birth (BWB), body weight on 1 wk (BW1), body weight on 3 wk (BW3), body weight on 5 wk (BW5) and body weight on 7 wk (BW7) using 4 statistical models. Comparative studies indicated that factored spectrally transformed linear mixed models (FaST-LMM) demonstrated the most superior efficiency, yielding more results with the minimal false positives. We discovered that PUS7, FBXO11, FOXN2, MSH6, and SLC4A4 were associated with BWB. RAG2, and TMEFF2 were candidate genes for BW1, and STARD13, Klotho, ZAR1L are likely candidates for BW3 and BW5. PLXNC1, ATP1A1, CD58, FRYL, OCIAD1, and OCIAD2 were linked to BW7. These findings provide a genetic reference for the selection and breeding of Jinling White ducks, while also deepened our understanding of Growth and development phenotypic in ducks.

Intragenic MFSD8 duplication and histopathological findings in a rabbit with neuronal ceroid lipofuscinosis

Neuronal ceroid lipofuscinoses (NCL) are among the most prevalent neurodegenerative disorders of early life in humans. Disease-causing variants have been described for 13 different NCL genes. In this study, a refined pathological characterization of a female rabbit with progressive neurological signs reminiscent of NCL was performed. Cytoplasmic pigment present in neurons was weakly positive with Sudan black B and autofluorescent. Immunohistology revealed astrogliosis, microgliosis and axonal degeneration. During the subsequent genetic investigation, the genome of the affected rabbit was sequenced and examined for private variants in NCL candidate genes. The analysis revealed a homozygous ~10.7 kb genomic duplication on chromosome 15 comprising parts of the MFSD8 gene, NC_013683.1:g.103,727,963_103,738,667dup. The duplication harbors two internal protein coding exons and is predicted to introduce a premature stop codon into the transcript, truncating ~50% of the wild-type MFSD8 open reading frame encoding the major facilitator superfamily domain containing protein 8, XP_002717309.2:p.(Glu235Leufs*23). Biallelic loss-of-function variants in MFSD8 have been described to cause NCL7 in human patients, dogs and a single cat. The available clinical and pathological data, together with current knowledge about MFSD8 variants and their functional impact in other species, point to the MFSD8 duplication as a likely causative defect for the observed phenotype in the affected rabbit.

A new Finnish flavor of feline coat coloration, "salmiak," is associated with a 95-kb deletion downstream of the KIT gene

Cats with a distinctive white hair pattern of unknown molecular cause have been discovered in the Finnish domestic cat population. Based on the unique appearance of these cats, we have named this phenotype salmiak ("salty licorice"). The use of a commercially available panel test to genotype four salmiak-colored cats revealed the absence of all known variants associated with white-haired phenotypic loci: full White (W), Spotting (Ws) and the Birman white Gloves associated (wg) allele of the KIT proto-oncogene (KIT) gene. Whole-genome sequencing on two salmiak-colored cats was conducted to search for candidate causal variants in the KIT gene. Despite a lack of coding variants, visual inspection of the short read alignments revealed a large ~95 kb deletion located ~65 kb downstream of the KIT gene in the salmiak cats. Additional PCR genotyping of 180 domestic cats and three salmiak-colored cats confirmed the homozygous derived variant genotype fully concordant with the salmiak phenotype. We suggest the newly identified variant be designated as wsal for "w salmiak".

Reovirus infection induces transcriptome-wide unique A-to-I editing changes in the murine fibroblasts

The conversion of Adenosine (A) to Inosine (I), by Adenosine Deaminases Acting on RNA or ADARs, is an essential post-transcriptional modification that contributes to proteome diversity and regulation in metazoans including humans. In addition to its transcriptome-regulating role, ADARs also play a major part in immune response to viral infection, where an interferon response activates interferon-stimulated genes, such as ADARp150, in turn dynamically regulating host-virus interactions. A previous report has shown that infection from reoviruses, despite strong activation of ADARp150, does not influence the editing of some of the major known editing targets, while likely editing others, suggesting a potentially nuanced editing pattern that may depend on different factors. However, the results were based on a handful of selected editing sites and did not cover the entire transcriptome. Thus, to determine whether and how reovirus infection specifically affects host ADAR editing patterns, we analyzed a publicly available deep-sequenced RNA-seq dataset, from murine fibroblasts infected with wild-type and mutant reovirus strains that allowed us to examine changes in editing patterns on a transcriptome-wide scale. To the best of our knowledge, this is the first transcriptome-wide report on host editing changes after reovirus infection. Our results demonstrate that reovirus infection induces unique nuanced editing changes in the host, including introducing sites uniquely edited in infected samples. Genes with edited sites are overrepresented in pathways related to immune regulation, cellular signaling, metabolism, and growth. Moreover, a shift in editing targets has also been observed, where the same genes are edited in infection and control conditions but at different sites, or where the editing rate is increased for some and decreased for other differential targets, supporting the hypothesis of dynamic and condition-specific editing by ADARs.

Site-directed genotype screening for elimination of antinutritional saponins in quinoa seeds identifies TSARL1 as a master controller of saponin biosynthesis selectively in seeds

Climate change may result in a drier climate and increased salinization, threatening agricultural productivity worldwide. Quinoa (Chenopodium quinoa) produces highly nutritious seeds and tolerates abiotic stresses such as drought and high salinity, making it a promising future food source. However, the presence of antinutritional saponins in their seeds is an undesirable trait. We mapped genes controlling seed saponin content to a genomic region that includes TSARL1. We isolated desired genetic variation in this gene by producing a large mutant library of a commercial quinoa cultivar and screening the library for specific nucleotide substitutions using droplet digital PCR. We were able to rapidly isolate two independent tsarl1 mutants, which retained saponins in the leaves and roots for defence, but saponins were undetectable in the seed coat. We further could show that TSARL1 specifically controls seed saponin biosynthesis in the committed step after 2,3-oxidosqualene. Our work provides new important knowledge on the function of TSARL1 and represents a breakthrough for quinoa breeding.

A Proof-of-Concept Study of Sequential Treatment with the HDAC Inhibitor Vorinostat following BRAF and MEK Inhibitors in BRAFV600-Mutated Melanoma

PURPOSE

The development of resistance limits the clinical benefit of BRAF and MEK inhibitors (BRAFi/MEKi) in BRAFV600-mutated melanoma. It has been shown that short-term treatment (14 days) with vorinostat was able to initiate apoptosis of resistant tumor cells. We aimed to assess the antitumor activity of sequential treatment with vorinostat following BRAFi/MEKi in patients with BRAFV600-mutated melanoma who progressed after initial response to BRAFi/MEKi.

PATIENTS AND METHODS

Patients with BRAFi/MEKi-resistant BRAFV600-mutated melanoma were treated with vorinostat 360 mg once daily for 14 days followed by BRAFi/MEKi. The primary endpoint was an objective response rate of progressive lesions of at least 30% according to Response Evaluation Criteria in Solid Tumors 1.1. Secondary endpoints included progression-free survival, overall survival, safety, pharmacokinetics of vorinostat, and translational molecular analyses using ctDNA and tumor biopsies.

RESULTS

Of the 26 patients with progressive BRAFi/MEKi-resistant BRAFV600-mutated melanoma receiving treatment with vorinostat, 22 patients were evaluable for response. The objective response rate was 9%, with one complete response for 31.2 months and one partial response for 14.9 months. Median progression-free survival and overall survival were 1.4 and 5.4 months, respectively. Common adverse events were fatigue (23%) and nausea (19%). ctDNA analysis showed emerging secondary mutations in NRAS and MEK in eight patients at the time of BRAFi/MEKi resistance. Elimination of these mutations by vorinostat treatment was observed in three patients.

CONCLUSIONS

Intermittent treatment with vorinostat in patients with BRAFi/MEKi-resistant BRAFV600-mutated melanoma is well tolerated. Although the primary endpoint of this study was not met, durable antitumor responses were observed in a minority of patients (9%).

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus)

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.