Region Of Locus Research Articles

Grapevine pangenome facilitates trait genetics and genomic breeding.

Grapevine breeding is hindered by a limited understanding of the genetic basis of complex agronomic traits. This study constructs a graph-based pangenome reference (Grapepan v.1.0) from 18 newly generated phased telomere-to-telomere assemblies and 11 published assemblies. Using Grapepan v.1.0, we build a variation map with 9,105,787 short variations and 236,449 structural variations (SVs) from the resequencing data of 466 grapevine cultivars. Integrating SVs into a genome-wide association study, we map 148 quantitative trait loci for 29 agronomic traits (50.7% newly identified), with 12 traits significantly contributed by SVs. The estimated heritability improves by 22.78% on average when including SVs. We discovered quantitative trait locus regions under divergent artificial selection in metabolism and berry development between wine and table grapes, respectively. Moreover, significant genetic correlations were detected among the 29 traits. Under a polygenic model, we conducted genomic predictions for each trait. In general, our study facilitates the breeding of superior cultivars via the genomic selection of multiple traits.

Evaluating the effects of silent genes on pairwise kinship testing

Short tandem repeat (STR) loci are frequently utilized in kinship testing, and mutations of a single base occurring in the primer-binding region of the STR locus can result in the failure of allelic amplification and the emergence of silent genes. Silent genes are not observable and, therefore, are excluded from the genotypes assessed. Pedigree likelihood ratios (LRs) are often employed in kinship testing to determine the likelihood of different kinship scenarios. LR values are derived from various types of genotypes. LRexact values are based on the exact or actual genotypes, which may include silent genes. Conversely, LRobserve values are based on observed genotypes that exclude silent genes, while LRadjust values incorporate all potential genotypes, including both observed and those with silent genes. Initially, the formulae for LRs in 1st degree, 2nd degree, and 3rd degree kinship testing are presented according to different genotype forms of pairwise individuals. The correctness of these formulae is then verified using the Familias software, and the results are compared with those from the GeneVisa software (www.genevisa.net). Lastly, the simulation modules of GeneVisa are used to assess the impact of silent genes on pairwise kinship testing. The findings indicate that the overall impact of silent genes is minimal, although in some cases, the effects can be relatively significant. The influence of silent genes generally decreases as the kinship relationship becomes more distant. In specific kinship tests, the effect of silent genes is reduced when the individuals are unrelated compared to when there is a kinship relationship. Utilizing the LRadjust value for 1st degree and 2nd degree kinship testing can substantially mitigate the effects of silent genes.

Pharmacogenomics of CRP response to statins: a GIST consortium study

Abstract Introduction Statins are first line treatments in the primary and secondary prevention of cardiovascular disease. Prior clinical studies have shown that statins act independently of lipid lowering mechanisms to decrease C-reactive protein (CRP), a marker of inflammation. Purpose To elucidate genetic loci associated with CRP response to statins. Methods CRP response was specified as the change from baseline in log CRP after at least 4 weeks of statin therapy. Cohort level Genome-wide Association Studies (GWAS) of this CRP response was performed by linear regression analysis adjusted for baseline CRP, age, sex and BMI covariates using genetic data imputed to 1000 Genomes, testing ~10 million single nucleotide polymorphisms (SNPs) of minor allele frequency (MAF) >2%. Following quality control with EasyQC, a meta-analysis was conducted with METAL to combine GWAS results from six different cohorts of European ancestry (CARDS, FHS, JUPITER, MESA, PARC, PROSPER) within the GIST consortium. 1Mb loci regions were defined, centred +/-500kb around the lead SNP. Conditional analysis was performed using GCTA. Results There were 11,075 statin-treated individuals. The meta-analysis results revealed two loci achieving genome wide significance (P<5e-8): APOE (chr 19) and HNF1A (chr 12) (Figure 1). A signal at the CRP locus was highly suggestive (P=2.3e-7). All three loci are known to be associated with CRP levels in the absence of statin. Conditional analysis did not reveal secondary signals. The most associated variant at the APOE locus was the missense SNP rs429358, which contributes to the APOE E4 haplotype and is a risk locus for both dyslipidaemia and Alzheimer’s dementia. It is a C>T missense variant of MAF 0.12 leading to a R/C amino acid change, more common in African ancestry populations (MAF 0.27), and least common in Asian populations (MAF 0.09). A nominally significant association (P=0.09) for interaction with randomized allocation to statin v. placebo within samples derived from clinical trials supports this pharmacogenetic effect at APOE beyond genetic effects on baseline CRP levels. The most associated variant at the HNF1A locus was the intronic SNP rs11065384, which is in strong LD with the HNF1A missense SNP (rs1169288). The HNF1A locus is associated with diabetes, cholesterol levels, and coronary artery disease. There was no between-study heterogeneity at the top three loci (P>0.1), and forest plots show contribution to these signals from all 6 studies within the meta-analysis (Figure 2 A,B,C). Conclusions APOE and HNF1A are associated with CRP response to statins at the level of genome-wide significance. TheAPOE E4 signal is also known to be associated with LDL-Cholesterol response to statins, whereas the HNF1A locus is identified here for the first time as a pharmacogenetic determinant of statin response. The effect of this interindividual variability in CRP response on non-cardiovascular clinical outcomes should be elucidated.Figure 1Figure 2 A,B,C

