A comprehensive genetic map of cytokine responses in Lyme borreliosis.

The incidence of Lyme borreliosis has risen, accompanied by persistent symptoms. The innate immune system and related cytokines are crucial in the host response and symptom development. We characterized cytokine production capacity before and after antibiotic treatment in 1,060 Lyme borreliosis patients. We observed a negative correlation between antibody production and IL-10 responses, as well as increased IL-1Ra responses in patients with disseminated disease. Genome-wide mapping the cytokine production allowed us to identify 34 cytokine quantitative trait loci (cQTLs), with 31 novel ones. We pinpointed the causal variant at the TLR1-6-10 locus and validated the regulation of IL-1Ra responses at transcritpome level using an independent cohort. We found that cQTLs contribute to Lyme borreliosis susceptibility and are relevant to other immune-mediated diseases. Our findings improve the understanding of cytokine responses in Lyme borreliosis and provide a genetic map of immune function as an expanded resource.

Transcriptome-wide association mapping provides insights into the genetic basis and candidate genes governing flowering, maturity and seed weight in rice bean (Vigna umbellata).

BACKGROUND: Rice bean (Vigna umbellata), an underrated legume, adapts to diverse climatic conditions with the potential to support food and nutritional security worldwide. It is used as a vegetable, minor food crop and a fodder crop, being a rich source of proteins, minerals, and essential fatty acids. However, little effort has been made to decipher the genetic and molecular basis of various useful traits in this crop. Therefore, we considered three economically important traits i.e., flowering, maturity and seed weight of rice bean and identified the associated candidate genes employing an associative transcriptomics approach on 100 diverse genotypes out of 1800 evaluated rice bean accessions from the Indian National Genebank. RESULTS: The transcriptomics-based genotyping of one-hundred diverse rice bean cultivars followed by pre-processing of genotypic data resulted in 49,271 filtered markers. The STRUCTURE, PCA and Neighbor-Joining clustering of 100 genotypes revealed three putative sub-populations. The marker-trait association analysis involving various genome-wide association study (GWAS) models revealed significant association of 82 markers on 48 transcripts for flowering, 26 markers on 22 transcripts for maturity and 22 markers on 21 transcripts for seed weight. The transcript annotation provided information on the putative candidate genes for the considered traits. The candidate genes identified for flowering include HSC80, P-II PsbX, phospholipid-transporting-ATPase-9, pectin-acetylesterase-8 and E3-ubiquitin-protein-ligase-RHG1A. Further, the WRKY1 and DEAD-box-RH27 were found to be associated with seed weight. Furthermore, the associations of PIF3 and pentatricopeptide-repeat-containing-gene with maturity and seed weight, and aldo-keto-reductase with flowering and maturity were revealed. CONCLUSION: This study offers insights into the genetic basis of key agronomic traits in rice bean, including flowering, maturity, and seed weight. The identified markers and associated candidate genes provide valuable resources for future exploration and targeted breeding, aiming to enhance the agronomic performance of rice bean cultivars. Notably, this research represents the first transcriptome-wide association study in pulse crop, uncovering the candidate genes for agronomically useful traits.

Unraveling the immunogenetic landscape of autism spectrum disorder: a comprehensive bioinformatics approach.

