Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong.

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

Haplotype-resolved assembly of a tetraploid potato genome using long reads and low-depth offspring data.

Potato is one of the world's major staple crops, and like many important crop plants, it has a polyploid genome. Polyploid haplotype assembly poses a major computational challenge. We introduce a novel strategy for the assembly of polyploid genomes and present an assembly of the autotetraploid potato cultivar Altus. Our method uses low-depth sequencing data from an offspring population to achieve chromosomal clustering and haplotype phasing on the assembly graph. Our approach generates high-quality assemblies of individual chromosomes with haplotype-specific sequence resolution of whole chromosome arms and can be applied in common breeding scenarios where collections of offspring are available.

Association between polymorphism and haplotype of ATP2B1 gene and skeletal fluorosis in Han population.

To evaluate the association between ATP2B1 gene polymorphisms and skeletal fluorosis, a cross-sectional study was conducted. In China, 962 individuals were recruited, including 342 cases of skeletal fluorosis. Four TP2BA1 polymorphisms (rs2070759, rs12817819, rs17249754, and rs7136259) were analysed. The results suggested that rs17249754 and rs7136259 were associated with skeletal fluorosis. After controlling confounders, the protective effect of GG genotype in rs17249754 was apparent in individuals over 45 years old, female, with urine fluoride concentration below 1.6 mg/L, serum calcium above 2.25 mmol/L or serum phosphorus between 1.1 and 1.3. Heterozygote TC in rs7136259 increased the risk of skeletal fluorosis in subjects who are elderly, female, with urinary fluoride more than 1.6 mg/L, serum calcium more than 2.25 mmol/L and blood phosphorus between 1.1 and 1.3 mmol/L. Four loci were found to be tightly related by linkage disequilibrium analysis, and the frequency of distribution of haplotype GCGT was lower in the skeletal fluorosis group.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones.

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

Haplotype-resolved chromosomal-level genome assembly of Buzhaye (Microcos paniculata).

Microcos paniculata is a shrub used traditionally as folk medicine and to make herbal teas. Previous research into this species has mainly focused on its chemical composition and medicinal value. However, the lack of a reference genome limits the study of the molecular mechanisms of active compounds in this species. Here, we assembled a haplotype-resolved chromosome-level genome of M. paniculata based on PacBio HiFi and Hi-C data. The assembly contains two haploid genomes with sizes 399.43 Mb and 393.10 Mb, with contig N50 lengths of 43.44 Mb and 30.17 Mb, respectively. About 99.93% of the assembled sequences could be anchored to 18 pseudo-chromosomes. Additionally, a total of 482 Mb repeat sequences were identified, accounting for 60.76% of the genome. A total of 49,439 protein-coding genes were identified, of which 48,979 (99%) were functionally annotated. This haplotype-resolved chromosome-level assembly and annotation of M. paniculata will serve as a valuable resource for investigating the biosynthesis and genetic basis of active compounds in this species, as well as advancing evolutionary phylogenomic studies in Malvales.

Human genetic adaptation related to cellular zinc homeostasis.

SLC30A9 encodes a ubiquitously zinc transporter (ZnT9) and has been consistently suggested as a candidate for positive selection in humans. However, no direct adaptive molecular phenotype has been demonstrated. Our results provide evidence for directional selection operating in two major complementary haplotypes in Africa and East Asia. These haplotypes are associated with differential gene expression but also differ in the Met50Val substitution (rs1047626) in ZnT9, which we show is found in homozygosis in the Denisovan genome and displays accompanying signatures suggestive of archaic introgression. Although we found no significant differences in systemic zinc content between individuals with different rs1047626 genotypes, we demonstrate that the expression of the derived isoform (ZnT9 50Val) in HEK293 cells shows a gain of function when compared with the ancestral (ZnT9 50Met) variant. Notably, the ZnT9 50Val variant was found associated with differences in zinc handling by the mitochondria and endoplasmic reticulum, with an impact on mitochondrial metabolism. Given the essential role of the mitochondria in skeletal muscle and since the derived allele at rs1047626 is known to be associated with greater susceptibility to several neuropsychiatric traits, we propose that adaptation to cold may have driven this selection event, while also impacting predisposition to neuropsychiatric disorders in modern humans.