Integrated GWAS and transcriptome analysis reveals key genes associated with muscle fibre and fat traits in Gushi chicken

ABSTRACT 1. In the following experiment meat quality traits of a Gushi-Anka F2 resource population were measured, and their heritability estimated. Intramuscular fat (IMF) had medium heritability (0.35) but leg muscle fibre density (LMD), leg muscle fibre diameter (LMF), breast muscle fibre density (BMD), fresh fat content (FFA), and absolute dry fat content (AFC) had low heritability (0–0.2). The IMF presented the most important genetic additive effect among the poultry meat quality-related traits studied. 2. The phenotypic data of meat quality traits in the Gushi-Anka F2 resource population were combined with genotyping by sequencing (GBS) data to obtain genotype data. Six meat quality traits in 734 birds were analysed by GWAS. Based on these variants, 83 significant (–log10(p) > 4.42) single nucleotide polymorphisms and four quantitative trait loci (QTL) regions corresponding to 175 genes were identified. Further linkage disequilibrium (LD) analysis was conducted on chromosome 13 (Chr13) and chromosome 27 (Chr27) QTL regions. 3. Based on the transcriptome data and GWAS results, 12 shared genes - ITGB3, DNAJC27, ETV4, C7orf50, FKBP1B, G3BP1, IGF2BP1, KCNH6, LOC416263, SCARA5, SMIM5 and TBL1XR1 were identified as candidate genes influencing muscle fibre and fat traits.

Identification of genomic regions associated with fatty acid metabolism across blood, liver, backfat and muscle in pigs

BackgroundThe composition and distribution of fatty acids (FA) are important factors determining the quality, flavor, and nutrient value of meat. In addition, FAs synthesized in the body participate in energy metabolism and are involved in different regulatory pathways in the form of signaling molecules or by acting as agonist or antagonist ligands of different nuclear receptors. Finally, synthesis and catabolism of FAs affect adaptive immunity by regulating lymphocyte metabolism. The present study performed genome-wide association studies using FA profiles of blood, liver, backfat and muscle from 432 commercial Duroc pigs.ResultsTwenty-five genomic regions located on 15 Sus scrofa chromosomes (SSC) were detected. Annotation of the quantitative trait locus (QTL) regions identified 49 lipid metabolism-related candidate genes. Among these QTLs, four were identified in more than one tissue. The ratio of C20:4n-6/C20:3n-6 was associated with the region on SSC2 at 7.56–14.26 Mb for backfat, liver, and muscle. Members of the fatty acid desaturase gene cluster (FADS1, FADS2, and FADS3) are the most promising candidate genes in this region. Two QTL regions on SSC14 (103.81–115.64 Mb and 100.91–128.14 Mb) were identified for FA desaturation in backfat and muscle. In addition, two separate regions on SSC9 at 0 – 14.55 Mb and on SSC12 at 0–1.91 Mb were both associated with the same multiple FA traits for backfat, with candidate genes involved in de novo FA synthesis and triacylglycerol (TAG) metabolism, such as DGAT2 and FASN. The ratio C20:0/C18:0 was associated with the region on SSC5 at 64.84–78.32 Mb for backfat. Furthermore, the association of the C16:0 content with the region at 118.92–123.95 Mb on SSC4 was blood specific. Finally, candidate genes involved in de novo lipogenesis regulate T cell differentiation and promote the generation of palmitoleate, an adipokine that alleviates inflammation.ConclusionsSeveral SNPs and candidate genes were associated with lipid metabolism in blood, liver, backfat, and muscle. These results contribute to elucidating the molecular mechanisms implicated in the determination of the FA profile in different pig tissues and can be useful in selection programs that aim to improve health and energy metabolism in pigs.

