Whole-genome resequencing reveals new mutations in candidate genes for Beichuan-white goat prolificacya

In this study, we evaluated the copy number variation in the genomes of two groups of Beichuan-white goat populations with large differences in litter size by FST method, and identified 1739 genes and 485 missense mutations in the genes subject to positive selection. Through functional enrichment, ITGAV, LRP4, CDH23, TPRN, RYR2 and CELSR1 genes, involved in embryonic morphogenesis, were essential for litter size trait, which received intensive attention. In addition, some mutation sites of these genes have been proposed (ITGAV: c.38C > T; TPRN: c.133A > T, c.1192A > G, c.1250A > C; CELSR1: c.7640T > C), whose allele frequencies were significantly changed in the high fecundity goat group. Besides, we found that new mutations at these sites altered the hydrophilicity and 3D structure of the protein. Candidate genes related to litter size in this study and their missense mutation sites were identified. These candidate genes are helpful to understand the genetic mechanism of fecundity in Beichuan white goat, and have important significance for future goat breeding.

Genome-wide association study identified five quantitative trait loci and two candidate genes for digestive traits in Suhuai pigs

This work aimed to identify markers and candidate genes underlying porcine digestive traits. In total, 331 pigs were genotyped by 80 K Chip data or 50 K Chip data. For apparent neutral detergent fiber digestibility, a total of 19 and 21 candidate single nucleotide polymorphisms (SNP) were respectively identified using a genome-wide efficient mixed-model association algorithm and linkage-disequilibrium adjusted kinship. Among them, three quantitative trait locus (QTL) regions were identified. For apparent acid detergent fiber digestibility, a total of 16 and 17 SNPs were identified by these two methods, respectively. Of these, three QTL regions were also identified. Moreover, two candidate genes (MST1 and LATS1), which are functionally related to intestinal homeostasis and health, were detected near these significant SNPs. Taken together, our results could provide a basis for deeper research on digestive traits in pigs.

Circadian rhythm and redox homeostasis candidate genes showed association with shallow elevation in Norway spruce

The analysis of genetic variation underlying local adaptation in natural populations, together with the response to different external stimuli, is currently a hot topic in forest sciences, with the aim of identifying genetic markers controlling key phenotypic traits of interest for their inclusion in restoration and breeding programs. In Europe, one of the main tree species is Norway spruce (Picea abies (L.) H.Karst.). Using the MassARRAY® platform, 568 trees from North Rhine-Westphalia (Germany) were genotyped with 94 single nucleotide polymorphisms (SNPs) related to circadian and growth rhythms, and to stress response. The association analysis of the selected markers with health status and elevation was performed using three different methods, and those identified by at least two of these were considered as high confidence associated SNPs. While just five markers showed a weak association with health condition, 32 SNPs were correlated with elevation, six of which were considered as high confidence associated SNPs, as indicated by at least two different association methods. Among these genes, thioredoxin and pseudo response regulator 1 (PRR1) are involved in redox homeostasis and ROS detoxification, APETALA2-like 3 (AP2L3), a transcription factor, is involved in seasonal apical growth, and a RPS2-like is a disease resistance gene. The function of some of these genes in controlling light-dependent reactions and metabolic processes suggests signatures of adaptation to local photoperiod and the synchronization of the circadian rhythm. This work provides new insights into the genetic basis of local adaptation over a shallow elevation gradient in Norway spruce.

Screening candidate genes for fruit size based on QTL-seq in Chinese jujube

Introduction

Fruit size is an important economic trait affecting jujube fruit quality, which has always been the focus of marker-assisted breeding of jujube traits. However, despite a large number of studies have been carried out, the mechanism and key genes regulating jujube fruit size are mostly unknown.

Methods

In this study, we used a new analysis method Quantitative Trait Loci sequencing (QTL-seq) (bulked segregant analysis) to screen the parents 'Yuhong' and 'Jiaocheng 5' with significant phenotypic differences and mixed offspring group with extreme traits of large fruit and small fruit, respectively, and, then, DNA mixed pool sequencing was carried out to further shortening the QTL candidate interval for fruit size trait and excavated candidate genes for controlling fruit size.

Results

The candidate intervals related to jujube fruit size were mainly located on chromosomes 1, 5, and 10, and the frequency of chromosome 1 was the highest. Based on the QTL-seq results, the annotation results of ANNOVAR were extracted from 424 SNPs (single-nucleotide polymorphisms) and 164 InDels (insertion-deletion), from which 40 candidate genes were selected, and 37 annotated candidate genes were found in the jujube genome. Four genes (LOC107428904, LOC107415626, LOC125420708, and LOC107418290) that are associated with fruit size growth and development were identified by functional annotation of the genes in NCBI (National Center for Biotechnology Information). The genes can provide a basis for further exploration and identification on genes regulating jujube fruit size.

Discussion

In summary, the data obtained in this study revealed that QTL intervals and candidate genes for fruit size at the genomic level provide valuable resources for future functional studies and jujube breeding.

Identification of genomic regions associated with resistance to gastrointestinal parasites in an indigenous sheep by single- and multiple-locus methods

We investigated the association between 157 SNPs located in 75 candidate genes involved in the immune system and proxy traits for resistance to gastrointestinal nematodes in sheep. A total of 211 lambs from eight flocks were sampled. Nematode eggs per gram were counted and classified as: (i) Strongyles, (ii) Nematodirus spp., (iii) Trichuris spp. and (iv) Marshallagia marshalli. Single- and multiple-locus models were used to test the marker-trait associations. Seven significant SNPs were identified on chromosomes OAR6, 15, 16, and 19. These findings provide insights for breeding nemarode-resistant traits in low-input production systems. General linear model, fixed and random model circulating probability unification, and Bayesian-information and linkage-disequilibrium iteratively nested keyway analyses identified a significant association between the eggs per gram of Strongyles nematodes and a specific variant of the PRLR gene.

Identification of candidate genes affecting the tibia quality in Nonghua duck

The skeleton is a vital organ providing structural support in poultry. Weakness in bone structure can lead to deformities, osteoporosis, cage fatigue, and fractures, resulting in economic losses. Research has substantiated that genetic factors play a significant role in influencing bone quality. The discovery of genetic markers associated with bone quality holds paramount importance for enhancing genetic traits related to the skeletal system in poultry. This study analyzed nine phenotypic indicators of tibia quality in 120-day-old ducks. The phenotypic correlation revealed a high correlation among diameter, Perimeter, and weight (0.69-0.78), and a strong correlation was observed between toughness and breaking strength (0.62). Then, we conducted a genome-wide association analysis of the phenotypic indicators to elucidate the genetic basis of tibial quality in Nonghua ducks. Among the 11 candidate genes that were annotated, TAPT1, BST1, and STIM2 were related to the diameter indicator, ZNF652, IGF2BP1, CASK, and GREB1L were associated with the weight and toughness indicators. RFX8, GLP1R, and DNAAF5 were identified for ash, calcium, and phosphorus content, respectively. Finally, KEGG and GO analysis for annotated genes were performed. STIM2 and BST1 were enriched into the Calcium signalling pathway and Niacin and nicotinamide metabolic pathway, which may be key candidate genes affecting bone quality phenotypes. Gene expression analysis of the candidate genes, such as STIM2, BST1, TAPT1, and CASK showed higher expression levels in bones compared to other tissues. The obtained results can contribute to new insights into tibial quality and provide new genetic biomarkers that can be employed in duck breeding.

The genetic architecture of complete blood counts in lactating Holstein dairy cows

Complete blood counts (CBCs) measure the abundance of individual immune cells, red blood cells, and related measures such as platelets in circulating blood. These measures can indicate the health status of an animal; thus, baseline circulating levels in a healthy animal may be related to the productive life, resilience, and production efficiency of cattle. The objective of this study is to determine the heritability of CBC traits and identify genomic regions that are associated with CBC measurements in lactating Holstein dairy cattle. The heritability of CBCs was estimated using a Bayes C0 model. The study population consisted of 388 cows with genotypes at roughly 75,000 markers and 16 different CBC phenotypes taken at one to three time points (n = 33, 131, and 224 for 1, 2, and 3 time points, respectively). Heritabilities ranged from 0.00 ± 0.00 (red cell distribution width) to 0.68 ± 0.06 (lymphocytes). A total of 96 different 1-Mb windows were identified that explained more than 1% of the genetic variance for at least one CBC trait, with 10 windows explaining more than 1% of the genetic variance for two or more traits. Multiple genes in the identified regions have functions related to immune response, cell differentiation, anemia, and disease. Positional candidate genes include RAD52 motif-containing protein 1 (RDM1), which is correlated with the degree of immune infiltration of immune cells, and C-X-C motif chemokine ligand 12 (CXCL12), which is critically involved in neutrophil bone marrow storage and release regulation and enhances neutrophil migration. Since animal health directly impacts feed intake, understanding the genetics of CBCs may be useful in identifying more disease-resilient and feed-efficient dairy cattle. Identification of genes responsible for variation in CBCs will also help identify the variability in how dairy cattle defend against illness and injury.