Background: Autism spectrum disorder (ASD) is a disease characterized by social disorder. Recently, the population affected by ASD has gradually increased around the world. There are great difficulties in diagnosis and treatment at present. Methods: The ASD datasets were obtained from the Gene Expression Omnibus database and the immune-relevant genes were downloaded from a previously published compilation. Subsequently, we used WGCNA to screen the modules related to the ASD and immune. We also choose the best combination and screen out the core genes from Consensus Machine Learning Driven Signatures (CMLS). Subsequently, we evaluated the genetic correlation between immune cells and ASD used GNOVA. And pleiotropic regions identified by PLACO and CPASSOC between ASD and immune cells. FUMA was used to identify pleiotropic regions, and expression trait loci (EQTL) analysis was used to determine their expression in different tissues and cells. Finally, we use qPCR to detect the gene expression level of the core gene. Results: We found a close relationship between neutrophils and ASD, and subsequently, CMLS identified a total of 47 potential candidate genes. Secondly, GNOVA showed a significant genetic correlation between neutrophils and ASD, and PLACO and CPASSOC identified a total of 14 pleiotropic regions. We annotated the 14 regions mentioned above and identified a total of 6 potential candidate genes. Through EQTL, we found that the CFLAR gene has a specific expression pattern in neutrophils, suggesting that it may serve as a potential biomarker for ASD and is closely related to its pathogenesis. Conclusions: In conclusion, our study yields unprecedented insights into the molecular and genetic heterogeneity of ASD through a comprehensive bioinformatics analysis. These valuable findings hold significant implications for tailoring personalized ASD therapies.

Exploring the interconnected between type 2 diabetes mellitus and nonalcoholic fatty liver disease: Genetic correlation and Mendelian randomization analysis.

Epidemiological and clinical studies have indicated a higher risk of nonalcoholic fatty liver disease (NAFLD) and type 2 diabetes mellitus (T2DM), implying a potentially shared genetic etiology, which is still less explored. Genetic links between T2DM and NAFLD were assessed using linkage disequilibrium score regression and pleiotropic analysis under composite null hypothesis. European GWAS data have identified shared genes, whereas SNP-level pleiotropic analysis under composite null hypothesis has explored pleiotropic loci. generalized gene-set analysis of GWAS data determines pleiotropic pathways and tissue enrichment using eQTL mapping to identify associated genes. Mendelian randomization analysis was used to investigate the causal relationship between NAFLD and T2DM. Linkage disequilibrium score regression analysis revealed a strong genetic correlation between T2DM and NAFLD, and identified 24 pleiotropic loci. These single-nucleotide polymorphisms are primarily involved in biosynthetic regulation, RNA biosynthesis, and pancreatic development. generalized gene-set analysis of GWAS data analysis revealed significant enrichment in multiple brain tissues. Gene mapping using these 3 methods led to the identification of numerous pleiotropic genes, with differences observed in liver and kidney tissues. These genes were mainly enriched in pancreas, brain, and liver tissues. The Mendelian randomization method indicated a significantly positive unidirectional causal relationship between T2DM and NAFLD. Our study identified a shared genetic structure between NAFLD and T2DM, providing new insights into the genetic pathogenesis and mechanisms of NAFLD and T2DM comorbidities.

Different control of resistance to two Colletotrichum orbiculare pathogenic races 0 and 1 in cucumber (Cucumis sativus L.).

KEY MESSAGE: Quantitative trait locus analysis identified independent novel loci in cucumbers responsible for resistance to races 0 and 1 of the anthracnose fungal pathogen Colletotrichum orbiculare. Cucumbers have been reported to be vulnerable to Colletotrichum orbiculare, causing anthracnose disease with significant yield loss under favorable conditions. The deployment of a single recessive Cssgr gene in cucumber breeding for anthracnose resistance was effective until a recent report on high-virulent strains infecting cucumbers in Japan conquering the resistance. QTL mapping was conducted to identify the resistance loci in the cucumber accession Ban Kyuri (G100) against C. orbiculare strains 104-T and CcM-1 of pathogenic races 0 and 1, respectively. A single dominant locus An5 was detected in the disease resistance hotspot on chromosome 5 for resistance to 104-T. Resistance to CcM-1 was governed by three loci with additive effects located on chromosomes 2 (An2) and 1 (An1.1 and An1.2). Molecular markers were developed based on variant calling between the corresponding QTL regions in the de novo assembly of the G100 genome and the publicly available cucumber genomes. Multiple backcrossed populations were deployed to fine-map An5 locus and narrow the region to approximately 222 kbp. Accumulation of An2 and An1.1 alleles displayed an adequate resistance to CcM-1 strain. This study provides functional molecular markers for pyramiding resistance loci that confer sufficient resistance against anthracnose in cucumbers.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus.

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

