Identification Of Candidate Genes Research Articles

Identification of candidate single-nucleotide polymorphisms (SNPs) and genes associated with sugarcane leaf scald disease

Leaf scald, caused by Xanthomonas albilineans, is a severe disease affecting sugarcane worldwide. One of the most practical ways to control it is by developing resistant sugarcane cultivars. It is essential to identify genes associated with the response to leaf scald. A panel of 170 sugarcane genotypes was evaluated for resistance to leaf scald in field conditions for 2 years, followed by a 1-year greenhouse experiment. The phenotypic evaluation data showed a wide continuous distribution, with heritability values ranging from 0.58 to 0.84. Thirteen single nucleotide polymorphisms (SNPs) were identified, significantly associated with leaf scald resistance. Among these, eight were stable across multiple environments and association models. The candidate genes identified and validated based on RNA-seq and qRT-PCR included two genes that encode NB-ARC leucine-rich repeat (LRR)-containing domain disease-resistance protein. These findings provide a basis for developing marker-assisted selection strategies in sugarcane breeding programs.

Identification of candidate genes controlling cold tolerance at the early seedling stage from Dongxiang wild rice by QTL mapping, BSA-Seq and RNA-Seq

BackgroundThe cold tolerance of rice is closely related to its production and geographic distribution. The identification of cold tolerance-related genes is of important significance for developing cold-tolerant rice. Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) is well-adapted to the cold climate of northernmost-latitude habitats ever found in the world, and is one of the most valuable rice germplasms for cold tolerance improvement.ResultsTranscriptome analysis revealed genes differentially expressed between Xieqingzao B (XB; a cold sensitive variety) and 19H19 (derived from an interspecific cross between DXWR and XB) in the room temperature (RT), low temperature (LT), and recovery treatments. The results demonstrated that chloroplast genes might be involved in the regulation of cold tolerance in rice. A high-resolution SNP genetic map was constructed using 120 BC5F2 lines derived from a cross between 19H19 and XB based on the genotyping-by-sequencing (GBS) technique. Two quantitative trait loci (QTLs) for cold tolerance at the early seedling stage (CTS), qCTS12 and qCTS8, were detected. Moreover, a total of 112 candidate genes associated with cold tolerance were identified based on bulked segregant analysis sequencing (BSA-seq). These candidate genes were divided into eight functional categories, and the expression trend of candidate genes related to ‘oxidation-reduction process’ and ‘response to stress’ differed between XB and 19H19 in the RT, LT and recovery treatments. Among these candidate genes, the expression level of LOC_Os12g18729 in 19H19 (related to ‘response to stress’) decreased in the LT treatment but restored and enhanced during the recovery treatment whereas the expression level of LOC_Os12g18729 in XB declined during recovery treatment. Additionally, XB contained a 42-bp deletion in the third exon of LOC_Os12g18729, and the genotype of BC5F2 individuals with a survival percentage (SP) lower than 15% was consistent with that of XB. Weighted gene coexpression network analysis (WGCNA) and modular regulatory network learning with per gene information (MERLIN) algorithm revealed a gene interaction/coexpression network regulating cold tolerance in rice. In the network, differentially expressed genes (DEGs) related to ‘oxidation-reduction process’, ‘response to stress’ and ‘protein phosphorylation’ interacted with LOC_Os12g18729. Moreover, the knockout mutant of LOC_Os12g18729 decreased cold tolerance in early rice seedling stage signifcantly compared with that of wild type.ConclusionsIn general, study of the genetic basis of cold tolerance of rice is important for the development of cold-tolerant rice varieties. In the present study, QTL mapping, BSA-seq and RNA-seq were integrated to identify two CTS QTLs qCTS8 and qCTS12. Furthermore, qRT-PCR, genotype sequencing and knockout analysis indicated that LOC_Os12g18729 could be the candidate gene of qCTS12. These results are expected to further exploration of the genetic mechanism of CTS in rice and improve cold tolerance of cultivated rice by introducing the cold tolerant genes from DXWR through marker-assisted selection.

