Identification of candidate genes associated with peanut pod length by combined analysis of QTL-seq and RNA-seq.

Pod length (PL) is one of the major traits determining pod size and yield of peanut. Discovering the quantitative trait loci (QTL) and identifying candidate genes associated with PL are essential for breeding high-yield peanut. In this study, quantitative trait loci sequencing (QTL-seq) was performed using the F2 population constructed by a short-pod variety Tifrunner (Tif) and a long-pod line Lps, and a 0.77 Mb genomic region on chromosome 07 was identified as the candidate region for PL. Then, the candidate region was narrowed to a 265.93 kb region by traditional QTL approach. RNA-seq analysis showed that there were four differentially expressed genes (DEGs) in the candidate region, among which Arahy.PF2L6F (AhCDC48) and Arahy.P4LK2T (AhTAA1) were speculated to be PL-related candidate genes. These results were informative for the elucidation of the underlying regulatory mechanism in peanut pod length and would facilitate further identification of valuable target genes.

Deepening genomic sequences of 1081 Gossypium hirsutum accessions reveals novel SNPs and haplotypes relevant for practical breeding utility.

Fiber quality is a major breeding goal in cotton, but phenotypically direct selection is often hindered. In this study, we identified fiber quality and yield related loci using GWAS based on 2.97 million SNPs obtained from 10.65× resequencing data of 1081 accessions. The results showed that 585 novel fiber loci, including two novel stable SNP peaks associated with fiber length on chromosomes At12 and Dt05 and one novel genome regions linked with fiber strength on chromosome Dt12 were identified. Furthermore, by means of gene expression analysis, GhM_A12G0090, GhM_D05G1692, GhM_D12G3135 were identified and GhM_D11G2208 function was identified in Arabidopsis. Additionally, 14 consistent and stable superior haplotypes were identified, and 25 accessions were detected as possessing these 14 superior haplotype in breeding. This study providing fundamental insight relevant to identification of genes associated with fiber quality and yield will enhance future efforts toward improvement of upland cotton.

Multi-environment association study highlights candidate genes for robust agronomic quantitative trait loci in a novel worldwide Capsicum core collection.

Investigating crop diversity through genome-wide association studies (GWAS) on core collections helps in deciphering the genetic determinants of complex quantitative traits. Using the G2P-SOL project world collection of 10 038 wild and cultivated Capsicum accessions from 10 major genebanks, we assembled a core collection of 423 accessions representing the known genetic diversity. Since complex traits are often highly dependent upon environmental variables and genotype-by-environment (G × E) interactions, multi-environment GWAS with a 10 195-marker genotypic matrix were conducted on a highly diverse subset of 350 Capsicum annuum accessions, extensively phenotyped in up to six independent trials from five climatically differing countries. Environment-specific and multi-environment quantitative trait loci (QTLs) were detected for 23 diverse agronomic traits. We identified 97 candidate genes potentially implicated in 53 of the most robust and high-confidence QTLs for fruit flavor, color, size, and shape traits, and for plant productivity, vigor, and earliness traits. Investigating the genetic architecture of agronomic traits in this way will assist the development of genetic markers and pave the way for marker-assisted selection. The G2P-SOL pepper core collection will be available upon request as a unique and universal resource for further exploitation in future gene discovery and marker-assisted breeding efforts by the pepper community.

Genome-wide identification of quantitative trait loci and candidate genes for seven carcass traits in a four-way intercross porcine population.

BACKGROUND: Carcass traits are essential economic traits in the commercial pig industry. However, the genetic mechanism of carcass traits is still unclear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to study seven carcass traits on 223 four-way intercross pigs, including dressing percentage (DP), number of ribs (RIB), skin thinkness (ST), carcass straight length (CSL), carcass diagonal length (CDL), loin eye width (LEW), and loin eye thickness (LET). RESULTS: A total of 227,921 high-quality single nucleotide polymorphisms (SNPs) were detected to perform GWAS. A total of 30 SNPs were identified for seven carcass traits using the mixed linear model (MLM) (p < 1.0 × 10- 5), of which 9 SNPs were located in previously reported quantitative trait loci (QTL) regions. The phenotypic variation explained (PVE) by the significant SNPs was from 2.43 to 16.32%. Furthermore, 11 candidate genes (LYPLAL1, EPC1, MATN2, ZFAT, ZBTB10, ZNF704, INHBA, SMYD3, PAK1, SPTBN2, and ACTN3) were found for carcass traits in pigs. CONCLUSIONS: The GWAS results will improve our understanding of the genetic basis of carcass traits. We hypothesized that the candidate genes associated with these discovered SNPs would offer a biological basis for enhancing the carcass quality of pigs in swine breeding.

