High-throughput tagging of endogenous loci for rapid characterization of protein function.

To facilitate the interrogation of protein function at scale, we have developed high-throughput insertion of tags across the genome (HITAG). HITAG enables users to rapidly produce libraries of cells, each with a different protein of interest C-terminally tagged. HITAG is based on a modified strategy for performing Cas9-based targeted insertions, coupled with an improved approach for selecting properly tagged lines. Analysis of the resulting clones generated by HITAG reveals high tagging specificity, with most successful tagging events being indel free. Using HITAG, we fuse mCherry to a set of 167 stress granule-associated proteins and elucidate the features that drive a subset of proteins to strongly accumulate within these transient RNA-protein granules.

Floxed Il1rl2 Locus with mCherry Reporter Element Reveals Distinct Expression Patterns of the IL-36 Receptor in Barrier Tissues.

IL-36 cytokines are emerging as beneficial in immunity against pathogens and cancers but can also be detrimental when dysregulated in autoimmune and autoinflammatory conditions. Interest in targeting IL-36 activity for therapeutic purposes is rapidly growing, yet many unknowns about the functions of these cytokines remain. Thus, the availability of robust research tools is essential for both fundamental basic science and pre-clinical studies to fully access outcomes of any manipulation of the system. For this purpose, a floxed Il1rl2, the gene encoding the IL-36 receptor, mouse strain was developed to facilitate the generation of conditional knockout mice. The targeted locus was engineered to contain an inverted mCherry reporter sequence that upon Cre-mediated recombination will be flipped and expressed under the control of the endogenous Il1rl2 promoter. This feature can be used to confirm knockout in individual cells but also as a reporter to determine which cells express the IL-36 receptor IL-1RL2. The locus was confirmed to function as intended and further used to demonstrate the expression of IL-1RL2 in barrier tissues. Il1rl2 expression was detected in leukocytes in all barrier tissues. Interestingly, strong expression was observed in epithelial cells at locations in direct contact with the environment such as the skin, oral mucosa, the esophagus, and the upper airways, but almost absent from epithelial cells at more inward facing sites, including lung alveoli, the small intestine, and the colon. These findings suggest specialized functions of IL-1RL2 in outward facing epithelial tissues and cells. The generated mouse model should prove valuable in defining such functions and may also facilitate basic and translational research.

Genome-wide association analysis explores the genetic loci of amino acid content in duck's breast muscle.

BACKGROUND: Amino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat. RESULT: In the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55-56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2-21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways. CONCLUSION: In this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

Association of Novel Loci With Keratoconus Susceptibility in a Chinese Genome-Wide Association Study.

Purpose: Keratoconus (KC) is a progressive corneal disease that can lead to corneal blindness if not properly managed. The purpose of this study was to identify genetic associations with KC in China and to investigate whether these genetic variants are associated with corneal thickness and corneal curvature in KC cases. Methods: A genome-wide association study was conducted on 853 patients with KC and 6248 controls. The KC cases were genotyped with the Illumina Infinium Human Asian Screening Array BeadChip, and the controls were genotyped with the Illumina Infinium Human Global Screening Array BeadChip. Genetic associations with KC, as well as correlations between the positive variants and corneal parameters including central corneal thickness (CCT) and mean keratometry (Km), were compared using PLINK version 1.90. Results: Our present study identified four single-nucleotide polymorphisms (SNPs) within four risk loci (PTGER3: rs2300163, EYA1: rs1077435, ASS1: rs141365191, and CHTF8: rs3743680) associated with KC in Chinese patients that reached genome-wide significance. Among the identified SNPs with P < 1.00 × 10-4, seven SNPs (FOSL2-PLB1: rs12622211, RXRA-COL5A1: rs3118515, rs3132306, rs1536482, rs3118520, KAT6B: rs192187772, RAP2A-IPO5: rs41361245) were observed to be associated with CCT, and one SNP (USP13: rs6767552) was found to be associated with Km. Conclusions: In the first genome-wide association study of KC with a relatively large study population in China, we identified four SNPs in four risk loci associated with the disease. The findings enriched the understanding of genetic susceptibility to KC and provided new insights into the genetic etiology of the disease.

