A fine-scale genetic map of the Japanese population

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.

High-Density Genetic Map Construction and Quantitative Trait Locus Analysis of Fruit- and Oil-Related Traits in Camellia oleifera Based on Double Digest Restriction Site-Associated DNA Sequencing

Camellia oleifera, an important tree species and source of edible oil in China, has received significant attention owing to the oil's high unsaturated fatty acid content, which has benefits for human health. However, the mechanisms underlying C. oleifera yield and oil quality are largely unknown. In this study, 180 F1 progenies were obtained from two parents with obvious differences in fruit- and oil-related traits. We constructed a high-density genetic map using a double digest restriction site-associated DNA sequencing (ddRAD-Seq) strategy in C. oleifera. This map spanned 3327 cM and anchored 2780 markers in 15 linkage groups (LGs), with an average marker interval of 1.20 cM. A total of 221 quantitative trait loci (QTLs) associated with fruit- and oil-related traits were identified across three years' worth of phenotypic data. Nine QTLs were detected simultaneously in at least two different years, located on LG02, LG04, LG05, LG06, and LG11, and explained 8.5-16.6% of the phenotypic variation in the corresponding traits, respectively. Seventeen major QTLs were obtained that explained 13.0-16.6% of the phenotypic variance. Eleven and five flanking SNPs of major QTLs for fruit- and oil-related traits were detected which could be used for marker-assisted selection in C. oleifera breeding programs. Furthermore, 202 potential candidate genes in QTL regions were identified based on the collinearity of the genetic map and the C. oleifera "CON" genome. A potential regulatory network controlling fruit development and oil biosynthesis was constructed to dissect the complex mechanism of oil accumulation. The dissection of these QTLs will facilitate the gene cloning underlying lipid synthesis and increase our understanding in order to enhance C. oleifera oil yield and quality.

Quantitative Trait Loci Mapping and Comparative Transcriptome Analysis of Fruit Weight (FW) in Watermelon (Citrullus lanatus L.)

The watermelon (Citrullus lanatus L.) holds substantial economic value as a globally cultivated horticultural crop. However, the genetic architecture of watermelon fruit weight (FW) remains poorly understood. In this study, we used sh14-11 with small fruit and N14 with big fruit to construct 100 recombinant inbred lines (RILs). Based on whole-genome resequencing (WGR), 218,127 single nucleotide polymorphisms (SNPs) were detected to construct a high-quality genetic map. After quantitative trait loci (QTL) mapping, a candidate interval of 31-38 Mb on chromosome 2 was identified for FW. Simultaneously, the bulked segregant analysis (BSA) in the F2 population corroborated the identification of the same interval, encompassing the homologous gene linked to the known FW-related gene fas. Additionally, RNA-seq was carried out across 11 tissues from sh14-11 and N14, revealing expression profiles that identified 1695 new genes and corrected the annotation of 2941 genes. Subsequent differential expression analysis unveiled 8969 differentially expressed genes (DEGs), with 354 of these genes exhibiting significant differences across four key developmental stages. The integration of QTL mapping and differential expression analysis facilitated the identification of 14 FW-related genes, including annotated TGA and NAC transcription factors implicated in fruit development. This combined approach offers valuable insights into the genetic basis of FW, providing crucial resources for enhancing watermelon cultivation.

High-density genetic map and quantitative trait loci map of skin color in hawthorn (Crataegus pinnatifida bge. Var. major N.E.Br.)

Fruit skin color is an important trait of the hawthorn tree, which has an important influence on fruit quality. Crataegus pinnatifida Bge. var. Major N.E.Br. Is one of the most widely cultivated varieties in China and has a long history of medicinal use. In recent years, it has attracted the attention of the world due to its nutritional and medicinal values. Skin color is the focus of breeders and food processors. At present, skin color-related genes have still not been mapped. In this study, "Shandong Da Mianqiu" (♀, red skin color), "Da Huang Mianzha" (♂, yellow skin color) and 131 F1 hybrids were used to construct genetic map of hawthorn by RAD-seq, and QTL mapping was performed by combining these features with the hue angle and the observed color. In this study, 13,260 SNP was assigned to 17 linkage groups, with an integrated map covering 2,297.75 cM was constructed. A total of 5 QTLs related to hawthorn skin color were detected on LG1, LG3 and LG15. Whether hue angle or pericarp color acts as phenotype for QTL mapping, the candidate genes include bHLH086, WD repeat regions and Myb-like. bHLH, WD and Myb play an important role in the color regulation of Hawthorn skin color. These results lay a solid foundation for QTL mapping and molecular marker-assisted breeding of hawthorn.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genome-wide QTL mapping for agronomic traits in the winter wheat cultivar Pindong 34 based on 90K SNP array