Genetic assessment and candidate genes identification for breed-specific characteristics of Qingyuan partridge chicken based on runs of homozygosity.

BACKGROUND: Several core breeding and supporting lines of the Qingyuan partridge chicken, a representative local chicken breed in China, have been developed over 20 years. Consequently, its economic traits related to growth and reproduction have been significantly improved by breeding selection and commercial utilization, but some characteristic traits, such as partridge feathers, high meat quality and sufficient flavor, have always been retained. However, effective methods for genetic assessment and functional gene exploration of similar trait groups are lacking. The presence of identical haplotype fragments transmitted from parent to offspring results in runs of homozygosity (ROH), which offer an efficient solution. In this study, genomes of 134 Qingyuan partridge chickens representing two breeding populations and one preserved population were re-sequenced to evaluate the genetic diversity and explore functional genes by analyzing the diversity, distribution, and frequency of ROH. RESULTS: The results showed a low level of genomic linkage and degree of inbreeding within both the bred and preserved populations, suggesting abundant genetic diversity and an adequate genetic potential of the Qingyuan partridge chicken. Throughout the long-term selection process, 21 genes, including GLI3, ANO5, BLVRA, EFNB2, SLC5A12, and SVIP, associated with breed-specific characteristics were accumulated within three ROH islands, whereas another 21 genes associated with growth traits including IRX1, IRX2, EGFR, TPK1, NOVA1, BDNF and so on were accumulated within five ROH islands. CONCLUSIONS: These findings provide new insights into the genetic assessment and identification of genes with breed-specific and selective characteristics, offering a solid genetic basis for breeding and protection of Qingyuan partridge chickens.

Genome-wide association study for cashmere traits in Inner Mongolia cashmere goat population reveals new candidate genes and haplotypes.

BACKGROUND: The cashmere goat industry is one of the main pillars of animal husbandry in Inner Mongolia Autonomous Region, and plays an irreplaceable role in local economic development. With the change in feeding methods and environment, the cashmere produced by Inner Mongolia cashmere goats shows a tendency of coarser, and the cashmere yield can not meet the consumption demand of people. However, the genetic basis behind these changes is not fully understood. We measured cashmere traits, including cashmere yield (CY), cashmere diameter (CD), cashmere thickness (CT), and fleece length (FL) traits for four consecutive years, and utilized Genome-wide association study of four cashmere traits in Inner Mongolia cashmere goats was carried out using new genomics tools to infer genomic regions and functional loci associated with cashmere traits and to construct haplotypes that significantly affect cashmere traits. RESULTS: We estimated the genetic parameters of cashmere traits in Inner Mongolia cashmere goats. The heritability of cashmere yield, cashmere diameter, and fleece length traits of Inner Mongolia cashmere goats were 0.229, 0.359, and 0.250, which belonged to the medium heritability traits (0.2 ~ 0.4). The cashmere thickness trait has a low heritability of 0.053. We detected 151 genome-wide significantly associated SNPs with four cashmere traits on different chromosomes, which were very close to the chromosomes of 392 genes (located within the gene or within ± 500 kb). Notch3, BMPR1B, and CCNA2 have direct functional associations with fibroblasts and follicle stem cells, which play important roles in hair follicle growth and development. Based on GO functional annotation and KEGG enrichment analysis, potential candidate genes were associated with pathways of hair follicle genesis and development (Notch, P13K-Akt, TGF-beta, Cell cycle, Wnt, MAPK). We calculated the effective allele number of the Inner Mongolia cashmere goat population to be 1.109-1.998, the dominant genotypes of most SNPs were wild-type, the polymorphic information content of 57 SNPs were low polymorphism (0 < PIC < 0.25), and the polymorphic information content of 79 SNPs were moderate polymorphism (0.25 < PIC < 0.50). We analyzed the association of SNPs with phenotypes and found that the homozygous mutant type of SNP1 and SNP3 was associated with the highest cashmere yield, the heterozygous mutant type of SNP30 was associated with the lowest cashmere thickness, the wild type of SNP76, SNP77, SNP78, SNP80, and SNP81 was associated with the highest cashmere thickness, and the wild type type of SNP137 was associated with the highest fleece length. 21 haplotype blocks and 68 haplotype combinations were constructed. Haplotypes A2A2, B2B2, C2C2, and D4D4 were associated with increased cashmere yield, haplotypes E2E2, F1F1, G5G5, and G1G5 were associated with decreased cashmere fineness, haplotypes H2H2 was associated with increased cashmere thickness, haplotypes I1I1, I1I2, J1J4, L5L3, N3N2, N3N3, O2O1, P2P2, and Q3Q3 were associated with increased cashmere length. We verified the polymorphism of 8 SNPs by KASP, and found that chr7_g.102631194A > G, chr10_g.82715068 T > C, chr1_g.124483769C > T, chr24_g.12811352C > T, chr6_g.114111249A > G, and chr6_g.115606026 T > C were significantly genotyped in verified populations (P < 0.05). CONCLUSIONS: In conclusion, the genetic effect of single SNP on phenotypes is small, and SNPs are more inclined to be inherited as a whole. By constructing haplotypes from SNPs that are significantly associated with cashmere traits, it will help to reveal the complex and potential causal variations in cashmere traits of Inner Mongolia cashmere goats. This will be a valuable resource for genomics and breeding of the cashmere goat.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