Single nucleotide polymorphisms of ANKK1, DDR4, and GRIN2B genes predict behavior in a prospective cohort of Mexican children and adolescents

Numerous studies have established associations between single nucleotide polymorphisms (SNPs) and various behavioral and neurodevelopmental conditions. This study explores the links between SNPs in candidate genes involved in central nervous system (CNS) physiology and their implications for the behavioral and emotional aspects in children and teenagers. A total of 590 participants, aged 7-15 years, from the Early Life Exposures In Mexico To Environmental Toxicants (ELEMENT) cohort study in Mexico City, underwent genotyping for at least one of 15 CNS gene-related SNPs at different timepoints. We employed multiple linear regression models to assess the potential impact of genetic variations on behavioral and cognitive traits, as measured by the Behavioral Assessment System for Children (BASC) and Conners parent rating scales. Significant associations were observed, including the rs1800497 TC genotype (ANKK1) with the Cognitive Problems/Inattention variable (p value = 0.003), the rs1800955 CT genotype (DDR4) with the Emotional Lability Global index variable (p value = 0.01), and the rs10492138 GA and rs7970177 TC genotypes (GRIN2B) with the Depression variable (p values 0.007 and 0.012, respectively). These finds suggest potential genetic profiles associated with "risk" and "protective" behaviors for these SNPs. Our results provide valuable insights into the role of genetic variations in neurobehavior and highlight the need for further research in the early identification and intervention in individuals at risk for these conditions.

Novel loci associated with resistance to downy and powdery mildew in grapevine

Among the main challenges in current viticulture, there is the increasing demand for sustainability in the protection from fungal diseases, such as downy mildew (DM) and powdery mildew (PM). Breeding disease-resistant grapevine varieties is a key strategy for better managing fungicide inputs. This study explores the diversity of grapevine germplasm (cultivated and wild) from Caucasus and neighboring areas to identify genotypes resistant to DM and PM, based on 13 Simple Sequence Repeat (SSR) loci and phenotypical (artificial pathogen inoculation) analysis, and to identify loci associated with DM and PM resistance, via Genome-Wide Association Analysis (GWAS) on Single Nucleotide Polymorphism (SNP) profiles. SSR analysis revealed resistant alleles for 16 out of 88 genotypes. Phenotypic data identified seven DM and 31 PM resistant genotypes. GWAS identified two new loci associated with DM resistance, located on chromosome 15 and 16 (designated as Rpv36 and Rpv37), and two with PM resistance, located on chromosome 6 and 17 (designated as Ren14 and Ren15). The four novel loci identified genomic regions rich in genes related to biotic stress response, such as genes involved in pathogen recognition, signal transduction and resistance response. This study highlights potential candidate genes associated with resistance to DM and PM, providing valuable insights for breeding programs for resistant varieties. To optimize their utilization, further functional characterization studies are recommended.

Genetic diversity and candidate genes for transient waterlogging tolerance in mungbean at the germination and seedling stages

Mungbean [Vigna radiata var. radiata (L.) Wilczek] production in Asia is detrimentally affected by transient soil waterlogging caused by unseasonal and increasingly frequent extreme precipitation events. While mungbean exhibits sensitivity to waterlogging, there has been insufficient exploration of germplasm for waterlogging tolerance, as well as limited investigation into the genetic basis for tolerance to identify valuable loci. This research investigated the diversity of transient waterlogging tolerance in a mini-core germplasm collection of mungbean and identified candidate genes for adaptive traits of interest using genome-wide association studies (GWAS) at two critical stages of growth: germination and seedling stage (i.e., once the first trifoliate leaf had fully-expanded). In a temperature-controlled glasshouse, 292 genotypes were screened for tolerance after (i) 4 days of waterlogging followed by 7 days of recovery at the germination stage and (ii) 8 days of waterlogging followed by 7 days of recovery at the seedling stage. Tolerance was measured against drained controls. GWAS was conducted using 3,522 high-quality DArTseq-derived SNPs, revealing five significant associations with five phenotypic traits indicating improved tolerance. Waterlogging tolerance was positively correlated with the formation of adventitious roots and higher dry masses. FGGY carbohydrate kinase domain-containing protein was identified as a candidate gene for adventitious rooting and mRNA-uncharacterized LOC111241851, Caffeoyl-CoA O-methyltransferase At4g26220 and MORC family CW-type zinc finger protein 3 and zinc finger protein 2B genes for shoot, root, and total dry matter production. Moderate to high broad-sense heritability was exhibited for all phenotypic traits, including seed emergence (81%), adventitious rooting (56%), shoot dry mass (81%), root dry mass (79%) and SPAD chlorophyll content (70%). The heritability estimates, marker-trait associations, and identification of sources of waterlogging tolerant germplasm from this study demonstrate high potential for marker-assisted selection of tolerance traits to accelerate breeding of climate-resilient mungbean varieties.

Genome-wide association study and genomic prediction for resistance to brown planthopper in rice

The brown planthopper (BPH) is the most destructive insect pest that threatens rice production globally. Developing rice varieties incorporating BPH-resistant genes has proven to be an effective control measure against BPH. In this study, we assessed the resistance of a core collection consisting of 502 rice germplasms by evaluating resistance scores, weight gain rates and honeydew excretions. A total of 117 rice varieties (23.31%) exhibited resistance to BPH. Genome-wide association studies (GWAS) were performed on both the entire panel of 502 rice varieties and its subspecies, and 6 loci were significantly associated with resistance scores (P value < 1.0e-8). Within these loci, we identified eight candidate genes encoding receptor-like protein kinase (RLK), nucleotide-binding and leucine-rich repeat (NB-LRR), or LRR proteins. Two loci had not been detected in previous study and were entirely novel. Furthermore, we evaluated the predictive ability of genomic selection for resistance to BPH. The results revealed that the highest prediction accuracy for BPH resistance reached 0.633. As expected, the prediction accuracy increased progressively with an increasing number of SNPs, and a total of 6.7K SNPs displayed comparable accuracy to 268K SNPs. Among various statistical models tested, the random forest model exhibited superior predictive accuracy. Moreover, increasing the size of training population improved prediction accuracy; however, there was no significant difference in prediction accuracy between a training population size of 737 and 1179. Additionally, when there existed close genetic relatedness between the training and validation populations, higher prediction accuracies were observed compared to scenarios when they were genetically distant. These findings provide valuable resistance candidate genes and germplasm resources and are crucial for the application of genomic selection for breeding durable BPH-resistant rice varieties.

Transcriptome analysis of genes involved in the pathogenesis mechanism of potato virus Y in potato cultivar YouJin

Introduction

Potatoes (Solanum tuberosum L.) can be infected by various viruses, but out of all of viruses, the potato virus Y (PVY) is the most detrimental. Research shows that the potato cultivar YouJin is especially vulnerable to PVY and displays severe symptoms, including leaf vein chlorosis, curled leaf margins, large necrotic spots on the leaf blades, and the growth of small new leaves.

Methods

PVY infection in potato cultivar YouJin was confirmed through symptom observation, RT-PCR, and Western blot analysis. Transcriptome sequencing was used to analyze the genes associated with PVY pathogenesis in this cultivar.

Result

Transcriptome analysis of differential genes was conducted in this study to examine the pathogenesis of PVY on YouJin. The results showed that 1,949 genes were differentially regulated, including 853 upregulated genes and 1,096 downregulated genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that carbohydrate synthesis and metabolism pathways were suppressed, and electron transferase and hydrolase activities were reduced. Moreover, there were increased expression levels of protein kinase genes. By focusing on plant-pathogen interaction pathways, six core genes all upregulating the WARK family of transcription factors were obtained. Additionally, a constructed PPI network revealed the identification of key modular differential genes, such as downregulated photosynthesis-related protein genes and upregulated AP2/ERF-ERF transcription factors. Functional network enrichment analysis revealed that PVY infection limited RNA metabolism, glutathionylation, and peroxiredoxin activity while triggering the expression of associated defense genes in YouJin. After analyzing the above, 26 DEGs were screened and 12 DEGs were confirmed via RT-qPCR.

Conclusion

These results establish a hypothetical framework for clarifying the pathogenesis of PVY in the YouJin variety of potatoes, which will help design the disease resistance of YouJin.