Genome scans reveal signals of selection associated with pollution in fish populations of Basilichthys microlepidotus, an endemic species of Chile

The Maipo River catchment is one of Chile’s most polluted basins. In recent decades, discharges of untreated sewage and organic matter have caused eutrophication and water quality degradation. We employed the indigenous silverfish species Basilichthys microlepidotus as a model organism to investigate the process of adaptation and selection on genes influenced by pollution. Using variation at single nucleotide polymorphisms (SNPs), we determined the temporal stability of the population structure patterns previously identified in this species by varying SNPs. We also examined local adaptation to pollution-selected genes. Using the genotypes of 7684 loci in 180 individuals, we identified 429 and 700 loci that may be undergoing selection. We detected these loci using the FSTHET and ARLEQUIN outlier detection software, respectively. Both software packages simultaneously identified a total of 250 loci. B. microlepidotus' population structure did not change over time at contaminated or unpolluted sites. In addition, our analysis found: (i) selection of genes associated with pollution, consistent with observations in other organisms; (ii) identification of candidate genes that are functionally linked to the same biological processes, molecular functions and/or cellular components that previously showed differential expression in the same populations; and (iii) a candidate gene with differential expression and a non-synonymous substitution.

Identification of SNP Markers and Candidate Genes Associated with Major Agronomic Traits in Coffea arabica

Genome-wide association studies (GWASs) allow for inferences about the relationships between genomic variants and phenotypic traits in natural or breeding populations. However, few have used this methodology in Coffea arabica. We aimed to identify chromosomal regions with significant associations between SNP markers and agronomic traits in C. arabica. We used a coffee panel consisting of 195 plants derived from 13 families in F2 generations and backcrosses of crosses between leaf rust-susceptible and -resistant genotypes. The plants were phenotyped for 18 agronomic markers and genotyped for 21,211 SNP markers. A GWAS enabled the identification of 110 SNPs with significant associations (p < 0.05) for several agronomic traits in C. arabica: plant height, plagiotropic branch length, number of vegetative nodes, canopy diameter, fruit size, cercosporiosis incidence, and rust incidence. The effects of each SNP marker associated with the traits were analyzed, such that they can be used for molecular marker-assisted selection. For the first time, a GWAS was used for these important agronomic traits in C. arabica, enabling applications in accelerated coffee breeding through marker-assisted selection and ensuring greater efficiency and time reduction. Furthermore, our findings provide preliminary knowledge to further confirm the genomic loci and potential candidate genes contributing to various structural and disease-related traits of C. arabica.

Identification of SNP Markers and Candidate Genes Associated with Major Agronomic Traits in Coffea arabica.

Genome-wide association studies (GWASs) allow for inferences about the relationships between genomic variants and phenotypic traits in natural or breeding populations. However, few have used this methodology in Coffea arabica. We aimed to identify chromosomal regions with significant associations between SNP markers and agronomic traits in C. arabica. We used a coffee panel consisting of 195 plants derived from 13 families in F2 generations and backcrosses of crosses between leaf rust-susceptible and -resistant genotypes. The plants were phenotyped for 18 agronomic markers and genotyped for 21,211 SNP markers. A GWAS enabled the identification of 110 SNPs with significant associations (p < 0.05) for several agronomic traits in C. arabica: plant height, plagiotropic branch length, number of vegetative nodes, canopy diameter, fruit size, cercosporiosis incidence, and rust incidence. The effects of each SNP marker associated with the traits were analyzed, such that they can be used for molecular marker-assisted selection. For the first time, a GWAS was used for these important agronomic traits in C. arabica, enabling applications in accelerated coffee breeding through marker-assisted selection and ensuring greater efficiency and time reduction. Furthermore, our findings provide preliminary knowledge to further confirm the genomic loci and potential candidate genes contributing to various structural and disease-related traits of C. arabica.

Heritability estimates and genome-wide association study of methane emission traits in Nellore cattle.