A study of genetic variants associated with skin traits in the Vietnamese population.

BACKGROUND: Most skin-related traits have been studied in Caucasian genetic backgrounds. A comprehensive study on skin-associated genetic effects on underrepresented populations such as Vietnam is needed to fill the gaps in the field. OBJECTIVES: We aimed to develop a computational pipeline to predict the effect of genetic factors on skin traits using public data (GWAS catalogs and whole-genome sequencing (WGS) data from the 1000 Genomes Project-1KGP) and in-house Vietnamese data (WGS and genotyping by SNP array). Also, we compared the genetic predispositions of 25 skin-related traits of Vietnamese population to others to acquire population-specific insights regarding skin health. METHODS: Vietnamese cohorts of whole-genome sequencing (WGS) of 1008 healthy individuals for the reference and 96 genotyping samples (which do not have any skin cutaneous issues) by Infinium Asian Screening Array-24 v1.0 BeadChip were employed to predict skin-associated genetic variants of 25 skin-related and micronutrient requirement traits in population analysis and correlation analysis. Simultaneously, we compared the landscape of cutaneous issues of Vietnamese people with other populations by assessing their genetic profiles. RESULTS: The skin-related genetic profile of Vietnamese cohorts was similar at most to East Asian cohorts (JPT: Fst = 0.036, CHB: Fst = 0.031, CHS: Fst = 0.027, CDX: Fst = 0.025) in the population study. In addition, we identified pairs of skin traits at high risk of frequent co-occurrence (such as skin aging and wrinkles (r = 0.45, p = 1.50e-5) or collagen degradation and moisturizing (r = 0.35, p = 1.1e-3)). CONCLUSION: This is the first investigation in Vietnam to explore genetic variants of facial skin. These findings could improve inadequate skin-related genetic diversity in the currently published database.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur.

BACKGROUND: Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change. RESULTS: Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change. CONCLUSIONS: The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

Genome-wide association study of agronomic traits related to nitrogen use efficiency in Henan wheat.