Copy number variation and elevated genetic diversity at immune trait loci in Atlantic and Pacific herring.

BACKGROUND: Genome-wide comparisons of populations are widely used to explore the patterns of nucleotide diversity and sequence divergence to provide knowledge on how natural selection and genetic drift affect the genome. In this study we have compared whole-genome sequencing data from Atlantic and Pacific herring, two sister species that diverged about 2 million years ago, to explore the pattern of genetic differentiation between the two species. RESULTS: The genome comparison of the two species revealed high genome-wide differentiation but with islands of remarkably low genetic differentiation, as measured by an FST analysis. However, the low FST observed in these islands is not caused by low interspecies sequence divergence (dxy) but rather by exceptionally high estimated intraspecies nucleotide diversity (π). These regions of low differentiation and elevated nucleotide diversity, termed high-diversity regions in this study, are not enriched for repeats but are highly enriched for immune-related genes. This enrichment includes genes from both the adaptive immune system, such as immunoglobulin, T-cell receptor and major histocompatibility complex genes, as well as a substantial number of genes with a role in the innate immune system, e.g. novel immune-type receptor, tripartite motif and tumor necrosis factor receptor genes. Analysis of long-read based assemblies from two Atlantic herring individuals revealed extensive copy number variation in these genomic regions, indicating that the elevated intraspecies nucleotide diversities were partially due to the cross-mapping of short reads. CONCLUSIONS: This study demonstrates that copy number variation is a characteristic feature of immune trait loci in herring. Another important implication is that these loci are blind spots in classical genome-wide screens for genetic differentiation using short-read data, not only in herring, likely also in other species harboring qualitatively similar variation at immune trait loci. These loci stood out in this study because of the relatively high genome-wide baseline for FST values between Atlantic and Pacific herring.

RNA-Seq Bulked Segregant Analysis of an Exotic B. napus ssp. napobrassica (Rutabaga) F2 Population Reveals Novel QTLs for Breeding Clubroot-Resistant Canola.

In this study, a rutabaga (Brassica napus ssp. napobrassica) donor parent FGRA106, which exhibited broad-spectrum resistance to 17 isolates representing 16 pathotypes of Plasmodiophora brassicae, was used in genetic crosses with the susceptible spring-type canola (B. napus ssp. napus) accession FG769. The F2 plants derived from a clubroot-resistant F1 plant were screened against three P. brassicae isolates representing pathotypes 3A, 3D, and 3H. Chi-square (χ2) goodness-of-fit tests indicated that the F2 plants inherited two major clubroot resistance genes from the CR donor FGRA106. The total RNA from plants resistant (R) and susceptible (S) to each pathotype were pooled and subjected to bulked segregant RNA-sequencing (BSR-Seq). The analysis of gene expression profiles identified 431, 67, and 98 differentially expressed genes (DEGs) between the R and S bulks. The variant calling method indicated a total of 12 (7 major + 5 minor) QTLs across seven chromosomes. The seven major QTLs included: BnaA5P3A.CRX1.1, BnaC1P3H.CRX1.2, and BnaC7P3A.CRX1.1 on chromosomes A05, C01, and C07, respectively; and BnaA8P3D.CRX1.1, BnaA8P3D.RCr91.2/BnaA8P3H.RCr91.2, BnaA8P3H.Crr11.3/BnaA8P3D.Crr11.3, and BnaA8P3D.qBrCR381.4 on chromosome A08. A total of 16 of the DEGs were located in the major QTL regions, 13 of which were on chromosome C07. The molecular data suggested that clubroot resistance in FGRA106 may be controlled by major and minor genes on both the A and C genomes, which are deployed in different combinations to confer resistance to the different isolates. This study provides valuable germplasm for the breeding of clubroot-resistant B. napus cultivars in Western Canada.

