Lycorine induces apoptosis of acute myeloid leukemia cells and inhibits triglyceride production via binding and targeting FABP5.

Acute myeloid leukemia (AML) is the most common hematopoietic malignancy with abnormal lipid metabolism. However, currently available information on the involvement of the alterations in lipid metabolism in AML development is limited. In this study, we demonstrate that FABP5 expression facilitates AML cell viability, protects AML cells from apoptosis, and maintains triglyceride production. Our bioinformatics analysis revealed that FABP5 expression was upregulated and correlated with unfavorable overall survival of AML patients. FABP5 expression may be used to distinguish normal and AML with high accuracy. FABP5-based risk score was an independent risk factor for AML patients. AML patients with highly expressed FABP5 predicted resistance to drugs. In vitro study showed that FABP5 expression was remarkably elevated in primary AML blasts and an AML cell line. Silencing FABP5 expression attenuated AML cell viability, reduced triglyceride production and lipid droplet accumulation, and induced apoptosis. We utilized AutoDock online tool to identify lycorine as an FABP5 inhibitor by binding FABP5 at amino acid residues Ile54, Thr56, Thr63, and Arg109. Lycorine treatment downregulated the expression levels of FABP5 and its target PPARγ, impaired AML cell viability, triggered apoptosis, and reduced triglyceride production in AML cells. These results demonstrate that FABP5 is critical for AML cell survival and highlight a novel metabolic vulnerability for AML.

Area postrema syndrome with linear enhancement along the surface of the brainstem and fourth ventricle in autoimmune GFAP astrocytopathy.

BACKGROUND: Glial fibrillary acidic protein (GFAP) astrocytopathy, a novel autoimmune disease of the nervous system, was first defined in 2016. To our knowledge, area postrema syndrome (APS) with linear enhancement along the surface of the brainstem and fourth ventricle is extremely rare in this disorder. CASE PRESENTATION: A Chinese woman presented with intractable nausea and vomiting after onset of flu-like symptoms. Brain magnetic resonance imaging (MRI) disclosed abnormal signal intensities in the dorsal medulla oblongata including area postrema. Besides, linear enhancement surrounding the surface of the brainstem and fourth ventricle was visualized after gadolinium injection. Cerebrospinal fluid (CSF) analysis showed increased cell count and protein. A cell-based assay was positive for anti-GFAP IgG in CSF. She was diagnosed with autoimmune GFAP astrocytopathy and treated with high-dose glucocorticoid. The patient received a quick recovery with entire resolution of the initial abnormalities. CONCLUSIONS: Isolated APS can be the initial manifestation of autoimmune GFAP astrocytopathy. Linear enhancement surrounding the surface of the brainstem and fourth ventricle is another neuroradiological hallmark.

The stroke mechanism, clinical presentation, and radiological feature of bilateral medial medullary infarction.

BACKGROUND: Bilateral medial medullary infarction (BMMI) is a rare type of posterior circulation stroke. The aim of this study is to characterize its stroke mechanisms, clinical manifestations, neuroradiological features, and prognosis. METHODS: From January 2015 to June 2021, a retrospective review of 15 patients diagnosed with BMMI was conducted. The clinical and neuroradiological features were summarized by our experienced neurologists. RESULTS: Fifteen patients (12 male, 3 female), ranging in age from 48 to 72 years, satisfied the inclusion criteria. The common clinical presentations included motor weakness (100%), deep sensory disturbance (93.3%), vertigo/dizziness (80%), dysarthria (93.3%), and dysphagia (66.7%). Vertically, infarct lesions in the rostral medulla were observed in all included patients. Horizontally, "heart appearance," "Y appearance," and "fan appearance" infarcts occurred in 9 cases (60%), 5 cases (33.3%), and 1 (6.7%) case, respectively. Patients (53.3%) had severe stenosis or occlusion in unilateral vertebral artery (VA), and 33.3% had normal findings in the vertebrobasilar artery. Patients (93.3%) achieved poor prognosis. CONCLUSION: BMMI is more frequently located in the rostral medulla and comprises three forms of infarction. The two main stroke etiologies of BMMI are large-artery atherosclerosis (LAA) and small vessel disease (SVD). BMMI is always associated with bad clinical outcome.