Introduction

Agronomic traits are key components of wheat yield. Exploitation of the major underlying quantitative trait loci (QTLs) can improve the yield potential in wheat breeding.

Methods

In this study, we constructed a recombinant inbred line (RIL) population from Mingxian 169 (MX169) and Pindong 34 (PD34) to determine the QTLs for grain length (GL), grain width (GW), grain length-to-width ratio (LWR), plant height (PH), spike length (SL), grain number per spike (GNS), and the thousand grain weight (TGW) across four environments using wheat 90K SNP array.

Results

A QTL associated with TGW, i.e., QTGWpd.swust-6BS, was identified on chromosome 6B, which explained approximately 14.1%-16.2% of the phenotypic variation. In addition, eight QTLs associated with GL were detected across six chromosomes in four different test environments. These were QGLpd.swust-1BL, QGLpd.swust-2BL, QGLpd.swust-3BL.1, QGLpd.swust-3BL.2, QGLpd.swust-5DL, QGLpd.swust-6AL, QGLpd.swust-6DL.1, and QGLpd.swust-6DL.2. They accounted for 9.0%-21.3% of the phenotypic variation. Two QTLs, namely, QGWpd.swust-3BS and QGWpd.swust-6DL, were detected for GW on chromosomes 3B and 6D, respectively. These QTLs explained 12.8%-14.6% and 10.8%-15.2% of the phenotypic variation, respectively. In addition, two QTLs, i.e., QLWRpd.swust-7AS.1 and QLWRpd.swust-7AS.2, were detected on chromosome 7A for the grain LWR, which explained 10.9%-11.6% and 11.6%-11.2% of the phenotypic variation, respectively. Another QTL, named QGNSpd-swust-6DS, was discovered on chromosome 6D, which determines the GNS and which accounted for 11.4%-13.8% of the phenotypic variation. Furthermore, five QTLs associated with PH were mapped on chromosomes 2D, 3A, 5A, 6B, and 7B. These QTLs were QPHpd.swust-2DL, QPHpd.swust-3AL, QPHpd.swust-5AL, QPHpd.swust-6BL, and QPHpd.swust-7BS, which accounted for 11.3%-19.3% of the phenotypic variation. Lastly, a QTL named QSLpd.swust-3AL, conferring SL, was detected on chromosome 3A and explained 16.1%-17.6% of the phenotypic variation. All of these QTLs were defined within the physical interval of the Chinese spring reference genome.

Discussion

The findings of this study have significant implications for the development of fine genetic maps, for genomic breeding, and for marker-assisted selection to enhance wheat grain yield.

A SLAF-based high-density genetic map construction and genetic architecture of thermotolerant traits in maize (Zea mays L.)

The leaf scorching trait at flowering is a crucial thermosensitive phenotype in maize under high temperature stress (HS), yet the genetic basis of this trait remains poorly understood. In this study, we genotyped a 254 RIL-F2:8 population, derived from the leaf scorch-free parental inbred line Abe2 and the leaf scorching maternal inbred line B73, using the specific-locus amplified fragment sequencing (SLAF-seq) method. A total of 10,112 polymorphic SLAF markers were developed, and a high-density genetic map with a total length of 1,475.88 cM was constructed. The average sequencing depth of the parents was 55.23X, and that of the progeny was 12.53X. Then, we identified a total of 16 QTLs associated with thermotolerant traits at flowering, of which four QTLs of leaf scorching damage (LS) were distributed on chromosomes 1 (qLS1), 2 (qLS2.1, qLS2.2) and 3 (qLS3), which could explain 19.73% of phenotypic variation. Combining one qLS1 locus with QTL-seq results led to the identification of 6 candidate genes. Expression experiments and sequence variation indicated that Zm00001d033328, encoding N-acetyl-gamma-glutamyl-phosphate reductase, was the most likely candidate gene controlling thermotolerant traits at flowering. In summary, the high-density genetic map and genetic basis of thermotolerant traits lay a critical foundation for mapping other complex traits and identifying the genes associated with thermotolerant traits in maize.

Genetic Dissection of Panicle Morphology Traits in Super High-Yield Hybrid Rice Chaoyou 1000