New Insights into the Genetic Basis of Lysine Accumulation in Rice Revealed by Multi-Model GWAS.

Lysine is an essential amino acid that cannot be synthesized in humans. Rice is a global staple food for humans but has a rather low lysine content. Identification of the quantitative trait nucleotides (QTNs) and genes underlying lysine content is crucial to increase lysine accumulation. In this study, five grain and three leaf lysine content datasets and 4,630,367 single nucleotide polymorphisms (SNPs) of 387 rice accessions were used to perform a genome-wide association study (GWAS) by ten statistical models. A total of 248 and 71 common QTNs associated with grain/leaf lysine content were identified. The accuracy of genomic selection/prediction RR-BLUP models was up to 0.85, and the significant correlation between the number of favorable alleles per accession and lysine content was up to 0.71, which validated the reliability and additive effects of these QTNs. Several key genes were uncovered for fine-tuning lysine accumulation. Additionally, 20 and 30 QTN-by-environment interactions (QEIs) were detected in grains/leaves. The QEI-sf0111954416 candidate gene LOC_Os01g21380 putatively accounted for gene-by-environment interaction was identified in grains. These findings suggested the application of multi-model GWAS facilitates a better understanding of lysine accumulation in rice. The identified QTNs and genes hold the potential for lysine-rich rice with a normal phenotype.

Identification of Quantitative Trait Loci Controlling Root Morphological Traits in an Interspecific Soybean Population Using 2D Imagery Data.

Roots are the hidden and most important part of plants. They serve as stabilizers and channels for uptaking water and nutrients and play a crucial role in the growth and development of plants. Here, two-dimensional image data were used to identify quantitative trait loci (QTL) controlling root traits in an interspecific mapping population derived from a cross between wild soybean 'PI366121' and cultivar 'Williams 82'. A total of 2830 single-nucleotide polymorphisms were used for genotyping, constructing genetic linkage maps, and analyzing QTLs. Forty-two QTLs were identified on twelve chromosomes, twelve of which were identified as major QTLs, with a phenotypic variation range of 36.12% to 39.11% and a logarithm of odds value range of 12.01 to 17.35. Two significant QTL regions for the average diameter, root volume, and link average diameter root traits were detected on chromosomes 3 and 13, and both wild and cultivated soybeans contributed positive alleles. Six candidate genes, Glyma.03G027500 (transketolase/glycoaldehyde transferase), Glyma.03G014500 (dehydrogenases), Glyma.13G341500 (leucine-rich repeat receptor-like protein kinase), Glyma.13G341400 (AGC kinase family protein), Glyma.13G331900 (60S ribosomal protein), and Glyma.13G333100 (aquaporin transporter) showed higher expression in root tissues based on publicly available transcriptome data. These results will help breeders improve soybean genetic components and enhance soybean root morphological traits using desirable alleles from wild soybeans.

African ancestry-derived APOL1 risk genotypes show proximal epigenetic associations.

Apolipoprotein L1 (APOL1) coding variants, termed G1 and G2, are established genetic risk factors for a growing spectrum of diseases, including kidney disease, in individuals of African ancestry. Evidence suggests that the risk variants, which show a recessive mode of inheritance, lead to toxic gain-of-function changes of the APOL1 protein. Disease occurrence and presentation vary, likely due to modifiers or second hits. To understand the role of the epigenetic landscape in relation to APOL1 risk variants, we performed methylation quantitative trait locus (meQTL) analysis to identify differentially methylated CpGs influenced by APOL1 risk variants in 611 African American individuals. We identified five CpGs that were significantly associated with APOL1 risk alleles in discovery and replication studies, and one CpG-APOL1 association was independent of other genomic variants. Our study highlights proximal DNA methylation alterations that may help explain the variable disease risk and clinical manifestation of APOL1 variants.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle.