Cytarabine-induced TNFα promotes the expansion and suppressive functions of myeloid-derived suppressor cells in acute myeloid leukaemia.

Acute myeloid leukaemia (AML) is very common haematopoietic malignancies with poor prognosis. Chemotherapy is still a mainstay therapy for AML patients. AML microenvironment plays critical roles in therapy response. However, the role of chemotherapy in AML microenvironment is poorly understood. In this study, we report that cytarabine (AraC)-triggered TNFα from AML cells expanded myeloid-derived suppressor cells (MDSCs) and enhanced MDSC functions and survival through activating IL-6/STAT3 and NFκB pathways. Blockade of TNFα in conditioned medium-derived AraC-treated AML cells (AraC_CM) impaired MDSC expansion and functions, reduced IL-6 secretion and the level of activated STAT3. Inhibiting IL6 or STAT3 abrogated AraC_CM-mediated MDSC suppressive function. Additionally, inhibiting TNFα also impaired AraC_CM-mediated NFκB activation. Blocking NFκB activation reduced MDSC viability induced by AraC_CM. Together, these results provided a role of AraC-induced TNFα in MDSC expansion and functions and suggest that targeting TNFα may benefit AML patients to current anticancer strategies by blocking MDSC-mediated immunosuppression.

Two divergent haplotypes from a highly heterozygous lychee genome suggest independent domestication events for early and late-maturing cultivars.

Lychee is an exotic tropical fruit with a distinct flavor. The genome of cultivar 'Feizixiao' was assembled into 15 pseudochromosomes, totaling ~470 Mb. High heterozygosity (2.27%) resulted in two complete haplotypic assemblies. A total of 13,517 allelic genes (42.4%) were differentially expressed in diverse tissues. Analyses of 72 resequenced lychee accessions revealed two independent domestication events. The extremely early maturing cultivars preferentially aligned to one haplotype were domesticated from a wild population in Yunnan, whereas the late-maturing cultivars that mapped mostly to the second haplotype were domesticated independently from a wild population in Hainan. Early maturing cultivars were probably developed in Guangdong via hybridization between extremely early maturing cultivar and late-maturing cultivar individuals. Variable deletions of a 3.7 kb region encompassed by a pair of CONSTANS-like genes probably regulate fruit maturation differences among lychee cultivars. These genomic resources provide insights into the natural history of lychee domestication and will accelerate the improvement of lychee and related crops.

Robust Benchmark Structural Variant Calls of An Asian Using State-of-the-art Long-read Sequencing Technologies.

The importance of structural variants (SVs) for human phenotypes and diseases is now recognized. Although a variety of SV detection platforms and strategies that vary in sensitivity and specificity have been developed, few benchmarking procedures are available to confidently assess their performances in biological and clinical research. To facilitate the validation and application of these SV detection approaches, we established an Asian reference material by characterizing the genome of an Epstein-Barr virus (EBV)-immortalized B lymphocyte line along with identified benchmark regions and high-confidence SV calls. We established a high-confidence SV callset with 8938 SVs by integrating four alignment-based SV callers, including 109× Pacific Biosciences (PacBio) continuous long reads (CLRs), 22× PacBio circular consensus sequencing (CCS) reads, 104× Oxford Nanopore Technologies (ONT) long reads, and 114× Bionano optical mapping platform, and one de novo assembly-based SV caller using CCS reads. A total of 544 randomly selected SVs were validated by PCR amplification and Sanger sequencing, demonstrating the robustness of our SV calls. Combining trio-binning-based haplotype assemblies, we established an SV benchmark for identifying false negatives and false positives by constructing the continuous high-confidence regions (CHCRs), which covered 1.46 gigabase pairs (Gb) and 6882 SVs supported by at least one diploid haplotype assembly. Establishing high-confidence SV calls for a benchmark sample that has been characterized by multiple technologies provides a valuable resource for investigating SVs in human biology, disease, and clinical research.