Haplotype-Resolution Transcriptome Analysis Reveals Important Responsive Gene Modules and Allele-Specific Expression Contributions under Continuous Salt and Drought in Camellia sinensis.

The tea plant, Camellia sinensis (L.) O. Kuntze, is one of the most important beverage crops with significant economic and cultural value. Global climate change and population growth have led to increased salt and drought stress, negatively affecting tea yield and quality. The response mechanism of tea plants to these stresses remains poorly understood due to the lack of reference genome-based transcriptional descriptions. This study presents a high-quality genome-based transcriptome dynamic analysis of C. sinensis' response to salt and drought stress. A total of 2244 upregulated and 2164 downregulated genes were identified under salt and drought stress compared to the control sample. Most of the differentially expression genes (DEGs) were found to involve divergent regulation processes at different time points under stress. Some shared up- and downregulated DEGs related to secondary metabolic and photosynthetic processes, respectively. Weighted gene co-expression network analysis (WGCNA) revealed six co-expression modules significantly positively correlated with C. sinensis' response to salt or drought stress. The MEpurple module indicated crosstalk between the two stresses related to ubiquitination and the phenylpropanoid metabolic regulation process. We identified 1969 salt-responsive and 1887 drought-responsive allele-specific expression (ASE) genes in C. sinensis. Further comparison between these ASE genes and tea plant heterosis-related genes suggests that heterosis likely contributes to the adversity and stress resistance of C. sinensis. This work offers new insight into the underlying mechanisms of C. sinensis' response to salt and drought stress and supports the improved breeding of tea plants with enhanced salt and drought tolerance.

Potato psyllids mount distinct gut responses against two different 'Candidatus Liberibacter solanacearum' haplotypes.

'Candidatus Liberibacter solanacearum' (Lso) is a bacterial pathogen infecting several crops and causing damaging diseases. Several Lso haplotypes have been identified. Among the seven haplotypes present in North America, LsoA and LsoB are transmitted by the potato psyllid, Bactericera cockerelli (Sulc), in a circulative and persistent manner. The gut, which is the first organ pathogen encounters, could be a barrier for Lso transmission. However, the molecular interactions between Lso and the psyllid vector at the gut interface remain largely unknown. In this study, we investigated the global transcriptional responses of the adult psyllid gut upon infection with two Lso haplotypes (LsoA and LsoB) using Illumina sequencing. The results showed that each haplotype triggers a unique transcriptional response, with most of the distinct genes elicited by the highly virulent LsoB. The differentially expressed genes were mainly associated with digestion and metabolism, stress response, immunity, detoxification as well as cell proliferation and epithelium renewal. Importantly, distinct immune pathways were triggered by LsoA and LsoB in the gut of the potato psyllid. The information in this study will provide an understanding of the molecular basis of the interactions between the potato psyllid gut and Lso, which may lead to the discovery of novel molecular targets for the control of these pathogens.

Using haplotype and QTL analysis to fix favorable alleles in diploid potato breeding.

At present, the potato (Solanum tuberosum L.) of international commerce is autotetraploid, and the complexity of this genetic system creates limitations for breeding. Diploid potato breeding has long been used for population improvement, and because of an improved understanding of the genetics of gametophytic self-incompatibility, there is now sustained interest in the development of uniform F1 hybrid varieties based on inbred parents. We report here on the use of haplotype and quantitative trait locus (QTL) analysis in a modified backcrossing (BC) scheme, using primary dihaploids of S. tuberosum as the recurrent parental background. In Cycle 1, we selected XD3-36, a self-fertile F2 individual homozygous for the self-compatibility gene Sli (S-locus inhibitor). Signatures of gametic and zygotic selection were observed at multiple loci in the F2 generation, including Sli. In the BC1 cycle, an F1 population derived from XD3-36 showed a bimodal response for vine maturity, which led to the identification of late versus early alleles in XD3-36 for the gene CDF1 (Cycling DOF Factor 1). Greenhouse phenotypes and haplotype analysis were used to select a vigorous and self-fertile F2 individual with 43% homozygosity, including for Sli and the early-maturing allele CDF1.3. Partially inbred lines from the BC1 and BC2 cycles have been used to initiate new cycles of selection, with the goal of reaching higher homozygosity while maintaining plant vigor, fertility, and yield.

