Biases in ARG-Based Inference of Historical Population Size in Populations Experiencing Selection.

Inferring the demographic history of populations provides fundamental insights into species dynamics and is essential for developing a null model to accurately study selective processes. However, background selection and selective sweeps can produce genomic signatures at linked sites that mimic or mask signals associated with historical population size change. While the theoretical biases introduced by the linked effects of selection have been well established, it is unclear whether ancestral recombination graph (ARG)-based approaches to demographic inference in typical empirical analyses are susceptible to misinference due to these effects. To address this, we developed highly realistic forward simulations of human and Drosophila melanogaster populations, including empirically estimated variability of gene density, mutation rates, recombination rates, purifying, and positive selection, across different historical demographic scenarios, to broadly assess the impact of selection on demographic inference using a genealogy-based approach. Our results indicate that the linked effects of selection minimally impact demographic inference for human populations, although it could cause misinference in populations with similar genome architecture and population parameters experiencing more frequent recurrent sweeps. We found that accurate demographic inference of D. melanogaster populations by ARG-based methods is compromised by the presence of pervasive background selection alone, leading to spurious inferences of recent population expansion, which may be further worsened by recurrent sweeps, depending on the proportion and strength of beneficial mutations. Caution and additional testing with species-specific simulations are needed when inferring population history with non-human populations using ARG-based approaches to avoid misinference due to the linked effects of selection.

Cell autonomous microglia defects in a stem cell model of frontotemporal dementia.

Neuronal dysfunction has been extensively studied as a central feature of neurodegenerative tauopathies. However, across neurodegenerative diseases, there is strong evidence for active involvement of immune cells like microglia in driving disease pathophysiology. Here, we demonstrate that tau mRNA and protein are expressed in microglia in human brains and in human induced pluripotent stem cell (iPSC)-derived microglia like cells (iMGLs). Using iMGLs harboring the MAPT IVS10+16 mutation and isogenic controls, we demonstrate that a tau mutation is sufficient to alter microglial transcriptional states. We discovered that MAPT IVS10+16 microglia exhibit cytoskeletal abnormalities, stalled phagocytosis, disrupted TREM2/TYROBP networks, and altered metabolism. Additionally, we found that secretory factors from MAPT IVS10+16 iMGLs impact neuronal health, reducing synaptic density in neurons. Key features observed in vitro were recapitulated in human brain tissue and cerebrospinal fluid from MAPT mutations carriers. Together, our findings that MAPT IVS10+16 drives cell-intrinsic dysfunction in microglia that impacts neuronal health has major implications for development of therapeutic strategies.

Long non-coding RNA SNHG8 drives stress granule formation in tauopathies.

Tauopathies are a heterogenous group of neurodegenerative disorders characterized by tau aggregation in the brain. In a subset of tauopathies, rare mutations in the MAPT gene, which encodes the tau protein, are sufficient to cause disease; however, the events downstream of MAPT mutations are poorly understood. Here, we investigate the role of long non-coding RNAs (lncRNAs), transcripts >200 nucleotides with low/no coding potential that regulate transcription and translation, and their role in tauopathy. Using stem cell derived neurons from patients carrying a MAPT p.P301L, IVS10 + 16, or p.R406W mutation and CRISPR-corrected isogenic controls, we identified transcriptomic changes that occur as a function of the MAPT mutant allele. We identified 15 lncRNAs that were commonly differentially expressed across the three MAPT mutations. The commonly differentially expressed lncRNAs interact with RNA-binding proteins that regulate stress granule formation. Among these lncRNAs, SNHG8 was significantly reduced in a mouse model of tauopathy and in FTLD-tau, progressive supranuclear palsy, and Alzheimer's disease brains. We show that SNHG8 interacts with tau and stress granule-associated RNA-binding protein TIA1. Overexpression of mutant tau in vitro is sufficient to reduce SNHG8 expression and induce stress granule formation. Rescuing SNHG8 expression leads to reduced stress granule formation and reduced TIA1 levels in immortalized cells and in MAPT mutant neurons, suggesting that dysregulation of this non-coding RNA is a causal factor driving stress granule formation via TIA1 in tauopathies.