The morphological characteristics of the rice panicle play a pivotal role in influencing yield. In our research, we employed F2 and F2:3 populations derived from the high-yielding hybrid rice variety Chaoyou 1000. We screened 123 pairs of molecular markers, which were available, to construct the genetic linkage map. Subsequently, we assessed the panicle morphology traits of F2 populations in Lingshui County, Hainan Province, in 2017, and F2:3 populations in Hangzhou City, Zhejiang Province, in 2018. These two locations represent two types of ecology. Hangzhou's climate is characterized by high temperatures and humidity, while Lingshui's climate is characterized by a tropical monsoon climate. In total, 33 QTLs were identified, with eight of these being newly discovered, and two of them were consistently detected in two distinct environments. We identified fourteen QTL-by-environment interactions (QEs), which collectively explained 4.93% to 59.95% of the phenotypic variation. While most of the detected QTLs are consistent with the results of previous tests, the novel-detected QTLs will lay the foundation for rice yield increase and molecular breeding.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3-45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

QTLs and Candidate Genes for Seed Protein Content in Two Recombinant Inbred Line Populations of Soybean

This study aimed to discover the quantitative trait loci (QTL) associated with a high seed protein content in soybean and unravel the potential candidate genes. We developed two recombinant inbred line populations: YS and SI, by crossing Saedanbaek (high protein) with YS2035-B-91-1-B-1 (low protein) and Saedanbaek with Ilmi (low protein), respectively, and evaluated the protein content for three consecutive years. Using single-nucleotide polymorphism (SNP)-marker-based linkage maps, four QTLs were located on chromosomes 15, 18, and 20 with high logarithm of odds values (5.9-55.0), contributing 5.5-66.0% phenotypic variance. In all three experimental years, qPSD20-1 and qPSD20-2 were stable and identified in overlapping positions in the YS and SI populations, respectively. Additionally, novel QTLs were identified on chromosomes 15 and 18. Considering the allelic sequence variation between parental lines, 28 annotated genes related to soybean seed protein-including starch, lipid, and fatty acid biosynthesis-related genes-were identified within the QTL regions. These genes could potentially affect protein accumulation during seed development, as well as sucrose and oil metabolism. Overall, this study offers insights into the genetic mechanisms underlying a high soybean protein content. The identified potential candidate genes can aid marker-assisted selection for developing soybean lines with an increased protein content.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

QTL Mapping and Transcriptome Analysis Reveal Candidate Genes Regulating Seed Color in Brassica napus

Yellow seeds are desirable in rapeseed breeding because of their higher oil content and better nutritional quality than black seeds. However, the underlying genes and formation mechanism of yellow seeds remain unclear. Here, a novel yellow-seeded rapeseed line (Huangaizao, HAZ) was crossed with a black-seeded rapeseed line (Zhongshuang11, ZS11) to construct a mapping population of 196 F₂ individuals, based on which, a high-density genetic linkage map was constructed. This map, comprising 4174 bin markers, was 1618.33 cM in length and had an average distance of 0.39 cM between its adjacent markers. To assess the seed color of the F₂ population, three methods (imaging, spectrophotometry, and visual scoring) were used and a common major quantitative trait locus (QTL) on chromosome A09, explaining 10.91-21.83% of the phenotypic variance, was detected. Another minor QTL, accounting for 6.19-6.69% of the phenotypic variance, was detected on chromosome C03, only by means of imaging and spectrophotometry. Furthermore, a dynamic analysis of the differential expressions between the parental lines showed that flavonoid biosynthesis-related genes were down-regulated in the yellow seed coats at 25 and 35 days after flowering. A coexpression network between the differentially expressed genes identified 17 candidate genes for the QTL intervals, including a flavonoid structure gene, novel4557 (BnaC03.TT4), and two transcription factor genes, namely, BnaA09G0616800ZS (BnaA09.NFYA8) and BnaC03G0060200ZS (BnaC03.NAC083), that may regulate flavonoid biosynthesis. Our study lays a foundation for further identifying the genes responsible for and understanding the regulatory mechanism of yellow seed formation in Brassica napus.

Construction of a high-resolution genetic map and identification of single nucleotide polymorphism markers relevant to flower stalk height in onion

INTRODUCTION

Onion (Allium cepa L., 2n=16) is an economically and nutritionally important vegetable crop worldwide. Construction of a high-resolution genetic map and map-based gene mining in onion have lagged behind other vegetable crops such as tomato and pepper.

METHODS