The objectives of the present study were to estimate the heritability for daily methane emission (CH4) and residual daily methane emission (CH4res) in Nellore cattle, as well as to perform genome-wide association studies (GWAS) to identify genomic regions and candidate genes influencing the genetic variation of CH4 and CH4res. Methane emission phenotypes of 743 Nellore animals belonging to three breeding programs were evaluated. CH4 was measured using the sulfur hexafluoride (SF6) tracer technique (which involves an SF6 permeation tube introduced into the rumen, and an appropriate apparatus on each animal), and CH4res was obtained as the difference between observed CH4 and CH4 adjusted for dry matter intake. A total of 6,252 genotyped individuals were used for genomic analyses. Data were analyzed with a univariate animal model by the single-step GBLUP method using the average information restricted maximum likelihood (AIREML) algorithm. The effects of single nucleotide polymorphisms (SNPs) were obtained using a single-step GWAS approach. Candidate genes were identified based on genomic windows associated with quantitative trait loci (QTLs) related to the two traits. Annotation of QTLs and identification of candidate genes were based on the initial and final coordinates of each genomic window considering the bovine genome ARS-UCD1.2 assembly. Heritability estimates were of moderate to high magnitude, being 0.42 ± 0.09 for CH4 and 0.21 ± 0.09 for CH4res, indicating that these traits will respond rapidly to genetic selection. GWAS revealed 11 and 15 SNPs that were significantly associated (P < 10-6) with genetic variation of CH4 and CH4res, respectively. QTLs associated with feed efficiency, residual feed intake, body weight and height overlapped with significant markers for the traits evaluated. Ten candidate genes were present in the regions of significant SNPs; three were associated with CH4 and seven with CH4res. The identified genes are related to different functions such as modulation of the rumen microbiota, fatty acid production, and lipid metabolism. CH4 and CH4res presented sufficient genetic variation and may respond rapidly to selection. Therefore, these traits can be included in animal breeding programs aimed at reducing enteric methane emissions across generations.

Exploring the genetic architecture of powdery mildew resistance in wheat through QTL meta-analysis.

Powdery mildew (PM), caused by Blumeria graminis f. sp. tritici, poses a significant threat to wheat production, necessitating the development of genetically resistant varieties for long-term control. Therefore, exploring genetic architecture of PM in wheat to uncover important genomic regions is an important area of wheat research. In recent years, the utilization of meta-QTL (MQTL) analysis has gained prominence as an essential tool for unraveling the complex genetic architecture underlying complex quantitative traits. The aim of this research was to conduct a QTL meta-analysis to pinpoint the specific genomic regions in wheat responsible for governing PM resistance. This study integrated 222 QTLs from 33 linkage-based studies using a consensus map with 54,672 markers. The analysis revealed 39 MQTLs, refined to 9 high-confidence MQTLs (hcMQTLs) with confidence intervals of 0.49 to 12.94 cM. The MQTLs had an average physical interval of 41.00 Mb, ranging from 0.000048 Mb to 380.71 Mb per MQTL. Importantly, 18 MQTLs co-localized with known resistance genes like Pm2, Pm3, Pm8, Pm21, Pm38, and Pm41. The study identified 256 gene models within hcMQTLs, providing potential targets for marker-assisted breeding and genomic prediction programs to enhance PM resistance. These MQTLs would serve as a foundation for fine mapping, gene isolation, and functional genomics studies, facilitating a deeper understanding of molecular mechanisms. The identification of candidate genes opens up exciting possibilities for the development of PM-resistant wheat varieties after validation.

Identification and functional characterization of candidate genes involved in biosynthesis of ginsenoside Rg_1