Risk loci involved in giant cell arteritis susceptibility: a genome-wide association study.

BACKGROUND: Giant cell arteritis is an age-related vasculitis that mainly affects the aorta and its branches in individuals aged 50 years and older. Current options for diagnosis and treatment are scarce, highlighting the need to better understand its underlying pathogenesis. Genome-wide association studies (GWAS) have emerged as a powerful tool for unravelling the pathogenic mechanisms involved in complex diseases. We aimed to characterise the genetic basis of giant cell arteritis by performing the largest GWAS of this vasculitis to date and to assess the functional consequences and clinical implications of identified risk loci. METHODS: We collected and meta-analysed genomic data from patients with giant cell arteritis and healthy controls of European ancestry from ten cohorts across Europe and North America. Eligible patients required confirmation of giant cell arteritis diagnosis by positive temporal artery biopsy, positive temporal artery doppler ultrasonography, or imaging techniques confirming large-vessel vasculitis. We assessed the functional consequences of loci associated with giant cell arteritis using cell enrichment analysis, fine-mapping, and causal gene prioritisation. We also performed a drug repurposing analysis and developed a polygenic risk score to explore the clinical implications of our findings. FINDINGS: We included a total of 3498 patients with giant cell arteritis and 15 550 controls. We identified three novel loci associated with risk of giant cell arteritis. Two loci, MFGE8 (rs8029053; p=4·96 × 10-8; OR 1·19 [95% CI 1·12-1·26]) and VTN (rs704; p=2·75 × 10-9; OR 0·84 [0·79-0·89]), were related to angiogenesis pathways and the third locus, CCDC25 (rs11782624; p=1·28 × 10-8; OR 1·18 [1·12-1·25]), was related to neutrophil extracellular traps (NETs). We also found an association between this vasculitis and HLA region and PLG. Variants associated with giant cell arteritis seemed to fulfil a specific regulatory role in crucial immune cell types. Furthermore, we identified several drugs that could represent promising candidates for treatment of this disease. The polygenic risk score model was able to identify individuals at increased risk of developing giant cell arteritis (90th percentile OR 2·87 [95% CI 2·15-3·82]; p=1·73 × 10-13). INTERPRETATION: We have found several additional loci associated with giant cell arteritis, highlighting the crucial role of angiogenesis in disease susceptibility. Our study represents a step forward in the translation of genomic findings to clinical practice in giant cell arteritis, proposing new treatments and a method to measure genetic predisposition to this vasculitis. FUNDING: Institute of Health Carlos III, Spanish Ministry of Science and Innovation, UK Medical Research Council, and National Institute for Health and Care Research.

Loci cg06256735 and cg15815843 in the MFAP5 gene regulatory regions are hypomethylated in varicose veins apparently due to active demethylation.

Varicose vein disease (VVD) is a common health problem worldwide. Microfibril-associated protein 5 (MFAP5) is one of the potential key players in its pathogenesis. Our previous microarray analysis revealed the cg06256735 and cg15815843 loci in the regulatory regions of the MFAP5 gene as hypomethylated in varicose veins which correlated with its up-regulation. The aim of this work was to validate preliminary microarray data, estimate the level of 5-hydroxymethylcytosine (5hmC) at these loci, and determine the methylation status of one of them in different layers of the venous wall. For this, methyl- and hydroxymethyl-sensitive restriction techniques were used followed by real-time PCR and droplet digital PCR, correspondingly, as well as bisulfite pyrosequencing of +/- oxidized DNA. Our microarray data on hypomethylation at the cg06256735 and cg15815843 loci in whole varicose vein segments were confirmed and it was also demonstrated that the level of 5hmC at these loci is increased in VVD. Specifically, among other layers of the venous wall, tunica (t.) intima is the main contributor to hypomethylation at the cg06256735 locus in varicose veins. Thus, it was shown that hypomethylation at the cg06256735 and cg15815843 loci takes place in VVD, with evidence to suggest that it happens through their active demethylation leading to up-regulation of the MFAP5 gene, and t. intima is most involved in this biochemical process.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci.