Haplotype-resolved genomes of two buckwheat crops provide insights into their contrasted rutin concentrations and reproductive systems.

BACKGROUND: Two widely cultivated annual buckwheat crops, Fagopyrum esculentum and F. tataricum, differ from each other in both rutin concentration and reproductive system. However, the underlying genetic mechanisms remain poorly elucidated. RESULTS: Here, we report the first haplotype-resolved chromosome-level genome assemblies of the two species. Two haplotype genomes of F. esculentum were assembled as 1.23 and 1.19 Gb with N50 = 9.8 and 12.4 Mb, respectively; the two haplotype genomes of F. tataricum were 453.7 and 446.2 Mb with N50 = 50 and 30 Mb, respectively. We further annotated protein-coding genes of each haplotype genome based on available gene sets and 48 newly sequenced transcriptomes. We found that more repetitive sequences, especially expansion of long terminal repeat retrotransposons (LTR-RTs), contributed to the large genome size of F. esculentum. Based on the well-annotated sequences, gene expressions, and luciferase experiments, we identified the sequence mutations of the promoter regions of two key genes that are likely to have greatly contributed to the high rutin concentration and selfing reproduction in F. tartaricum. CONCLUSIONS: Our results highlight the importance of high-quality genomes to identify genetic mutations underlying phenotypic differences between closely related species. F. tataricum may have been experienced stronger selection than F. esculentum through choosing these two non-coding alleles for the desired cultivation traits. These findings further suggest that genetic manipulation of the non-coding promoter regions could be widely employed for breeding buckwheat and other crops.

Intralymphatic glutamic acid decarboxylase administration in type 1 diabetes patients induced a distinctive early immune response in patients with DR3DQ2 haplotype.

GAD-alum given into lymph nodes to Type 1 diabetes (T1D) patients participating in a multicenter, randomized, placebo-controlled double-blind study seemed to have a positive effect for patients with DR3DQ2 haplotype, who showed better preservation of C-peptide than the placebo group. Here we compared the immunomodulatory effect of GAD-alum administered into lymph nodes of patients with T1D versus placebo with focus on patients with DR3DQ2 haplotype. Methods: GAD autoantibodies, GADA subclasses, GAD65-induced cytokine secretion (Luminex panel) and proliferation of peripheral mononuclear cells were analyzed in T1D patients (n=109) who received either three intra-lymphatic injections (one month apart) with 4 µg GAD-alum and oral vitamin D supplementation (2000 IE daily for 120 days), or placebo. Results: Higher GADA, GADA subclasses, GAD65-induced proliferation and cytokine secretion was observed in actively treated patients after the second injection of GAD-alum compared to the placebo group. Following the second injection of GAD-alum, actively treated subjects with DR3DQ2 haplotype had higher GAD65-induced secretion of several cytokine (IL4, IL5, IL7, IL10, IL13, IFNγ, GM-CSF and MIP1ß) and proliferation compared to treated individuals without DR3DQ2. Stratification of samples from GAD-alum treated patients according to C-peptide preservation at 15 months revealed that "good responder" individuals with better preservation of C-peptide secretion, independently of the HLA haplotype, had increased GAD65-induced proliferation and IL13 secretion at 3 months, and a 2,5-fold increase of IL5 and IL10 as compared to "poor responders". The second dose of GAD-alum also induced a more pronounced cytokine secretion in "good responders" with DR3DQ2, compared to few "good responders" without DR3DQ2 haplotype. Conclusion: Patients with DR3DQ2 haplotype had a distinct early cellular immune response to GAD-alum injections into the lymph node, and predominant GAD65-induced IL13 secretion and proliferation that seems to be associated with a better clinical outcome. If confirmed in the ongoing larger randomized double-blind placebo-controlled clinical trial (DIAGNODE-3), including only patients carrying DR3DQ2 haplotype, these results might be used as early surrogate markers for clinical efficacy.

Insights into the Genetic Determination of the Autotetraploid Potato Plant Height.

