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

Role of U11/U12 minor spliceosome gene ZCRB1 in Ciliogenesis and WNT Signaling

Despite the fact that 0.5% of human introns are processed by the U11/U12 minor spliceosome, the latter influences gene expression across multiple cellular processes. The ZCRB1 protein is a recently described core component of the U12 mono-snRNP minor spliceosome, but its functional significance to minor splicing, gene regulation, and biological signaling cascades is poorly understood. Using CRISPR-Cas9 and siRNA targeted knockout and knockdown strategies, we show that human cell lines with a partial reduction in ZCRB1 expression exhibit significant dysregulation of the splicing and expression of U12-type genes, primarily due to dysregulation of U12 mono-snRNA. RNA-Seq and targeted analyses of minor intron-containing genes indicate a downregulation in the expression of genes involved in ciliogenesis, and consequentially an upregulation in WNT signaling. Additionally, zcrb1 CRISPR-Cas12a knockdown in zebrafish embryos led to gross developmental and body axis abnormalities, disrupted ciliogenesis, and upregulated WNT signaling, complementing our human cell studies. This work highlights a conserved and essential biological role of the minor spliceosome in general, and the ZCRB1 protein specifically in cellular and developmental processes across species, shedding light on the multifaceted relationship between splicing regulation, ciliogenesis, and WNT signaling.

Frequency of pathogenic germline variants in pediatric medulloblastoma survivors

Background

Medulloblastoma is the most common malignant brain tumor in children. Most cases are sporadic, but well characterized germline alterations in APC, ELP1, GPR161, PTCH1, SUFU, and TP53 predispose to medulloblastoma. However, knowledge about pathogenic/likely pathogenic (P/LP) variants that predispose to medulloblastoma vary based on genes evaluated, patient demographics, and pathogenicity definitions.

Methods

Germline exome sequencing was conducted on 160 childhood survivors of medulloblastoma. Analyses focused on rare variants in 239 known cancer susceptibility genes (CSGs). P/LP variants were identified using ClinVar and InterVar. Variants of unknown significance in known medulloblastoma predisposing genes (APC, ELP1, GPR161, PTCH1, SUFU, TP53) were further classified for loss of function variants. We compared the frequency of P/LP variants in cases to that in 1,259 cancer-free adult controls.

Results

Twenty cases (12.5%) had a P/LP variant in an autosomal dominant CSG versus 5% in controls (p=1.0 x10-3), and 10 (6.3%) of these were P/LP variants in a known medulloblastoma gene, significantly greater than 0.2% observed in controls (p=1.4x10-8). The CSGs with the most P/LP variants in cases, and significantly higher than controls, were ELP1 (p=3.0x10-4) and SUFU (p=1.4x10-3).

Conclusion

Approximately one in eight pediatric medulloblastoma survivors had an autosomal dominant P/LP CSG variant. We confirm several known associated genes and identify novel genes that may be important in medulloblastoma.

Genomic landscape of hepatocellular carcinoma in Egyptian patients by whole exome sequencing

BACKGROUND

Hepatocellular carcinoma (HCC) is the most common primary liver cancer. Chronic hepatitis and liver cirrhosis lead to accumulation of genetic alterations driving HCC pathogenesis. This study is designed to explore genomic landscape of HCC in Egyptian patients by whole exome sequencing.

METHODS

Whole exome sequencing using Ion Torrent was done on 13 HCC patients, who underwent surgical intervention (7 patients underwent living donor liver transplantation (LDLT) and 6 patients had surgical resection}.

RESULTS