BACKGROUND: Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle. RESULTS: We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits. CONCLUSIONS: This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

Selection of rice breeding lines for resistance to biotic and abiotic stresses.

Rice (Oryza sativa L.) grown in many countries around the world with different climatic conditions and a huge number of environmental stresses, both biotic (fungi, bacteria, viruses, insects) and abiotic (cold, drought, salinity) limit rice productivity. In this regard, breeders and scientists are trying to create rice lines that are resistant to multiple stresses. The aim of this work was to screen and select cold and blast resistant rice breeding lines (RBLs) using molecular markers. Molecular screening of RBLs and parental varieties to cold tolerance was carried out using markers RM24545, RM1377, RM231 and RM569 associated with QTLs (qPSST-3, qPSST-7, qPSST-9). It was discovered that the presence of three QTLs characterizes the cold resistance of studied genotypes, and the absence of one of them leads to cold sensitivity. As a result, 21 cold-resistant out of the 28 studied RBLs were identified. These cold resistant 21 RBLs were further tested to blast resistance using markers Pi-ta, Pita3, Z56592, 195R-1, NMSMPi9-1, TRS26, Pikh MAS, MSM6, 9871.T7E2b, RM224 and RM1233. It was revealed that 16 RBLs from 21 studied lines contain 5-6 blast resistance genes. In accordance with the blast resistance strategy, the presence of 5 or more genes ensures the formation of stable resistance to Magnaporthe oryzae. Thus, 16 lines resistant to multiple stresses, such as cold and blast disease were developed. It should be noted that 6 of these selected lines are high-yielding, which is very important in rice breeding program. These RBLs can be used in breeding process as starting lines, germplasm exchange as a source of resistant genes for the development of new rice varieties resistant to multiple stress factors.

Genome-wide association study and genomic selection of spike-related traits in bread wheat.

KEY MESSAGE: Identification of 337 stable MTAs for wheat spike-related traits improved model accuracy, and favorable alleles of MTA259 and MTA64 increased grain weight and yield per plant. Wheat (Triticum aestivum L.) is one of the three primary global, staple crops. Improving spike-related traits in wheat is crucial for optimizing spike and plant morphology, ultimately leading to increased grain yield. Here, we performed a genome-wide association study using a dataset of 24,889 high-quality unique single-nucleotide polymorphisms (SNPs) and phenotypic data from 314 wheat accessions across eight diverse environments. In total, 337 stable and significant marker-trait associations (MTAs) related to spike-related traits were identified. MTA259 and MTA64 were consistently detected in seven and six environments, respectively. The presence of favorable alleles associated with MTA259 and MTA64 significantly reduced wheat spike exsertion length and spike length, while enhancing thousand kernel weight and yield per plant. Combined gene expression and network analyses identified TraesCS6D03G0692300 and TraesCS6D03G0692700 as candidate genes for MTA259 and TraesCS2D03G0111700 and TraesCS2D03G0112500 for MTA64. The identified MTAs significantly improved the prediction accuracy of each model compared with using all the SNPs, and the random forest model was optimal for genome selection. Additionally, the eight stable and major MTAs, including MTA259, MTA64, MTA66, MTA94, MTA110, MTA165, MTA180, and MTA164, were converted into cost-effective and efficient detection markers. This study provided valuable genetic resources and reliable molecular markers for wheat breeding programs.

Genetic tapestry of Capsicum fruit colors: a comparative analysis of four cultivated species.