Panax ginseng is a perennial herb with the main active compounds of ginsenosides. Among the reported ginsenosides, ginsenoside Rg_1 not only has a wide range of medicinal functions and abundant content but also is one of the major ginsenoside for the quality evaluation of this herb in the Chinese Pharmacopoeia. The main biosynthesis pathway of ginsenoside Rg_1 in P. ginseng has been clarified, which lays a foundation for the comprehensive and in-depth analysis of the biosynthesis and regulatory mechanism of ginseno-side Rg_1. However, the biosynthesis of ginsenoside Rg_1 is associated with other complex processes involving a variety of regulatory genes and catalyzing enzyme genes, which remain to be studied comprehensively. With the transcriptome data of 344 root samples from 4-year-old P. ginseng plants and their corresponding ginsenoside Rg_1 content obtained in the previous study, this study screened out 217 differentially expressed genes(DEGs) with Rg_1 content changes by DEseq2 analysis in R language. Furthermore, the weighted gene co-expression network analysis(WGCNA) revealed 40 hub genes among the DEGs.Pearsoncorrelation analysis was further perforned to yield 20 candidate genes significantly correlated with ginsenoside Rg_1 content, and these genes were annotated to multiple metabolic processes including primary metabolism and secondary metabolism. Finally, the treatment of P. ginseng adventitious roots with methyl jasmonate indicated that 16 of these genes promoted the biosynthesis of ginsenoside Rg_1 in response to methyl jasmonate induction. Finally, one of the 16 genes was randomly selected to verify the function of the gene by genetic transformation and qRT-PCR and to confirm the rationality of the methodology of this study. The above results lay a foundation for studying the mechanism for regulation on the synthesis of ginsenoside Rg_1 and provide genetic resources for the industrial production of ginsenoside Rg_1.

Identification of candidate gene associated with maize northern leaf blight resistance in a multi-parent population

Key messageQTL mapping combined with genome-wide association studies, revealed a potential candidate gene for resistance to northern leaf blight in the tropical CATETO-related maize line YML226, providing a basis for marker-assisted selection of maize varietiesNorthern leaf blight (NLB) is a foliar disease that can cause severe yield losses in maize. Identifying and utilizing NLB-resistant genes is the most effective way to prevent and control this disease. In this study, five important inbred lines of maize were used as parental lines to construct a multi-parent population for the identification of NLB-resistant loci. QTL mapping and GWAS analysis revealed that QTL qtl_YML226_1, which had the largest phenotypic variance explanation (PVE) of 9.28%, and SNP 5-49,193,921 were co-located in the CATETO-related line YML226. This locus was associated with the candidate gene Zm00001d014471, which encodes a pentatricopeptide repeat (PPR) protein. In the coding region of Zm00001d014471, YML226 had more specific SNPs than the other parental lines. qRT-PCR showed that the relative expressions of Zm00001d014471 in inoculated and uninoculated leaves of YML226 were significantly higher, indicating that the expression of the candidate gene was correlated with NLB resistance. The analysis showed that the higher expression level in YML226 might be caused by SNP mutations. This study identified NLB resistance candidate loci and genes in the tropical maize inbred line YML226 derived from the CATETO germplasm, thereby providing a theoretical basis for using modern marker-assisted breeding techniques to select genetic resources resistant to NLB.

Metabolism-Based Herbicide Resistance to Mesosulfuron-methyl and Identification of Candidate Genes in Bromus japonicus.

The evolved resistance of Bromus japonicus Houtt. to ALS-inhibiting herbicides is well established. Previous studies have primarily focused on target-site resistance; however, non-target-site resistance has not been well characterized. This investigation demonstrated that ALS gene sequencing did not detect any previously known resistance mutations in a mesosulfuron-methyl-resistant (MR) population, and notably, treatment with the P450 monooxygenase (P450) inhibitor malathion markedly heightened susceptibility to mesosulfuron-methyl. Utilizing UPLC-MS/MS analysis confirmed elevated mesosulfuron-methyl metabolism in MR plants. The integration of Isoform Sequencing (Iso-Seq) and RNA Sequencing (RNA-Seq) facilitated the identification of candidate genes associated with non-target sites in a subpopulation with two generations of herbicide selection. Through qRT-PCR analysis, 21 differentially expressed genes were characterized, and among these, 10 genes (comprising three P450s, two glutathione S-transferases, one glycosyltransferase, two ATP-binding cassette transporters, one oxidase, and one hydrolase) exhibited constitutive upregulation in resistant plants. Our findings substantiated that increased herbicide metabolism is a driving force behind mesosulfuron-methyl resistance in this B. japonicus population.

Considerations for reporting variants in novel candidate genes identified during clinical genomic testing.