High Clinical Exome Sequencing Diagnostic Rates and Novel Phenotypic Expansions for Nonisolated Microphthalmia, Anophthalmia, and Coloboma

Purpose

A molecular diagnosis is only made in a subset of individuals with nonisolated microphthalmia, anophthalmia, and coloboma (MAC). This may be due to underutilization of clinical (whole) exome sequencing (cES) and an incomplete understanding of the genes that cause MAC. The purpose of this study is to determine the efficacy of cES in cases of nonisolated MAC and to identify new MAC phenotypic expansions.

Methods

We determined the efficacy of cES in 189 individuals with nonisolated MAC. We then used cES data, a validated machine learning algorithm, and previously published expression data, case reports, and animal models to determine which candidate genes were most likely to contribute to the development of MAC.

Results

We found the efficacy of cES in nonisolated MAC to be between 32.3% (61/189) and 48.1% (91/189). Most genes affected in our cohort were not among genes currently screened in clinically available ophthalmologic gene panels. A subset of the genes implicated in our cohort had not been clearly associated with MAC. Our analyses revealed sufficient evidence to support low-penetrance MAC phenotypic expansions involving nine of these human disease genes.

Conclusions

We conclude that cES is an effective means of identifying a molecular diagnosis in individuals with nonisolated MAC and may identify putatively damaging variants that would be missed if only a clinically available ophthalmologic gene panel was obtained. Our data also suggest that deleterious variants in BRCA2, BRIP1, KAT6A, KAT6B, NSF, RAC1, SMARCA4, SMC1A, and TUBA1A can contribute to the development of MAC.

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.

Deciphering the genetic landscape of seedling drought stress tolerance in wheat (Triticum aestivum L.) through genome-wide association studies

Wheat is an important cereal crop constrained by several biotic and abiotic stresses including drought stress. Understating the effect of drought stress and the genetic basis of stress tolerance is important to develop drought resilient, high-yielding wheat cultivars. In this study, we investigated the effects of drought stress on seedling characteristics in an association panel consisting of 198 germplasm lines. Our findings revealed that drought stress had a detrimental effect on all the seedling characteristics under investigation with a maximum effect on shoot length (50.94% reduction) and the minimum effect on germination percentage (7.9% reduction). To gain a deeper understanding, we conducted a genome-wide association analysis using 12,511 single nucleotide polymorphisms (SNPs), which led to the identification of 39 marker-trait associations (MTAs). Of these 39 MTAs, 13 were particularly noteworthy as they accounted for >10% of the phenotypic variance with a LOD score >5. These high-confidence MTAs were further utilized to extract 216 candidate gene (CGs) models within 1 Mb regions. Gene annotation and functional characterization identified 83 CGs with functional relevance to drought stress. These genes encoded the WD40 repeat domain, Myb/SANT-like domain, WSD1-like domain, BTB/POZ domain, Protein kinase domain, Cytochrome P450, Leucine-rich repeat domain superfamily, BURP domain, Calmodulin-binding protein60, Ubiquitin-like domain, etc. Findings from this study hold significant promise for wheat breeders as they provide direct assistance in selecting lines harboring favorable alleles for improved drought stress tolerance. Additionally, the identified SNPs and CGs will enable marker-assisted selection of potential genomic regions associated with enhanced drought stress tolerance in wheat.

Revisiting Exome Data Identified Missed Splice Site Variant of the Asparagine Synthetase ( ASNS ) Gene

Next-generation sequencing, such as whole-exome sequencing (WES), is increasingly used in the study of Mendelian disorders, yet many are reported as "negative." Inappropriate variant annotation and filtering steps are reasons for missing the molecular diagnosis. Noncoding variants, including splicing mutations, are examples of variants that can be overlooked. Herein, we report a family of four affected newborns, and all presented with severe congenital microcephaly. Initial research WES analysis identified a damaging homozygous variant in NME1 gene as a possible cause of primary microcephaly phenotype in these patients. However, reanalysis of the exome data uncovered a biallelic splice site variant in asparagine synthetase gene which seems to be the possible cause of the phenotype in these patients. This study highlights the importance of revisiting the exome data and the issue of "negative" exome and the afterward approaches to identify and prove new candidate genes.

Whole-Genome Scanning for Selection Signatures Reveals Candidate Genes Associated with Growth and Tail Length in Sheep

Compared to Chinese indigenous sheep, Western sheep have rapid growth rate, larger physique, and higher meat yield. These excellent Western sheep were introduced into China for crossbreeding to expedite the enhancement of production performance and mutton quality in local breeds. Here, we investigated population genetic structure and genome-wide selection signatures among the Chinese indigenous sheep and the introduced sheep based on whole-genome resequencing data. The PCA, N-J tree and ADMIXTURE results showed significant genetic difference between Chinese indigenous sheep and introduced sheep. The nucleotide diversity (π) and linkage disequilibrium (LD) decay results indicated that the genomic diversity of introduced breeds were lower. Then, Fst & π ratio, XP-EHH, and de-correlated composite of multiple signals (DCMS) methods were used to detect the selection signals. The results showed that we identified important candidate genes related to growth rate and body size in the introduced breeds. Selected genes with stronger selection signatures are associated with growth rate (CRADD), embryonic development (BVES, LIN28B, and WNT11), body size (HMGA2, MSRB3, and PTCH1), muscle development and fat metabolism (MSTN, PDE3A, LGALS12, GGPS1, and SAR1B), wool color (ASIP), and hair development (KRT71, KRT74, and IRF2BP2). Thus, these genes have the potential to serve as candidate genes for enhancing the growth traits of Chinese indigenous sheep. We also identified tail-length trait-related candidate genes (HOXB13, LIN28A, PAX3, and VEGFA) in Chinese long-tailed breeds. Among these genes, HOXB13 is the main candidate gene for sheep tail length phenotype. LIN28A, PAX3, and VEGFA are related to embryonic development and angiogenesis, so these genes may be candidate genes for sheep tail type traits. This study will serve as a foundation for further genetic improvement of Chinese indigenous sheep and as a reference for studies related to growth and development of sheep.

Identification of Selection Signatures and Candidate Genes Related to Environmental Adaptation and Economic Traits in Tibetan Pigs

Tibetan pigs are indigenous to the Qinghai-Tibet Plateau and have been the subject of extensive genomic research primarily focused on their adaptation to high altitudes. However, genetic modifications associated with their response to low-altitude acclimation have not been thoroughly explored. To investigate the genetic basis underlying the low-altitude acclimation of Tibetan pigs, we generated and analyzed genotyping data of Tibetan pigs that inhabit high-altitude regions (average altitude 4000 m) and Tibetan pigs that have inhabited nearby low-altitude regions (average altitude 500 m) for approximately 20 generations. We found that the highland and lowland Tibetan pigs have distinguishable genotype and phenotype variations. We identified 46 and 126 potentially selected SNPs associated with 29 and 56 candidate genes in highland and lowland Tibetan pigs, respectively. Candidate genes in the highland Tibetan pigs were involved in immune response (NFYC and STAT1) and radiation (NABP1), whereas candidate genes in the lowland Tibetan pigs were related to reproduction (ESR2, DMRTA1, and ZNF366), growth and development (NTRK3, FGF18, and MAP1B), and blood pressure regulation (CARTPT). These findings will help to understand the mechanisms of environmental adaptation in Tibetan pigs and offer valuable information into the genetic improvement of Tibetan pigs pertaining to low-altitude acclimation and economic traits.

Dissection of the E8 locus in two early maturing Canadian soybean populations

Soybean [Glycine max (L.) Merr.] is a short-day crop for which breeders want to expand the cultivation range to more northern agro-environments by introgressing alleles involved in early reproductive traits. To do so, we investigated quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) regions comprised within the E8 locus, a large undeciphered region (~7.0 Mbp to 44.5 Mbp) associated with early maturity located on chromosome GM04. We used a combination of two mapping algorithms, (i) inclusive composite interval mapping (ICIM) and (ii) genome-wide composite interval mapping (GCIM), to identify major and minor regions in two soybean populations (QS15524F2:F3 and QS15544RIL) having fixed E1, E2, E3, and E4 alleles. Using this approach, we identified three main QTL regions with high logarithm of the odds (LODs), phenotypic variation explained (PVE), and additive effects for maturity and pod-filling within the E8 region: GM04:16,974,874-17,152,230 (E8-r1); GM04:35,168,111-37,664,017 (E8-r2); and GM04:41,808,599-42,376,237 (E8-r3). Using a five-step variant analysis pipeline, we identified Protein far-red elongated hypocotyl 3 (Glyma.04G124300; E8-r1), E1-like-a (Glyma.04G156400; E8-r2), Light-harvesting chlorophyll-protein complex I subunit A4 (Glyma.04G167900; E8-r3), and Cycling dof factor 3 (Glyma.04G168300; E8-r3) as the most promising candidate genes for these regions. A combinatorial eQTL mapping approach identified significant regulatory interactions for 13 expression traits (e-traits), including Glyma.04G050200 (Early flowering 3/E6 locus), with the E8-r3 region. Four other important QTL regions close to or encompassing major flowering genes were also detected on chromosomes GM07, GM08, and GM16. In GM07:5,256,305-5,404,971, a missense polymorphism was detected in the candidate gene Glyma.07G058200 (Protein suppressor of PHYA-105). These findings demonstrate that the locus known as E8 is regulated by at least three distinct genomic regions, all of which comprise major flowering genes.

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

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