Comprehensive preimplantation genetic testing by massively parallel sequencing.

STUDY QUESTION: Can whole genome sequencing (WGS) offer a relatively cost-effective approach for embryonic genome-wide haplotyping and preimplantation genetic testing (PGT) for monogenic disorders (PGT-M), aneuploidy (PGT-A) and structural rearrangements (PGT-SR)? SUMMARY ANSWER: Reliable genome-wide haplotyping, PGT-M, PGT-A and PGT-SR could be performed by WGS with 10× depth of parental and 4× depth of embryonic sequencing data. WHAT IS KNOWN ALREADY: Reduced representation genome sequencing with a genome-wide next-generation sequencing haplarithmisis-based solution has been verified as a generic approach for automated haplotyping and comprehensive PGT. Several low-depth massively parallel sequencing (MPS)-based methods for haplotyping and comprehensive PGT have been developed. However, an additional family member, such as a sibling, or a proband, is required for PGT-M haplotyping using low-depth MPS methods. STUDY DESIGN, SIZE, DURATION: In this study, 10 families that had undergone traditional IVF-PGT and 53 embryos, including 13 embryos from two PGT-SR families and 40 embryos from eight PGT-M families, were included to evaluate a WGS-based method. There were 24 blastomeres and 29 blastocysts in total. All embryos were used for PGT-A. Karyomapping validated the WGS results. Clinical outcomes of the 10 families were evaluated. PARTICIPANTS/MATERIALS, SETTING, METHODS: A blastomere or a few trophectoderm cells from the blastocyst were biopsied, and multiple displacement amplification (MDA) was performed. MDA DNA and bulk DNA of family members were used for library construction. Libraries were sequenced, and data analysis, including haplotype inheritance deduction for PGT-M and PGT-SR and read-count analysis for PGT-A, was performed using an in-house pipeline. Haplotyping with a proband and parent-only haplotyping without additional family members were performed to assess the WGS methodology. Concordance analysis between the WGS results and traditional PGT methods was performed. MAIN RESULTS AND THE ROLE OF CHANCE: For the 40 PGT-M and 53 PGT-A embryos, 100% concordance between the WGS and single-nucleotide polymorphism (SNP)-array results was observed, regardless of whether additional family members or a proband was included for PGT-M haplotyping. For the 13 embryos from the two PGT-SR families, the embryonic balanced translocation was detected and 100% concordance between WGS and MicroSeq with PCR-seq was demonstrated. LIMITATIONS, REASONS FOR CAUTION: The number of samples in this study was limited. In some cases, the reference embryo for PGT-M or PGT-SR parent-only haplotyping was not available owing to failed direct genotyping. WIDER IMPLICATIONS OF THE FINDINGS: WGS-based PGT-A, PGT-M and PGT-SR offered a comprehensive PGT approach for haplotyping without the requirement for additional family members. It provided an improved complementary method to PGT methodologies, such as low-depth MPS- and SNP array-based methods. STUDY FUNDING/COMPETING INTEREST(S): This research was supported by the research grant from the National Key R&D Program of China (2018YFC0910201 and 2018YFC1004900), the Guangdong province science and technology project of China (2019B020226001), the Shenzhen Birth Defect Screening Project Lab (JZF No. [2016] 750) and the Shenzhen Municipal Government of China (JCYJ20170412152854656). This work was also supported by the National Natural Science Foundation of China (81771638, 81901495 and 81971344), the National Key R&D Program of China (2018YFC1004901 and 2016YFC0905103), the Shanghai Sailing Program (18YF1424800), the Shanghai Municipal Commission of Science and Technology Program (15411964000) and the Shanghai 'Rising Stars of Medical Talent' Youth Development Program Clinical Laboratory Practitioners Program (201972). The authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Performance of copy number variants detection based on whole-genome sequencing by DNBSEQ platforms.