Since the first novel gene discovery for a Mendelian condition was made via exome sequencing (ES), the rapid increase in the number of genes known to underlie Mendelian conditions coupled with the adoption of exome (and more recently, genome) sequencing by diagnostic testing labs has changed the landscape of genomic testing for rare disease. Specifically, many individuals suspected to have a Mendelian condition are now routinely offered clinical ES. This commonly results in a precise genetic diagnosis but frequently overlooks the identification of novel candidate genes. Such candidates are also less likely to be identified in the absence of large-scale gene discovery research programs. Accordingly, clinical laboratories have both the opportunity, and some might argue a responsibility, to contribute to novel gene discovery which should in turn increase the diagnostic yield for many conditions. However, clinical diagnostic laboratories must necessarily balance priorities for throughput, turnaround time, cost efficiency, clinician preferences, and regulatory constraints, and often do not have the infrastructure or resources to effectively participate in either clinical translational or basic genome science research efforts. For these and other reasons, many laboratories have historically refrained from broadly sharing potentially pathogenic variants in novel genes via networks like Matchmaker Exchange, much less reporting such results to ordering providers. Efforts to report such results are further complicated by a lack of guidelines for clinical reporting and interpretation of variants in novel candidate genes. Nevertheless, there are myriad benefits for many stakeholders, including patients/families, clinicians, researchers, if clinical laboratories systematically and routinely identify, share, and report novel candidate genes. To facilitate this change in practice, we developed criteria for triaging, sharing, and reporting novel candidate genes that are most likely to be promptly validated as underlying a Mendelian condition and translated to use in clinical settings.

Linkage Mapping and Genome-Wide Association Study Identified Two Peanut Late Leaf Spot Resistance Loci, PLLSR-1 and PLLSR-2, Using Nested Association Mapping.

Identification of candidate genes and molecular markers for late leaf spot (LLS) disease resistance in peanut (Arachis hypogaea) has been a focus of molecular breeding for the U.S. industry-funded peanut genome project. Efforts have been hindered by limited mapping resolution due to low levels of genetic recombination and marker density available in traditional biparental mapping populations. To address this, a multi-parental nested association mapping population has been genotyped with the peanut 58K single-nucleotide polymorphism (SNP) array and phenotyped for LLS severity in the field for 3 years. Joint linkage-based quantitative trait locus (QTL) mapping identified nine QTLs for LLS resistance with significant phenotypic variance explained up to 47.7%. A genome-wide association study identified 13 SNPs consistently associated with LLS resistance. Two genomic regions harboring the consistent QTLs and SNPs were identified from 1,336 to 1,520 kb (184 kb) on chromosome B02 and from 1,026.9 to 1,793.2 kb (767 kb) on chromosome B03, designated as peanut LLS resistance loci, PLLSR-1 and PLLSR-2, respectively. PLLSR-1 contains 10 nucleotide-binding site leucine-rich repeat disease resistance genes. A nucleotide-binding site leucine-rich repeat disease resistance gene, Arahy.VKVT6A, was also identified on homoeologous chromosome A02. PLLSR-2 contains five significant SNPs associated with five different genes encoding callose synthase, pollen defective in guidance protein, pentatricopeptide repeat, acyl-activating enzyme, and C2 GRAM domains-containing protein. This study highlights the power of multi-parent populations such as nested association mapping for genetic mapping and marker-trait association studies in peanuts. Validation of these two LLS resistance loci will be needed for marker-assisted breeding.

Transcriptome analysis reveals the molecular mechanism of differences in growth between photoautotrophy and heterotrophy in Chlamydomonas reinhardtii.