Genome-wide association studies for milk production traits in two autochthonous Aosta cattle breeds

Genome-wide association studies (GWASs) are used to identify quantitative trait loci for phenotypic traits of interest. The use of multilocus mixed models allows to correct for population stratification and account for long-range linkage disequilibrium. In this study, GWASs were conducted to identify the genetic bases of milk production (milk yield, protein and fat composition, and yield) in two autochthonous dual-purpose cattle breeds from the Aosta Valley. Using either the breeding values or the deregressed proofs, common significative single nucleotide polymorphisms have been identified for milk yield, protein percentage, and fat percentage. Two major quantitative trait loci regions have been identified on the chromosomes 5 and 14 for the fat percentage, harbouring the MGST1, CYHR1, VPS28, and CPSF1 genes. For the protein percentage, a candidate region has been identified on BTA 6; in this region, the CSN1S1, CSN2, HSTN, CSN3, and RUFY3 genes are annotated. Most of the identified genes have already been associated with milk composition in other studies on cosmopolitan and local cattle. These results show that the genes involved in milk composition quantitative traits in the Aosta cattle are common also in other cattle breeds and they can be further investigated with the use of whole genome sequencing data.

A GWAS study highlights significant associations between a series of indels in a FLOWERING LOCUS T gene promoter and flowering time in white lupin (Lupinus albus L.)

BackgroundWhite lupin (Lupinus albus L.) is a high-protein Old World grain legume with remarkable food and feed production interest. It is sown in autumn or early spring, depending on the local agroclimatic conditions. This study aimed to identify allelic variants associated with vernalization responsiveness, in order to improve our knowledge of legume flowering regulatory pathways and develop molecular selection tools for the desired phenology as required for current breeding and adaptation to the changing climate.ResultsSome 120 white lupin accessions originating from a wide range of environments of Europe, Africa, and Asia were phenotyped under field conditions in three environments with different intensities of vernalization, namely, a Mediterranean and a subcontinental climate sites of Italy under autumn sowing, and a suboceanic climate site of France under spring sowing. Two hundred sixty-two individual genotypes extracted from them were phenotyped in a greenhouse under long-day photoperiod without vernalization. Phenology data, and marker data generated by Diversity Arrays Technology sequencing (DArT-seq) and by PCR-based screening targeting published quantitative trait loci (QTLs) from linkage map and newly identified insertion/deletion polymorphisms in the promoter region of the FLOWERING LOCUS T homolog, LalbFTc1 gene (Lalb_Chr14g0364281), were subjected to a genome-wide association study (GWAS). Population structure followed differences in phenology and isolation by distance pattern. The GWAS highlighted numerous loci significantly associated with flowering time, including four LalbFTc1 gene promoter deletions: 2388 bp and 2126 bp deletions at the 5’ end, a 264 bp deletion in the middle and a 28 bp deletion at the 3’ end of the promoter. Besides LalbFTc1 deletions, this set contained DArT-seq markers that matched previously published major QTLs in chromosomes Lalb_Chr02, Lalb_Chr13 and Lalb_Chr16, and newly discovered QTLs in other chromosomes.ConclusionsThis study highlighted novel QTLs for flowering time and validated those already published, thereby providing novel evidence on the convergence of FTc1 gene functional evolution into the vernalization pathway in Old World lupin species. Moreover, this research provided the set of loci specific for extreme phenotypes (the earliest or the latest) awaiting further implementation in marker-assisted selection for spring- or winter sowing.