BACKGROUND: DNBSEQ™ platforms are new massively parallel sequencing (MPS) platforms that use DNA nanoball technology. Use of data generated from DNBSEQ™ platforms to detect single nucleotide variants (SNVs) and small insertions and deletions (indels) has proven to be quite effective, while the feasibility of copy number variants (CNVs) detection is unclear. RESULTS: Here, we first benchmarked different CNV detection tools based on Illumina whole-genome sequencing (WGS) data of NA12878 and then assessed these tools in CNV detection based on DNBSEQ™ sequencing data from the same sample. When the same tool was used, the CNVs detected based on DNBSEQ™ and Illumina data were similar in quantity, length and distribution, while great differences existed within results from different tools and even based on data from a single platform. We further estimated the CNV detection power based on available CNV benchmarks of NA12878 and found similar precision and sensitivity between the DNBSEQ™ and Illumina platforms. We also found higher precision of CNVs shorter than 1 kbp based on DNBSEQ™ platforms than those based on Illumina platforms by using Pindel, DELLY and LUMPY. We carefully compared these two available benchmarks and found a large proportion of specific CNVs between them. Thus, we constructed a more complete CNV benchmark of NA12878 containing 3512 CNV regions. CONCLUSIONS: We assessed and benchmarked CNV detections based on WGS with DNBSEQ™ platforms and provide guidelines for future studies.

An integrated Asian human SNV and indel benchmark established using multiple sequencing methods.

Sequencing technologies have been rapidly developed recently, leading to the breakthrough of sequencing-based clinical diagnosis, but accurate and complete genome variation benchmark would be required for further assessment of precision medicine applications. Despite the human cell line of NA12878 has been successfully developed to be a variation benchmark, population-specific variation benchmark is still lacking. Here, we established an Asian human variation benchmark by constructing and sequencing a stabilized cell line of a Chinese Han volunteer. By using seven different sequencing strategies, we obtained ~3.88 Tb clean data from different laboratories, hoping to reach the point of high sequencing depth and accurate variation detection. Through the combination of variations identified from different sequencing strategies and different analysis pipelines, we identified 3.35 million SNVs and 348.65 thousand indels, which were well supported by our sequencing data and passed our strict quality control, thus should be high confidence variation benchmark. Besides, we also detected 5,913 high-quality SNVs which had 969 sites were novel and  located in the high homologous regions supported by long-range information in both the co-barcoding single tube Long Fragment Read (stLFR) data and PacBio HiFi CCS data. Furthermore, by using the long reads data (stLFR and HiFi CCS), we were able to phase more than 99% heterozygous SNVs, which helps to improve the benchmark to be haplotype level. Our study provided comprehensive sequencing data as well as the integrated variation benchmark of an Asian derived cell line, which would be valuable for future sequencing-based clinical development.

Impact of sequencing depth and technology on de novo RNA-Seq assembly.

BACKGROUND: RNA-Seq data is inherently nonuniform for different transcripts because of differences in gene expression. This makes it challenging to decide how much data should be generated from each sample. How much should one spend to recover the less expressed transcripts? The sequencing technology used is another consideration, as there are inevitably always biases against certain sequences. To investigate these effects, we first looked at high-depth libraries from a set of well-annotated organisms to ascertain the impact of sequencing depth on de novo assembly. We then looked at libraries sequenced from the Universal Human Reference RNA (UHRR) to compare the performance of Illumina HiSeq and MGI DNBseq™ technologies. RESULTS: On the issue of sequencing depth, the amount of exomic sequence assembled plateaued using data sets of approximately 2 to 8 Gbp. However, the amount of genomic sequence assembled did not plateau for many of the analyzed organisms. Most of the unannotated genomic sequences are single-exon transcripts whose biological significance will be questionable for some users. On the issue of sequencing technology, both of the analyzed platforms recovered a similar number of full-length transcripts. The missing "gap" regions in the HiSeq assemblies were often attributed to higher GC contents, but this may be an artefact of library preparation and not of sequencing technology. CONCLUSIONS: Increasing sequencing depth beyond modest data sets of less than 10 Gbp recovers a plethora of single-exon transcripts undocumented in genome annotations. DNBseq™ is a viable alternative to HiSeq for de novo RNA-Seq assembly.

Reliable multiplex sequencing with rare index mis-assignment on DNB-based NGS platform.