Genome-wide association study of Nelore and Angus heifers with low and high ovarian follicle counts

The number of antral follicles is considered an important fertility trait because animals with a high follicle count (HFC) produce more oocytes and embryos per cycle. Identification of these animals by genetic markers such as single nucleotide polymorphisms (SNPs) can accelerate selection of future generations. The aim of this study was to perform a genome wide association study (GWAS) on Nelore and Angus heifers with HFC and low (LFC) antral follicle counts. The groups HFC and LFC for genotyping were formed based on the average of total follicles (≥ 3 mm) counted in each breed consistently ± standard deviation. A total of 72 Nelore heifers (32 HFC and 40 LFC) and 48 Angus heifers (21 HFC and 27 LFC) were selected and the DNA was extracted from blood and hair bulb. Genotyping was done using the Illumina Bovine HD 770K BeadChip. The GWAS analysis showed 181 and 201 SNPs with genotype/phenotype association (P ≤ 0.01) in Nelore and Angus heifers, respectively. Functional enrichment analysis was performed on candidate genes that were associated with SNPs. A total of 97 genes were associated to the 181 SNPs in the Nelore heifers and the functional analysis identified genes (ROBO1 and SLIT3) in the ROBO-SLIT pathway that can be involved in the control of germ cell migration in the ovary as it is involved in lutheal cell migration and fetal ovary development. In the Angus heifers, 57 genes were associated with the 201 SNPs, highlighting Fribilin 1 (FBN1) gene, involved in regulation of growth factors directly involved in follicle activation and development. In summary, GWAS for Nelore and Angus heifers showed SNPs associated with higher follicle count phenotype. Furthermore, these findings offer valuable insights for the further investigation of potential mechanism involved in follicle formation and development, important for breeding programs for both breeds.

Advancements and Prospects of Genome-Wide Association Studies (GWAS) in Maize

Genome-wide association studies (GWAS) have emerged as a powerful tool for unraveling intricate genotype-phenotype association across various species. Maize (Zea mays L.), renowned for its extensive genetic diversity and rapid linkage disequilibrium (LD), stands as an exemplary candidate for GWAS. In maize, GWAS has made significant advancements by pinpointing numerous genetic loci and potential genes associated with complex traits, including responses to both abiotic and biotic stress. These discoveries hold the promise of enhancing adaptability and yield through effective breeding strategies. Nevertheless, the impact of environmental stress on crop growth and yield is evident in various agronomic traits. Therefore, understanding the complex genetic basis of these traits becomes paramount. This review delves into current and future prospectives aimed at yield, quality, and environmental stress resilience in maize and also addresses the challenges encountered during genomic selection and molecular breeding, all facilitated by the utilization of GWAS. Furthermore, the integration of omics, including genomics, transcriptomics, proteomics, metabolomics, epigenomics, and phenomics has enriched our understanding of intricate traits in maize, thereby enhancing environmental stress tolerance and boosting maize production. Collectively, these insights not only advance our understanding of the genetic mechanism regulating complex traits but also propel the utilization of marker-assisted selection in maize molecular breeding programs, where GWAS plays a pivotal role. Therefore, GWAS provides robust support for delving into the genetic mechanism underlying complex traits in maize and enhancing breeding strategies.

Genome-Wide Association Study Reveals a QTL Region for Ashy Stem Blight Resistance in PRA154 Andean Common Bean

Ashy stem blight (ASB) caused by Macrophomina phaseolina (Tassi) Goidanich affects the common bean (Phaseolus vulgaris L.) at all growing stages. Higher levels of resistance were observed in Andean common beans, but specific resistant quantitative trait loci (QTLs) conferring resistance to this pathogen have not been reported in this gene pool. The objectives of this research were to: (i) conduct a genome-wide association study (GWAS) and QTL mapping for resistance in the Andean breeding line PRA154; and (ii) identify single nucleotide polymorphism (SNP) markers and candidate genes for ASB resistance. Phenotyping was conducted under greenhouse conditions by inoculating the 107 F6:7 recombinant inbred lines (RILs) derived from the cross between the susceptible cultivar 'Verano' and the partial-resistant breeding line PRA154 twice with the M. phaseolina isolate PRI21. Genotyping was performed with 109,040 SNPs distributed across all 11 P. vulgaris chromosomes. A novel major QTL was located between 28,761,668 and 31,263,845 bp, extending 2.5 Mbp on chromosome Pv07, and the highest significant SNP markers were Chr07_28761668_A_G, Chr07_29131720_G_A, and Chr07_31263845_C_T with the highest LOD (more than 10 in most of the cases) and R-squared values, explaining 40% of the phenotypic variance of the PRI21 isolate. The gene Phvul.007G173900 (methylcrotonyl-CoA carboxylase alpha chain and mitochondrial 3-methylcrotonyl-CoA carboxylase 1 [MCCA]) with a size of 10,891 bp, located between 29,131,591 and 29,142,481 bp on Pv07, was identified as one candidate for ASB resistance in PRA154, and it contained Chr07_29131720_G_A. The QTL and genetic marker information could be used to assist common bean breeders to develop germplasm and cultivars with ASB resistance through molecular breeding.

Genome-wide analyses based on a novel donkey 40K liquid chip reveal the gene responsible for coat color diversity in Chinese Dezhou donkey

Dezhou donkey is one of the representative local breeds in China, which is mainly divided into two strains: Sanfen and Wutou. There are obvious differences in coat color between the two strains. The former shows light points around the eyes, around the muzzle and under the belly, while the latter is completely solid black. In this study, genome-wide association analysis was performed for the differences in coat color traits between the Sanfen (n = 97) and Wutou (n = 108) strains using a novel donkey 40K liquid chip developed based on GenoBaits technology, to identify genomic regions and causal genes that could explain this variation. We also used FST and The cross-population composite likelihood ratio test (XPCLR) analyses to explore selected regions related to coat color differences. We identified one significant region on chromosome 15, with the most significant SNP located within the agouti signaling protein (ASIP) gene. At the same time, both FST and XPCLR methods detected the same selected region on chromosome 15, and ASIP was the gene with the strongest signal. ASIP and melanocortin 1 receptor (MC1R) control the ratio of eumelanin to pheomelanin through their protein activity. They are deeply involved in the process of melanosome organation and melanogenesis, thus affecting mammals' coat color variation. We used a range of genome-wide approach to identify the genetic basis of coat color variation in Dezhou donkeys. The results provide a supplement to the color variation study in Chinese donkeys at the genome-wide level, and preliminarily verified the reliability of the Molbreeding Donkey No. 1 40K liquid chip.

Genome-wide genetic architecture for plant maturity and drought tolerance in diploid potatoes

Cultivated potato (Solanum tuberosum) is known to be highly susceptible to drought. With climate change and its frequent episodes of drought, potato growers will face increased challenges to achieving their yield goals. Currently, a high proportion of untapped potato germplasm remains within the diploid potato relatives, and the genetic architecture of the drought tolerance and maturity traits of diploid potatoes is still unknown. As such, a panel of 384 ethyl methanesulfonate-mutagenized diploid potato clones were evaluated for drought tolerance and plant maturity under field conditions. Genome-wide association studies (GWAS) were conducted to dissect the genetic architecture of the traits. The results obtained from the genetic structure analysis of the panel showed five main groups and seven subgroups. Using the Genome Association and Prediction Integrated Tool-mixed linear model GWAS statistical model, 34 and 17 significant quantitative trait nucleotides (QTNs) were found associated with maturity and drought traits, respectively. Chromosome 5 carried most of the QTNs, some of which were also detected by using the restricted two-stage multi-locus multi-allele-GWAS haploblock-based model, and two QTNs were found to be pleiotropic for both maturity and drought traits. Using the non-parametric U-test, one and three QTNs, with 5.13%-7.4% phenotypic variations explained, showed favorable allelic effects that increase the maturity and drought trait values. The quantitaive trait loci (QTLs)/QTNs associated with maturity and drought trait were found co-located in narrow (0.5-1 kb) genomic regions with 56 candidate genes playing roles in plant development and senescence and in abiotic stress responses. A total of 127 potato clones were found to be late maturing and tolerant to drought, while nine were early to moderate-late maturing and tolerant to drought. Taken together, the data show that the studied germplasm panel and the identified candidate genes are prime genetic resources for breeders and biologists in conventional breeding and targeted gene editing as climate adaptation tools.