Comparative transcriptome profiling to untie the drought tolerance molecular mechanism of backcross rice (Oryza sativa L.) inbred

AbstractRice production is severely hampered by drought stress, which causes enormous economic losses. The issue of global climate change is gaining importance, and hence development of rice genotypes tolerant to drought stress is becoming more critical. To address this issue, backcross inbred lines (BILs) developed were subjected to drought stress, and their molecular mechanism was studied. The drought‐tolerant parent Apo and drought‐susceptible, high‐yielding IR64 along with two BILs, namely, CB 229 (qDTY2.2 + qDTY3.1 + qDTY8.1) and CB 193‐3 (qDTY3.1 + qDTY8.1) were tested in a greenhouse for their response to drought. In this study, CB 229 showed better performance under water stress irrigated conditions; it was on par with IR64. Drought‐responsive transcriptome profiling was carried out in both the parents and the superior BIL CB 229 through the RNA‐Seq approach. About 3050 differentially expressed genes (DEGs) (2021 upregulated and 1029 downregulated) were detected in tolerant BIL CB 229 in drought stress. Most of the DEGs were involved in carbohydrate metabolism and the formation of cell walls, as well as genes associated with metabolite adaptability, ROS homeostasis and post‐transcriptional regulation. Genes such as chaperone protein, senescence‐induced receptor‐like serine, mannose‐6‐phosphate isomerase, aquaporin and heat shock proteins (LOC_Os02g26840, LOC_Os02g25720, LOC_Os07g35570, LOC_Os01g64970, etc.) were upregulated in the tolerant Apo and CB 229. It was observed that the BIL CB 229 yielded higher grains than both parents under moisture stress. Ninety‐four genes in the quantitative trait loci (QTLs) region were found to be differentially regulated in CB 229, Apo and IR64. Out of 94, nine genes co‐localized within the QTL qDTY2.2, 12 genes within qDTY3.1 and four genes within qDTY8.1 were differentially upregulated in CB 229 and downregulated in the susceptible genotype. The study revealed that the QTLs qDTY2.2, qDTY3.1 and qDTY8.1 are found to have complementary effects and offer enhanced levels of tolerance against drought due to complementation. Additionally, this analysis discovered new DEGs that may be involved in functions related to drought tolerance but lack function annotations. Future research can focus on a few of the significant DEGs found in this study. Taken together, our findings provide insight into drought tolerance mechanisms and could assist the development of novel strategies for improving drought tolerance in rice.

Kdm4a is an activity downregulated barrier to generate engrams for memory separation

Memory engrams are a subset of learning activated neurons critical for memory recall, consolidation, extinction and separation. While the transcriptional profile of engrams after learning suggests profound neural changes underlying plasticity and memory formation, little is known about how memory engrams are selected and allocated. As epigenetic factors suppress memory formation, we developed a CRISPR screening in the hippocampus to search for factors controlling engram formation. We identified histone lysine-specific demethylase 4a (Kdm4a) as a negative regulator for engram formation. Kdm4a is downregulated after neural activation and controls the volume of mossy fiber boutons. Mechanistically, Kdm4a anchors to the exonic region of Trpm7 gene loci, causing the stalling of nascent RNAs and allowing burst transcription of Trpm7 upon the dismissal of Kdm4a. Furthermore, the YTH domain containing protein 2 (Ythdc2) recruits Kdm4a to the Trpm7 gene and stabilizes nascent RNAs. Reducing the expression of Kdm4a in the hippocampus via genetic manipulation or artificial neural activation facilitated the ability of pattern separation in rodents. Our work indicates that Kdm4a is a negative regulator of engram formation and suggests a priming state to generate a separate memory.

GEGA (Gallus Enriched Gene Annotation): an online tool providing genomics and functional information across 47 tissues for a chicken gene-enriched atlas gathering Ensembl and Refseq genome annotations.