BACKGROUND: Massively-parallel-sequencing, coupled with sample multiplexing, has made genetic tests broadly affordable. However, intractable index mis-assignments (commonly exceeds 1%) were repeatedly reported on some widely used sequencing platforms. RESULTS: Here, we investigated this quality issue on BGI sequencers using three library preparation methods: whole genome sequencing (WGS) with PCR, PCR-free WGS, and two-step targeted PCR. BGI's sequencers utilize a unique DNA nanoball (DNB) technology which uses rolling circle replication for DNA-nanoball preparation; this linear amplification is PCR free and can avoid error accumulation. We demonstrated that single index mis-assignment from free indexed oligos occurs at a rate of one in 36 million reads, suggesting virtually no index hopping during DNB creation and arraying. Furthermore, the DNB-based NGS libraries have achieved an unprecedentedly low sample-to-sample mis-assignment rate of 0.0001 to 0.0004% under recommended procedures. CONCLUSIONS: Single indexing with DNB technology provides a simple but effective method for sensitive genetic assays with large sample numbers.

Germline and somatic variant identification using BGISEQ-500 and HiSeq X Ten whole genome sequencing.

Technological innovation and increased affordability have contributed to the widespread adoption of genome sequencing technologies in biomedical research. In particular large cancer research consortia have embraced next generation sequencing, and have used the technology to define the somatic mutation landscape of multiple cancer types. These studies have primarily utilised the Illumina HiSeq platforms. In this study we performed whole genome sequencing of three malignant pleural mesothelioma and matched normal samples using a new platform, the BGISEQ-500, and compared the results obtained with Illumina HiSeq X Ten. Germline and somatic, single nucleotide variants and small insertions or deletions were independently identified from data aligned human genome reference. The BGISEQ-500 and HiSeq X Ten platforms showed high concordance for germline calls with genotypes from SNP arrays (>99%). The germline and somatic single nucleotide variants identified in both sequencing platforms were highly concordant (86% and 72% respectively). These results indicate the potential applicability of the BGISEQ-500 platform for the identification of somatic and germline single nucleotide variants by whole genome sequencing. The BGISEQ-500 datasets described here represent the first publicly-available cancer genome sequencing performed using this platform.

Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments.

Pearl millet [Cenchrus americanus (L.) Morrone] is a staple food for more than 90 million farmers in arid and semi-arid regions of sub-Saharan Africa, India and South Asia. We report the ∼1.79 Gb draft whole genome sequence of reference genotype Tift 23D2B1-P1-P5, which contains an estimated 38,579 genes. We highlight the substantial enrichment for wax biosynthesis genes, which may contribute to heat and drought tolerance in this crop. We resequenced and analyzed 994 pearl millet lines, enabling insights into population structure, genetic diversity and domestication. We use these resequencing data to establish marker trait associations for genomic selection, to define heterotic pools, and to predict hybrid performance. We believe that these resources should empower researchers and breeders to improve this important staple crop.

The pearl oyster Pinctada fucata martensii genome and multi-omic analyses provide insights into biomineralization.