Mutational signature was mostly S1, S5, S6, and S12 in HCC. Analysis of highly mutated genes in both HCC and Non-HCC revealed the presence of highly mutated genes in HCC (AHNAK2, MUC6, MUC16, TTN, ZNF17, FLG, MUC12, OBSCN, PDE4DIP, MUC5b, and HYDIN). Among the 26 significantly mutated HCC genes-identified across 10 genome sequencing studies-in addition to TCGA, APOB and RP1L1 showed the highest number of mutations in both HCC and Non-HCC tissues. Tier 1, Tier 2 variants in TCGA SMGs in HCC and Non-HCC (TP53, PIK3CA, CDKN2A, and BAP1). Cancer Genome Landscape analysis revealed Tier 1 and Tier 2 variants in HCC (MSH2) and in Non-HCC (KMT2D and ATM). For KEGG analysis, the significantly annotated clusters in HCC were Notch signaling, Wnt signaling, PI3K-AKT pathway, Hippo signaling, Apelin signaling, Hedgehog (Hh) signaling, and MAPK signaling, in addition to ECM-receptor interaction, focal adhesion, and calcium signaling. Tier 1 and Tier 2 variants KIT, KMT2D, NOTCH1, KMT2C, PIK3CA, KIT, SMARCA4, ATM, PTEN, MSH2, and PTCH1 were low frequency variants in both HCC and Non-HCC.

CONCLUSION

Our results are in accordance with previous studies in HCC regarding highly mutated genes, TCGA and specifically enriched pathways in HCC. Analysis for clinical interpretation of variants revealed the presence of Tier 1 and Tier 2 variants that represent potential clinically actionable targets. The use of sequencing techniques to detect structural variants and novel techniques as single cell sequencing together with multiomics transcriptomics, metagenomics will integrate the molecular pathogenesis of HCC in Egyptian patients.

A synthetic peptide mimic kills Candida albicans and synergistically prevents infection

More than two million people worldwide are affected by life-threatening, invasive fungal infections annually. Candida species are the most common cause of nosocomial, invasive fungal infections and are associated with mortality rates above 40%. Despite the increasing incidence of drug-resistance, the development of novel antifungal formulations has been limited. Here we investigate the antifungal mode of action and therapeutic potential of positively charged, synthetic peptide mimics to combat Candida albicans infections. Our data indicates that these synthetic polymers cause endoplasmic reticulum stress and affect protein glycosylation, a mode of action distinct from currently approved antifungal drugs. The most promising polymer composition damaged the mannan layer of the cell wall, with additional membrane-disrupting activity. The synergistic combination of the polymer with caspofungin prevented infection of human epithelial cells in vitro, improved fungal clearance by human macrophages, and significantly increased host survival in a Galleria mellonella model of systemic candidiasis. Additionally, prolonged exposure of C. albicans to the synergistic combination of polymer and caspofungin did not lead to the evolution of tolerant strains in vitro. Together, this work highlights the enormous potential of these synthetic peptide mimics to be used as novel antifungal formulations as well as adjunctive antifungal therapy.

Clinical and Genomic Perspective of SARS CoV-2 Infection in Liver Disease Patients: A Single-Centre Retrospective Study

The limited literature on the clinical course of COVID-19 among patients with underlying liver disease (LD) is available from India. The present study aimed to evaluate the clinical and mutational profile of SARS-CoV-2 among LD cases. This was a retrospective study including admitted LD cases in whom SARS-CoV-2 RT-PCR testing was performed. Complete demographic and clinical details were retrieved from Hospital Information System. Detailed mutational analysis was performed by comparing LD COVID-19 positive study group, i.e. LD-CoV(+) with COVID-19 positive outpatients without any underlying LD as control, i.e. NLD-CoV(+). Out of 232 enrolled LD cases, 137 (59.1%) were LD-CoV(+). LD cases with existing co-morbidities were affected more (P = 0.002) and had 2.29 times (OR 2.29, CI 95%, 1.25-4.29) higher odds of succumbing to COVID-19 (P = 0.006). On multivariate regression analysis, ascites (P = 0.05), severe COVID-19 pneumonia (P = 0.046), and an increased levels of bilirubin (P = 0.005) and alkaline phosphatase (P = 0.003) were found to be associated with adverse outcome in LD-CoV(+).On mutational analysis, we found certain differences between LD- and NLD-CoV(+) infected with Delta [LD- and NLD-CoV (+ /D)] and Omicron [LD- and NLD-CoV(+/O)]. More mutations were shared between LD- and NLD-CoV(+/O) compared to LD- and NLD-CoV(+/D). There were differences in prevalence of indel mutations specific to LD-CoV ( +) for both Delta and Omicron. Moreover, we also reported an interesting genic bias between LD- and NLD-CoV( +) in harbouring deleterious/tolerated mutations. To conclude, LD cases with comorbidities were affected more and had higher odds of mortality due to COVID-19. The definite difference between LD- and NLD-CoV(+) groups with respect to frequency of harboured mutations and an inherent genic bias between them is of noteworthy importance.