Sequence Analysis of Six Candidate Genes in Miniature Schnauzers with Primary Hypertriglyceridemia

Miniature Schnauzers are predisposed to primary hypertriglyceridemia (HTG). In this study, we performed whole genome sequencing (WGS) of eight Miniature Schnauzers with primary HTG and screened for risk variants in six HTG candidate genes: LPL, APOC2, APOA5, GPIHBP1, LMF1, and APOE. Variants were filtered to identify those present in ≥2 Miniature Schnauzers with primary HTG and uncommon (<10% allele frequency) in a WGS variant database including 613 dogs from 61 other breeds. Three variants passed filtering: an APOE TATA box deletion, an LMF1 intronic SNP, and a GPIHBP1 missense variant. The APOE and GPIHBP1 variants were genotyped in a cohort of 108 Miniature Schnauzers, including 68 with primary HTG and 40 controls. A multivariable regression model, including age and sex, did not identify an effect of APOE (estimate = 0.18, std. error = 0.14; p = 0.20) or GPIHBP1 genotypes (estimate = -0.26, std. error = 0.42; p = 0.54) on triglyceride concentration. In conclusion, we did not identify a monogenic cause for primary HTG in Miniature Schnauzers in the six genes evaluated. However, if HTG in Miniature Schnauzers is a complex disease resulting from the cumulative effects of multiple variants and environment, the identified variants cannot be ruled out as contributing factors.

Genomic Regions Associated with Resistance to Gastrointestinal Nematode Parasites in Sheep-A Review

Gastrointestinal nematodes (GINs) can be a major constraint and global challenge to the sheep industry. These nematodes infect the small intestine and abomasum of grazing sheep, causing symptoms such as weight loss, diarrhea, hypoproteinemia, and anemia, which can lead to death. The use of anthelmintics to treat infected animals has led to GIN resistance, and excessive use of these drugs has resulted in residue traced in food and the environment. Resistance to GINs can be measured using multiple traits, including fecal egg count (FEC), Faffa Malan Chart scores, hematocrit, packed cell volume, eosinophilia, immunoglobulin (Ig), and dagginess scores. Genetic variation among animals exists, and understanding these differences can help identify genomic regions associated with resistance to GINs in sheep. Genes playing important roles in the immune system were identified in several studies in this review, such as the CFI and MUC15 genes. Results from several studies showed overlapping quantitative trait loci (QTLs) associated with multiple traits measuring resistance to GINs, mainly FEC. The discovery of genomic regions, positional candidate genes, and QTLs associated with resistance to GINs can help increase and accelerate genetic gains in sheep breeding programs and reveal the genetic basis and biological mechanisms underlying this trait.

A Genome-Wide Association Study Reveals the Genetic Mechanisms of Nutrient Accumulation in Spinach

Spinach is a significant source of vitamins, minerals, and antioxidants. These nutrients make it delicious and beneficial for human health. However, the genetic mechanism underlying the accumulation of nutrients in spinach remains unclear. In this study, we analyzed the content of chlorophyll a, chlorophyll b, oxalate, nitrate, crude fiber, soluble sugars, manganese, copper, and iron in 62 different spinach accessions. Additionally, 3,356,182 high-quality, single-nucleotide polymorphisms were found using resequencing and used in a genome-wide association study. A total of 2077 loci were discovered that significantly correlated with the concentrations of the nutritional elements. Data mining identified key genes in these intervals for four traits: chlorophyll, oxalate, soluble sugar, and Fe. Our study provides insights into the genetic architecture of nutrient variation and facilitates spinach breeding for good nutrition.

[Mapping of QTL associated with rice cooking quality and candidate gene analysis]

Excavating the quantitative trait locus (QTL) associated with rice cooking quality, analyzing candidate genes, and improving cooking quality-associated traits of rice varieties by genetic breeding can effectively improve the taste of rice. In this study, we used the indica rice HZ, the japonica rice Nekken2 and 120 recombinant inbred lines (RILs) populations constructed from them as experimental materials to measure the gelatinization temperature (GT), gel consistency (GC) and amylose content (AC) of rice at the maturity stage. We combined the high-density genetic map for QTL mapping. A total of 26 QTLs associated with rice cooking quality (1 QTL associated with GT, 13 QTLs associated with GC, and 12 QTLs associated with AC) were detected, among which the highest likelihood of odd (LOD) value reached 30.24. The expression levels of candidate genes in the localization interval were analyzed by quantitative real-time polymerase chain reaction (qRT-PCR), and it was found that the expression levels of six genes were significantly different from that in parents. It was speculated that the high expression of LOC_Os04g20270 and LOC_Os11g40100 may greatly increase the GC of rice, while the high expression of LOC_Os01g04920 and LOC_Os02g17500 and the low expression of LOC_Os03g02650 and LOC_Os05g25840 may reduce the AC. The results lay a molecular foundation for the cultivation of new high-quality rice varieties, and provide important genetic resources for revealing the molecular regulation mechanism of rice cooking quality.

Genetic Analysis of Partially Resistant and Susceptible Chickpea Cultivars in Response to Ascochyta rabiei Infection

The molecular mechanism involved in chickpea (Cicer arietinum L.) resistance to the necrotrophic fungal pathogen Ascochyta rabiei is not well documented. A. rabiei infection can cause severe damage in chickpea, resulting in significant economic losses. Understanding the resistance mechanism against ascochyta blight can help to define strategies to develop resistant cultivars. In this study, differentially expressed genes from two partially resistant cultivars (CDC Corinne and CDC Luna) and a susceptible cultivar (ICCV 96029) to ascochyta blight were identified in the early stages (24, 48 and 72 h) of A. rabiei infection using RNA-seq. Altogether, 3073 genes were differentially expressed in response to A. rabiei infection across different time points and cultivars. A larger number of differentially expressed genes (DEGs) were found in CDC Corinne and CDC Luna than in ICCV 96029. Various transcription factors including ERF, WRKY, bHLH and MYB were differentially expressed in response to A. rabiei infection. Genes involved in pathogen detection and immune signalings such as receptor-like kinases (RLKs), Leucine-Rich Repeat (LRR)-RLKs, and genes associated with the post-infection defence response were differentially expressed among the cultivars. GO functional enrichment and pathway analysis of the DEGs suggested that the biological processes such as metabolic process, response to stimulus and catalytic activity were overrepresented in both resistant and susceptible chickpea cultivars. The expression patterns of eight randomly selected genes revealed by RNA-seq were confirmed by quantitative PCR (qPCR) analysis. The results provide insights into the complex molecular mechanism of the chickpea defence in response to the A. rabiei infection.

Genome-wide survey reveals the genetic background of Xinjiang Brown cattle in China

Introduction: Xinjiang Brown cattle are a famous dual-purpose (dairy-beef) cultivated breed in China that occupy a pivotal position within the cattle breeding industry in Xinjiang, China. However, little information is available on the genetic background of this breed. To fill this research gap, we conducted a whole-genome screen using specific-locus amplified fragment sequencing to examine the genetic structure and diversity of 130 Xinjiang Brown cattle-grazing type (XBG, traditional type) cattle. Methods: A subsequent joint analysis incorporating two ancestral breeds, specifically 19 Brown Swiss (BS) foreign and nine Kazakh (KZ) Chinese cattle, as well as 20 Xinjiang Brown cattle-housing type (XBH) cattle, was used to explore the genetic background of the Xinjiang Brown cattle. Results: The results showed that, after nearly a century of crossbreeding, XBG cattle formed a single population with a stable genetic performance. The genetic structure, genetic diversity, and selection signature analysis of the two ancestral types showed highly different results compared to that of XBH cattle. Local ancestry inference showed that the average proportions of XGB cattle within the BS and KZ cattle lineages were 37.22% and 62.78%, respectively, whereas the average proportions of XBH cattle within the BS and KZ cattle lineages were 95.14% and 4.86%, respectively. Thus, XGB cattle are more representative of all Xinjiang Brown cattle, in line with their breeding history, which involves crossbreeding. Two complementary approaches, fixation index and mean nucleotide diversity, were used to detect selection signals in the four aforementioned cattle breeds. Finally, the analysis of 26 candidate genes in Xinjiang Brown cattle revealed significant enrichment in 19 Gene Ontology terms, and seven candidate genes were enriched in three pathways related to disease resistance (CDH4, SIRPB1, and SIRPα) and the endocrine system (ADCY5, ABCC8, KCNJ11, and KCNMA1). Finally, development of the core SNPs in XBG cattle yielded 8,379 loci. Conclusion: The results of this study detail the evolutionary process of crossbreeding in Xinjiang Brown cattle and provide guidance for selecting and breeding new strains of this species.