In this study, we constructed a high-resolution genetic map of onion using 321 F2 individuals from a cross between two double haploid lines DH-1×DH-17 and employing specific length amplified fragment (SLAF)-seq technology. The genetic map containing 10,584 polymorphic SLAFs with 21,250 single nucleotide polymorphism (SNP) markers and 8 linkage groups was developed for onion, which spanned 928.32 cM, with an average distance of 0.09 cM between adjacent markers.

RESULTS

Using this map, we carried out QTL mapping of Ms locus related to the male-fertile trait and reproduced previous mapping results, which proved that this map was of good quality. Then, four QTLs (located on LG2, LG5, and LG8) were detected for flower stalk height, explaining 26.60% of the phenotypic variance. Among them, we proposed that 20 SLAF markers (in three QTLs) of flower stalk height trait were effective favorable allelic variant markers associated with heterosis.

DISCUSSION

Overall, the genetic map was structured using SLAF-seq based on DH lines, and it is the highest-quality and highest-resolution linkage map of onion to date. It lays a foundation for the fine mapping and candidate gene identification of flower stalk height, and provides new insights into the developmental genetic mechanisms in onion breeding.

QTL Mapping and a Transcriptome Integrative Analysis Uncover the Candidate Genes That Control the Cold Tolerance of Maize Introgression Lines at the Seedling Stage

Chilling injury owing to low temperatures severely affects the growth and development of maize (Zea mays.L) seedlings during the early and late spring seasons. The existing maize germplasm is deficient in the resources required to improve maize's ability to tolerate cold injury. Therefore, it is crucial to introduce and identify excellent gene/QTLs that confer cold tolerance to maize for sustainable crop production. Wild relatives of maize, such as Z. perennis and Tripsacum dactyloides, are strongly tolerant to cold and can be used to improve the cold tolerance of maize. In a previous study, a genetic bridge among maize that utilized Z. perennis and T. dactyloides was created and used to obtain a highly cold-tolerant maize introgression line (MIL)-IB030 by backcross breeding. In this study, two candidate genes that control relative electrical conductivity were located on MIL-IB030 by forward genetics combined with a weighted gene co-expression network analysis. The results of the phenotypic, genotypic, gene expression, and functional verification suggest that two candidate genes positively regulate cold tolerance in MIL-IB030 and could be used to improve the cold tolerance of cultivated maize. This study provides a workable route to introduce and mine excellent genes/QTLs to improve the cold tolerance of maize and also lays a theoretical and practical foundation to improve cultivated maize against low-temperature stress.

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.

Understanding the molecular mechanism of leaf morphogenesis in vegetable crops conduces to breeding process

Leaf morphology can affect the development and yield of plants by regulating plant architecture and photosynthesis. Several factors can determine the final leaf morphology, including the leaf complexity, size, shape, and margin type, which suggests that leaf morphogenesis is a complex regulation network. The formation of diverse leaf morphology is precisely controlled by gene regulation on translation and transcription levels. To further reveal this, more and more genome data has been published for different kinds of vegetable crops and advanced genotyping approaches have also been applied to identify the causal genes for the target traits. Therefore, the studies on the molecular regulation of leaf morphogenesis in vegetable crops have also been largely improved. This review will summarize the progress on identified genes or regulatory mechanisms of leaf morphogenesis and development in vegetable crops. These identified markers can be applied for further molecular-assisted selection (MAS) in vegetable crops. Overall, the review will contribute to understanding the leaf morphology of different crops from the perspective of molecular regulation and shortening the breeding cycle for vegetable crops.

Analysis of Genetic Diversity and Resistance to Foliar Pathogens Based on Genotyping-by-Sequencing of a Para Rubber Diversity Panel and Progeny of an Interspecific Cross

Para rubber trees (Hevea brasiliensis) are the largest major source of natural rubber in the world. Its major pathogens are Phytophthora spp., Corynespora cassiicola, and Colletotrichum spp. A rubber diversity panel of 116 clones using over 12,000 single nucleotide polymorphisms (SNPs) from DArTSeq genotyping revealed clear phylogenetic differences in clones that originated from different geographical regions of the world. An integrated linkage map constructed with an F₁ progeny of 86 from an interspecific cross between H. brasiliensis and H. benthamiana using 23,978 markers [10,323 SNPs and 13,655 SilicoDArTs] spanned 3947.83 cM with 0.83 cM average marker-interval. The genome scaffolds that were anchored to the linkage map, covering 1.44 Gb of H. brasiliensis reference genome, revealed a high level of collinearity between the genetic map and reference genome. Association analysis identified 12 SNPs significantly associated with the resistance against Phytophthora, Corynespora, and Colletotrichum in six linkage groups: 2, 6, 12, 14, 17, and 18. Kompetitive Allele-Specific PCR marker assays were developed for those 12 SNPs, screened with 178 individuals, and detected clear separation between two genotypes. Within the proximity to those SNPs, 41 potentially key genes that have previously been reported to associate with plant disease resistance were predicted with high confidence.