Genomic ascertainment of CHEK2 -related cancer predisposition

Purpose

There is clear evidence that deleterious germline variants in CHEK2 increases risk for breast and prostate cancers; there is limited or conflicting evidence for other cancers. Genomic ascertainment was used to quantify cancer risk in CHEK2 germline pathogenic variant heterozygotes.

Patients and Methods

Germline CHEK2 variants were extracted from two exome-sequenced biobanks linked to the electronic health record: UK Biobank (n= 469,765 ) and Geisinger MyCode (n=170,503 ) . Variants were classified as per American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) criteria. Heterozygotes harbored a CHEK2 pathogenic/likely pathogenic (P/LP) variant; controls harbored benign/likely benign CHEK2 variation or wildtype CHEK2 . Tumor phenotype and demographic data were retrieved; to adjust for relatedness, association analysis was performed with SAIGE-GENE+ with Bonferroni correction.

Results

In CHEK2 heterozygotes in both MyCode and UK Biobank, there was a significant excess risk of all cancers tested, including breast cancer (C50; OR=1.54 and 1.84, respectively), male genital organ cancer (C60-C63; OR=1.61 and 1.77 respectively), urinary tract cancer (C64-C68; OR=1.56 and 1.75, respectively) and lymphoid, hematopoietic, and related tissue cancer (C81-C96; OR=1.42 and 2.11, respectively). Compared to controls, age-dependent cancer penetrance in CHEK2 heterozygotes was significantly younger in both cohorts; no significant difference was observed between the penetrance of truncating and missense variants for cancer in either cohort. Overall survival was significantly decreased in CHEK2 heterozygotes in UK Biobank but there was no statistical difference in MyCode.

Conclusion

Using genomic ascertainment in two population-scale cohorts, this investigation quantified the prevalence, penetrance, cancer phenotype and survival in CHEK2 heterozygotes. Tailored treatment options and surveillance strategies to manage those risks are warranted.

Establishment and application of a SNP molecular identification system in Grifola frondosa

Germplasm resources of edible mushrooms are essential for the breeding of varieties with improved traits. Analysis of the genetic diversity of Grifola frondosa germplasm resources and clarification of the genetic relationships among strains can provide valuable information for the selection of breeding parents. A total of 829,488 high-quality SNP loci were screened from 2,125,382 SNPs obtained by sequencing 60 G. frondose. Phylogenetic analysis, PCA, and population structure analysis based on the high-quality SNPs showed that the 60 strains could be divided into five subgroups, and the clustering results were consistent with the geographical distributions of these strains. Based on high-quality SNP loci, a core collection containing 18 representative germplasm resources was constructed, and 1,473 Kompetitive Allele-Specific PCR markers were obtained. A total of 722 SNP markers in the exonic regions were screened using KASP-genotyping experiments, and 50 candidate SNP markers and 12 core SNP markers were obtained. Genetic fingerprints of G. frondosa germplasm resources were constructed based on the selected SNP markers; these fingerprints provide an accurate, rapid, convenient, and efficient method for the identification of G. frondosa germplasm resources. The results of this study have important implications for the preservation and utilization of G. frondosa germplasm resources and the identification of varieties.