BACKGROUND: Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space. RESULTS: To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4. CONCLUSIONS: Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene.

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation - a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

Two genome-wide interaction loci modify the association of nonsteroidal anti-inflammatory drugs with colorectal cancer.

Regular, long-term aspirin use may act synergistically with genetic variants, particularly those in mechanistically relevant pathways, to confer a protective effect on colorectal cancer (CRC) risk. We leveraged pooled data from 52 clinical trial, cohort, and case-control studies that included 30,806 CRC cases and 41,861 controls of European ancestry to conduct a genome-wide interaction scan between regular aspirin/nonsteroidal anti-inflammatory drug (NSAID) use and imputed genetic variants. After adjusting for multiple comparisons, we identified statistically significant interactions between regular aspirin/NSAID use and variants in 6q24.1 (top hit rs72833769), which has evidence of influencing expression of TBC1D7 (a subunit of the TSC1-TSC2 complex, a key regulator of MTOR activity), and variants in 5p13.1 (top hit rs350047), which is associated with expression of PTGER4 (codes a cell surface receptor directly involved in the mode of action of aspirin). Genetic variants with functional impact may modulate the chemopreventive effect of regular aspirin use, and our study identifies putative previously unidentified targets for additional mechanistic interrogation.

Genome-wide association analyses of breast cancer in women of African ancestry identify new susceptibility loci and improve risk prediction.

We performed genome-wide association studies of breast cancer including 18,034 cases and 22,104 controls of African ancestry. Genetic variants at 12 loci were associated with breast cancer risk (P < 5 × 10-8), including associations of a low-frequency missense variant rs61751053 in ARHGEF38 with overall breast cancer (odds ratio (OR) = 1.48) and a common variant rs76664032 at chromosome 2q14.2 with triple-negative breast cancer (TNBC) (OR = 1.30). Approximately 15.4% of cases with TNBC carried six risk alleles in three genome-wide association study-identified TNBC risk variants, with an OR of 4.21 (95% confidence interval = 2.66-7.03) compared with those carrying fewer than two risk alleles. A polygenic risk score (PRS) showed an area under the receiver operating characteristic curve of 0.60 for the prediction of breast cancer risk, which outperformed PRS derived using data from females of European ancestry. Our study markedly increases the population diversity in genetic studies for breast cancer and demonstrates the utility of PRS for risk prediction in females of African ancestry.

Deep learning of left atrial structure and function provides link to atrial fibrillation risk.

Increased left atrial volume and decreased left atrial function have long been associated with atrial fibrillation. The availability of large-scale cardiac magnetic resonance imaging data paired with genetic data provides a unique opportunity to assess the genetic contributions to left atrial structure and function, and understand their relationship with risk for atrial fibrillation. Here, we use deep learning and surface reconstruction models to measure left atrial minimum volume, maximum volume, stroke volume, and emptying fraction in 40,558 UK Biobank participants. In a genome-wide association study of 35,049 participants without pre-existing cardiovascular disease, we identify 20 common genetic loci associated with left atrial structure and function. We find that polygenic contributions to increased left atrial volume are associated with atrial fibrillation and its downstream consequences, including stroke. Through Mendelian randomization, we find evidence supporting a causal role for left atrial enlargement and dysfunction on atrial fibrillation risk.

Genome-wide association study reveals a locus in ADARB2 for complete freedom from headache in Danish Blood Donors.