BACKGROUND: Nitrogen use efficiency (NUE) is closely related to crop yield and nitrogen fertilizer application rate. Although NUE is susceptible to environments, quantitative trait nucleotides (QTNs) for NUE in wheat germplasm populations have been rarely reported in genome-wide associated study. RESULTS: In this study, 244 wheat accessions were phenotyped by three NUE-related traits in three environments and genotyped by 203,224 SNPs. All the phenotypes for each trait were used to associate with all the genotypes of these SNP markers for identifying QTNs and QTN-by-environment interactions via 3VmrMLM. Among 279 QTNs and one QTN-by-environment interaction for low nitrogen tolerance, 33 were stably identified, especially, one large QTN (r2 > 10%), qPHR3A.2, was newly identified for plant height ratio in one environment and multi-environment joint analysis. Among 52 genes around qPHR3A.2, four genes (TraesCS3A01G101900, TraesCS3A01G102200, TraesCS3A01G104100, and TraesCS3A01G105400) were found to be differentially expressed in low-nitrogen-tolerant wheat genotypes, while TaCLH2 (TraesCS3A01G101900) was putatively involved in porphyrin metabolism in KEGG enrichment analyses. CONCLUSIONS: This study identified valuable candidate gene for low-N-tolerant wheat breeding and provides new insights into the genetic basis of low N tolerance in wheat.

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice.

BACKGROUND: Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice. RESULTS: In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation. CONCLUSION: This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Exploring the gene expression network involved in the heat stress response of a thermotolerant tomato genotype.

BACKGROUND: The increase in temperatures due to the current climate change dramatically affects crop cultivation, resulting in yield losses and altered fruit quality. Tomato is one of the most extensively grown and consumed horticultural products, and although it can withstand a wide range of climatic conditions, heat stress can affect plant growth and development specially on the reproductive stage, severely influencing the final yield. In the present work, the heat stress response mechanisms of one thermotolerant genotype (E42) were investigated by exploring its regulatory gene network. This was achieved through a promoter analysis based on the identification of the heat stress elements (HSEs) mapping in the promoters, combined with a gene co-expression network analysis aimed at identifying interactions among heat-related genes. RESULTS: Results highlighted 82 genes presenting HSEs in the promoter and belonging to one of the 52 gene networks obtained by the GCN analysis; 61 of these also interact with heat shock factors (Hsfs). Finally, a list of 13 candidate genes including two Hsfs, nine heat shock proteins (Hsps) and two GDSL esterase/lipase (GELPs) were retrieved by focusing on those E42 genes exhibiting HSEs in the promoters, interacting with Hsfs and showing variants, compared to Heinz reference genome, with HIGH and/or MODERATE impact on the translated protein. Among these, the Gene Ontology annotation analysis evidenced that only LeHsp100 (Solyc02g088610) belongs to a network specifically involved in the response to heat stress. CONCLUSIONS: As a whole, the combination of bioinformatic analyses carried out on genomic and trascriptomic data available for tomato, together with polymorphisms detected in HS-related genes of the thermotolerant E42 allowed to determine a subset of candidate genes involved in the HS response in tomato. This study provides a novel approach in the investigation of abiotic stress response mechanisms and further studies will be conducted to validate the role of the highlighted genes.

Diversity of Expression Patterns of Lr34, Lr67, and Candidate Genes towards Lr46 with Analysis of Associated miRNAs in Common Wheat Hybrids in Response to Puccinia triticina Fungus.

Leaf rust caused by Puccinia triticina (Pt) is one of the most dangerous diseases causing significant losses in common wheat crops. In adult plants resistant to rust, a horizontal adult plant resistance (APR) type is observed, which protects the plant against multiple pathogen races and is distinguished by greater persistence under production conditions. Crucial pleiotropic slow-rust genes such as Lr34, Lr46, Lr67, and Lr68, in combination with other genes of lesser influence, continue to increase durable resistance to rust diseases. Based on our previous results, we selected four candidate genes for Lr46 out of ten candidates and analysed them for expression before and after inoculation by P. triticina. As part of our study, we also investigated the expression patterns of miRNA molecules complementary to Lr34 and the candidate genes. The aim of the study was to analyse the expression profiles of candidate genes for the Lr46 gene and the Lr34 and Lr67 genes responsible for the differential leaf-rust resistance of hybrid forms of the F1 generation resulting from crosses between the Glenlea cultivar and cultivars from Polish breeding companies. In addition, the expression of five miRNAs (tae-miR9653b, tae-miR5384-3p, tae-miR9780, tae-miR9775 and tae-miR164), complementary to Lr34, and selected candidate genes were analysed using stem-loop RT-PCR and ddPCR. Biotic stress was induced in adult plants by inoculation with Pt fungal spores, under controlled conditions. Plant material was collected before and 6, 12, 24, and 48 h after inoculation (hpi). Differences in expression patterns of Lr34, Lr67, and candidate genes (for Lr46) were analysed by qRT-PCR and showed that gene expression changed at the analysed time points. Identification of molecular markers coupled to the Lr genes studied was also carried out to confirm the presence of these genes in wheat hybrids. qRT-PCR was used to examine the expression levels of the resistance genes. The highest expression of Lr46/Yr29 genes (Lr46-Glu2, Lr46-RLK1, Lr46-RLK2, and Lr46-RLK3) occurred at 12 and 24 hpi, and such expression profiles were obtained for only one candidate gene among the four genes analysed (Lr46-Glu2), indicating that it may be involved in resistance mechanisms of response to Pt infection.