Genomic data reveal a north-south split and introgression history of blood fluke (Schistosoma haematobium) populations from across Africa

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the extent of hybridization in nature is unclear. We analyzed 34.6 million single nucleotide variants in 162 samples from 18 African countries, revealing a sharp genetic discontinuity between northern and southern S. haematobium. We found no evidence for recent hybridization. Instead the data reveal admixture events that occurred 257-879 generations ago in northern S. haematobium populations. Fifteen introgressed S. bovis genes are approaching fixation in northern S. haematobium with four genes potentially driving adaptation. We identified 19 regions that were resistant to introgression; these were enriched on the sex chromosomes. These results (i) demonstrate strong barriers to gene flow between these species, (ii) indicate that hybridization may be less common than currently envisaged, but (iii) reveal profound genomic consequences of interspecific hybridization between schistosomes of medical and veterinary importance.

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.

Identification of stripe rust resistance gene YrBDT in Chinese landrace wheat Baidatou using BSE-seq and BSR-seq

KEY MESSAGE

A new stripe rust resistance gene YrBDT in Chinese landrace wheat Baidatou was mapped to a 943.6-kb interval on chromosome arm 6DS and co-segregated with a marker CAPS3 developed from candidate gene TraesCS6D03G0027300. Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is a devastating foliar disease of wheat. Chinese landrace wheat Baidatou has shown high resistance to a broad spectrum of Pst races at both the seedling and adult-plant stages for decades in the Longnan region of Gansu province, a hot spot for stripe rust epidemics. Here, we report fine mapping and candidate gene analysis of stripe rust resistance gene YrBDT in Baidatou. Analysis of F1, F2 plants and F2:3 lines indicated that resistance in Baidatou to Pst race CYR31 was conferred by a single dominant gene, temporarily designated YrBDT. Bulked segregant exome capture sequencing (BSE-seq) analysis revealed 61 high-confidence polymorphic SNPs concentrated in a 5.4-Mb interval at the distal of chromosome arm 6DS. Several SNPs and InDels were also identified by genome mining of DNA sampled from the parents and contrasting bulks. The YrBDT locus was mapped to a 943.6-kb (4,658,322-5,601,880 bp) genomic region spanned by markers STS2 and STS3 based on IWGSC RefSeq v2.1, including five putative disease resistance genes. There was high collinearity of the target interval among Chinese Spring RefSeq v2.1, Ae. tauschii AL8/78 and Fielder genomes. The expression level of TraesCS6D03G0027300 showed significant association with Pst infection, and a gene-specific marker CAPS3 developed from TraesCS6D03G0027300 co-segregated with YrBDT suggesting this gene as a candidate of YrBDT. The resistance gene and flanking markers can be used in marker-assisted selection for improvement of stripe rust resistance.

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.

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

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.

A Novel Pathogenic TUBA1A Variant in a Croatian Infant Is Linked to a Severe Tubulinopathy with Walker-Warburg-like Features

Tubulinopathies are associated with malformations of cortical development but not Walker-Warburg Syndrome. Intensive monitoring of a Croatian infant presenting as Walker-Warburg Syndrome in utero began at 21 weeks due to increased growth of cerebral ventricles and foetal biparietal diameter. Monitoring continued until Caesarean delivery at 34 weeks where the infant was eutrophic. Clinical assessment of a progressive neurological disorder of unknown aetiology found a macrocephalic head and markedly hypoplastic genitalia with a micropenis. Neurological examination showed generalized hypotonia with very rare spontaneous movements, hypotonia-induced respiratory insufficiency and ventilator dependence, and generalized myoclonus intensifying during manipulation. With clinical features of hypotonia, lissencephaly, and brain malformations, Walker-Warburg Syndrome was suspected; however, eye anomalies were absent. Genetic trio analysis via whole-exome sequencing only identified a novel de novo mutation in the TUBA1A gene (NM_006009.4:c.848A>G; NP_006000.2:p.His283Arg) in the infant, who died at 2 months of age, as the likely cause. We report a previously unpublished, very rare heterozygous TUBA1A mutation with clinical features of macrocephaly and hypoplastic genitalia which have not previously been associated with the gene. The absence of eye phenotypes or mutations in Walker-Warburg-associated genes confirm this as not a new presentation of Walker-Warburg Syndrome but a novel TUBA1A tubulinopathy for neonatologists to be aware of.

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.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

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.

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.

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.

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.

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.

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.

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.

