Search

Abstract Background Single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) have led to groundbreaking advancements in life sciences. To develop bioinformatics tools for scRNA-seq and SRT data and perform unbiased benchmarks, data simulation has been widely adopted by providing explicit ground truth and generating customized datasets. However, the performance of simulation methods under multiple scenarios has not been comprehensively assessed, making it challenging to choose suitable methods without practical guidelines. Results We systematically evaluated 49 simulation methods developed for scRNA-seq and/or SRT data in terms of accuracy, functionality, scalability, and usability using 152 reference datasets derived from 24 platforms. SRTsim, scDesign3, ZINB-WaVE, and scDesign2 have the best accuracy performance across various platforms. Unexpectedly, some methods tailored to scRNA-seq data have potential compatibility for simulating SRT data. Lun, SPARSim, and scDesign3-tree outperform other methods under corresponding simulation scenarios. Phenopath, Lun, Simple, and MFA yield high scalability scores but they cannot generate realistic simulated data. Users should consider the trade-offs between method accuracy and scalability (or functionality) when making decisions. Additionally, execution errors are mainly caused by failed parameter estimations and appearance of missing or infinite values in calculations. We provide practical guidelines for method selection, a standard pipeline Simpipe ( https://github.com/duohongrui/simpipe ; https://doi.org/10.5281/zenodo.11178409 ), and an online tool Simsite ( https://www.ciblab.net/software/simshiny/ ) for data simulation. Conclusions No method performs best on all criteria, thus a good-yet-not-the-best method is recommended if it solves problems effectively and reasonably. Our comprehensive work provides crucial insights for developers on modeling gene expression data and fosters the simulation process for users.

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets.

Erythritol has shown excellent insecticidal performance against a wide range of insect species, but the molecular mechanism by which it causes insect mortality and sterility is not fully understood. The mortality and sterility of Drosophila melanogaster were assessed after feeding with 1M erythritol for 72 h and 96 h, and gene expression profiles were further compared through RNA sequencing. Enrichment analysis of GO and KEGG revealed that expressions of the adipokinetic hormone gene (Akh), amylase gene (Amyrel), α-glucosidase gene (Mal-B1/2, Mal-A1-4, Mal-A7/8), and triglyceride lipase gene (Bmm) were significantly up-regulated, while insulin-like peptide genes (Dilp2, Dilp3 and Dilp5) were dramatically down-regulated. Seventeen genes associated with eggshell assembly, including Dec-1 (down 315-fold), Vm26Ab (down 2014-fold) and Vm34Ca (down 6034-fold), were significantly down-regulated or even showed no expression. However, there were no significant differences in the expression of three diuretic hormone genes (DH44, DH31, CAPA) and eight aquaporin genes (Drip, Big brain, AQP, Eglp1, Eglp2, Eglp3, Eglp4 and Prip) involved in osmolality regulation (all p value > 0.05). We concluded that erythritol, a competitive inhibitor of α-glucosidase, severely reduced substrates and enzyme binding, inhibiting effective carbohydrate hydrolysis in the midgut and eventually causing death due to energy deprivation. It was clear that Drosophila melanogaster did not die from the osmolality of the hemolymph. Our findings elucidate the molecular mechanism underlying the mortality and sterility in Drosophila melanogaster induced by erythritol feeding. It also provides an important theoretical basis for the application of erythritol as an environmentally friendly pesticide.

Tumor heterogeneity remains a major challenge for disease subtyping, risk stratification, and accurate clinical management. Exosome-based liquid biopsy can effectively overcome the limitations of tissue biopsy, achieving minimal invasion, multi-point dynamic monitoring, and good prognosis assessment, and has broad clinical prospects. However, there is still lacking comprehensive analysis of tumor-derived exosome (TDE)-based stratification of risk patients and prognostic assessment for breast cancer with systematic dissection of biological heterogeneity. In this study, the robust corroborative analysis for biomarker discovery (RCABD) strategy was used for the identification of exosome molecules, differential expression verification, risk prediction modeling, heterogenous dissection with multi-ome (6101 molecules), our ExoBCD database (306 molecules), and 53 independent studies (481 molecules). Our results showed that a 10-molecule exosome-derived signature (exoSIG) could successfully fulfill breast cancer risk stratification, making it a novel and accurate exosome prognostic indicator (Cox P = 9.9E-04, HR = 3.3, 95% CI 1.6–6.8). Interestingly, HLA-DQB2 and COL17A1, closely related to tumor metastasis, achieved high performance in prognosis prediction (86.35% contribution) and accuracy (Log-rank P = 0.028, AUC = 85.42%). With the combined information of patient age and tumor stage, they formed a bimolecular risk signature (Clinmin-exoSIG) and a convenient nomogram as operable tools for clinical applications. In conclusion, as an extension of ExoBCD, this study conducted systematic analyses to identify prognostic multi-molecular panel and risk signature, stratify patients and dissect biological heterogeneity based on breast cancer exosomes from a multi-omics perspective. Our results provide an important reference for in-depth exploration of the ''biological heterogeneity - risk stratification - prognosis prediction''.

Abstract Background The olfactory system plays a crucial role in regulating insect behaviors. The detection of odorants is mainly mediated by various odorant receptors (ORs) that are expressed in the dendrites of olfactory neurons of chemosensilla. Anopheles sinensis is a major malaria vector in Eastern Asia and its genome has recently been successfully sequenced and annotated. In this study, we present genome-wide identification and expression profiling of OR genes in different chemosensory tissues of An. sinensis. Methods The OR genes were identified using the available genome sequences of An. sinensis. A series of bioinformatics analyses were conducted to investigate the structure, genome distribution, selective pressure and phylogenetic relationships of OR genes, the conserved domains and specific functional sites in the OR amino acid sequences. The expression levels of OR genes were analyzed from transcriptomic data from An. sinensis antennae, proboscis and maxillary palps of both sexes. Results A total of 59 putative OR genes have been identified and characterized in An. sinensis. This number is significantly less than that in An. gambiae. Whether this difference is caused by the contraction or expansion of OR genes after divergence of the two species remains unknown. The RNA-seq analysis showed that AsORs have obvious tissue- and sex-specific expression patterns. Most AsORs are highly expressed in the antennae and the expression pattern and number of AsORs expressed in antennae are similar in males and females. However, the relative levels of AsOR transcripts are much higher in female antennae than in male antennae, which indicates that the odor sensitivity is likely to be increased in female mosquitoes. Based on the expression patterns and previous studies, we have speculated on the functions of some OR genes but this needs to be validated by further behavioral, molecular and electrophysiological studies. Further studies are necessary to compare the olfactory-driven behaviors and identify receptors that respond strongly to components of human odors that may act in the process of human recognition. Conclusions This is the first genome-wide analysis of the entire repertoire of OR genes in An. sinensis. Characterized features and profiled expression patterns of ORs suggest their involvement in the odorous reception of this species. Our findings provide a basis for further research on the functions of OR genes and additional genetic and behavioral targets for more sustainable management of An. sinensis in the future. Graphical Abstract

Phosphatidylethanolamine-binding proteins (PEBP) family plays important roles in regulating plant flowering time and morphogenesis. However, geneme-wide identification and functional analysis of PEBP genes in the rigorous short-day plant Perilla frutescens (PfPEBP) have not been studied. In this study, 10 PfPEBP were identified and divided into three subfamilies based on their phylogenetic relationships: FT-like, TFL1-like and MFT-like. Gene structure analysis showed that all PfPEBP genes contain 4 exons and 3 introns. Motifs DPDxP and GIHR essential for anion-binding activity are highly conserved in PfPEBP. A large number of light-responsive elements were detected in promoter regions of PfPEBP. Gene expression of PfFT1 exhibited a diurnal rhythm. It was highly expressed in leaves under the short-day photoperiod, but higher in flowers and seeds under the long-day photoperiod. Overexpression of PfFT1 in Arabidopsis thaliana not only promoted early flowering of Col-0 or Ler, but also rescued the late flowering phenotype of ft-1 mutant. We concluded that PfFT1 promotes early flowering by regulating the expression of flowering-related genes AtAP1, AtLFY, AtFUL and AtSOC1. In conclusion, our results provided valuable information for elucidating the functions of PfPEBP in P. frutescens and shed light on the promoting effect of PfFT1 on flowering.

BackgroundMelanoma is a highly aggressive skin cancer with a poor prognosis and mortality. Immune checkpoint blockade (ICB) therapy (e.g., anti-PD-1 therapy) has opened a new horizon in melanoma treatment, but some patients present a non-responsive state. Cancer-associated fibroblasts (CAFs) make up the majority of stromal cells in the tumor microenvironment (TME) and have an important impact on the response to immunotherapy. There is still a lack of identification of CAFs-related predictors for anti-PD-1 therapy, although the establishment of immunotherapy biomarkers is well underway. This study aims to explore the potential CAFs-related gene panel for predicting the response to anti-PD-1 therapy in melanoma patients and elucidating their potential effect on TME.MethodsThree gene expression datasets from melanoma patients without anti-PD-1 treatment, in a total of 87 samples, were downloaded from Gene Expression Omnibus (GEO) as the discovery sets (GSE91061) and validation sets (GSE78220 and GSE122220). The CAFs-related module genes were identified from the discovery sets by weighted gene co-expression network analysis (WGCNA). Concurrently, we utilized differential gene analysis on the discovery set to obtain differentially expressed genes (DEGs). Then, CAFs-related key genes were screened with the intersection of CAFs-related module genes and DEGs, succeeded by supervised machine learning-based identification. As a consequence of expression analysis, gene set enrichment analysis, survival analysis, staging analysis, TME analysis, and correlation analysis, the multidimensional systematic characterizations of the key genes were uncovered. The diagnostic performance of the CAFs-related gene panel was assessed by receiver operating characteristic (ROC) curves in the validation sets. Eventually, the CAFs-related gene panel was verified by the expression from the single-cell analysis.ResultsThe six-gene panel associated with CAFs were finally identified for predicting the response to anti-PD-1 therapy, including CDK14, SYNPO2, TCF4, GJA1, CPXM1, and TFPI. The multigene panel demonstrated excellent combined diagnostic performance with the area under the curve of ROC reaching 90.5 and 75.4% ~100% in the discovery and validation sets, respectively.ConclusionConfirmed by clinical treatment outcomes, the identified CAFs-related genes can be used as a promising biomarker panel for prediction to anti-PD-1 therapy response, which may serve as new immunotherapeutic targets to improve survival outcomes of melanoma patients.

Insect-associated microorganisms play important roles in the health and development of insects. This study aimed to investigate the similarities and differences in bacterial community structure and composition between the larval gut of Osmia excavata, nest soil, and brood provision from the nest tube. We sequenced larvae gut and their environments’ microorganisms of O. excavata from four locations based on full-length 16S rRNA gene amplicons. The results showed 156, 280, and 366 bacterial OTUs from gut, brood provision, and nest soil, respectively, and three groups shared 131 bacterial OTUs. In the gut, the top two dominant bacteria were Sodalis praecaptivus (68.99%), Lactobacillus micheneri (17.95%). In the brood provision, the top two dominant bacteria were S. praecaptivus (26.66%), Acinetobacter nectaris (13.05%), and in the nest soil, the two most abundant bacteria were Gaiella occulta (4.33%), Vicinamibacter silvestris (3.88%). There were significant differences in diversity between the brood provision groups and the nest soil groups, respectively. Three of the four locations did not differ for gut microbial diversity. Bacteria similar to other solitary bees also existed in the gut of the larvae. Results indicated when the habitat environments were similar, the bacterial community diversity of the gut of O. excavata was similar, despite significant differences among brood provisions and soils, respectively.

To meet the pollination need of economic crops, Osmia excavata has been successfully used to improve the pollination efficiency of Rosaceae and Brassicaceae plants. As a widely used pollinator of economic crops, a systematic study of flower-visiting species and diversities of O. excavata stocked in China was not found. To investigate the foraging pollen species and diversities of O. excavata, beebread from 20 experimental plots in China was collected by the trap-nesting method and analyzed by DNA metabarcoding technology. A total of 26 pollen plants in 14 genera and nine families were identified. A further analysis showed that the richness and abundance of the wild flowering plants in orchards and farmlands were lower than those in the nearby semi-natural habitats. The favorite pollen comes from economic crops apple and rape and wild flowering plants Juncus interior, Rosa gymnocarpa, and Rosa laevigata. Through a diversity index analysis, it was found that the Anhui region has the highest pollen plant diversity, while the Liaoning region has the lowest. Our results can provide a basis for flower-visiting species and diversities of O. excavata.

Abstract Background Perilla frutescens is well known for its high α-linolenic acid (ALA) accumulation in seeds and medicinal values as well as a source of edible and general-purpose oils. However, the regulatory mechanisms of the biosynthesis of fatty acid in its seeds remain poorly understood due to the lacking of sequenced genome. For better understanding the regulation of lipid metabolism and further increase its oil content or modify oil composition, time-course transcriptome and lipid composition analyses were performed. Results Analysis of fatty acid content and composition showed that the α-linolenic acid and oleic acid accumulated rapidly from 5 DAF to 15 DAF and then kept relatively stable. However, the amount of palmitic acid and linoleic acid decreased quickly from 5 DAF to 15DAF. No significant variation of stearic acid content was observed from 5 DAF to 25DAF. Our transcriptome data analyses revealed that 110,176 unigenes were generated from six seed libraries at 5, 10, 20 DAF. Of these, 53 (31 up, 22 down) and 653 (259 up, 394 down) genes showed temporal and differentially expression during the seed development in 5 DAF vs 10 DAF, 20 vs 10 DAF, respectively. The differentially expressed genes were annotated and found to be involved in distinct functional categories and metabolic pathways. Deep mining of transcriptome data led to the identification of key genes involved in fatty acid and triacylglycerol biosynthesis and metabolism. Thirty seven members of transcription factor family AP2, B3 and NFYB putatively involved in oil synthesis and deposition were differentially expressed during seed development. The results of qRT-PCR for selected genes showed a strong positive correlation with the expression abundance measured in RNA-seq analysis. Conclusions The present study provides valuable genomic resources for characterizing Perilla seed gene expression at the transcriptional level and will extend our understanding of the complex molecular and cellular events of oil biosynthesis and accumulation in oilseed crops.

The complete mitogenome of Lasioglossum affine (Hymenoptera: Halictidae) was sequenced and analyzed. The whole mitogenome is 17,352 bp (AT%=84.1%) and encodes 37 typical eukaryotic mitochondrial genes, including 13 protein-coding genes (PCGs), 22 tRNAs, two rRNAs, and an AT-rich region. Further analysis found three gene rearrangements, where trn I-Q-M → trn M-I-Q, trn W-C-Y → trn C-W-Y, and trn K-D → trn D-K were shuffled. The phylogenetic relationships of 19 species of Hymenoptera were established using maximum-likelihood method based on 13 concatenated PCGs. The result showed that Lasioglossum affine is a sister of Lasioglossum sp. SJW-2017.

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.

For accurate gene expression quantification, normalization of gene expression data against reliable reference genes is required. It is known that the expression levels of commonly used reference genes vary considerably under different experimental conditions, and therefore, their use for data normalization is limited. In this study, an unbiased identification of reference genes in Caenorhabditis elegans was performed based on 145 microarray datasets (2296 gene array samples) covering different developmental stages, different tissues, drug treatments, lifestyle, and various stresses. As a result, thirteen housekeeping genes (rps-23, rps-26, rps-27, rps-16, rps-2, rps-4, rps-17, rpl-24.1, rpl-27, rpl-33, rpl-36, rpl-35, and rpl-15) with enhanced stability were comprehensively identified by using six popular normalization algorithms and RankAggreg method. Functional enrichment analysis revealed that these genes were significantly overrepresented in GO terms or KEGG pathways related to ribosomes. Validation analysis using recently published datasets revealed that the expressions of newly identified candidate reference genes were more stable than the commonly used reference genes. Based on the results, we recommended using rpl-33 and rps-26 as the optimal reference genes for microarray and rps-2 and rps-4 for RNA-sequencing data validation. More importantly, the most stable rps-23 should be a promising reference gene for both data types. This study, for the first time, successfully displays a large-scale microarray data driven genome-wide identification of stable reference genes for normalizing gene expression data and provides a potential guideline on the selection of universal internal reference genes in C. elegans, for quantitative gene expression analysis.

Asthma is a common chronic airway disease worldwide. Due to its clinical and genetic heterogeneity, the cellular and molecular processes in asthma are highly complex and relatively unknown. To discover novel biomarkers and the molecular mechanisms underlying asthma, several studies have been conducted by focusing on gene expression patterns in epithelium through microarray analysis. However, few robust specific biomarkers were identified and some inconsistent results were observed. Therefore, it is imperative to conduct a robust analysis to solve these problems. Herein, an integrated gene expression analysis of ten independent, publicly available microarray data of bronchial epithelial cells from 348 asthmatic patients and 208 healthy controls was performed. As a result, 78 up- and 75 down-regulated genes were identified in bronchial epithelium of asthmatics. Comprehensive functional enrichment and pathway analysis revealed that response to chemical stimulus, extracellular region, pathways in cancer, and arachidonic acid metabolism were the four most significantly enriched terms. In the protein-protein interaction network, three main communities associated with cytoskeleton, response to lipid, and regulation of response to stimulus were established, and the most highly ranked 6 hub genes (up-regulated CD44, KRT6A, CEACAM5, SERPINB2, and down-regulated LTF and MUC5B) were identified and should be considered as new biomarkers. Pathway cross-talk analysis highlights that signaling pathways mediated by IL-4/13 and transcription factor HIF-1α and FOXA1 play crucial roles in the pathogenesis of asthma. Interestingly, three chemicals, polyphenol catechin, antibiotic lomefloxacin, and natural alkaloid boldine, were predicted and may be potential drugs for asthma treatment. Taken together, our findings shed new light on the common molecular pathogenesis mechanisms of asthma and provide theoretical support for further clinical therapeutic studies.

Microsporidia have attracted considerable attention because they infect a wide range of hosts, from invertebrates to vertebrates, and cause serious human diseases and major economic losses in the livestock industry. There are no prospective drugs to counteract this pathogen. Eukaryotic protein kinases (ePKs) play a central role in regulating many essential cellular processes and are therefore potential drug targets. In this study, a comprehensive summary and comparative analysis of the protein kinases in four microsporidia—Enterocytozoon bieneusi, Encephalitozoon cuniculi, Nosema bombycis and Nosema ceranae—was performed. The results show that there are 34 ePKs and 4 atypical protein kinases (aPKs) in E. bieneusi, 29 ePKs and 6 aPKs in E. cuniculi, 41 ePKs and 5 aPKs in N. bombycis, and 27 ePKs and 4 aPKs in N. ceranae. These data support the previous conclusion that the microsporidian kinome is the smallest eukaryotic kinome. Microsporidian kinomes contain only serine-threonine kinases and do not contain receptor-like and tyrosine kinases. Many of the kinases related to nutrient and energy signaling and the stress response have been lost in microsporidian kinomes. However, cell cycle-, development- and growth-related kinases, which are important to parasites, are well conserved. This reduction of the microsporidian kinome is in good agreement with genome compaction, but kinome density is negatively correlated with proteome size. Furthermore, the protein kinases in each microsporidian genome are under strong purifying selection pressure. No remarkable differences in kinase family classification, domain features, gain and/or loss, and selective pressure were observed in these four species. Although microsporidia adapt to different host types, the coevolution of microsporidia and their hosts was not clearly reflected in the protein kinases. Overall, this study enriches and updates the microsporidian protein kinase database and may provide valuable information and candidate targets for the design of treatments for pathogenic diseases.

Steinernema carpocapsae is an entomopathogenic nematode widely used for the control of insect pests due to its virulence, which is mainly attributed to the ability the parasitic stage has to overcome insect defences. To identify the mechanisms underlying such a characteristic, we studied a novel serpin-like inhibitor (sc-srp-6) that was detected in a transcriptome analysis. Recombinant Sc-SRP-6 produced in Escherichia coli had a native fold of serpins belonging to the α-1-peptidase family and exhibited inhibitory activity against trypsin and α-chymotrypsin with Ki of 0.42 × 10(-7) M and 1.22 × 10(-7) M, respectively. Functional analysis revealed that Sc-SRP-6 inhibits insect digestive enzymes, thus preventing the hydrolysis of ingested particles. Moreover, Sc-SRP-6 impaired the formation of hard clots at the injury site, a major insect defence mechanism against invasive pathogens. Sc-SRP-6 does not prevent the formation of clot fibres and the activation of prophenoloxidases but impairs the incorporation of the melanin into the clot. Binding assays showed a complex formation between Sc-SRP-6 and three proteins in the hemolymph of lepidopteran required for clotting, apolipophorin, hexamerin and trypsin-like, although the catalytic inhibition occurred exclusively in trypsin-like. This data allowed the conclusion that Sc-SRP-6 promotes nematode virulence by inhibiting insect gut juices and by impairing immune clot reaction.

Melanoma is the most malignant form of skin cancer, which seriously threatens human life and health. Anti-PD-1 immunotherapy has shown clinical benefits in improving patients’ overall survival, but some melanoma patients failed to respond. Effective therapeutic biomarkers are vital to evaluate and optimize benefits from anti-PD-1 treatment. Although the establishment of immunotherapy biomarkers is well underway, studies that identify predictors by gene network-based approaches are lacking. Here, we retrieved the existing datasets (GSE91061, GSE78220 and GSE93157, 79 samples in total) on anti-PD-1 therapy to explore potential therapeutic biomarkers in melanoma using weighted correlation network analysis (WGCNA), function validation and clinical corroboration. As a result, 13 hub genes as critical nodes were traced from the key module associated with clinical features. After receiver operating characteristic (ROC) curve validation by an independent dataset (GSE78220), six hub genes with diagnostic significance were further recovered. Moreover, these six genes were revealed to be closely associated not only with the immune system regulation, immune infiltration, and validated immunotherapy biomarkers, but also with excellent prognostic value and significant expression level in melanoma. The random forest prediction model constructed using these six genes presented a great diagnostic ability for anti-PD-1 immunotherapy response. Taken together, IRF1, JAK2, CD8A, IRF8, STAT5B, and SELL may serve as predictive therapeutic biomarkers for melanoma and could facilitate future anti-PD-1 therapy. [ABSTRACT FROM AUTHOR]

Transforming growth factor β1 (TGF-β1) is a cytokine that plays an important role in the control of cell proliferation and differentiation in breast cancer. The -509C/T polymorphism in the TGF-β1 gene has been implicated in breast cancer risk. However, studies on the association between this polymorphism and breast cancer risk have produced conflicting results. To derive a more precise estimation of the relationship, a meta-analysis of the -509C/T polymorphism (5,825 cases and 7,953 controls) from seven published case-control studies was performed. Our analysis suggests that -509C/T has no association with breast cancer risk when using either dominant [odds ratio (OR) = 1.01, 95% confidence intervals (CI): 0.82-1.24], or recessive models (OR = 0.91, 95% CI: 0.66-1.27), or other genetic models to analyze the data. In ethnic subgroups analysis, -509C/T also did not appear to be a risk factor for breast cancer. However, larger scale primary studies are still required to further evaluate the interaction of TGF-β1 -509C/T polymorphism and breast cancer risk in specific populations. [ABSTRACT FROM AUTHOR]