Nacre, the iridescent material found in pearls and shells of molluscs, is formed through an extraordinary process of matrix-assisted biomineralization. Despite recent advances, many aspects of the biomineralization process and its evolutionary origin remain unknown. The pearl oyster Pinctada fucata martensii is a well-known master of biomineralization, but the molecular mechanisms that underlie its production of shells and pearls are not fully understood. We sequenced the highly polymorphic genome of the pearl oyster and conducted multi-omic and biochemical studies to probe nacre formation. We identified a large set of novel proteins participating in matrix-framework formation, many in expanded families, including components similar to that found in vertebrate bones such as collagen-related VWA-containing proteins, chondroitin sulfotransferases, and regulatory elements. Considering that there are only collagen-based matrices in vertebrate bones and chitin-based matrices in most invertebrate skeletons, the presence of both chitin and elements of collagen-based matrices in nacre suggests that elements of chitin- and collagen-based matrices have deep roots and might be part of an ancient biomineralizing matrix. Our results expand the current shell matrix-framework model and provide new insights into the evolution of diverse biomineralization systems.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype.

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (Ar ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with Ar , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

A reference human genome dataset of the BGISEQ-500 sequencer.

Background: BGISEQ-500 is a new desktop sequencer developed by BGI. Using DNA nanoball and combinational probe anchor synthesis developed from Complete Genomics™ sequencing technologies, it generates short reads at a large scale. Here, we present the first human whole-genome sequencing dataset of BGISEQ-500. The dataset was generated by sequencing the widely used cell line HG001 (NA12878) in two sequencing runs of paired-end 50 bp (PE50) and two sequencing runs of paired-end 100 bp (PE100). We also include examples of the raw images from the sequencer for reference. Finally, we identified variations using this dataset, estimated the accuracy of the variations, and compared to that of the variations identified from similar amounts of publicly available HiSeq2500 data. We found similar single nucleotide polymorphism (SNP) detection accuracy for the BGISEQ-500 PE100 data (false positive rate [FPR] = 0.00020%, sensitivity = 96.20%) compared to the PE150 HiSeq2500 data (FPR = 0.00017%, sensitivity = 96.60%) better SNP detection accuracy than the PE50 data (FPR = 0.0006%, sensitivity = 94.15%). But for insertions and deletions (indels), we found lower accuracy for BGISEQ-500 data (FPR = 0.00069% and 0.00067% for PE100 and PE50 respectively, sensitivity = 88.52% and 70.93%) than the HiSeq2500 data (FPR = 0.00032%, sensitivity = 96.28%). Our dataset can serve as the reference dataset, providing basic information not just for future development, but also for all research and applications based on the new sequencing platform.

Draft genome of the living fossil Ginkgo biloba.

BACKGROUND: Ginkgo biloba L. (Ginkgoaceae) is one of the most distinctive plants. It possesses a suite of fascinating characteristics including a large genome, outstanding resistance/tolerance to abiotic and biotic stresses, and dioecious reproduction, making it an ideal model species for biological studies. However, the lack of a high-quality genome sequence has been an impediment to our understanding of its biology and evolution. FINDINGS: The 10.61 Gb genome sequence containing 41,840 annotated genes was assembled in the present study. Repetitive sequences account for 76.58% of the assembled sequence, and long terminal repeat retrotransposons (LTR-RTs) are particularly prevalent. The diversity and abundance of LTR-RTs is due to their gradual accumulation and a remarkable amplification between 16 and 24 million years ago, and they contribute to the long introns and large genome. Whole genome duplication (WGD) may have occurred twice, with an ancient WGD consistent with that shown to occur in other seed plants, and a more recent event specific to ginkgo. Abundant gene clusters from tandem duplication were also evident, and enrichment of expanded gene families indicates a remarkable array of chemical and antibacterial defense pathways. CONCLUSIONS: The ginkgo genome consists mainly of LTR-RTs resulting from ancient gradual accumulation and two WGD events. The multiple defense mechanisms underlying the characteristic resilience of ginkgo are fostered by a remarkable enrichment in ancient duplicated and ginkgo-specific gene clusters. The present study sheds light on sequencing large genomes, and opens an avenue for further genetic and evolutionary research.

The genetic basis for ecological adaptation of the Atlantic herring revealed by genome sequencing.

Ecological adaptation is of major relevance to speciation and sustainable population management, but the underlying genetic factors are typically hard to study in natural populations due to genetic differentiation caused by natural selection being confounded with genetic drift in subdivided populations. Here, we use whole genome population sequencing of Atlantic and Baltic herring to reveal the underlying genetic architecture at an unprecedented detailed resolution for both adaptation to a new niche environment and timing of reproduction. We identify almost 500 independent loci associated with a recent niche expansion from marine (Atlantic Ocean) to brackish waters (Baltic Sea), and more than 100 independent loci showing genetic differentiation between spring- and autumn-spawning populations irrespective of geographic origin. Our results show that both coding and non-coding changes contribute to adaptation. Haplotype blocks, often spanning multiple genes and maintained by selection, are associated with genetic differentiation.