QTL mapping for the flag leaf-related traits using RILs derived from Trititrigia germplasm line SN304 and wheat cultivar Yannong15 in multiple environments.

BACKGROUND: Developing and enriching genetic resources plays important role in the crop improvement. The flag leaf affects plant architecture and contributes to the grain yield of wheat (Triticum aestivum L.). The genetic improvement of flag leaf traits faces problems such as a limited genetic basis. Among the various genetic resources of wheat, Thinopyrum intermedium has been utilized as a valuable resource in genetic improvement due to its disease resistance, large spikes, large leaves, and multiple flowers. In this study, a recombinant inbred line (RIL) population was derived from common wheat Yannong15 and wheat-Th. intermedium introgression line SN304 was used to identify the quantitative trait loci (QTL) for flag leaf-related traits. RESULTS: QTL mapping was performed for flag leaf length (FLL), flag leaf width (FLW) and flag leaf area (FLA). A total of 77 QTLs were detected, and among these, 51 QTLs with positive alleles were contributed by SN304. Fourteen major QTLs for flag leaf traits were detected on chromosomes 2B, 3B, 4B, and 2D. Additionally, 28 QTLs and 8 QTLs for flag leaf-related traits were detected in low-phosphorus and drought environments, respectively. Based on major QTLs of positive alleles from SN304, we identified a pair of double-ended anchor primers mapped on chromosome 2B and amplified a specific band of Th. intermedium in SN304. Moreover, there was a major colocated QTL on chromosome 2B, called QFll/Flw/Fla-2B, which was delimited to a physical interval of approximately 2.9 Mb and contained 20 candidate genes. Through gene sequence and expression analysis, four candidate genes associated with flag leaf formation and growth in the QTL interval were identified. CONCLUSION: These results promote the fine mapping of QFll/Flw/Fla-2B, which have pleiotropic effects, and will facilitate the identification of candidate genes for flag leaf-related traits. Additionally, this work provides a theoretical basis for the application of Th. intermedium in wheat breeding.

High-density genetic map construction and QTL mapping of a zigzag-shaped stem trait in tea plant (Camellia sinensis).

The highly unique zigzag-shaped stem phenotype in tea plants boasts significant ornamental value and is exceptionally rare. To investigate the genetic mechanism behind this trait, we developed BC1 artificial hybrid populations. Our genetic analysis revealed the zigzag-shaped trait as a qualitative trait. Utilizing whole-genome resequencing, we constructed a high-density genetic map from the BC1 population, incorporating 5,250 SNP markers across 15 linkage groups, covering 3,328.51 cM with an average marker interval distance of 0.68 cM. A quantitative trait locus (QTL) for the zigzag-shaped trait was identified on chromosome 4, within a 61.2 to 97.2 Mb range, accounting for a phenotypic variation explained (PVE) value of 13.62%. Within this QTL, six candidate genes were pinpointed. To better understand their roles, we analyzed gene expression in various tissues and individuals with erect and zigzag-shaped stems. The results implicated CsXTH (CSS0035625) and CsCIPK14 (CSS0044366) as potential key contributors to the zigzag-shaped stem formation. These discoveries lay a robust foundation for future functional genetic mapping and tea plant genetic enhancement.