KEY MESSAGE: Genome-wide association study of color spaces across the four cultivated Capsicum spp. revealed a shared set of genes influencing fruit color, suggesting mechanisms and pathways across Capsicum species are conserved during the speciation. Notably, Cytochrome P450 of the carotenoid pathway, MYB transcription factor, and pentatricopeptide repeat-containing protein are the major genes responsible for fruit color variation across the Capsicum species. Peppers (Capsicum spp.) rank among the most widely consumed spices globally. Fruit color, serving as a determinant for use in food colorants and cosmeceuticals and an indicator of nutritional contents, significantly influences market quality and price. Cultivated Capsicum species display extensive phenotypic diversity, especially in fruit coloration. Our study leveraged the genetic variance within four Capsicum species (Capsicum baccatum, Capsicum chinense, Capsicum frutescens, and Capsicum annuum) to elucidate the genetic mechanisms driving color variation in peppers and related Solanaceae species. We analyzed color metrics and chromatic attributes (Red, Green, Blue, L*, a*, b*, Luminosity, Hue, and Chroma) on samples cultivated over six years (2015-2021). We resolved genomic regions associated with fruit color diversity through the sets of SNPs obtained from Genotyping by Sequencing (GBS) and genome-wide association study (GWAS) with a Multi-Locus Mixed Linear Model (MLMM). Significant SNPs with FDR correction were identified, within the Cytochrome P450, MYB-related genes, Pentatricopeptide repeat proteins, and ABC transporter family were the most common among the four species, indicating comparative evolution of fruit colors. We further validated the role of a pentatricopeptide repeat-containing protein (Chr01:31,205,460) and a cytochrome P450 enzyme (Chr08:45,351,919) via competitive allele-specific PCR (KASP) genotyping. Our findings advance the understanding of the genetic underpinnings of Capsicum fruit coloration, with developed KASP assays holding potential for applications in crop breeding and aligning with consumer preferences. This study provides a cornerstone for future research into exploiting Capsicum's diverse fruit color variation.

Natural variations of heterosis-related allele-specific expression genes in promoter regions lead to allele-specific expression in maize.

BACKGROUND: Heterosis has successfully enhanced maize productivity and quality. Although significant progress has been made in delineating the genetic basis of heterosis, the molecular mechanisms underlying its genetic components remain less explored. Allele-specific expression (ASE), the imbalanced expression between two parental alleles in hybrids, is increasingly being recognized as a factor contributing to heterosis. ASE is a complex process regulated by both epigenetic and genetic variations in response to developmental and environmental conditions. RESULTS: In this study, we explored the differential characteristics of ASE by analyzing the transcriptome data of two maize hybrids and their parents under four light conditions. On the basis of allele expression patterns in different hybrids under various conditions, ASE genes were divided into three categories: bias-consistent genes involved in basal metabolic processes in a functionally complementary manner, bias-reversal genes adapting to the light environment, and bias-specific genes maintaining cell homeostasis. We observed that 758 ASE genes (ASEGs) were significantly overlapped with heterosis quantitative trait loci (QTLs), and high-frequency variations in the promoter regions of heterosis-related ASEGs were identified between parents. In addition, 10 heterosis-related ASEGs participating in yield heterosis were selected during domestication. CONCLUSIONS: The comprehensive analysis of ASEGs offers a distinctive perspective on how light quality influences gene expression patterns and gene-environment interactions, with implications for the identification of heterosis-related ASEGs to enhance maize yield.

Common variation in a long non-coding RNA gene modulates variation of circulating TGF-ß2 levels in metastatic colorectal cancer patients (Alliance).

