Region Of Locus Research Articles

Advancing the Indian cattle pangenome: characterizing non-reference sequences in Bos indicus

BackgroundIndia harbors the world’s largest cattle population, encompassing over 50 distinct Bos indicus breeds. This rich genetic diversity underscores the inadequacy of a single reference genome to fully capture the genomic landscape of Indian cattle. To comprehensively characterize the genomic variation within Bos indicus and, specifically, dairy breeds, we aim to identify non-reference sequences and construct a comprehensive pangenome.ResultsFive representative genomes of prominent dairy breeds, including Gir, Kankrej, Tharparkar, Sahiwal, and Red Sindhi, were sequenced using 10X Genomics ‘linked-read’ technology. Assemblies generated from these linked-reads ranged from 2.70 Gb to 2.77 Gb, comparable to the Bos indicus Brahman reference genome. A pangenome of Bos indicus cattle was constructed by comparing the newly assembled genomes with the reference using alignment and graph-based methods, revealing 8 Mb and 17.7 Mb of novel sequence respectively. A confident set of 6,844 Non-reference Unique Insertions (NUIs) spanning 7.57 Mb was identified through both methods, representing the pangenome of Indian Bos indicus breeds. Comparative analysis with previously published pangenomes unveiled 2.8 Mb (37%) commonality with the Chinese indicine pangenome and only 1% commonality with the Bos taurus pangenome. Among these, 2,312 NUIs encompassing ~ 2 Mb, were commonly found in 98 samples of the 5 breeds and designated as Bos indicus Common Insertions (BICIs) in the population. Furthermore, 926 BICIs were identified within 682 protein-coding genes, 54 long non-coding RNAs (lncRNA), and 18 pseudogenes. These protein-coding genes were enriched for functions such as chemical synaptic transmission, cell junction organization, cell-cell adhesion, and cell morphogenesis. The protein-coding genes were found in various prominent quantitative trait locus (QTL) regions, suggesting potential roles of BICIs in traits related to milk production, reproduction, exterior, health, meat, and carcass. Notably, 63.21% of the bases within the BICIs call set contained interspersed repeats, predominantly Long Interspersed Nuclear Elements (LINEs). Additionally, 70.28% of BICIs are shared with other domesticated and wild species, highlighting their evolutionary significance.ConclusionsThis is the first report unveiling a robust set of NUIs defining the pangenome of Bos indicus breeds of India. The analyses contribute valuable insights into the genomic landscape of desi cattle breeds.

The Construction of Heterothallic Strains of Komagataella kurtzmanii Using the I-SceI Meganuclease.

The methylotrophic yeast Komagataella kurtzmanii belongs to the group of homothallic fungi that are able to spontaneously change their mating type by inversion of chromosomal DNA in the MAT locus region. As a result, natural and genetically engineered cultures of these yeasts typically contain a mixture of sexually dimorphic cells that are prone to self-diploidisation and spore formation accompanied by genetic rearrangements. These characteristics pose a significant challenge to the development of genetically stable producers for industrial use. In the present study, we constructed heterothallic strains of K. kurtzmanii, ensuring a constant mating type by unifying the genetic sequences in the active and silent MAT loci. To obtain such strains, we performed site-directed inactivation of one of the two yeast MAT loci, replacing its sequence with a selective HIS4 gene surrounded by I-SceI meganuclease recognition sites. We then used transient expression of the SCE1 gene, encoding a recombinant I-SceI meganuclease, to induce site-specific cleavage of HIS4, followed by damage repair by homologous recombination in mutant cells. As a result, heterothallic strains designated 'Y-727-2(alpha)' and 'Y-727-9(a)', which correspond to the α and a mating type, respectively, were obtained. The strains demonstrated a loss of the ability to self-diploidize. The results of PCR and whole genome analysis confirmed the identity of the contents of the MAT loci. Analysis of the genomes of the final strains, however, revealed a fusion of chromosome 3 and chromosome 4 in strain Y-727-2(alpha)-1. This finding was subsequently confirmed by pulsed-field gel electrophoresis of yeast chromosomes. However, the ability of the Y-727-2(alpha)-derived producers to efficiently secrete recombinant β-galactosidase was unaffected by this genomic rearrangement.

Lower expression of BIN1's neuronal isoform in vulnerable excitatory neurons increases risk in Alzheimer's disease

Background: Neurons in Alzheimer's disease (AD) experience elevated DNA damage, with DNA repair sites enriched at enhancer regions of genes essential for neuronal survival. Excitatory neurons in the cortical superficial layers expressing CUX2 and RORB (Cux2+/Rorb+), are selectively vulnerable in AD, but their relationship to single nucleotide polymorphisms (SNPs) in AD genome-wide association studies (GWAS) is unclear. Objective: This study aimed to identify and characterize functional AD-GWAS SNPs using single-nucleus RNA sequencing data, focusing on selectively vulnerable neurons and DNA repair hotspots. Methods: Filters were applied to identify candidate SNPs based on overlap with repair hotspots, RNA expression, transcription factor binding, AD association, and epigenetic significance. In vitro assays and analyses of large datasets from bulk RNA-seq (n = 1894), proteomics (n = 400), and single-nucleus RNA-seq (n = 424, 1.6 M cells) were conducted. Results: BIN1 SNP, rs78710909, met multiple criteria - located in an AD-GWAS locus, repair hotspot, and promoter region. rs78710909C exhibits 1.52× higher AD risk and 5.4× differential transcription factor binding. In vitro, rs78710909C shows greater enhancer activity and weaker p53 but stronger E2F1 binding. BIN1's neuronal isoform is neuroprotective, but its AD expression is lower (p < 0.01). Moreover, only in AD and Cux2+/Rorb + neurons, rs78710909C is associated with a lower average BIN1 neuronal isoform ratio (p < 0.01). The genes upregulated in neurons with lower neuronal isoform ratio were associated with the hallmarks of AD pathology. Conclusions: In a disease-relevant mechanism, the BIN1 SNP rs78710909C is associated with a lower ratio of BIN1's neuronal isoform which increases the vulnerability of specific excitatory neurons in AD patients.