Genome-wide association study of stem structural characteristics that extracted by a high-throughput phenotypic analysis "LabelmeP rice" in rice

Stem is important for assimilating transport and plant strength; however, less is known about the genetic basis of its structural characteristics. In this study, a high-throughput method, "LabelmeP rice" was developed to generate 14 traits related to stem regions and vascular bundles, which allows the establishment of a stem cross-section phenotype dataset containing anatomical information of 1738 images from hand-cut transections of stems collected from 387 rice germplasm accessions grown over two successive seasons. Then, the phenotypic diversity of the rice accessions was evaluated. Genome-wide association studies identified 94, 83, and 66 significant single nucleotide polymorphisms (SNPs) for the assayed traits in 2 years and their best linear unbiased estimates, respectively. These SNPs can be integrated into 29 quantitative trait loci (QTL), and 11 of them were common in 2 years, while correlated traits shared 19. In addition, 173 candidate genes were identified, and six located at significant SNPs were repeatedly detected and annotated with a potential function in stem development. By using three introgression lines (chromosome segment substitution lines), four of the 29 QTLs were validated. LOC_Os01g70200, located on the QTL uq1.4, is detected for the area of small vascular bundles (SVB) and the rate of large vascular bundles number to SVB number. Besides, the CRISPR/Cas9 editing approach has elucidated the function of the candidate gene LOC_Os06g46340 in stem development. In conclusion, the results present a time- and cost-effective method that provides convenience for extracting rice stem anatomical traits and the candidate genes/QTL, which would help improve rice.

A machine learning enhanced EMS mutagenesis probability map for efficient identification of causal mutations in Caenorhabditis elegans

Chemical mutagenesis-driven forward genetic screens are pivotal in unveiling gene functions, yet identifying causal mutations behind phenotypes remains laborious, hindering their high-throughput application. Here, we reveal a non-uniform mutation rate caused by Ethyl Methane Sulfonate (EMS) mutagenesis in the C. elegans genome, indicating that mutation frequency is influenced by proximate sequence context and chromatin status. Leveraging these factors, we developed a machine learning enhanced pipeline to create a comprehensive EMS mutagenesis probability map for the C. elegans genome. This map operates on the principle that causative mutations are enriched in genetic screens targeting specific phenotypes among random mutations. Applying this map to Whole Genome Sequencing (WGS) data of genetic suppressors that rescue a C. elegans ciliary kinesin mutant, we successfully pinpointed causal mutations without generating recombinant inbred lines. This method can be adapted in other species, offering a scalable approach for identifying causal genes and revitalizing the effectiveness of forward genetic screens.

Implementation of long-read sequencing for routine molecular diagnosis of familial mediterranean fever

Background

Long-read sequencing technology, widely used in research, is proving useful in clinical diagnosis, especially for infectious diseases. Despite recent advances, it hasn't been routinely applied to constitutional human diseases. Long-read sequencing detects intronic variants and phases variants, crucial for identifying recessive diseases.

Methods

We integrated long-read sequencing into the clinical diagnostic workflow for the MEFV gene, responsible for familial Mediterranean fever (FMF), using a Nanopore-based workflow. This involved long-range PCR amplification, native barcoding kit library preparation, GridION sequencing, and in-house bioinformatics. We compared this new workflow against our validated method using 39 patient samples and 3 samples from an external quality assessment scheme to ensure compliance with ISO15189 standards.

Results

Our evaluation demonstrated excellent performance, meeting ISO15189 requirements for reproducibility, repeatability, sensitivity, and specificity. Since October 2022, 150 patient samples were successfully analyzed with no failures. Among these samples, we identified 13 heterozygous carriers of likely pathogenic (LP) or pathogenic (P) variants, 1 patient with a homozygous LP/P variant in MEFV, and 4 patients with compound heterozygous variants.