The green alga Chlamydomonas reinhardtii can grow photoautotrophically utilizing light and CO2, and heterotrophically utilizing acetate. The physiological and biochemical responses of autotrophy and heterotrophy are different in C. reinhardtii. However, there is no complete understanding of the molecular physiology between autotrophy and heterotrophy. Therefore, we performed biochemical, molecular and transcriptome analysis of C. reinhardtii between autotrophy and heterotrophy. The cell growth characterization demonstrated that heterotrophic cell had enhanced growth rates, and autotrophic cell accumulated more chlorophyll. The transcriptome data showed that a total of 2,970 differentially expressed genes (DEGs) were identified from photoautotrophy 12h (P12h) to heterotrophy 12h (H12h). The DEGs were involved in photosynthesis, the tricarboxylic acid cycle (TCA), pyruvate and oxidative phosphorylation metabolisms. Moreover, the results of qRT-PCR revealed that the relative expression levels of malate dehydrogenase (MDH), succinate dehydrogenase (SDH), ATP synthase (ATPase), and starch synthase (SSS) were increased significantly from P12h and H12h. The protein activity of NAD-malate dehydrogenase (NAD-MDH) and succinate dehydrogenase (SDH) were significantly higher in the H12h group. The above results indicated that the high growth rate observed in heterotrophic cell may be the effects of environmental or genetic regulation of photosynthesis. Therefore, the identification of novel candidate genes in heterotrophy will contribute to the development of microalga strains with higher growth capacity and better performance for biomass production.

Identification of ethanol analgesia quantitative trait loci and candidate genes in BXD recombinant inbred mouse lines.

Alcohol consumption produces acute analgesic effects, and people experiencing pain conditions may drink alcohol to alleviate discomfort. However, tolerance to the analgesic properties of alcohol could prompt escalating consumption and dependence. Both nociception and alcohol-induced analgesia are under significant genetic control. Understanding the genetic architecture of these processes could inform better treatment options for people with pain conditions. This study aims to identify quantitative trait loci (QTL) driving variation in ethanol-induced analgesia across BXD recombinant inbred mouse lines. Male and female mice from 62 BXD strains received ethanol or saline oral gavage for five days and were tested for hot plate (HP) latency at baseline, Day 1, and Day 5. QTL mapping of HP phenotypes identified a significant provisional QTL on chromosome 17 for Day 1 HP latency in mice receiving ethanol. An additional highly suggestive QTL was present on chromosome 9 for the difference in pre- and post-ethanol thermal nociception. Candidate genes within QTL support intervals were provisionally identified using HP phenotypic correlations to transcriptomic database, expression QTL analysis, and other bioinformatics inquiries. The combined behavioral and bioinformatic analyses yielded strong ethanol analgesia candidate genes, specifically Myo6 . Thus, the results of this genetic study of ethanol-induced analgesia in BXD mouse strains may contribute significantly to our understanding of the molecular basis for individual variation in the analgesic response to acute ethanol.

Coloration in Equine: Overview of Candidate Genes Associated with Coat Color Phenotypes.

Variation in coat color among equids has attracted significant interest in genetics and breeding research. The range of colors is primarily determined by the type, concentration, and distribution of melanin pigments, with the balance between eumelanin and pheomelanin influenced by numerous genetic factors. Advances in genomic and sequencing technologies have enabled the identification of several candidate genes that influence coat color, thereby clarifying the genetic basis of these diverse phenotypes. In this review, we concisely categorize coat coloration in horses and donkeys, focusing on the biosynthesis and types of melanin involved in pigmentation. Moreover, we highlight the regulatory roles of some key candidate genes, such as MC1R, TYR, MITF, ASIP, and KIT, in coat color variation. Moreover, the review explores how coat color relates to selective breeding and specific equine diseases, offering valuable insights for developing breeding strategies that enhance both the esthetic and health aspects of equine species.

Identification of candidate genes associated with primary feathers of tianfu nonghua ducks based on Genome-wide association studies