Multi-model genome-wide association studies for appearance quality in rice

Improving the quality of the appearance of rice is critical to meet market acceptance. Mining putative quality-related genes has been geared towards the development of effective breeding approaches for rice. In the present study, two SL-GWAS (CMLM and MLM) and three ML-GWAS (FASTmrEMMA, mrMLM, and FASTmrMLM) genome-wide association studies were conducted in a subset of 3K-RGP consisting of 198 rice accessions with 553,831 SNP markers. A total of 594 SNP markers were identified using the mixed linear model method for grain quality traits. Additionally, 70 quantitative trait nucleotides (QTNs) detected by the ML-GWAS models were strongly associated with grain aroma (AR), head rice recovery (HRR, %), and percentage of grains with chalkiness (PGC, %). Finally, 39 QTNs were identified using single- and multi-locus GWAS methods. Among the 39 reliable QTNs, 20 novel QTNs were identified for the above-mentioned three quality-related traits. Based on annotation and previous studies, four functional candidate genes (LOC_Os01g66110, LOC_Os01g66140, LOC_Os07g44910, and LOC_Os02g14120) were found to influence AR, HRR (%), and PGC (%), which could be utilized in rice breeding to improve grain quality traits.

Uncovering Novel QTLs and Candidate Genes for Salt Tolerance at the Bud Burst Stage in Rice through Genome-Wide Association Study

Salt stress is one of the most important factors limiting rice growth and yield increase. Salt tolerance of rice at the bud burst (STB) stage determines whether germinated seeds can grow normally under salt stress, which is very important for direct seeding. However, reports on quantitative trait loci (QTLs) and candidate genes for STB in rice are very limited. In this study, a natural population of 130 indica and 81 japonica rice accessions was used to identify STB-related QTLs and candidate genes using a genome-wide association study (GWAS). Nine QTLs, including five for relative shoot length (RSL), two for relative root length (RRL), and two for relative root number (RRN), were identified. Five of these STB-related QTLs are located at the same site as the characterized salt tolerance genes, such as OsMDH1, OsSRFP1, and OsCDPK7. However, an important QTL related to RSL, qRSL1-2, has not been previously identified and was detected on chromosome 1. The candidate region for qRSL1-2 was identified by linkage disequilibrium analysis, 18 genes were found to have altered expression levels under salt stress through the RNA-seq database, and 10 of them were found to be highly expressed in the shoot. It was also found that, eight candidate genes (LOC_Os01g62980, LOC_Os01g63190, LOC_Os01g63230, LOC_Os01g63280, LOC_Os01g63400, LOC_Os01g63460, and LOC_Os01g63580) for qRSL1-2 carry different haplotypes between indica and japonica rice, which exactly corresponds to the significant difference in RSL values between indica and japonica rice in this study. Most of the accessions with elite haplotypes were indica rice, which had higher RSL values. These genes with indica-japonica specific haplotypes were identified as candidate genes. Rice accessions with elite haplotypes could be used as important resources for direct seeding. This study also provides new insights into the genetic mechanism of STB.

Identification and Analysis of Genes Related to Testicular Size in 14-Day-Old Piglets

The RNA-Seq technology was used to screen the key genes that affect the early development of the testes of Duroc × Landrace × Yorkshire piglets, to determine the regulatory pathway and provide reference for subsequent reproductive performance research, breeding, and other production practices. This study selected 14-day-old Duroc × Landrace × Yorkshire piglets as the trial animals. Testes from piglets with similar weights and no pathological changes were divided into small testis (ST) and large testis (LT) groups, and the RNA-Seq screening of differentially expressed genes (DEGs) was performed to find candidate genes and regulatory pathways related to early testicular development. The results show that 570 DEGs were found in the ST and LT groups, with 281 upregulated and 289 downregulated. The DEGs were mainly enriched on 47 gene ontology (GO) functional items. The Kyoto encyclopedia of genes and genotypes (KEGG) enrichment analysis found that there were 44 significantly enriched KEGG signal pathways, and the regulation of testicular development mainly focused on the arachidonic acid metabolism, Wnt signaling pathway and GnRH secretion pathways. The PTGES, SFRP1, SPP1, PLA2G4E, KCNJ5, PTGS2, and HCN1 genes were found to be as closely related to the testicular development of these Duroc × Landrace × Yorkshire piglets, and the differential gene expression was consistent with the real-time quantitative reverse transcription PCR (real-time qRT-PCR) validation results. This study was validated by high-throughput sequencing analysis and real-time qRT-PCR, and showed that the PTGES, SFRP1, SPP1, PLA2G4E, KCNJ5, PTGS2, and HCN1 genes may be involved in the regulation of germ cell development, spermatogenesis and semen traits. These should be further studied as candidate genes for early testicular development and reproductive trait regulation in boars.

QTL Detection and Candidate Gene Identification for Eating and Cooking Quality Traits in Rice (Oryza sativa L.) via a Genome-Wide Association Study

The eating and cooking quality (ECQ) directly affects the taste of rice, being closely related to factors such as gelatinization temperature (GT), gel consistency (GC) and amylose content (AC). Mining the quantitative trait loci (QTLs), and gene loci controlling ECQ-related traits is vital. A genome-wide association study on ECQ-related traits was conducted, combining 1.2 million single nucleotide polymorphisms (SNPs) with the phenotypic data of 173 rice accessions. Two QTLs for GT, one for GC and five for AC were identified, of which two were found in previously reported genes, and six were newly found. There were 28 positional candidate genes in the region of qAC11. Based on a linkage disequilibrium (LD) analysis, three candidate genes were screened within the LD region associated with AC. There were significant differences between the haplotypes of LOC_Os11g10170, but no significant differences were found for the other two genes. The qRT-PCR results showed that the gene expression levels in the accessions with high ACs were significantly larger than those in the accessions with low ACs at 35d and 42d after flowering. Hap 2 and Hap 3 of LOC_Os11g10170 reduced the AC by 13.09% and 10.77%, respectively. These results provide a theoretical and material basis for improving the ECQ of rice.

Ascertaining the genetic background of the Celtic-Iberian pig strain: A signatures of selection approach

Celtic-Iberian pig breeds were majority in Spain and Portugal until the first half of the 20th century. In the 1990s, they were nearly extinct as a result of the introduction of foreign improved pig breeds. Despite its historical importance, the genetic background of the Celtic-Iberian pig strain is poorly documented. In this study, we have identified genomic regions that might contain signatures of selection peculiar of the Celtic-Iberian genetic lineage. A total of 153 DNA samples of Celtic-Iberian pigs (Spanish Gochu Asturcelta and Portuguese Bísara breeds), Iberian pigs (Spanish Iberian and Portuguese Alentejano breeds), Cinta Senese pig, Korean local pig and Cosmopolitan pig (Hampshire, Landrace and Large White individuals) were analysed. A pairwise-comparison approach was applied: the Gochu Asturcelta and the Bísara samples as test populations and the five other pig populations as reference populations. Three different statistics (XP-EHH, FST and ΔDAF) were computed on each comparison. Strict criteria were used to identify selection sweeps in order to reduce the noise brought on by the Gochu Asturcelta and Bísara breeds' severe population bottlenecks. Within test population, SNPs used to construct potential candidate genomic areas under selection were only considered if they were identified in four of ten two-by-two pairwise comparisons and in at least two of three statistics. Genomic regions under selection constructed within test population were subsequently overlapped to construct candidate regions under selection putatively unique to the Celtic-Iberian pig strain. These genomic regions were finally used for enrichment analyses. A total of 39 candidate regions, mainly located on SSC5 and SSC9 and covering 3130.5 kb, were identified and could be considered representative of the ancient genomic background of the Celtic-Iberian strain. Enrichment analysis allowed to identify a total of seven candidate genes (NOL12, LGALS1, PDXP, SH3BP1, GGA1, WIF1, and LYPD6). Other studies reported that the WIF1 gene is associated with ear size, one of the characteristic traits of the Celtic-Iberian pig strain. The function of the other candidate genes could be related to reproduction, adaptation and immunity traits, indirectly fitting with the rusticity of a non-improved pig strain traditionally exploited under semi-extensive conditions.

Transcriptome Profiling Identifies Differentially Expressed Genes in Skeletal Muscle Development in Native Chinese Ducks