Integrated transcriptomic and WGCNA analyses reveal candidate genes regulating mainly flavonoid biosynthesis in Litsea coreana var. sinensis.

Litsea coreana Levl. var. sinensis (Allen) Yang et P. H. Huang is a popular ethnic herb and beverage plant known for its high flavonoid content, which has been linked to a variety of pharmacological benefits and crucial health-promoting impacts in humans. The progress in understanding the molecular mechanisms of flavonoid accumulation in this plant has been hindered due to the deficiency of genomic and transcriptomic resources. We utilized a combination of Illumina and Oxford Nanopore Technology (ONT) sequencing to generate a de novo hybrid transcriptome assembly. In total, 126,977 unigenes were characterized, out of which 107,977 were successfully annotated in seven public databases. Within the annotated unigenes, 3,781 were categorized into 58 transcription factor families. Furthermore, we investigated the presence of four valuable flavonoids-quercetin-3-O-ß-D-galactoside, quercetin-3-O-ß-D-glucoside, kaempferol-3-O-ß-D-galactoside, and kaempferol-3-O-ß-D-glucoside in 98 samples, using high-performance liquid chromatography. A weighted gene co-expression network analysis identified two co-expression modules, MEpink and MEturquoise, that showed strong positive correlation with flavonoid content. Within these modules, four transcription factor genes (R2R3-MYB, NAC, WD40, and ARF) and four key enzyme-encoding genes (CHI, F3H, PAL, and C4H) emerged as potential hub genes. Among them, the R2R3-MYB (LcsMYB123) as a homologous gene to AtMYB123/TT2, was speculated to play a significant role in flavonol biosynthesis based on phylogenetic analysis. Our findings provided a theoretical foundation for further research into the molecular mechanisms of flavonoid biosynthesis. Additionally, The hybrid transcriptome sequences will serve as a valuable molecular resource for the transcriptional annotation of L. coreana var. sinensis, which will contribute to the improvement of high-flavonoid materials.

Insights into the genetic architecture of Phytophthora capsici root rot resistance in chile pepper (Capsicum spp.) from multi-locus genome-wide association study.

BACKGROUND: Phytophthora root rot, a major constraint in chile pepper production worldwide, is caused by the soil-borne oomycete, Phytophthora capsici. This study aimed to detect significant regions in the Capsicum genome linked to Phytophthora root rot resistance using a panel consisting of 157 Capsicum spp. genotypes. Multi-locus genome wide association study (GWAS) was conducted using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS). Individual plants were separately inoculated with P. capsici isolates, 'PWB-185', 'PWB-186', and '6347', at the 4-8 leaf stage and were scored for disease symptoms up to 14-days post-inoculation. Disease scores were used to calculate disease parameters including disease severity index percentage, percent of resistant plants, area under disease progress curve, and estimated marginal means for each genotype. RESULTS: Most of the genotypes displayed root rot symptoms, whereas five accessions were completely resistant to all the isolates and displayed no symptoms of infection. A total of 55,117 SNP markers derived from GBS were used to perform multi-locus GWAS which identified 330 significant SNP markers associated with disease resistance. Of these, 56 SNP markers distributed across all the 12 chromosomes were common across the isolates, indicating association with more durable resistance. Candidate genes including nucleotide-binding site leucine-rich repeat (NBS-LRR), systemic acquired resistance (SAR8.2), and receptor-like kinase (RLKs), were identified within 0.5 Mb of the associated markers. CONCLUSIONS: Results will be used to improve resistance to Phytophthora root rot in chile pepper by the development of Kompetitive allele-specific markers (KASP®) for marker validation, genomewide selection, and marker-assisted breeding.

New insights into QTNs and potential candidate genes governing rice yield via a multi-model genome-wide association study.