BACKGROUND: Herein, we report results from a genome-wide study conducted to identify protein quantitative trait loci (pQTL) for circulating angiogenic and inflammatory protein markers in patients with metastatic colorectal cancer (mCRC). The study was conducted using genotype, protein marker, and baseline clinical and demographic data from CALGB/SWOG 80405 (Alliance), a randomized phase III study designed to assess outcomes of adding VEGF or EGFR inhibitors to systemic chemotherapy in mCRC patients. Germline DNA derived from blood was genotyped on whole-genome array platforms. The abundance of protein markers was quantified using a multiplex enzyme-linked immunosorbent assay from plasma derived from peripheral venous blood collected at baseline. A robust rank-based method was used to assess the statistical significance of each variant and protein pair against a strict genome-wide level. A given pQTL was tested for validation in two external datasets of prostate (CALGB 90401) and pancreatic cancer (CALGB 80303) patients. Bioinformatics analyses were conducted to further establish biological bases for these findings. RESULTS: The final analysis was carried out based on data from 540,021 common typed genetic variants and 23 protein markers from 869 genetically estimated European patients with mCRC. Correcting for multiple testing, the analysis discovered a novel cis-pQTL in LINC02869, a long non-coding RNA gene, for circulating TGF-ß2 levels (rs11118119; AAF = 0.11; P-value < 1.4e-14). This finding was validated in a cohort of 538 prostate cancer patients from CALGB 90401 (AAF = 0.10, P-value < 3.3e-25). The analysis also validated a cis-pQTL we had previously reported for VEGF-A in advanced pancreatic cancer, and additionally identified trans-pQTLs for VEGF-R3, and cis-pQTLs for CD73. CONCLUSIONS: This study has provided evidence of a novel cis germline genetic variant that regulates circulating TGF-ß2 levels in plasma of patients with advanced mCRC and prostate cancer. Moreover, the validation of previously identified pQTLs for VEGF-A, CD73, and VEGF-R3, potentiates the validity of these associations.

The impact of exercise on gene regulation in association with complex trait genetics.

Endurance exercise training is known to reduce risk for a range of complex diseases. However, the molecular basis of this effect has been challenging to study and largely restricted to analyses of either few or easily biopsied tissues. Extensive transcriptome data collected across 15 tissues during exercise training in rats as part of the Molecular Transducers of Physical Activity Consortium has provided a unique opportunity to clarify how exercise can affect tissue-specific gene expression and further suggest how exercise adaptation may impact complex disease-associated genes. To build this map, we integrate this multi-tissue atlas of gene expression changes with gene-disease targets, genetic regulation of expression, and trait relationship data in humans. Consensus from multiple approaches prioritizes specific tissues and genes where endurance exercise impacts disease-relevant gene expression. Specifically, we identify a total of 5523 trait-tissue-gene triplets to serve as a valuable starting point for future investigations [Exercise; Transcription; Human Phenotypic Variation].

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

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

Complex genetic architecture of the chicken Growth1 QTL region.

The genetic complexity of polygenic traits represents a captivating and intricate facet of biological inheritance. Unlike Mendelian traits controlled by a single gene, polygenic traits are influenced by multiple genetic loci, each exerting a modest effect on the trait. This cumulative impact of numerous genes, interactions among them, environmental factors, and epigenetic modifications results in a multifaceted architecture of genetic contributions to complex traits. Given the well-characterized genome, diverse traits, and range of genetic resources, chicken (Gallus gallus) was employed as a model organism to dissect the intricate genetic makeup of a previously identified major Quantitative Trait Loci (QTL) for body weight on chromosome 1. A multigenerational advanced intercross line (AIL) of 3215 chickens whose genomes had been sequenced to an average of 0.4x was analyzed using genome-wide association study (GWAS) and variance-heterogeneity GWAS (vGWAS) to identify markers associated with 8-week body weight. Additionally, epistatic interactions were studied using the natural and orthogonal interaction (NOIA) model. Six genetic modules, two from GWAS and four from vGWAS, were strongly associated with the studied trait. We found evidence of both additive- and non-additive interactions between these modules and constructed a putative local epistasis network for the region. Our screens for functional alleles revealed a missense variant in the gene ribonuclease H2 subunit B (RNASEH2B), which has previously been associated with growth-related traits in chickens and Darwin's finches. In addition, one of the most strongly associated SNPs identified is located in a non-coding region upstream of the long non-coding RNA, ENSGALG00000053256, previously suggested as a candidate gene for regulating chicken body weight. By studying large numbers of individuals from a family material using approaches to capture both additive and non-additive effects, this study advances our understanding of genetic complexities in a highly polygenic trait and has practical implications for poultry breeding and agriculture.