crosshap: R package for local haplotype visualization for trait association analysis.

SUMMARY: Genome-wide association studies (GWAS) excels at harnessing dense genomic variant datasets to identify candidate regions responsible for producing a given phenotype. However, GWAS and traditional fine-mapping methods do not provide insight into the complex local landscape of linkage that contains and has been shaped by the causal variant(s). Here, we present crosshap, an R package that performs robust density-based clustering of variants based on their linkage profiles to capture haplotype structures in a local genomic region of interest. Following this, crosshap is equipped with visualization tools for choosing optimal clustering parameters (ɛ) before producing an intuitive figure that provides an overview of the complex relationships between linked variants, haplotype combinations, phenotype, and metadata traits. AVAILABILITY AND IMPLEMENTATION: The crosshap package is freely available under the MIT license and can be downloaded directly from CRAN with R >4.0.0. The development version is available on GitHub alongside issue support (https://github.com/jacobimarsh/crosshap). Tutorial vignettes and documentation are available (https://jacobimarsh.github.io/crosshap/).

Diversity of fatty acid biosynthesis genes across the soybean pangenome.

Soybean (Glycine max) is a major crop that contributes more than half of global oilseed production. Much research has been directed towards improvement of the fatty acid profile of soybean seeds through marker assisted breeding. Recently published soybean pangenomes, based on thousands of soybean lines, provide an opportunity to identify new alleles that may be involved in fatty acid biosynthesis. In this study, we identify fatty acid biosynthesis genes in soybean pangenomes based on sequence identity with known genes and examine their sequence diversity across diverse soybean collections. We find three possible instances of a gene missing in wild soybean, including FAD8 and FAD2-2D, which may be involved in oleic and linoleic acid desaturation, respectively, although we recommend follow-up research to verify the absence of these genes. More than half of the 53 fatty acid biosynthesis genes identified contained missense variants, including one linked with a previously identified QTL for oil quality. These variants were present in multiple studies based on either short read mappings or alignment of reference grade genomes. Missense variants were found in previously characterized genes including FAD2-1A and FAD2-1B, both of which are involved in desaturation of oleic acid, as well as uncharacterized candidate fatty acid biosynthesis genes. We find that the frequency of missense alleles in fatty acid biosynthesis genes has been reduced significantly more than the global average frequency of missense mutations during domestication, and missense variation in some genes is near absent in modern cultivars. This could be due to the selection for fatty acid profiles in seed, though future work should be conducted towards understanding the phenotypic impacts of these variants.

Defects in lysosomal function and lipid metabolism in human microglia harboring a TREM2 loss of function mutation.

TREM2 is an innate immune receptor expressed by microglia in the adult brain. Genetic variation in the TREM2 gene has been implicated in risk for Alzheimer's disease and frontotemporal dementia, while homozygous TREM2 mutations cause a rare leukodystrophy, Nasu-Hakola disease (NHD). Despite extensive investigation, the role of TREM2 in NHD pathogenesis remains poorly understood. Here, we investigate the mechanisms by which a homozygous stop-gain TREM2 mutation (p.Q33X) contributes to NHD. Induced pluripotent stem cell (iPSC)-derived microglia (iMGLs) were generated from two NHD families: three homozygous TREM2 p.Q33X mutation carriers (termed NHD), two heterozygous mutation carriers, one related non-carrier, and two unrelated non-carriers. Transcriptomic and biochemical analyses revealed that iMGLs from NHD patients exhibited lysosomal dysfunction, downregulation of cholesterol genes, and reduced lipid droplets compared to controls. Also, NHD iMGLs displayed defective activation and HLA antigen presentation. This defective activation and lipid droplet content were restored by enhancing lysosomal biogenesis through mTOR-dependent and independent pathways. Alteration in lysosomal gene expression, such as decreased expression of genes implicated in lysosomal acidification (ATP6AP2) and chaperone mediated autophagy (LAMP2), together with reduction in lipid droplets were also observed in post-mortem brain tissues from NHD patients, thus closely recapitulating in vivo the phenotype observed in iMGLs in vitro. Our study provides the first cellular and molecular evidence that the TREM2 p.Q33X mutation in microglia leads to defects in lysosomal function and that compounds targeting lysosomal biogenesis restore a number of NHD microglial defects. A better understanding of how microglial lipid metabolism and lysosomal machinery are altered in NHD and how these defects impact microglia activation may provide new insights into mechanisms underlying NHD and other neurodegenerative diseases.

Long non-coding RNA SNHG8 drives stress granule formation in tauopathies.

Tauopathies are a heterogenous group of neurodegenerative disorders characterized by tau aggregation in the brain. In a subset of tauopathies, rare mutations in the MAPT gene, which encodes the tau protein, are sufficient to cause disease; however, the events downstream of MAPT mutations are poorly understood. Here, we investigate the role of long non-coding RNAs (lncRNAs), transcripts >200 nucleotides with low/no coding potential that regulate transcription and translation, and their role in tauopathy. Using stem cell derived neurons from patients carrying a MAPT p.P301L, IVS10+16, or p.R406W mutation, and CRISPR-corrected isogenic controls, we identified transcriptomic changes that occur as a function of the MAPT mutant allele. We identified 15 lncRNAs that were commonly differentially expressed across the three MAPT mutations. The commonly differentially expressed lncRNAs interact with RNA-binding proteins that regulate stress granule formation. Among these lncRNAs, SNHG8 was significantly reduced in a mouse model of tauopathy and in FTLD-tau, progressive supranuclear palsy, and Alzheimer’s disease brains. We show that SNHG8 interacts with tau and stress granule-associated RNA-binding protein TIA1. Overexpression of mutant tau in vitro is sufficient to reduce SNHG8 expression and induce stress granule formation. Rescuing SNHG8 expression leads to reduced stress granule formation and reduced TIA1 levels, suggesting that dysregulation of this non-coding RNA is a causal factor driving stress granule formation via TIA1 in tauopathies.

Legume-wide comparative analysis of pod shatter locus PDH1 reveals phaseoloid specificity, high cowpea expression, and stress responsive genomic context.

Pod dehiscence is a major source of yield loss in legumes, which is exacerbated by aridity. Disruptive mutations in "Pod indehiscent 1" (PDH1), a pod sclerenchyma-specific lignin biosynthesis gene, has been linked to significant reductions in dehiscence in several legume species. We compared syntenic PDH1 regions across 12 legumes and two outgroups to uncover key historical evolutionary trends at this important locus. Our results clarified the extent to which PDH1 orthologs are present in legumes, showing the typical genomic context surrounding PDH1 has only arisen relatively recently in certain phaseoloid species (Vigna, Phaseolus, Glycine). The notable absence of PDH1 in Cajanus cajan may be a major contributor to its indehiscent phenotype compared with other phaseoloids. In addition, we identified a novel PDH1 ortholog in Vigna angularis and detected remarkable increases in PDH1 transcript abundance during Vigna unguiculata pod development. Investigation of the shared genomic context of PDH1 revealed it lies in a hotspot of transcription factors and signaling gene families that respond to abscisic acid and drought stress, which we hypothesize may be an additional factor influencing expression of PDH1 under specific environmental conditions. Our findings provide key insights into the evolutionary history of PDH1 and lay the foundation for optimizing the pod dehiscence role of PDH1 in major and understudied legume species.

Application of crop wild relatives in modern breeding: An overview of resources, experimental and computational methodologies.

Global agricultural industries are under pressure to meet the future food demand; however, the existing crop genetic diversity might not be sufficient to meet this expectation. Advances in genome sequencing technologies and availability of reference genomes for over 300 plant species reveals the hidden genetic diversity in crop wild relatives (CWRs), which could have significant impacts in crop improvement. There are many ex-situ and in-situ resources around the world holding rare and valuable wild species, of which many carry agronomically important traits and it is crucial for users to be aware of their availability. Here we aim to explore the available ex-/in- situ resources such as genebanks, botanical gardens, national parks, conservation hotspots and inventories holding CWR accessions. In addition we highlight the advances in availability and use of CWR genomic resources, such as their contribution in pangenome construction and introducing novel genes into crops. We also discuss the potential and challenges of modern breeding experimental approaches (e.g. de novo domestication, genome editing and speed breeding) used in CWRs and the use of computational (e.g. machine learning) approaches that could speed up utilization of CWR species in breeding programs towards crop adaptability and yield improvement.

Isolation and Sequencing of Chromosome Arm 7RS of Rye, Secale cereale.

Rye (Secale cereale) is a climate-resilient cereal grown extensively as grain or forage crop in Northern and Eastern Europe. In addition to being an important crop, it has been used to improve wheat through introgression of genomic regions for improved yield and disease resistance. Understanding the genomic diversity of rye will assist both the improvement of this crop and facilitate the introgression of more valuable traits into wheat. Here, we isolated and sequenced the short arm of rye chromosome 7 (7RS) from Triticale 380SD using flow cytometry and compared it to the public Lo7 rye whole genome reference assembly. We identify 2747 Lo7 genes present on the isolated chromosome arm and two clusters containing seven and sixty-five genes that are present on Triticale 380SD 7RS, but absent from Lo7 7RS. We identified 29 genes that are not assigned to chromosomal locations in the Lo7 assembly but are present on Triticale 380SD 7RS, suggesting a chromosome arm location for these genes. Our study supports the Lo7 reference assembly and provides a repertoire of genes on Triticale 7RS.

An SGSGeneloss-Based Method for Constructing a Gene Presence-Absence Table Using Mosdepth.

Presence-absence variants (PAV) are genomic regions present in some individuals of a species, but not others. PAVs have been shown to contribute to genomic diversity, especially in bacteria and plants. These structural variations have been linked to traits and can be used to track a species' evolutionary history. PAVs are usually called by aligning short read sequence data from one or more individuals to a reference genome or pangenome assembly, and then comparing coverage. Regions where reads do not align define absence in that individual, and the regions are classified as PAVs. The method below details how to align sequence reads to a reference and how to use the sequencing-coverage calculator Mosdepth to identify PAVs and construct a PAV table for use in downstream comparative genome analysis.

Genetic and Genomic Resources for Soybean Breeding Research.

Soybean (Glycine max) is a legume species of significant economic and nutritional value. The yield of soybean continues to increase with the breeding of improved varieties, and this is likely to continue with the application of advanced genetic and genomic approaches for breeding. Genome technologies continue to advance rapidly, with an increasing number of high-quality genome assemblies becoming available. With accumulating data from marker arrays and whole-genome resequencing, studying variations between individuals and populations is becoming increasingly accessible. Furthermore, the recent development of soybean pangenomes has highlighted the significant structural variation between individuals, together with knowledge of what has been selected for or lost during domestication and breeding, information that can be applied for the breeding of improved cultivars. Because of this, resources such as genome assemblies, SNP datasets, pangenomes and associated databases are becoming increasingly important for research underlying soybean crop improvement.

Haplotype mapping uncovers unexplored variation in wild and domesticated soybean at the major protein locus cqProt-003.

KEY MESSAGE: The major soy protein QTL, cqProt-003, was analysed for haplotype diversity and global distribution, and results indicate 304 bp deletion and variable tandem repeats in protein coding regions are likely causal candidates. Here, we present association and linkage analysis of 985 wild, landrace and cultivar soybean accessions in a pan genomic dataset to characterize the major high-protein/low-oil associated locus cqProt-003 located on chromosome 20. A significant trait-associated region within a 173 kb linkage block was identified, and variants in the region were characterized, identifying 34 high confidence SNPs, 4 insertions, 1 deletion and a larger 304 bp structural variant in the high-protein haplotype. Trinucleotide tandem repeats of variable length present in the second exon of gene Glyma.20G085100 are strongly correlated with the high-protein phenotype and likely represent causal variation. Structural variation has previously been found in the same gene, for which we report the global distribution of the 304 bp deletion and have identified additional nested variation present in high-protein individuals. Mapping variation at the cqProt-003 locus across demographic groups suggests that the high-protein haplotype is common in wild accessions (94.7%), rare in landraces (10.6%) and near absent in cultivated breeding pools (4.1%), suggesting its decrease in frequency primarily correlates with domestication and continued during subsequent improvement. However, the variation that has persisted in under-utilized wild and landrace populations holds high breeding potential for breeders willing to forego seed oil to maximize protein content. The results of this study include the identification of distinct haplotype structures within the high-protein population, and a broad characterization of the genomic context and linkage patterns of cqProt-003 across global populations, supporting future functional characterization and modification.

Expanding Gene-Editing Potential in Crop Improvement with Pangenomes.

Pangenomes aim to represent the complete repertoire of the genome diversity present within a species or cohort of species, capturing the genomic structural variance between individuals. This genomic information coupled with phenotypic data can be applied to identify genes and alleles involved with abiotic stress tolerance, disease resistance, and other desirable traits. The characterisation of novel structural variants from pangenomes can support genome editing approaches such as Clustered Regularly Interspaced Short Palindromic Repeats and CRISPR associated protein Cas (CRISPR-Cas), providing functional information on gene sequences and new target sites in variant-specific genes with increased efficiency. This review discusses the application of pangenomes in genome editing and crop improvement, focusing on the potential of pangenomes to accurately identify target genes for CRISPR-Cas editing of plant genomes while avoiding adverse off-target effects. We consider the limitations of applying CRISPR-Cas editing with pangenome references and potential solutions to overcome these limitations.

Linkage Disequilibrium Statistics and Block Visualization.

Linkage disequilibrium analysis enables researchers to interrogate the genome for patterns of coinheritance between genetic markers. Visualizing these patterns, and the characteristic haplotype "blocks" of linked variants can be challenging; however, advancements are being made through the development of bioinformatics software. Here, we introduce methods for producing linkage disequilibrium statistics using the widely applicable population genomics tool PLINK, before plotting linkage blocks generated in R and utilizing visualization software LDBlockShow to compare different measures of linkage and definitions of blocks.

Sequencing the USDA core soybean collection reveals gene loss during domestication and breeding.

The gene content of plants varies between individuals of the same species due to gene presence/absence variation, and selection can alter the frequency of specific genes in a population. Selection during domestication and breeding will modify the genomic landscape, though the nature of these modifications is only understood for specific genes or on a more general level (e.g., by a loss of genetic diversity). Here we have assembled and analyzed a soybean (Glycine spp.) pangenome representing more than 1,000 soybean accessions derived from the USDA Soybean Germplasm Collection, including both wild and cultivated lineages, to assess genomewide changes in gene and allele frequency during domestication and breeding. We identified 3,765 genes that are absent from the Lee reference genome assembly and assessed the presence/absence of all genes across this population. In addition to a loss of genetic diversity, we found a significant reduction in the average number of protein-coding genes per individual during domestication and subsequent breeding, though with some genes and allelic variants increasing in frequency associated with selection for agronomic traits. This analysis provides a genomic perspective of domestication and breeding in this important oilseed crop.

ELAVL4, splicing, and glutamatergic dysfunction precede neuron loss in MAPT mutation cerebral organoids.

Frontotemporal dementia (FTD) because of MAPT mutation causes pathological accumulation of tau and glutamatergic cortical neuronal death by unknown mechanisms. We used human induced pluripotent stem cell (iPSC)-derived cerebral organoids expressing tau-V337M and isogenic corrected controls to discover early alterations because of the mutation that precede neurodegeneration. At 2 months, mutant organoids show upregulated expression of MAPT, glutamatergic signaling pathways, and regulators, including the RNA-binding protein ELAVL4, and increased stress granules. Over the following 4 months, mutant organoids accumulate splicing changes, disruption of autophagy function, and build-up of tau and P-tau-S396. By 6 months, tau-V337M organoids show specific loss of glutamatergic neurons as seen in individuals with FTD. Mutant neurons are susceptible to glutamate toxicity, which can be rescued pharmacologically by the PIKFYVE kinase inhibitor apilimod. Our results demonstrate a sequence of events that precede neurodegeneration, revealing molecular pathways associated with glutamate signaling as potential targets for therapeutic intervention in FTD.

Crop breeding for a changing climate: integrating phenomics and genomics with bioinformatics.

KEY MESSAGE: Safeguarding crop yields in a changing climate requires bioinformatics advances in harnessing data from vast phenomics and genomics datasets to translate research findings into climate smart crops in the field. Climate change and an additional 3 billion mouths to feed by 2050 raise serious concerns over global food security. Crop breeding and land management strategies will need to evolve to maximize the utilization of finite resources in coming years. High-throughput phenotyping and genomics technologies are providing researchers with the information required to guide and inform the breeding of climate smart crops adapted to the environment. Bioinformatics has a fundamental role to play in integrating and exploiting this fast accumulating wealth of data, through association studies to detect genomic targets underlying key adaptive climate-resilient traits. These data provide tools for breeders to tailor crops to their environment and can be introduced using advanced selection or genome editing methods. To effectively translate research into the field, genomic and phenomic information will need to be integrated into comprehensive clade-specific databases and platforms alongside accessible tools that can be used by breeders to inform the selection of climate adaptive traits. Here we discuss the role of bioinformatics in extracting, analysing, integrating and managing genomic and phenomic data to improve climate resilience in crops, including current, emerging and potential approaches, applications and bottlenecks in the research and breeding pipeline.

TFEB regulates lysosomal exocytosis of tau and its loss of function exacerbates tau pathology and spreading.

Neurofibrillary tangles (NFTs) composed of hyperphosphorylated and misfolded tau protein are a pathological hallmark of Alzheimer's disease and other tauopathy conditions. Tau is predominantly an intraneuronal protein but is also secreted in physiological and pathological conditions. The extracellular tau has been implicated in the seeding and propagation of tau pathology and is the prime target of the current tau immunotherapy. However, truncated tau species lacking the microtubule-binding repeat (MTBR) domains essential for seeding have been shown to undergo active secretion and the mechanisms and functional consequences of the various extracellular tau are poorly understood. We report here that the transcription factor EB (TFEB), a master regulator of lysosomal biogenesis, plays an essential role in the lysosomal exocytosis of selected tau species. TFEB loss of function significantly reduced the levels of interstitial fluid (ISF) tau in PS19 mice expressing P301S mutant tau and in conditioned media of mutant tau expressing primary neurons, while the secretion of endogenous wild-type tau was not affected. Mechanistically we found that TFEB regulates the secretion of truncated mutant tau lacking MTBR and this process is dependent on the lysosomal calcium channel TRPML1. Consistent with the seeding-incompetent nature of the truncated tau and supporting the concept that TFEB-mediated lysosomal exocytosis promotes cellular clearance, we show that reduced ISF tau in the absence of TFEB is associated with enhanced intraneuronal pathology and accelerated spreading. Our results support the idea that TFEB-mediated tau exocytosis serves as a clearance mechanism to reduce intracellular tau under pathological conditions and that effective tau immunotherapy should devoid targeting these extracellular tau species.