BACKGROUND: Rice (Oryza sativa L.) is one of the globally important staple food crops, and yield-related traits are prerequisites for improved breeding efficiency in rice. Here, we used six different genome-wide association study (GWAS) models for 198 accessions, with 553,229 single nucleotide markers (SNPs) to identify the quantitative trait nucleotides (QTNs) and candidate genes (CGs) governing rice yield. RESULTS: Amongst the 73 different QTNs in total, 24 were co-localized with already reported QTLs or loci in previous mapping studies. We obtained fifteen significant QTNs, pathway analysis revealed 10 potential candidates within 100kb of these QTNs that are predicted to govern plant height, days to flowering, and plot yield in rice. Based on their superior allelic information in 20 elite and 6 inferior genotypes, we found a higher percentage of superior alleles in the elite genotypes in comparison to inferior genotypes. Further, we implemented expression analysis and enrichment analysis enabling the identification of 73 candidate genes and 25 homologues of Arabidopsis, 19 of which might regulate rice yield traits. Of these candidate genes, 40 CGs were found to be enriched in 60 GO terms of the studied traits for instance, positive regulator metabolic process (GO:0010929), intracellular part (GO:0031090), and nucleic acid binding (GO:0090079). Haplotype and phenotypic variation analysis confirmed that LOC_OS09G15770, LOC_OS02G36710 and LOC_OS02G17520 are key candidates associated with rice yield. CONCLUSIONS: Overall, we foresee that the QTNs, putative candidates elucidated in the study could summarize the polygenic regulatory networks controlling rice yield and be useful for breeding high-yielding varieties.

The utility of exome sequencing in diagnosing pediatric neurodevelopmental disorders in a highly consanguineous population.

Exome sequencing (ES) has been utilized in diagnosing children with neurodevelopmental manifestations, this study aimed to investigate the utility of ES in children within a highly consanguineous population that presented with neurodevelopmental complaints. A retrospective chart review was performed for 405 children with neurodevelopmental complaints who have had ES and were evaluated in multiple centers in the United Arab Emirates over a four-year period. Within the cohort of 405 children, consanguinity was reported in 35% (144/405). The primary clinical presentations were developmental delay/cognitive impairment, distinctive facial features, hypotonia, seizures, and weakness. The diagnostic yield was 57% (231/405). Novel variants were identified in 54% (125/231) of positive cases. Within the positive cases, specific treatment was available in 6% (13/231) and copy number variants (CNV) were reported in 3% (8/231) of cases. In eight children, variants in genes that have not yet been linked to human disease that could potentially be the cause of the observed phenotype "candidate genes" were identified. ES was utilized effectively within this cohort with a high diagnostic yield and through the identification of novel gene variants, CNVs, candidate genes and secondary findings as well as the alteration of the treatment plan in cases where treatment was available.

Whole-exome sequencing reveals candidate high-risk susceptibility genes for endometriosis.

BACKGROUND: Endometriosis is a common, chronic disease among fertile-aged women. Disease course may be highly invasive, requiring extensive surgery. The etiology of endometriosis remains elusive, though a high level of heritability is well established. Several low-penetrance predisposing loci have been identified, but high-risk susceptibility remains undetermined. Endometriosis is known to increase the risk of epithelial ovarian cancers, especially of endometrioid and clear cell types. Here, we have analyzed a Finnish family where four women have been diagnosed with surgically verified, severely symptomatic endometriosis and two of the patients also with high-grade serous carcinoma. RESULTS: Whole-exome sequencing revealed three rare candidate predisposing variants segregating with endometriosis. The variants were c.1238C>T, p.(Pro413Leu) in FGFR4, c.5065C>T, p.(Arg1689Trp) in NALCN, and c.2086G>A, p.(Val696Met) in NAV2. The only variant predicted deleterious by in silico tools was the one in FGFR4. Further screening of the variants in 92 Finnish endometriosis and in 19 endometriosis-ovarian cancer patients did not reveal additional carriers. Histopathology, positive p53 immunostaining, and genetic analysis supported the high-grade serous subtype of the two tumors in the family. CONCLUSIONS: Here, we provide FGFR4, NALCN, and NAV2 as novel high-risk candidate genes for familial endometriosis. Our results also support the association of endometriosis with high-grade serous carcinoma. Further studies are required to validate the findings and to reveal the exact pathogenesis mechanisms of endometriosis. Elucidating the genetic background of endometriosis defines the etiology of the disease and provides opportunities for expedited diagnostics and personalized treatments.