The primary feathers of ducks have important economic value in the poultry industry. This study quantified the primary feather phenotype of Nonghua ducks, including the primary feathers' length, area, distribution of black spots, and feather symmetry. And genome-wide association analysis was used to screen candidate genes that affect the primary feather traits. The genome-wide association study (GWAS) results identified the genetic region related to feather length (FL) on chromosome 2. Through Linkage disequilibrium (LD) analysis, candidate regions (chr2: 115,246,393-116,501,448 bp) were identified and were further annotated to 5 genes: MRS2, GPLD1, ALDH5A1, KIAA0319, and ATP9B. Secondly, candidate regions related to feather black spots were identified on chromosome 21. Through LD analysis, the candidate regions (chr21: 163,552-2,183,853 bp) were screened and further annotated to 47 genes. Among them, STK4, CCN5, and YWHAB genes were related to melanin-related pathways or pigment deposition, which may be key genes affecting the distribution of black spots on feathers. In addition, we also screened 125 genes on multiple chromosomes that may be related to feather symmetry. Among them, significant SNPs on chromosome 1 were further identified as candidate regions (chr1: 142,118,209-142,223,605 bp) through LD analysis and annotated into 2 genes, TGFBRAP1 and LOC113839965. These results reported the genetic basis of the primary feather from multiple phenotypes, and offered valuable insights into the genetic basis for the growth and development of duck feathers and feather color pattern.

Genome-wide analysis reveals genomic diversity and signatures of selection in Qinchuan beef cattle

BackgroundIndigenous Chinese cattle have abundant genetic diversity and a long history of artificial selection, giving local breeds advantages in adaptability, forage tolerance and resistance. The detection of selective sweeps and comparative genome analysis of selected breeds and ancestral populations provide a basis for understanding differences among breeds and for the identification and utilization of candidate genes. We investigated genetic diversity, population structure, and signatures of selection using genome-wide sequencing data for a new breed of Qinchuan cattle (QNC, n = 21), ancestral Qinchuan cattle (QCC, n = 20), and Zaosheng cattle (ZSC, n = 19).ResultsA population structure analysis showed that the ancestry components of QNC and ZSC were similar. In addition, the QNC and ZSC groups showed higher proportions of European taurine ancestry than that of QCC, and this may explain the larger body size of QNC, approaching that of European cattle under long-term domestication and selection. A neighbor-joining tree revealed that QCC individuals were closely related, whereas QNC formed a distinct group. To search for signatures of selection in the QNC genome, we evaluated nucleotide diversity (θπ), the fixation index (FST) and Tajima’s D. Overlapping selective sweeps were enriched for one KEGG pathway, the apelin signaling pathway, and included five candidate genes (MEF2A, SMAD2, CAMK4, RPS6, and PIK3CG). We performed a comprehensive review of genomic variants in QNC, QCC, and ZSC using whole-genome sequencing data. QCC was rich in novel genetic diversity, while diversity in QNC and ZSC cattle was reduced due to strong artificial selection, with divergence from the original cattle.ConclusionsWe identified candidate genes associated with production traits. These results support the success of selective breeding and can guide further breeding and resource conservation of Qinchuan cattle.

Rhytidome- and cork-type barks of holm oak, cork oak and their hybrids highlight processes leading to cork formation

BackgroundThe periderm is basic for land plants due to its protective role during radial growth, which is achieved by the polymers deposited in the cell walls. In most trees, like holm oak, the first periderm is frequently replaced by subsequent internal periderms yielding a heterogeneous outer bark made of a mixture of periderms and phloem tissues, known as rhytidome. Exceptionally, cork oak forms a persistent or long-lived periderm which results in a homogeneous outer bark of thick phellem cell layers known as cork. Cork oak and holm oak distribution ranges overlap to a great extent, and they often share stands, where they can hybridize and produce offspring showing a rhytidome-type bark.ResultsHere we use the outer bark of cork oak, holm oak, and their natural hybrids to analyse the chemical composition, the anatomy and the transcriptome, and further understand the mechanisms underlying periderm development. We also include a unique natural hybrid individual corresponding to a backcross with cork oak that, interestingly, shows a cork-type bark. The inclusion of hybrid samples showing rhytidome-type and cork-type barks is valuable to approach cork and rhytidome development, allowing an accurate identification of candidate genes and processes. The present study underscores that abiotic stress and cell death are enhanced in rhytidome-type barks whereas lipid metabolism and cell cycle are enriched in cork-type barks. Development-related DEGs showing the highest expression, highlight cell division, cell expansion, and cell differentiation as key processes leading to cork or rhytidome-type barks.ConclusionTranscriptome results, in agreement with anatomical and chemical analyses, show that rhytidome and cork-type barks are active in periderm development, and suberin and lignin deposition. Development and cell wall-related DEGs suggest that cell division and expansion are upregulated in cork-type barks whereas cell differentiation is enhanced in rhytidome-type barks.