Headache disorders are the most common disorders of the nervous system. The lifetime prevalence of headache disorders show that some individuals never experience headache. The etiology of complete freedom from headache is not known. To assess genetic variants associated with complete freedom from headache, we performed a genome-wide association study of individuals who have never experienced a headache. We included 63,992 individuals (2,998 individuals with complete freedom from headache and 60,994 controls) from the Danish Blood Donor Study Genomic Cohort. Participants were included in two rounds, from 2015 to 2018 and in 2020. We discovered a genome-wide significant association, with the lead variant rs7904615[G] in ADARB2 (EAF = 27%, OR = 1.20 [1.13-1.27], p = 3.92 × 10-9). The genomic locus was replicated in a non-overlapping cohort of 13,032 individuals (539 individuals with complete freedom from headache and 12,493 controls) from the Danish Blood Donor Study Genomic Cohort (p < 0.05, two-sided). Participants for the replication were included from 2015 to 2020. In conclusion, we show that complete freedom from headache has a genetic component, and we suggest that ADARB2 is involved in complete freedom from headache. The genomic locus was specific for complete freedom from headache and was not associated with any primary headache disorders.

Epidemiologic association and shared genetic architecture between cataract and hearing difficulties among middle-aged and older adults.

Age-related cataract and hearing difficulties are major sensory disorders that often co-exist in the global-wide elderly and have a tangible influence on the quality of life. However, the epidemiologic association between cataract and hearing difficulties remains unexplored, while little is known about whether the two share their genetic etiology. We first investigated the clinical association between cataract and hearing difficulties using the UK Biobank covering 502,543 individuals. Both unmatched analysis (adjusted for confounders) and a matched analysis (one control matched for each patient with cataract according to confounding factors) were undertaken and confirmed that cataract was associated with hearing difficulties (OR, 2.12; 95% CI, 1.98-2.27; OR, 2.03; 95% CI, 1.86-2.23, respectively). Furthermore, we explored and quantified the shared genetic architecture of these two complex sensory disorders at the common variant level using the bivariate causal mixture model (MiXeR) and conditional/conjunctional false discovery rate method based on the largest available genome-wide association studies of cataract (N = 585,243) and hearing difficulties (N = 323,978). Despite detecting only a negligible genetic correlation, we observe polygenic overlap between cataract and hearing difficulties and identify 6 shared loci with mixed directions of effects. Follow-up analysis of the shared loci implicates candidate genes QKI, STK17A, TYR, NSF, and TCF4 likely contribute to the pathophysiology of cataracts and hearing difficulties. In conclusion, this study demonstrates the presence of epidemiologic association between cataract and hearing difficulties and provides new insights into the shared genetic architecture of these two disorders at the common variant level.

Genetic insights into Cyphocharax magdalenae (Characiformes: Curimatidae): Microsatellite loci development and population analysis in the Cauca River, Colombia.

Cyphocharax magdalenae, a Colombian freshwater fish species, plays a vital role in nutrients distribution and serves as a significant food source for other fish species and local fishing communities. Considered a short-distance migratory species, C. magdalenae populations face substantial extinction risk due to human activities impacting their habitats. To address the lack of knowledge on genetic diversity and population structure, this study used next-generation sequencing technology to develop species-specific microsatellite loci and conducted a population genetics analysis of C. magdalenae in the middle and lower sections of the Cauca River, Colombia. Out of 30 pairs of microsatellite primers evaluated in 324 individuals, 14 loci were found to be polymorphic, at linkage equilibrium and, in at least one population, their genotypic frequencies were in Hardy-Weinberg equilibrium. Results showed high genetic diversity levels compared to other neotropical Characiformes, with inbreeding coefficients similar to those reported for phylogenetically related species. Moreover, C. magdalenae exhibits seasonal population structure (rainy-dry) consisting of two genetic stocks showing bottleneck signals and high effective population sizes. This information is essential for understanding the current species genetics and developing future management programs for this fishery resource.

Genome composition and GC content influence loci distribution in reduced representation genomic studies.

BACKGROUND: Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored. RESULTS: Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected. CONCLUSIONS: Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

GWAS determined genetic loci associated with callus induction in oil palm tissue culture.