Chromosome-level reference genome for European flat oyster (Ostrea edulis L.)

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome-level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold-N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD-Seq technology, comprising 4016 SNPs from four full-sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab-initio gene prediction identified 35,699 protein-coding genes. Chromosome-level synteny was demonstrated against multiple high-quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high-resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

High-density linkage map construction in an autotetraploid blueberry population and detection of quantitative trait loci for anthocyanin content

Highbush blueberry (Vaccinium corymbosum, 2n = 4x = 48) is the most cultivated type of blueberry, both in New Zealand and overseas. Its perceived nutritional value is conferred by phytonutrients, particularly anthocyanins. Identifying the genetic mechanisms that control the biosynthesis of these metabolites would enable faster development of cultivars with improved fruit qualities. Here, we used recently released tools for genetic mapping in autotetraploids to build a high-density linkage map in highbush blueberry and to detect quantitative trait loci (QTLs) for fruit anthocyanin content. Genotyping was performed by target sequencing, with ∼18,000 single nucleotide polymorphism (SNP) markers being mapped into 12 phased linkage groups (LGs). Fruits were harvested when ripe for two seasons and analyzed with high-performance liquid chromatography-mass spectrometry (HPLC-MS): 25 different anthocyanin compounds were identified and quantified. Two major QTLs that were stable across years were discovered, one on LG2 and one on LG4, and the underlying candidate genes were identified. Interestingly, the presence of anthocyanins containing acylated sugars appeared to be under strong genetic control. Information gained in this study will enable the design of molecular markers for marker-assisted selection and will help build a better understanding of the genetic control of anthocyanin biosynthesis in this crop.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Quantitative trait loci identification and breeding value estimation of grain weight-related traits based on a new wheat 50K single nucleotide polymorphism array-derived genetic map

Mining novel and less utilized thousand grain weight (TGW) related genes are useful for improving wheat yield. In this study, a recombinant inbred line population from a cross between Zhongkemai 138 (ZKM138, high TGW) and Chuanmai 44 (CM44, low TGW) was used to construct a new Wheat 50K SNP array-derived genetic map that spanned 1,936.59 cM and contained 4, 139 markers. Based on this map, ninety-one quantitative trait loci (QTL) were detected for eight grain-related traits in six environments. Among 58 QTLs, whose superior alleles were contributed by ZKM138, QTgw.cib-6A was a noticeable major stable QTL and was also highlighted by bulked segregant analysis with RNA sequencing (BSR-Seq). It had a pyramiding effect on TGW enhancement but no significant trade-off effect on grain number per spike or tiller number, with two other QTLs (QTgw.cib-2A.2 and QTgw.cib-6D), possibly explaining the excellent grain performance of ZKM138. After comparison with known loci, QTgw.cib-6A was deduced to be a novel locus that differed from nearby TaGW2 and TaBT1. Seven simple sequence repeat (SSR) and thirty-nine kompetitive allele-specific PCR (KASP) markers were finally developed to narrow the candidate interval of QTgw.cib-6A to 4.1 Mb. Only six genes in this interval were regarded as the most likely candidate genes. QTgw.cib-6A was further validated in different genetic backgrounds and presented 88.6% transmissibility of the ZKM138-genotype and a 16.4% increase of TGW in ZKM138 derivatives. And the geographic pattern of this locus revealed that its superior allele is present in only 6.47% of 433 Chinese modern wheat varieties, indicating its potential contribution to further high-yield breeding.

Quantitative trait loci mapping reveals an oligogenic architecture of a rapidly adapting trait during the European invasion of common ragweed