QTL mapping provides new insights into emamectin benzoate resistance in salmon lice, Lepeophtheirus salmonis

BackgroundThe salmon louse (Lepeophtheirus salmonis) is a parasite of wild and farmed salmonid fish, causing huge economic damage to the commercial farming of Atlantic salmon (Salmo salar) in the northern hemisphere. The avermectin emamectin benzoate (EMB) is widely used for salmon delousing. While resistance to EMB is widespread in Atlantic populations of L. salmonis, the molecular mechanisms of resistance remain to be elucidated. The aim of the present work was to obtain insights into potential EMB resistance mechanisms by identifying genetic and transcriptomic markers associated with EMB resistance.ResultsCrosses were performed between EMB-susceptible and -resistant L. salmonis, sourced from two parental strains isolated in Scotland, producing fully pedigreed families. The EMB susceptibility of individual parasites was characterised using time-to-response bioassays. Parasites of two families were subjected to double digest restriction site-associated DNA sequencing (ddRAD-seq) for simultaneous discovery of single nucleotide polymorphisms (SNPs) and genotyping. Data analysis revealed that EMB resistance is associated with one quantitative trait locus (QTL) region on L. salmonis chromosome 5. Marker-trait association was confirmed by genotyping assays for 7 SNPs in two additional families. Furthermore, the transcriptome of male parasites of the EMB-susceptible and -resistant L. salmonis parental strains was assessed. Among eighteen sequences showing higher transcript expression in EMB-resistant as compared to drug-susceptible lice, the most strongly up-regulated gene is located in the above QTL region and shows high homology to β spectrin, a cytoskeleton protein that has roles in neuron architecture and function. Further genes differentially regulated in EMB-resistant lice include a glutathione S-transferase (GST), and genes coding for proteins with predicted roles in mitochondrial function, intracellular signalling or transcription.ConclusionsMajor determinants of EMB resistance in L. salmonis are located on Chromosome 5. Resistance can be predicted using a limited number of genetic markers. Genes transcriptionally up-regulated in EMB resistant parasites include a β spectrin, a cytoskeletal protein with still incompletely understood roles in neuron structure and function, as well as glutathione S-transferase, an enzyme with potential roles in the biochemical defence against toxicants.

The TRIM37 variant rs57141087 contributes to triple-negative breast cancer outcomes in Black women

Triple-negative breast cancer (TNBC) disproportionately affects younger Black women, who show more aggressive phenotypes and poorer outcomes than women of other racial identities. While the impact of socioenvironmental inequities within and beyond health systems is well documented, the genetic influence in TNBC-associated racial disparities remains elusive. Here, we report that cancer-free breast tissue from Black women expresses TRIM37 at a significantly higher level relative to White women. A reporter-based screen for regulatory variants identifies a non-coding risk variant rs57141087 in the 5’ gene upstream region of the TRIM37 locus with enhancer activity. Mechanistically, rs57141087 increases enhancer–promoter interactions through NRF1, resulting in stronger TRIM37 promoter activity. Phenotypically, high TRIM37 levels drive neoplastic transformations in immortalized breast epithelial cells. Finally, context-dependent TRIM37 expression reveals that early-stage TRIM37 levels affect the initiation and trajectory of breast cancer progression. Together, our results indicate a genotype-informed association of oncogenic TRIM37 with TNBC risk in Black women and implicate TRIM37 as a predictive biomarker to better identify patients at risk of aggressive TNBC.

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.

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.

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.

Multi-model GWAS reveals key loci for horticultural traits in reconstructed garden strawberry.

The cultivated garden strawberry (Fragaria × ananassa) has a rich history, originating from the hybridization of two wild octoploid strawberry species in the 18th century. Two-step reconstruction of Fragaria × ananassa through controlled crossings between pre-improved selections of its parental species is a promising approach for enriching the breeding germplasm of strawberry for wider adaptability. We created a population of reconstructed strawberry by hybridizing elite selections of F. virginiana and F. chiloensis. A replicated field experiment was conducted to evaluate the population's performance for eleven horticulturally important traits, over multiple years. Population structure analyses based on Fana-50 k SNP array data confirmed pedigree-based grouping of the progenies into four distinct groups. As complex traits are often influenced by environmental variables, and population structure can lead to spurious associations, we tested multiple genome-wide association study (GWAS) models. GWAS uncovered 39 quantitative trait loci (QTL) regions for eight traits distributed across twenty chromosomes, including 11 consistent and 28 putative QTLs. Candidate genes for traits including winter survival, flowering time, runnering vigor, and hermaphrodism were identified within the QTL regions. To our knowledge, this study marks the first comprehensive investigation of adaptive and horticultural traits in a large, multi-familial reconstructed strawberry population using SNP markers.

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.