Conclusion

This study represents the first integration of long-read sequencing for FMF clinical diagnosis, achieving 100 % sensitivity and specificity. Our findings highlight its potential to identify pathogenic variants without parental segregation analysis, offering faster, cost-effective, and accurate clinical diagnosis. This successful implementation lays the groundwork for future applications in other constitutional human diseases, advancing precision medicine.

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.

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.

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.

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.

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.

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.

Genome-wide identification, characterization and expression analysis of the bZIP transcription factors in garlic (Allium sativum L.)

Introduction

The bZIP genes (bZIPs) are essential in numerous biological processes, including development and stress responses. Despite extensive research on bZIPs in many plants, a comprehensive genome-wide analysis of bZIPs in garlic has yet to be undertaken.

Methods

In this study, we identified and classified 64 AsbZIP genes (AsbZIPs) into 10 subfamilies. A systematic analysis of the evolutionary characteristics of these AsbZIPs, including chromosome location, gene structure, conserved motifs, and gene duplication, was conducted. Furthermore, we also examined the nucleotide diversity, cis-acting elements, and expression profiles of AsbZIPs in various tissues and under different abiotic stresses and hormone treatments.

Results and Discussion

Our findings revealed that gene replication plays a crucial role in the expansion of AsbZIPs, with a minor genetic bottleneck observed during domestication. Moreover, the identification of cis-acting elements suggested potential associations of AsbZIPs with garlic development, hormone, and stress responses. Several AsbZIPs exhibited tissue-preferential and stress/hormone-responsive expression patterns. Additionally, Asa7G01972 and Asa7G01379 were notably differentially expressed under various stresses and hormone treatments. Subsequent yeast two-hybridization and yeast induction experiments validated their interactions with Asa1G01577, a homologue of ABI5, reinforcing their importance in hormone and abiotic stress responses. This study unveiled the characteristics of the AsbZIP superfamily and lays a solid foundation for further functional analysis of AsbZIP in garlic.

Candida albicans' inorganic phosphate transport and evolutionary adaptation to phosphate scarcity

Phosphorus is essential in all cells' structural, metabolic and regulatory functions. For fungal cells that import inorganic phosphate (Pi) up a steep concentration gradient, surface Pi transporters are critical capacitators of growth. Fungi must deploy Pi transporters that enable optimal Pi uptake in pH and Pi concentration ranges prevalent in their environments. Single, triple and quadruple mutants were used to characterize the four Pi transporters we identified for the human fungal pathogen Candida albicans, which must adapt to alkaline conditions during invasion of the host bloodstream and deep organs. A high-affinity Pi transporter, Pho84, was most efficient across the widest pH range while another, Pho89, showed high-affinity characteristics only within one pH unit of neutral. Two low-affinity Pi transporters, Pho87 and Fgr2, were active only in acidic conditions. Only Pho84 among the Pi transporters was clearly required in previously identified Pi-related functions including Target of Rapamycin Complex 1 signaling, oxidative stress resistance and hyphal growth. We used in vitro evolution and whole genome sequencing as an unbiased forward genetic approach to probe adaptation to prolonged Pi scarcity of two quadruple mutant lineages lacking all 4 Pi transporters. Lineage-specific genomic changes corresponded to divergent success of the two lineages in fitness recovery during Pi limitation. Initial, large-scale genomic alterations like aneuploidies and loss of heterozygosity eventually resolved, as populations gained small-scale mutations. Severity of some phenotypes linked to Pi starvation, like cell wall stress hypersensitivity, decreased in parallel to evolving populations' fitness recovery in Pi scarcity, while severity of others like membrane stress responses diverged from Pi scarcity fitness. Among preliminary candidate genes for contributors to fitness recovery, those with links to TORC1 were overrepresented. Since Pi homeostasis differs substantially between fungi and humans, adaptive processes to Pi deprivation may harbor small-molecule targets that impact fungal growth, stress resistance and virulence.

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.