China boasts a rich diversity of indigenous duck species, some of which exhibit desirable economic traits. Here, we generated transcriptome sequencing datasets of breast muscle tissue samples from 1D of four groups: Pekin duck pure breeding group (P), Jinling White duck breeding group (J), P ♂ × J ♀ orthogonal group (PJ) and J ♂ × P ♀ reciprocal-cross group (JP) (n = 3), chosen based on the distinctive characteristics of duck muscle development during the embryonic period. We identified 5053 differentially expressed genes (DEGs) among the four groups. Network prediction analysis showed that ribosome and oxidative phosphorylation-related genes were the most enriched, and muscular protein-related genes were found in the 14-day-old embryonic group. We found that previously characterized functional genes, such as FN1, AGRN, ADNAMST3, APOB and FGF9, were potentially involved in muscle development in 14-day-old embryos. Functional enrichment analysis suggested that genes that participated in molecular function and cell component and key signaling pathways (e.g., hippo, ribosome, oxidative phosphorylation) were significantly enriched in the development of skeletal muscle at 14 days of embryonic age. These results indicate a possible role of muscle metabolism and myoglobin synthesis in skeletal muscle development in both duck parents and hybrids.

QTL Mapping of Trichome Traits and Analysis of Candidate Genes in Leaves of Wheat (Triticum aestivum L.)

Trichome plays an important role in heat dissipation, cold resistance, water absorption, protection of leaves from mechanical damage, and direct exposure to ultraviolet rays. It also plays an important role in the photosynthesis, transpiration, and respiration of plants. However, the genetic basis of trichome traits is not fully understood in wheat. In this study, wheat DH population (Hanxuan 10 × Lumai 14) was used to map quantitative trait loci (QTL) for trichome traits in different parts of flag leaf at 10 days after anther with growing in Zhao County, Hebei Province, and Taigu County, Shanxi Province, respectively. The results showed that trichome density (TD) was leaf center > leaf tip > leaf base and near vein > middle > edge, respectively, in both environments. The trichome length (TL) was leaf tip > leaf center > leaf base and edge > middle > near vein. Significant phenotypic positive correlations were observed between the trichome-related traits of different parts. A total of 83 QTLs for trichome-related traits were mapped onto 18 chromosomes, and each one accounted for 2.41 to 27.99% of the phenotypic variations. Two QTL hotspots were detected in two marker intervals: AX-95232910~AX-95658735 on 3A and AX-94850949~AX-109507404 on 7D. Six possible candidate genes (TraesCS3A02G406000, TraesCS3A02G414900, TraesCS3A02G440900, TraesCS7D02G145200, TraesCS7D02G149200, and TraesCS7D02G152400) for trichome-related traits of wheat leaves were screened out according to their predicted expression levels in wheat leaves. The expression of these genes may be induced by a variety of abiotic stresses. The results provide the basis for further validation and functional characterization of the candidate genes.

Molecular Genetic Aspects of Sporadic Multiglandular Primary Hyperparathyroidism

Multiglandular primary hyperparathyroidism (MGD) represents a rare form of primary hyperparathyroidism (PHPT). MGD is associated with hereditary PHPT, but the sporadic MGD is more common and affects a similar patient profile as single gland parathyroid disease (SGD). The distinction between SGD and MGD is of great clinical importance, especially for the strategy of parathyroidectomy. Based on the limited knowledge available, MGD is likely to be a genetically heterogeneous disease resulting from the interaction of germline and somatic DNA mutations together with epigenetic alterations. Furthermore, these events may combine and occur independently in parathyroid tumors within the same individual with MGD. Gene expression profiling has shown that SGD and MGD may represent distinct entities in parathyroid tumorigenesis. We are waiting for studies to analyze exactly which genes are different in SGD and MGD in order to identify potential biomarkers that can distinguish between the two forms of the disease.

Corrigendum: Unravelling the genetic framework associated with grain quality and yield-related traits in maize (Zea mays L.)

[This corrects the article DOI: 10.3389/fgene.2023.1248697.].

A Single-Step Genome-Wide Association Study for Semen Traits of Egyptian Buffalo Bulls

The present study aimed to contribute to the limited research on buffalo (Bubalus bubalis) semen traits by incorporating genomic data. A total of 8465 ejaculates were collected. The genotyping procedure was conducted using the Axiom® Buffalo Genotyping 90 K array designed by the Affymetrix Expert Design Program. After conducting a quality assessment, we utilized 67,282 SNPs genotyped in 192 animals. We identified several genomic loci explaining high genetic variance by employing single-step genomic evaluation. The aforementioned regions were located on buffalo chromosomes no. 3, 4, 6, 7, 14, 16, 20, 22, and the X-chromosome. The X-chromosome exhibited substantial influence, accounting for 4.18, 4.59, 5.16, 5.19, and 4.31% of the genomic variance for ejaculate volume, mass motility, livability, abnormality, and concentration, respectively. In the examined genomic regions, we identified five novel candidate genes linked to male fertility and spermatogenesis, four in the X-chromosome and one in chromosome no. 16. Additional extensive research with larger sample sizes and datasets is imperative to validate these findings and evaluate their applicability for genomic selection.

Genetic analysis of milk citrate predicted by milk mid-infrared spectra of Holstein cows in early lactation

Milk citrate is regarded as an early biomarker of negative energy balance (NEB) in dairy cows during early lactation and serves as a suitable candidate phenotype for genomic selection due to its wide availability across a large number of cows through milk mid-infrared spectra prediction. However, its genetic background is not well known. Therefore, the objectives of this study were to (1) analyze the genetic parameters of milk citrate; (2) identify genomic regions associated with milk citrate; (3) analyze the functional annotation of candidate genes and quantitative trait loci (QTL) related to milk citrate in Walloon Holstein cows. In total, 134,517 test-day milk citrate phenotypes (mmol/L) collected within the first 50 d in milk (DIM) on 52,198 Holstein cows were used. These milk citrate phenotypes, predicted by milk mid-infrared spectra, were divided into 3 traits according to the first (citrate1), second (citrate2), and third to fifth parity (citrate3+). Genomic information for 566,170 SNPs was available for 4,479 animals. A multiple-trait repeatability model was used to estimate genetic parameters. A single-step genome-wide association study (GWAS) was used to identify candidate genes for citrate and post-GWAS analysis was done to investigate relationship and function of the identified candidate genes. The heritabilities estimated for citrate1, citrate2 and citrate3+ were 0.40, 0.37 and 0.35, respectively. The genetic correlations between the 3 traits ranged from 0.98 to 0.99. The genomic correlations between the 3 traits were also nearly 1.00 across the genomic regions (1 Mb) in the whole genome, which means that citrate can be considered as a single trait in the first 5 parities. In total, 603 significant SNPs located on 3 genomic regions (chromosome7 68.569 - 68.575 Mb, 14 1.31 - 3.05 Mb, and 20 54.00 - 64.28 Mb), were identified to be associated with milk citrate. We identified 89 candidate genes including GPT, ANKH, PPP1R16A and 32 QTL reported in the literature related to the identified significant SNPs. These identified QTL were mainly reported associated with milk fatty acids and metabolic diseases in dairy cows. This study suggests that milk citrate in Holstein cows is highly heritable and has the potential to be used as an early proxy for the NEB of Holstein cows in a breeding objective.

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.

Multitrait meta-analyses identify potential candidate genes for growth-related traits in Holstein heifers

Understanding the underlying pleiotropic relationships among growth and body size traits is important for refining breeding strategies in dairy cattle for optimal body size and growth rate. Therefore, we performed single-trait GWAS for monthly-recorded body weight (BW), hip height, body length, and chest girth from birth to 12 mo of age in Holstein animals, followed by stepwise multiple regression of independent or lowly-linked markers from GWAS loci using conditional and joint association analyses (COJO). Subsequently, we conducted a multitrait meta-analysis to detect pleiotropic markers. Based on the single-trait GWAS, we identified 170 significant SNPs, in which 59 of them remained significant after the COJO analyses. The most significant SNP, located at BTA7:3,676,741, explained 2.93% of the total phenotypic variance for BW6 (BW at 6 mo of age). We identified 17 SNPs with potential pleiotropic effects based on the multitrait meta-analyses, which resulted in 3 additional SNPs in comparison to those detected based on the single-trait GWAS. The identified quantitative trait loci regions overlap with genes known to influence human growth-related traits. According to positional and functional analyses, we proposed HMGA2, HNF4G, MED13L, BHLHE40, FRZB, DMP1, TRIB3, and GATAD2A as important candidate genes influencing the studied traits. The combination of single-trait GWAS and meta-analyses of GWAS results improved the efficiency of detecting associated SNPs, and provided new insights into the genetic mechanisms of growth and development in Holstein cattle.