KEY MESSAGE: GWAS identified six loci at 25 kb downstream of WAK2, a crucial gene for cell wall and callus formation, enabling development of a SNP marker for enhanced callus induction potential. Efficient callus induction is vital for successful oil palm tissue culture, yet identifying genomic loci and markers for early detection of genotypes with high potential of callus induction remains unclear. In this study, immature male inflorescences from 198 oil palm accessions (dura, tenera and pisifera) were used as explants for tissue culture. Callus induction rates were collected at one-, two- and three-months after inoculation (C1, C2 and C3) as phenotypes. Resequencing generated 11,475,258 high quality single nucleotide polymorphisms (SNPs) as genotypes. GWAS was then performed, and correlation analysis revealed a positive association of C1 with both C2 (R = 0.81) and C3 (R = 0.50), indicating that C1 could be used as the major phenotype for callus induction rate. Therefore, only significant SNPs (P ≤ 0.05) in C1 were identified to develop markers for screening individuals with high potential of callus induction. Among 21 significant SNPs in C1, LD block analysis revealed six SNPs on chromosome 12 (Chr12) potentially linked to callus formation. Subsequently, 13 SNP markers were identified from these loci and electrophoresis results showed that marker C-12 at locus Chr12_12704856 can be used effectively to distinguish the GG allele, which showed the highest probability (69%) of callus induction. Furthermore, a rapid SNP variant detection method without electrophoresis was established via qPCR-based melting curve analysis. Our findings facilitated marker-assisted selection for specific palms with high potential of callus induction using immature male inflorescence as explant, aiding ortet palm selection in oil palm tissue culture.

GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis.

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Unsupervised deep representation learning enables phenotype discovery for genetic association studies of brain imaging.

Understanding the genetic architecture of brain structure is challenging, partly due to difficulties in designing robust, non-biased descriptors of brain morphology. Until recently, brain measures for genome-wide association studies (GWAS) consisted of traditionally expert-defined or software-derived image-derived phenotypes (IDPs) that are often based on theoretical preconceptions or computed from limited amounts of data. Here, we present an approach to derive brain imaging phenotypes using unsupervised deep representation learning. We train a 3-D convolutional autoencoder model with reconstruction loss on 6130 UK Biobank (UKBB) participants' T1 or T2-FLAIR (T2) brain MRIs to create a 128-dimensional representation known as Unsupervised Deep learning derived Imaging Phenotypes (UDIPs). GWAS of these UDIPs in held-out UKBB subjects (n = 22,880 discovery and n = 12,359/11,265 replication cohorts for T1/T2) identified 9457 significant SNPs organized into 97 independent genetic loci of which 60 loci were replicated. Twenty-six loci were not reported in earlier T1 and T2 IDP-based UK Biobank GWAS. We developed a perturbation-based decoder interpretation approach to show that these loci are associated with UDIPs mapped to multiple relevant brain regions. Our results established unsupervised deep learning can derive robust, unbiased, heritable, and interpretable brain imaging phenotypes.

Gap-free X and Y chromosome assemblies of Salix arbutifolia reveal an evolutionary change from male to female heterogamety in willows, without a change in the position of the sex-determining locus.

In the Vetrix clade of Salix, a genus of woody flowering plants, sex determination involves chromosome 15, but an XY system has changed to a ZW system. We studied the detailed genetic changes involved. We used genome sequencing, with chromosome conformation capture (Hi-C) and PacBio HiFi reads to assemble chromosome level gap-free X and Y of Salix arbutifolia, and distinguished the haplotypes in the 15X- and 15Y-linked regions, to study the evolutionary history of the sex-linked regions (SLRs). Our sequencing revealed heteromorphism of the X and Y haplotypes of the SLR, with the X-linked region being considerably larger than the corresponding Y region, mainly due to accumulated repetitive sequences and gene duplications. The phylogenies of single-copy orthogroups within the SLRs indicate that S. arbutifolia and Salix purpurea share an ancestral SLR within a repeat-rich region near the chromosome 15 centromere. During the change in heterogamety, the X-linked region changed to a W-linked one, while the Z was derived from the Y.