Plant height is an important characteristic, the modification of which can improve the ability of stress adaptation as well as the yield. In this study, genome-wide association analysis was performed for plant height traits in 370 potato cultivars using the tetraploid potato genome as a reference. A total of 92 significant single nucleotide polymorphism (SNP) loci for plant height were obtained, which were particularly significant in haplotypes A3 and A4 on chromosome 1 and A1, A2, and A4 on chromosome 5. Thirty-five candidate genes were identified that were mainly involved in the gibberellin and brassinolide signal transduction pathways, including the FAR1 gene, methyltransferase, ethylene response factor, and ubiquitin protein ligase. Among them, PIF3 and GID1a were only present on chromosome 1, with PIF3 in all four haplotypes and GID1a in haplotype A3. This could lead to more effective genetic loci for molecular marker-assisted selection breeding as well as more precise localization and cloning of genes for plant height traits in potatoes.

Acquisition and Transmission of 'Candidatus Liberibacter solanacearum' Differs Among Wolbachia-Infected and -Uninfected Haplotypes of Bactericera cockerelli.

'Candidatus Liberibacter solanacearum' (Lso) causes disease symptoms and economic losses in potato, tomato, and other solanaceous crops in North America. Lso is transmitted to plants by the potato psyllid, Bactericera cockerelli, which occurs as distinct haplotypes named western, central, and northwestern that differ in the presence or absence of the bacterial endosymbiont, Wolbachia. Previous work showed that all three vector haplotypes can transmit Lso, but it was not clear whether acquisition and transmission rates of Lso were equal among the haplotypes. The goal of our study was to compare Lso infection rates among psyllids of the western, central, and northwestern haplotypes. Using data collected from several years of periodic testing of Lso infection of laboratory-reared potato psyllid colonies, we showed that psyllids of the western and central haplotypes are more likely to harbor Lso than are psyllids of the northwestern haplotype. We then used greenhouse assays to demonstrate that psyllids of the northwestern haplotype are less likely to acquire and transmit Lso than those of the western haplotype. Lso infection rates corresponded with Wolbachia infection among the three psyllid haplotypes. The Wolbachia-infected central and western haplotypes were more likely to harbor and transmit Lso than the Wolbachia-free northwestern haplotype. Results demonstrate that potato psyllids of the western and central haplotypes pose a greater risk for spread of Lso in crops and suggest a pattern between infection with Lso and Wolbachia in potato psyllid.

Development of a Loop-Mediated Isothermal Amplification (LAMP) Method to Detect the Potato Zebra Chip Pathogen 'Candidatus Liberibacter solanacearum' (Lso) and Differentiate Haplotypes A and B.

'Candidatus Liberibacter solanacearum' (Lso) is the causal agent of zebra chip of potato (Solanum tuberosum), which can significantly reduce potato yield. In this study, a loop-mediated isothermal amplification (LAMP) method for the detection of Lso haplotypes A and B was developed and evaluated. Two sets of LAMP primers named LAMP-A and LAMP-B were designed and tested for specificity and sensitivity. Both LAMP-A and LAMP-B were specific to Lso in in silico analysis using the Primer-Blast tool. The LAMP-A and LAMP-B could only produce positive signals from DNA mixtures of Lso-infected tomato but not from the genomic DNA of 37 nontarget plant pathogens. The sensitivity of LAMP-A and LAMP-B on Lso haplotypes A and B were tested on gBlocks and genomic DNA from Lso-infected tomato. On the genomic DNA for LAMP-A, the lowest amount of template DNA for a positive LAMP reaction was 2 to 20 ng on four haplotype A strains and 20 to 80 ng on four haplotype B strains; for LAMP-B, the lowest amount of template DNA for a positive LAMP reaction was 0.02 to 2 ng on four haplotype B strains and 20 ng to no amplification on four haplotype A strains. On gBlocks for LAMP-A, the lowest number of copies for a positive LAMP reaction was 60 on haplotype A and 600 on haplotype B; for LAMP-B, the lowest number of copies for a positive LAMP reaction was 60 on haplotype B and 600 on haplotype A. Therefore, considering the convenience of the LAMP technique, as well as the high specificity and sensitivity, the LAMP-A and LAMP-B primers can be used together to test the probable Lso-infected plant or psyllid samples to rapidly, accurately, and directly differentiate haplotypes A and B. We highly recommend this LAMP system to plant pathology practitioners and diagnostic labs for routine detection of Lso and confirmation of zebra chip disease on potato or tomato.