Candidate signatures of positive selection for environmental adaptation in indigenous African cattle: A review

Environmental adaptation traits of indigenous African cattle are increasingly being investigated to respond to the need for sustainable livestock production in the context of unpredictable climatic changes. Several studies have highlighted genomic regions under positive selection probably associated with adaptation to environmental challenges (e.g. heat stress, trypanosomiasis, tick and tick-borne diseases). However, little attention has focused on pinpointing the candidate causative variant(s) controlling the traits. This review compiled information from 22 studies on signatures of positive selection in indigenous African cattle breeds to identify regions under positive selection. We highlight some key candidate genome regions and genes of relevance to the challenges of living in extreme environments (high temperature, high altitude, high infectious disease prevalence). They include candidate genes involved in biological pathways relating to innate and adaptive immunity (e.g. BoLAs, SPAG11, IL1RL2 and GFI1B), heat stress (e.g. HSPs, SOD1 and PRLH) and hypoxia responses (e.g. BDNF and INPP4A). Notably, the highest numbers of candidate regions are found on BTA3, BTA5 and BTA7. They overlap with genes playing roles in several biological functions and pathways. These include but are not limited to growth and feed intake, cell stability, protein stability and sweat gland development. This review may further guide targeted genome studies aiming to assess the importance of candidate causative mutations, within regulatory and protein-coding genome regions, to further understand the biological mechanisms underlying African cattle's unique adaption.

Genome-wide association study for wattles trait in the dairy goat breed

Dairy goats are significant livestock that provide high-quality milk sources in the world. The wattles trait is an evident phenotypic character on the neck of a dairy goat, which is considered to be under genetic control. We collected samples of 189 dairy goats, including 94 with wattles and 95 without wattles, from four different farms and multiple dairy goat breeds. The samples were genotyped with the GeneSeek Genomic Profiler Goat 70 K SNP chip. Genome-wide association studies (GWAS) in wattles have identified associations with single nucleotide polymorphisms (SNPs) at chromosome 10. In this area, an extremely strong association locus was assigned to FMN1 (Formin 1) belongs to the formin homology family and is associated with limb deformity, other candidate genes of interest confirmed for wattles were ARHGAP11A (Rho GTPase Activating Protein 11 A) and GJD2 (Gap Junction Protein Delta 2). Meanwhile, we found the presence or absence of wattles had no significant effect on milk yield. This research will provide genetic resources useful to explore genetic factors affecting the trait.

Genome-wide association analysis of time to heading and maturity in bread wheat using 55K microarrays

To investigate the genetic mechanisms underlying the reproductive traits (time to flowering and maturity) in wheat and identify candidate genes associated, a phenotypic analysis was conducted on 239 wheat accessions (lines) from around the world. A genome-wide association study (GWAS) of wheat heading and maturity phases was performed using the MLM (Q+K) model in the TASSLE software, combined with the Wheat 55K SNP array. The results revealed significant phenotypic variation in heading and maturity among the wheat accessions across different years, with coefficients of variation ranging from 0.96% to 1.97%. The phenotypic data from different years exhibited excellent correlation, with a genome-wide linkage disequilibrium (LD) attenuation distance of 3 Mb. Population structure analysis, evolutionary tree analysis, and principal component analysis indicated that the 239 wheat accessions formed a relatively homogeneous natural population, which could be divided into three subgroups. The GWAS results identified a total of 293 SNP marker loci that were significantly associated with wheat heading and maturity stages (P ≤ 0.001) in different environments. Among them, nine stable SNP marker loci were consistently detected in multiple environments. These marker loci were distributed on wheat chromosomes 1A、1B、2D、3A、5B、6D and 7A. Each individual locus explained 4.03%-16.06% of the phenotypic variation. Furthermore, through careful analysis of the associated loci with large phenotypic effect values and stable inheritance, a total of nine candidate genes related to wheat heading and maturity stages were identified. These findings have implications for molecular marker-assisted selection breeding programs targeting specific wheat traits at the heading and maturity stages. In summary, this study conducted a comprehensive GWAS of wheat heading and maturity phases, revealing significant associations between genetic markers and key developmental stages in wheat. The identification of candidate genes and marker loci provides valuable information for further studies on wheat breeding and genetic improvement targeted at enhancing heading and maturity traits.

Investigation of the molecular genetic causes of non-syndromic primary ovarian ınsufficiency by next generation sequencing analysis

Objective

The aim of this study is to investigate the molecular genetic causes of non-syndromic primary ovarian insufficiency (POI) cases with the gene panel basedon next generation sequencing analysis and to establish the relationship between genotype and phenotype.

Materials and methods

Twenty three cases aged 14-40 years followed up with POI were included. Patients with a karyotype of 46, XX, primary or secondary amenorrhea before the age of 40, with elevated FSH (>40 IU/mL) and low AMH levels (<0.03 ng/mL) were included in the study. Molecular genetic analyzes were performed by the next generation sequencing analysis method targeted with the TruSight TM Exome panel.

Results

Median age of the cases was 17.8 (14.0-24.3) years, and 12 (52%) cases admitted before the age of 18. Fifteen (65%) patients had consanguineous parents. In2 (8.6%) cases, variants detected were in genes that have been previously proven to cause POI. One was homozygous variant in FIGLA gene and the other was homozygous variant in PSMC3IP gene. Heterozygous variants were detected in PROK2, WDR11 and CHD7 associated with hypogonadotropic hypogonadism, but these variants are insufficient to contribute to the POI phenotype.

Conclusion

Genetic panels based on next generation sequencing analysis technologies can be used to determine the molecular genetic diagnosis of POI, which has a highly heterogeneous genetic basis.

Genome-Wide Association Study Reveals Novel Powdery Mildew Resistance Loci in Bread Wheat

Powdery mildew (PM), caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt), significantly threatens global bread wheat production. Although the use of resistant cultivars is an effective strategy for managing PM, currently available wheat cultivars lack sufficient levels of resistance. To tackle this challenge, we conducted a comprehensive genome-wide association study (GWAS) using a diverse panel of 286 bread wheat genotypes. Over three consecutive years (2020-2021, 2021-2022, and 2022-2023), these genotypes were extensively evaluated for PM severity under field conditions following inoculation with virulent Bgt isolates. The panel was previously genotyped using the Illumina 90K Infinium iSelect assay to obtain genome-wide single-nucleotide polymorphism (SNP) marker coverage. By applying FarmCPU, a multilocus mixed model, we identified a total of 113 marker-trait associations (MTAs) located on chromosomes 1A, 1B, 2B, 3A, 3B, 4A, 4B, 5A, 5B, 6B, 7A, and 7B at a significance level of p ≤ 0.001. Notably, four novel MTAs on chromosome 6B were consistently detected in 2020-2021 and 2021-2022. Furthermore, within the confidence intervals of the identified SNPs, we identified 96 candidate genes belonging to different proteins including 12 disease resistance/host-pathogen interaction-related protein families. Among these, protein kinases, leucine-rich repeats, and zinc finger proteins were of particular interest due to their potential roles in PM resistance. These identified loci can serve as targets for breeding programs aimed at developing disease-resistant wheat cultivars.

Fine mapping of qBK1.2, a major QTL governing resistance to bakanae disease in rice

Bakanae disease caused by Fusarium fujikuroi is an emerging disease of rice causing losses in all rice-growing regions around the world. A BC2F2 population was developed by backcrossing the recurrent parent Pusa Basmati 1121 (PB1121) with the recombinant inbred line RIL28, which harbors a major quantitative trait locus (QTL) governing resistance to bakanae, qBK1.2. MassARRAY-based single-nucleotide polymorphism (SNP) assays targeting the genomic region of qBK1.2 helped in fine mapping the QTL to a region of 130 kb between the SNP markers rs3164311 and rs3295562 using 24 recombinants. In-silico mining of the fine-mapped region identified 11 putative candidate genes with functions related to defense. The expression analysis identified two significantly differentially expressed genes, that is, LOC_Os01g06750 and LOC_Os01g06870, between the susceptible genotype PB1121 and the resistant genotypes Pusa1342 and R-NIL4. Furthermore, the SNPs identified in LOC_Os01g06750 produced minor substitutions of amino acids with no major effect on the resistance-related functional motifs. However, LOC_Os01g06870 had 21 amino acid substitutions, which led to the creation of the leucine-rich repeat (LRR) domain in the resistant genotype Pusa1342, thereby making it a potential candidate underlying the major bakanae-resistant QTL qBK1.2. The markers used in the fine mapping program are of immense utility in marker-assisted breeding for bakanae resistance in rice.