GEGA is a user-friendly tool designed to navigate through various genomic and functional information related to an enriched gene atlas in chicken that integrates the gene catalogues from the two reference databases, NCBI-RefSeq and EMBL-Ensembl/GENCODE, along with four additional rich resources such as FAANG and NONCODE. Using the latest GRCg7b genome assembly, GEGA encompasses a total of 78 323 genes, including 24 102 protein-coding genes (PCGs) and 44 428 long non-coding RNAs (lncRNAs), significantly increasing the number of genes provided by each resource independently. However, GEGA is more than just a gene database. It offers a range of features that allow us to go deeper into the functional aspects of these genes. Users can explore gene expression and co-expression profiles across 47 tissues from 36 datasets and 1400 samples, discover tissue-specific variations and their expression as a function of sex or age and extract orthologous genes or their genomic configuration relative to the closest gene. For the communities interested in a specific gene, a list of genes or a quantitative trait locus region in chicken, GEGA's user-friendly interface facilitates efficient gene analysis, easy downloading of results and a multitude of graphical representations, from genomic information to detailed visualization of expression levels.

Generation of New β-Conglycinin-Deficient Soybean Lines by Editing the lincRNA lincCG1 Using the CRISPR/Cas9 System.

Soybean β-conglycinin is a major allergen that adversely affects the nutritional properties of soybean. Soybean deficient in β-conglycinin is associated with low allergenicity and high nutritional value. Long intergenic noncoding RNAs (lincRNAs) regulate gene expression and are considered important regulators of essential biological processes. Despite increasing knowledge of the functions of lincRNAs, relatively little is known about the effects of lincRNAs on the accumulation of soybean β-conglycinin. The current study presents the identification of a lincRNA lincCG1 that was mapped to the intergenic noncoding region of the β-conglycinin α-subunit locus. The full-length lincCG1 sequence was cloned and found to regulate the expression of soybean seed storage protein (SSP) genes via both cis- and trans-acting regulatory mechanisms. Loss-of-function lincCG1 mutations generated using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) system led to the deficiency of the allergenic α'-, α-, and β-subunits of soybean β-conglycinin as well as higher content of proteins, sulfur-containing amino acids, and free arginine. The dominant null allele LincCG1, and consequently, the β-conglycinin-deficient phenotype associated with the lincCG1-gene-edited line was stably inherited by the progenies in a Mendelian fashion. The dominant null allele LincCG1 may therefore be exploited for engineering/developing novel hypoallergenic soybean varieties. Furthermore, Cas9-free and β-conglycinin-deficient homozygous mutant lines were obtained in the T1 generation. This study is the first to employ the CRISPR/Cas9 technology for editing a lincRNA gene associated with the soybean allergenic protein β-conglycinin. Moreover, this study reveals that lincCG1 plays a crucial role in regulating the expression of the β-conglycinin subunit gene cluster, besides highlighting the efficiency of employing the CRISPR/Cas9 system for modulating lincRNAs, and thereby regulating soybean seed components.

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

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

The Presence of esat-6 and cfp10 and Other Gene Orthologs of the RD 1 Region in Non-Tuberculous Mycobacteria, Mycolicibacteria, Mycobacteroides and Mycolicibacter as Possible Impediments for the Diagnosis of (Animal) Tuberculosis.