Biological invasions offer a unique opportunity to investigate evolution over contemporary timescales. Rapid adaptation to local climates during range expansion can be a major determinant of invasion success, yet fundamental questions remain about its genetic basis. This study sought to investigate the genetic basis of climate adaptation in invasive common ragweed (Ambrosia artemisiifolia). Flowering time adaptation is key to this annual species' invasion success, so much so that it has evolved repeated latitudinal clines in size and phenology across its native and introduced ranges despite high gene flow among populations. Here, we produced a high-density linkage map (4493 SNPs) and paired this with phenotypic data from an F2 mapping population (n = 336) to identify one major and two minor quantitative trait loci (QTL) underlying flowering time and height differentiation in this species. Within each QTL interval, several candidate flowering time genes were also identified. Notably, the major flowering time QTL detected in this study was found to overlap with a previously identified haploblock (putative inversion). Multiple genetic maps of this region identified evidence of suppressed recombination in specific genotypes, consistent with inversions. These discoveries support the expectation that a concentrated genetic architecture with fewer, larger, and more tightly linked alleles should underlie rapid local adaptation during invasion, particularly when divergently adapting populations experience high levels of gene flow.

Mapping Genetic Events of SARS-CoV-2 Variants

Genetic mutation and recombination are driving the evolution of SARS-CoV-2, leaving many genetic imprints which could be utilized to track the evolutionary pathway of SARS-CoV-2 and explore the relationships among variants. Here, we constructed a complete genetic map, showing the explicit evolutionary relationship among all SARS-CoV-2 variants including 58 groups and 46 recombination types identified from 3,392,553 sequences, which enables us to keep well informed of the evolution of SARS-CoV-2 and quickly determine the parents of novel variants. We found that the 5' and 3' of the spike and nucleoprotein genes have high frequencies to form the recombination junctions and that the RBD region in S gene is always exchanged as a whole. Although these recombinants did not show advantages in community transmission, it is necessary to keep a wary eye on the novel genetic events, in particular, the mutants with mutations on spike and recombinants with exchanged moieties on spike gene.

Identification of QTL for Stem Traits in Wheat (Triticum aestivum L.)

Lodging in wheat (Triticum aestivum L.) is a complicated phenomenon that is influenced by physiological, genetics, and external factors. It causes a great yield loss and reduces grain quality and mechanical harvesting efficiency. Lodging resistance is contributed by various traits, including increased stem strength. The aim of this study was to map quantitative trait loci (QTL) controlling stem strength-related features (the number of big vascular bundles, stem diameter, stem wall thickness) using a doubled haploid (DH) population derived from a cross between Baiqimai and Neixiang 5. Field experiments were conducted during 2020-2022, and glasshouse experiments were conducted during 2021-2022. Significant genetic variations were observed for all measured traits, and they were all highly heritable. Fifteen QTL for stem strength-related traits were identified on chromosomes 2D, 3A, 3B, 3D, 4B, 5A, 6B, 7A, and 7D, respectively, and 7 QTL for grain yield-related traits were identified on chromosomes 2B, 2D, 3D, 4B, 7A, and 7B, respectively. The superior allele of the major QTL for the number of big vascular bundle (VB) was independent of plant height (PH), making it possible to improve stem strength without a trade-off of PH, thus improving lodging resistance. VB also showed positive correlations with some of the yield components. The result will be useful for molecular marker-assisted selection (MAS) for high stem strength and high yield potential.

Dynamic Quantitative Trait Loci Mapping for Plant Height in Recombinant Inbred Line Population of Upland Cotton

Plant height (PH) is a key plant architecture trait for improving the biological productivity of cotton. Ideal PH of cotton is conducive to lodging resistance and mechanized harvesting. To detect quantitative trait loci (QTL) and candidate genes of PH in cotton, a genetic map was constructed with a recombinant inbred line (RIL) population of upland cotton. PH phenotype data under nine environments and three best linear unbiased predictions (BLUPs) were used for QTL analyses. Based on restriction-site-associated DNA sequence (RAD-seq), the genetic map contained 5,850 single-nucleotide polymorphism (SNP) markers, covering 2,747.12 cM with an average genetic distance of 0.47 cM. Thirty-seven unconditional QTL explaining 1.03-12.50% of phenotypic variance, including four major QTL and seven stable QTL, were identified. Twenty-eight conditional QTL explaining 3.27-28.87% of phenotypic variance, including 1 major QTL, were identified. Importantly, five QTL, including 4 stable QTL, were both unconditional and conditional QTL. Among the 60 PH QTL (including 39 newly identified), none of them were involved in the whole period of PH growth, indicating that QTL related to cotton PH development have dynamic expression characteristics. Based on the functional annotation of Arabidopsis homologous genes and transcriptome data of upland cotton TM-1, 14 candidate genes were predicted within 10 QTL. Our research provides valuable information for understanding the genetic mechanism of PH development, which also increases the economic production of cotton.