Identification of three new 'Candidatus Liberibacter solanacearum' haplotypes in four psyllid species (Hemiptera: Psylloidea).

Eleven haplotypes of the bacterium, 'Candidatus Liberibacter solanacearum', have been identified worldwide, several of which infect important agricultural crops. In the United States, haplotypes A and B are associated with yield and quality losses in potato, tomato, and other crops of the Solanaceae. Both haplotypes are vectored by potato psyllid, Bactericera cockerelli. Recently, a third haplotype, designated F, was identified in southern Oregon potato fields. To identify the vector of this haplotype, psyllids of multiple species were collected from yellow sticky cards placed near potato fields during two growing seasons. Over 2700 specimens were tested for 'Ca. L. solanacearum' by polymerase chain reaction. Forty-seven psyllids harbored the bacterium. The infected specimens comprised four psyllid species in two families, Aphalaridae and Triozidae (Hemiptera: Psylloidea). Nucleic acid and/or amino acid sequence analysis of the 'Ca. L. solanacearum' 16S ribosomal RNA, 50S ribosomal proteins L10/L12, and outer membrane protein identified three new haplotypes of the bacterium, designated as Aph1, Aph2 and Aph3, including two variants of Aph2 (Aph2a and Aph2b). The impact of these new haplotypes on solanaceous or other crops is not known. The vector of 'Ca. L. solanacearum' haplotype F was not detected in this study.

Phase III, randomised, double-blind, placebo-controlled, multicentre trial to evaluate the efficacy and safety of rhGAD65 to preserve endogenous beta cell function in adolescents and adults with recently diagnosed type 1 diabetes, carrying the genetic HLA DR3-DQ2 haplotype: the DIAGNODE-3 study protocol.

INTRODUCTION: Type 1 diabetes (T1D) is an autoimmune disease leading to the destruction of the insulin-producing beta cells resulting in insulin deficiency and hyperglycaemic. Today, no approved therapy exists to halt this detrimental immunologic process. In a recent phase 2b study, intralymphatic administration of recombinant human glutamic acid decarboxylase 65 kDa (rhGAD65) adsorbed to Alhydrogel adjuvant to individuals recently diagnosed with T1D and carrying the HLA DR3-DQ2 haplotype showed promising results in preserving endogenous insulin secretion, confirming the results of a large meta-analysis of three randomised placebo-controlled trials of subcutaneous rhGAD65. The aim of the current precision medicine phase 3 study is to determine whether intralymphatic administration of rhGAD65 preserves insulin secretion and improves glycaemic control in presumed responder individuals with recently diagnosed T1D carrying HLA DR3-DQ2. METHODS AND ANALYSIS: Individuals ≥12 and <29 years recently diagnosed with T1D (<6 months) will be screened for the HLA DR3-DQ2 haplotype, endogenous insulin production estimated by fasting C-peptide and presence of GAD65 antibodies. 330 patients are planned to be randomised to 3 monthly intralymphatic injections of rhGAD65 or placebo (both accompanied by oral vitamin D supplementation), followed by 22 months of follow-up. The study is powered to detect a treatment effect in the two coprimary endpoints; change from baseline in AUC(0-120min) C-peptide levels during a mixed meal tolerance test, and change from baseline in glycaemic control estimated by haemoglobin A1c at 24 months. Secondary endpoints include effects on glucose patterns collected by masked continuous glucose monitoring, proportion of patients in partial remission and number of episodes of severe hypoglycaemia and/or diabetic ketoacidosis. ETHICS AND DISSEMINATION: The trial is approved by Ethics Committees in Poland (124/2021), the Netherlands (R21.089), Sweden (2021-05063), Czech Republic (EK-1144/21), Germany (2021361) and Spain (21/2021). Results will be published in international peer-reviewed scientific journals and presented at national and international conferences. TRIAL REGISTRATION NUMBER: EudraCT identifier: 2021-002731-32, NCT identifier: NCT05018585.