The Esx-1 family proteins of the Type VII secretion systems of Mycobacterium bovis and Mycobacterium tuberculosis have been assessed and are frequently used as candidates for tuberculosis (TB) diagnosis in both humans and animals. The presence of ESAT-6 and CFP 10 proteins, which are the most immunogenic proteins of the Esx-1 system and have been widely investigated for the immunodiagnosis of tuberculosis, in some Mycobacteriaceae and in Mycobacterium leprae, poses limitations for their use in specific diagnoses of TB. As such, to improve the specificity of the ESAT-6/CFP 10-based cell-mediated immunity (CMI) assays, other proteins encoded by genes within and outside the RD 1 region of the esx-1 locus have been evaluated as candidate antigens for CMI, as well as to investigate humoral responses in combination with ESAT-6 and or CFP 10, with varying specificity and sensitivity results. Hence, in this study, we evaluated various non-tuberculous mycobacteria (NTM), Mycolicibacterium, Mycolicibacter and Mycobacteroides species genomes available on the NCBI database for the presence and composition of the RD1 region of the esx-1 locus. In addition, we also assayed by polymerase chain reaction (PCR) and sequencing of Mycobacteriaceae available in our culture collection for the presence and sequence diversity of esxA and esxB genes encoding ESAT-6 and CFP 10, respectively. Whole genome sequence (WGS) data analysis revealed the presence of RD 1 gene orthologs in 70 of the over 100 published genomes of pathogenic and non- pathogenic Mycobcteriaceae other than tuberculosis. Among species evaluated from our culture collection, in addition to earlier reports of the presence of esxA and esxB in certain Mycolicibacterium, Mycolicibacterium septicum/peregrinum, Mycolicibacterium porcinum and Mycobacterium sp. N845T were also found to harbour orthologs of both genes. Orthologs of esxA only were detected in Mycobacterium brasiliensis, Mycolicibacterium elephantis and Mycolicibacterium flouroantheinivorans, whereas in Mycolicibacter engbackii, Mycolicibacterium mageritense and Mycobacterium paraffinicum, only esxB orthologs were detected. A phylogenetic analysis based on esxA and esxB sequences separated slow-growing from rapidly growing bacteria. These findings strengthen previous suggestions that esxA and esxB may be encoded in the majority of Mycobacteriaceae. The role of the Esx-1 system in both pathogenic and non-pathogenic Mycobacteriaceae needs further investigation, as these species may pose limitations to immunological assays for TB.

Targeting dependency on a paralog pair of CBP/p300 against de-repression of KREMEN2 in SMARCB1-deficient cancers

SMARCB1, a subunit of the SWI/SNF chromatin remodeling complex, is the causative gene of rhabdoid tumors and epithelioid sarcomas. Here, we identify a paralog pair of CBP and p300 as a synthetic lethal target in SMARCB1-deficient cancers by using a dual siRNA screening method based on the “simultaneous inhibition of a paralog pair” concept. Treatment with CBP/p300 dual inhibitors suppresses growth of cell lines and tumor xenografts derived from SMARCB1-deficient cells but not from SMARCB1-proficient cells. SMARCB1-containing SWI/SNF complexes localize with H3K27me3 and its methyltransferase EZH2 at the promotor region of the KREMEN2 locus, resulting in transcriptional downregulation of KREMEN2. By contrast, SMARCB1 deficiency leads to localization of H3K27ac, and recruitment of its acetyltransferases CBP and p300, at the KREMEN2 locus, resulting in transcriptional upregulation of KREMEN2, which cooperates with the SMARCA1 chromatin remodeling complex. Simultaneous inhibition of CBP/p300 leads to transcriptional downregulation of KREMEN2, followed by apoptosis induction via monomerization of KREMEN1 due to a failure to interact with KREMEN2, which suppresses anti-apoptotic signaling pathways. Taken together, our findings indicate that simultaneous inhibitors of CBP/p300 could be promising therapeutic agents for SMARCB1-deficient cancers.

Runs of homozygosity analysis for selection signatures in the Yellow Korean native chicken.

Yellow Korean native chicken (KNC-Y) is one of the five pure Korean indigenous chicken breeds that were restored through a government project in 1992. KNC-Y is recognized for its superior egg production performance compared to other KNC lines. In this study, we performed runs of homozygosity (ROH) analysis to discover selection signatures associated with egg production traits in the KNC-Y population. A total of 675 DNA samples from KNC-Y were genotyped to generate single nucleotide polymorphism (SNP) data using custom 60K Affymetrix SNP chips. ROH analysis was performed using PLINK software, with predefined parameters set for the analysis. The threshold of ROH island was defined as the top 1% frequency of SNPs withing the ROH among the population. In the KNC-Y population, a total of 29,958 runs of homozygosity (ROH) fragments were identified. The average total length of ROH was 120.84 Mb, with each ROH fragment having an average length of 2.71 Mb. The calculated ROH-based inbreeding coefficient (FROH) was 0.13. Furthermore, we revealed the presence of ROH islands on chromosomes 1, 2, 4, 5, 7, 8, and 11. Within the identified regions, a total of 111 genes were annotated, and among them were genes related to economic traits, including PRMT3, ANO5, HDAC4, LSS, PLA2G4A, and PTGS2. Most of the overlapping quantitative trait locus regions with ROH islands were found to be associated with production traits. This study conducted a comprehensive analysis of ROH in the KNC-Y population. Notably, among the findings, the PTGS2 gene is believed to play a crucial role in influencing the laying performance of KNC-Y.