HLA Variants and Their Association with IgE-Mediated Banana Allergy: A Cross-sectional Study

BackgroundBanana allergy is on the rise in tropical regions. Advances in genomics and candidate gene identification have increased interest in genetic factors in food allergies. However, the genetic basis of IgE-mediated banana allergy is underexplored. ObjectiveTo characterize HLA variants and their association with IgE-mediated banana allergy. MethodsThis cross-sectional study recruited banana-allergic adults, confirmed by allergology tests, with non-allergic individuals as controls. Genomic DNA extraction and sequencing BAM files for HLA typing were conducted. Allele frequency was calculated using the direct counting method, and odds ratio (OR) with 95 % confidence interval (CI) were determined. Fisher's exact or chi-square tests were used to assess associations with Bonferroni's correction for multiple tests. The allele frequency of the Thai population from The Allele Frequency Net Database was used to compute the allele enrichment ratio (ER). ResultsA total of 59 cases and 64 controls were recruited. HLA genotyping indicated potential associations of HLA-B*15:25 (OR 11.872; p-value 0.027), HLA-C*04:03 (OR 7.636; p-value 0.033), and HLA-DQB1*06:09 (OR 11.558; p-value 0.039) with banana allergy. However, after Bonferroni correction, none of these associations reached statistical significance. Comparing allele frequency with the general population from The Allele Frequency Net Database, our ER analysis revealed a higher prevalence in the banana allergy group for B*15:25 (ER 1.849), C*04:03 (ER 1.332), and DQB1*06:09 (ER 6.602) alleles. ConclusionsThis study provides initial genetic insights into banana allergy, suggesting potential links with specific HLA alleles. Despite 12 initially identifying alleles, none were statistically significant after multiple testing correction. Larger studies are needed to detect possible significant correlations.

Comparative transcriptome analysis provides a glance into the regulatory of insulin-like peptide 1 gene in the Pacific white shrimp Litopenaeus vannamei

Insulin-like peptides (ILPs) have been reported to exert significant and diverse regulatory effects through the insulin signaling pathway in invertebrates. However, there is a large research gap in the functional analysis of ILP genes in crustaceans. In this study, we combined RNA interference (RNAi) technology with transcriptome sequencing to investigate the regulatory effects of the ILP1 gene in the Pacific white shrimp Litopenaeus vannamei. After two weeks of continuous RNAi of LvILP1, the shrimp in the dsLvILP1 group exhibited reduced body length gain, weight gain, molting rate, and hemolymph glucose concentration compared to the control dsEGFP group, indicating that knockdown of LvILP1 slowed shrimp growth, molting, and metabolism. Comparative transcriptome analysis of three tissues revealed that knockdown of LvILP1 inhibited muscle growth and development by modulating the mTOR and MAPK signaling pathways, and prevented the renewal of epidermal tissue through the Hippo signaling pathway. Meanwhile, it also suppressed certain growth signals in the hepatopancreas, such as juvenile hormone (JH) synthesis and 20-hydroxyecdysteroid (20E) signal activation. Additionally, it promoted the conversion pathway of some sugars into amino acid metabolism instead of energy metabolism after glycolysis in the hepatopancreas. Interestingly, knockdown of LvILP1 also affected carbohydrate metabolism by regulating the intracellular and extracellular transport of glucose and trehalose, as well as regulating the expression levels of key metabolic genes in three shrimp tissues. This disruption led to a decrease in blood glucose levels and suppressed energy production, ultimately providing less energy for shrimp growth and molting. This study provides valuable insights into the underlying genetic mechanism through which ILP1 regulates glucose metabolism and growth in shrimp. It enhances our understanding of the structure and function of the insulin pathway in crustaceans and serves as an important reference for the identification of candidate genes in molecular genetic breeding of shrimp.