Genetic and particle modelling approaches to assessing population connectivity in a deep sea lobster.

The emergence of high resolution population genetic techniques, such as genotyping-by-sequencing (GBS), in combination with recent advances in particle modelling of larval dispersal in marine organisms, can deliver powerful new insights to support fisheries conservation and management. In this study, we used this combination to investigate the population connectivity of a commercial deep sea lobster species, the New Zealand scampi, Metanephrops challengeri, which ranges across a vast area of seafloor around New Zealand. This species has limited dispersal capabilities, including larvae with weak swimming abilities and short pelagic duration, while the reptant juvenile/adult stages of the lifecycle are obligate burrow dwellers with limited home ranges. Ninety-one individuals, collected from five scampi fishery management areas around New Zealand, were genotyped using GBS. Using 983 haplotypic genomic loci, three genetically distinct groups were identified: eastern, southern and western. These groups showed significant genetic differentiation with clear source-sink dynamics. The direction of gene flow inferred from the genomic data largely reflected the hydrodynamic particle modelling of ocean current flow around New Zealand. The modelled dispersal during pelagic larval phase highlights the strong connectivity among eastern sampling locations and explains the low genetic differentiation detected among these sampled areas. Our results highlight the value of using a transdisciplinary approach in the inference of connectivity among populations for informing conservation and fishery management.

Haplotype-resolved genome assembly of Bletilla striata (Thunb.) Reichb.f. to elucidate medicinal value.

Bletilla striata, commonly known as baiji, is a species used in traditional Chinese medicine; it is highly regarded for its medicinal applications and therefore has high economic value. Here, we report a high-quality haplotype-resolved genome of B. striata, haplotype A (2.37 Gb, with a scaffold N50 of 146.39 Mb and a contig N50 of 1.65 Mb) and haplotype B (2.43 Gb, with a scaffold N50 of 150.22 Mb and a contig N50 of 1.66 Mb), assembled from high-fidelity (HiFi) reads and chromosome conformation capture (Hi-C) reads. We find evidence that B. striata has undergone two whole-genome duplication (WGD) events: an ancient WGD event shared by most monocots and a recent WGD event unique to all orchids. We also reconstructed the ancestral orchid karyotype (AOK) of 18 ancient chromosomes and the evolutionary trajectories of 16 modern B. striata chromosomes. Comparative genomic analysis suggests that the expanded gene families of B. striata might play important roles in secondary metabolite biosynthesis and environmental adaptation. By combining genomic and transcriptomic data, we identified the 10 core members from nine gene families that were probably involved in B. striata polysaccharide (BSP) biosynthesis. Based on virus-induced gene silencing (VIGS) and yeast two-hybrid experiments, we present an MYB transcription factor (TF), BsMYB2, that can regulate BSP biosynthesis by directly interacting with eight key BSP-related genes: sacA1, HK1, scrK1, scrK2, GPI1, manA1, GMPP1 and UGP2_1. Our study will enhance the understanding of orchid evolution and accelerate the molecular-assisted breeding of B. striata for improving traits of medicinal value.

Natural allelic variation of GmST05 controlling seed size and quality in soybean.

Seed size is one of the most important agronomic traits determining the yield of crops. Cloning the key genes controlling seed size and pyramiding their elite alleles will facilitate yield improvement. To date, few genes controlling seed size have been identified in soybean, a major crop that provides half of the plant oil and one quarter of the plant protein globally. Here, through a genome-wide association study of over 1800 soybean accessions, we determined that natural allelic variation at GmST05 (Seed Thickness 05) predominantly controlled seed thickness and size in soybean germplasm. Further analyses suggested that the two major haplotypes of GmST05 differed significantly at the transcriptional level. Transgenic experiments demonstrated that GmST05 positively regulated seed size and influenced oil and protein contents, possibly by regulating the transcription of GmSWEET10a. Population genetic diversity analysis suggested that allelic variations of GmST05 were selected during geographical differentiation but have not been fixed. In summary, natural variation in GmST05 determines transcription levels and influences seed size and quality in soybean, making it an important gene resource for soybean molecular breeding.