Variable effects of captivity on microbiomes in populations of IUCN-endangered Blanding's turtles (Emydoidea blandingii).

Microbiome composition is increasingly considered in species reintroduction efforts and may influence survival and reproductive success. Many turtle species are threatened by anthropogenic pressures and are frequently raised in captivity for reintroduction efforts, yet little is known about turtle microbiome composition in either wild or captive settings. Here, we investigated trends in microbiome composition of captive and wild IUCN-endangered Blanding's turtles (Emydoidea blandingii). We amplified and sequenced the V4 region of the 16S rDNA locus from plastron, cloaca, and water samples of wild E. blandingii adults and two populations of captive E. blandingii juveniles being raised for headstarting. Plastron, cloaca, and water-associated microbiomes differed strongly from each other and were highly variable among captive sites and between captive and wild sites. Across plastron, cloaca, and water-associated microbial communities, microbial diversity changed over time, but not in a predictable direction between captive sites. Plastron beta diversity correlated with growth rate in captive samples, indicating that external microbiomes may correlate with individual fitness. Our results indicate that external and internal microbiomes vary between captive and wild turtles and may reflect differences in fitness of captive-raised individuals.

A simple sequencing protocol for genotyping the &lt;i&gt;HLA-C&lt;/i&gt;locus by the Sanger method

The HLA-C gene, which belongs to the HLA superfamily (Human Leukocyte Antigen), codes for the Major Histocompatibility Complex (MHC), which plays crucial roles in the human immune system. This study aimed to develop a simple sequencing protocol by the Sanger method using fewer primers and reactions for genotyping the HLA-C gene than published protocols. The simple protocol with three primers includes one PCR reaction and two sequencing reactions. The primer set comprising SEQ ID1 and SEQ ID2 was used for the PCR reaction to specifically amplify the exon 2 – exon 3 region of the HLA-C locus, which contains the typical SNPs of each HLA-C allele. The PCR product was purified and used as a template for the sequencing reactions. Two forward primers, SEQ ID1 and SEQ ID3 were used for sequencing, in which, the SEQ ID1 forward primer is located in the intron 1 region and the SEQ ID3 forward primer is located in the intron 2 region of the HLA-C gene. Testing the simple sequencing protocol on four samples of known HLA-C genotypes showed 100% accurate results. The established Sanger sequencing protocol is simple to implement, and reduces cost and time. Thus, this protocol can be used for HLA-C sequencing for pharmacogenetic studies and applications.

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.

Molecular mapping of the broad bean wilt virus 2 resistance locus bwvr in Capsicum annuum using BSR-seq

Key messageBulked segregant RNA seq of pools of pepper accessions that are susceptible or resistant to Broad bean wilt virus 2 identifies a gene that might confer resistance to this devastating pathogen.The single-stranded positive-sense RNA virus Broad bean wilt virus 2 (BBWV2) causes substantial damage to pepper (Capsicum annuum) cultivation. Here, we describe mapping the BBWV2 resistance locus bwvr using a F7:8 recombinant inbred line (RIL) population constructed by crossing the BBWV2-resistant pepper accession ‘SNU-C’ with the susceptible pepper accession ‘ECW30R.’ All F1 plants infected with the BBWV2 strain PAP1 were susceptible to the virus, and the RIL population showed a 1:1 ratio of resistance to susceptibility, indicating that this trait is controlled by a single recessive gene. To map bwvr, we performed bulked segregant RNA-seq (BSR-seq). We sequenced pools of resistant and susceptible lines from the RILs and aligned the reads to the high-quality ‘Dempsey’ reference genome to identify variants between the pools. This analysis identified 519,887 variants and selected the region from 245.9–250.8 Mb of the Dempsey reference genome as the quantitative trait locus region for bwvr. To finely map bwvr, we used newly designed high-resolution melting (HRM) and Kompetitive allele specific PCR (KASP) markers based on variants obtained from the BSR-seq reads and the PepperSNP16K array. Comparative analysis identified 11 SNU-C-specific SNPs within the bwvr locus. Using markers derived from these variants, we mapped the candidate bwvr locus to the region from 246.833–246.949 kb. SNU-C-specific variants clustered near DEM.v1.00035533 within the bwvr locus. DEM.v1.00035533 encodes the nitrate transporter NPF1.2 and contains a SNP within its 5′ untranslated region. The bwvr locus, which contains four genes including DEM.v1.00035533, could represent a valuable resource for global pepper breeding programs.

Abstract 5592: HMGA1 modulates chromatin state and transcriptional networks involved in plasticity in refractory myeloid leukemia

Abstract Although High Mobility Group A1 (HMGA1) chromatin regulators are overexpressed in refractory acute myeloid leukemia (AML) and diverse solid tumors, targetable mechanisms underlying HMGA1 remain elusive. HMGA1 is highly expressed in stem cells and poorly differentiated cancers where it modulates chromatin structure to activate developmental gene networks. In AML arising from JAK2V617F myeloproliferative neoplasms (MPN), HMGA1 is required for leukemogenesis by inducing transcriptional networks that drive aberrant proliferation and cell fate. Like MPN AML, KMT2A-rearranged (KMT2A-r) AML is refractory to therapy and highly lethal. KMT2A-r AML is caused by rearrangements of the KMT2A gene, which encode abnormal fusion proteins that induce pro-leukemogenic genes, including HOX genes. While prior work identified KMT2A fusion partners and protein complexes, this has not led to better therapies for most patients. We therefore sought to examine HMGA1 function in KMT2A-r AML. We hypothesized that: 1) HMGA1 drives plasticity in KMT2A-r AML by altering chromatin state and transcriptional networks, and, 2) Targeting HMGA1 will enhance sensitivity to therapy. In two KMT2A-r AML patient cohorts, HMGA1 is highly overexpressed compared to CD34+ stem and progenitors from healthy controls, suggesting it drives leukemogenesis. To define HMGA1 function, we silenced HMGA1 using CRISPR or short hairpin RNA in KMT2A-r cell lines (MOLM-14, THP-1, MV4-11). Strikingly, silencing HMGA1 disrupts proliferation and clonogenicity in vitro, while prolonging survival in immunosuppressed mice transplanted with KMT2A-r AML cells. Intriguingly, leukemic cells that engraft in the bone marrow from the pool of cells with HMGA1 silencing express higher levels of HMGA1, suggesting that escape from gene silencing and a specific level of HMGA1 is required for leukemic engraftment and expansion. HMGA1 also increases in KMT2A-r AML cells that become resistant to cytarabine (AraC), suggesting that HMGA1 endows AML cells with the capacity to survive and expand after AraC. Because repression of tumor suppressors was recently defined as a mechanism of therapy resistance in AML, we also examined tumor suppressor pathways in both JAK2V617F and KMT2A-r AML. Intriguingly, HMGA1 occupies the promoter region of the CDKN1A locus by chromatin immunoprecipitation and represses its expression in KMT2A-r and JAK2V617F AML cells. To define additional mechanisms underlying HMGA1 in KMT2A-r AML, multi-omics studies are underway. In preliminary work thus far, we also found that HMGA1 induces HOX genes. Together, these findings implicate HMGA1 as a promising therapeutic target and novel epigenetic regulator in KMT2A-r AML by altering chromatin state and downstream gene expression. Citation Format: Bailey E. West, Jung-Hyun Kim, Audrey-Ann Supreme, Iliana Herrera, Zanshé Thompson, Li Z. Luo, Faiza Shaik, Joseph Kim, Hyunsung Woo, Soheil Meshinchi, Rhonda E. Ries, Linda M. Resar. HMGA1 modulates chromatin state and transcriptional networks involved in plasticity in refractory myeloid leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 5592.