Brassinosteroids repress the seed maturation program during the seed-to-seedling transition.

In flowering plants, repression of the seed maturation program is essential for the transition from the seed to the vegetative phase, but the underlying mechanisms remain poorly understood. The B3-domain protein VIVIPAROUS1/ABSCISIC ACID-INSENSITIVE3-LIKE 1 (VAL1) is involved in repressing the seed maturation program. Here we uncovered a molecular network triggered by the plant hormone brassinosteroid (BR) that inhibits the seed maturation program during the seed-to-seedling transition in Arabidopsis (Arabidopsis thaliana). val1-2 mutant seedlings treated with a BR biosynthesis inhibitor form embryonic structures, whereas BR signaling gain-of-function mutations rescue the embryonic structure trait. Furthermore, the BR-activated transcription factors BRI1-EMS-SUPPRESSOR 1 and BRASSINAZOLE-RESISTANT 1 bind directly to the promoter of AGAMOUS-LIKE15 (AGL15), which encodes a transcription factor involved in activating the seed maturation program, and suppress its expression. Genetic analysis indicated that BR signaling is epistatic to AGL15 and represses the seed maturation program by downregulating AGL15. Finally, we showed that the BR-mediated pathway functions synergistically with the VAL1/2-mediated pathway to ensure the full repression of the seed maturation program. Together, our work uncovered a mechanism underlying the suppression of the seed maturation program, shedding light on how BR promotes seedling growth.

AtPER1 enhances primary seed dormancy and reduces seed germination by suppressing the ABA catabolism and GA biosynthesis in Arabidopsis seeds.

Seed is vital to the conservation of germplasm and plant biodiversity. Seed dormancy is an adaptive trait in numerous seed-plant species, enabling plants to survive under stressful conditions. Seed dormancy is mainly controlled by abscisic acid (ABA) and gibberellin (GA) and can be classified as primary and secondary seed dormancy. The primary seed dormancy is induced by maternal ABA. Here we found that AtPER1, a seed-specific peroxiredoxin, is involved in enhancing primary seed dormancy. Two loss-of-function atper1 mutants, atper1-1 and atper1-2, displayed suppressed primary seed dormancy accompanied with reduced ABA and increased GA contents in seeds. Furthermore, atper1 mutant seeds were insensitive to abiotic stresses during seed germination. The expression of several ABA catabolism genes (CYP707A1, CYP707A2, and CYP707A3) and GA biosynthesis genes (GA20ox1, GA20ox3, and KAO3) in atper1 mutant seeds was increased compared to wild-type seeds. The suppressed primary seed dormancy of atper1-1 was completely reduced by deletion of CYP707A genes. Furthermore, loss-of-function of AtPER1 cannot enhance the seed germination ratio of aba2-1 or ga1-t, suggesting that AtPER1-enhanced primary seed dormancy is dependent on ABA and GA. Additionally, the level of reactive oxygen species (ROS) in atper1 mutant seeds was significantly higher than that in wild-type seeds. Taken together, our results demonstrate that AtPER1 eliminates ROS to suppress ABA catabolism and GA biosynthesis, and thus improves the primary seed dormancy and make the seeds less sensitive to adverse environmental conditions.

The H3K27me3 Demethylase RELATIVE OF EARLY FLOWERING6 Suppresses Seed Dormancy by Inducing Abscisic Acid Catabolism.

Seed dormancy is an adaptive trait that is crucial to plant survival. Abscisic acid (ABA) is the primary phytohormone that induces seed dormancy. However, little is known about how the level of ABA in seeds is determined. Here we show that the Arabidopsis (Arabidopsis thaliana) H3K27me3 demethylase RELATIVE OF EARLY FLOWERING6 (REF6) suppresses seed dormancy by inducing ABA catabolism in seeds. Seeds of the ref6 loss-of-function mutants displayed enhanced dormancy that was associated with increased endogenous ABA content. We further show that the transcripts of two genes key to ABA catabolism, CYP707A1 and CYP707A3, but not genes involved in ABA biosynthesis, were significantly reduced in ref6 mutants during seed development and germination. In developing siliques, REF6 bound directly to CYP707A1 and CYP707A3, and was responsible for reducing their H3K27me3 levels. Genetic analysis demonstrated that the enhanced seed dormancy and ABA concentration in ref6 depended mainly on the reduced expression of CYP707A1 and CYP707A3 Conversely, overexpression of CYP707A1 could offset the enhanced seed dormancy of ref6 Taken together, our results revealed an epigenetic regulation mechanism that is involved in the regulation of ABA content in seeds.

Single-molecule optical mapping enables quantitative measurement of D4Z4 repeats in facioscapulohumeral muscular dystrophy (FSHD).

PURPOSE: Facioscapulohumeral muscular dystrophy (FSHD) is a common adult muscular dystrophy. Over 95% of FSHD cases are associated with contraction of the D4Z4 tandem repeat (~3.3 kb per unit) at 4q35 with a specific genomic configuration (haplotype) called 4qA. Molecular diagnosis of FSHD typically requires pulsed-field gel electrophoresis with Southern blotting. We aim to develop novel genomic and computational methods for characterising D4Z4 repeat numbers in FSHD. METHODS: We leveraged a single-molecule optical mapping platform that maps locations of restriction enzyme sites on high molecular weight (>150 kb) DNA molecules. We developed bioinformatics methods to address several challenges, including the differentiation of 4qA with 4qB alleles, the differentiation of 4q35 and 10q26 segmental duplications, the quantification of repeat numbers with different enzymes that may or may not have recognition sites within D4Z4 repeats. We evaluated the method on 25 human subjects (13 patients, 3 individual control subjects, 9 control subjects from 3 families) labelled by the Nb.BssSI and/or Nt.BspQI enzymes. RESULTS: We demonstrated that the method gave a direct quantitative measurement of repeat numbers on D4Z4 repeats with 4qA allelic configuration and the levels of postzygotic mosaicism. Our method had high concordance with Southern blots from several cohorts on two platforms (Bionano Saphyr and Bionano Irys), but with improved quantification of repeat numbers. CONCLUSION: While the study is limited by small sample size, our results demonstrated that single-molecule optical mapping is a viable approach for more refined analysis on genotype-phenotype relationships in FSHD, especially when postzygotic mosaicism is present.

RNA polymerase II-independent recruitment of SPT6L at transcription start sites in Arabidopsis.

SPT6 is a conserved elongation factor that is associated with phosphorylated RNA polymerase II (RNAPII) during transcription. Recent transcriptome analysis in yeast mutants revealed its potential role in the control of transcription initiation at genic promoters. However, the mechanism by which this is achieved and how this is linked to elongation remains to be elucidated. Here, we present the genome-wide occupancy of Arabidopsis SPT6-like (SPT6L) and demonstrate its conserved role in facilitating RNAPII occupancy across transcribed genes. We also further demonstrate that SPT6L enrichment is unexpectedly shifted, from gene body to transcription start site (TSS), when its association with RNAPII is disrupted. Protein domains, required for proper function and enrichment of SPT6L on chromatin, are subsequently identified. Finally, our results suggest that recruitment of SPT6L at TSS is indispensable for its spreading along the gene body during transcription. These findings provide new insights into the mechanisms underlying SPT6L recruitment in transcription and shed light on the coordination between transcription initiation and elongation.

[Clinical practice guidelines for phenylketonuria].

Phenylketonuria (PKU) is an autosomal recessive hereditary disease and a common disorder of amino acid metabolism. The average incidence of PKU in China is approximately 1/11 000. It is characterized by lower incidence in the South and higher incidence in the North, particularly the Northwest. PKU is a treatable disease and has been listed in the national newborn screening program. Neonates with positive indication of screening can achieve satisfactory therapeutic effect by timely control of phenylalanine intake after the definite diagnosis. This guideline aims to summarize the knowledge of medical genetics and key points of clinical management of PKU, so as to improve the diagnostic level and standardize newborn screening and clinical treatment of patients.

[A consensus on the standardization of the next generation sequencing process for the diagnosis of genetic diseases (1) - Procedures prior to genetic testing].

Pre-testing preparation is the basis and starting point of genetic testing. The process includes collection of clinical information, formulation of testing scheme, genetic counseling before testing, and completion of informed consent and testing authorization. To effectively identify genetic diseases in clinics can greatly improve the diagnostic rate of next generation sequencing (NGS), thereby reducing medical cost and improving clinical efficacy. The analysis of NGS results relies, to a large extent, on the understanding of genotype-phenotype correlations, therefore it is particularly important to collect and evaluate clinical phenotypes and describe them in uniform standard terms. Different types of genetic diseases or mutations may require specific testing techniques, which can yield twice the result with half the effort. Pre-testing genetic counseling can help patients and their families to understand the significance of relevant genetic testing, formulate individualized testing strategies, and lay a foundation for follow-up.

[A consensus on the standardization of the next generation sequencing process for the diagnosis of genetic diseases (2) - Sample collection, processing and detection].

With high accuracy and precision, next generation sequencing (NGS) has provided a powerful tool for clinical testing of genetic diseases. To follow a standardized experimental procedure is the prerequisite to obtain stable, reliable, and effective NGS data for the assistance of diagnosis and/or screening of genetic diseases. At a conference of genetic testing industry held in Shanghai, May 2019, physicians engaged in the diagnosis and treatment of genetic diseases, experts engaged in clinical laboratory testing of genetic diseases and experts from third-party genetic testing companies have fully discussed the standardization of NGS procedures for the testing of genetic diseases. Experts from different backgrounds have provided opinions for the operation and implementation of NGS testing procedures including sample collection, reception, preservation, library construction, sequencing and data quality control. Based on the discussion, a consensus on the standardization of the testing procedures in NGS laboratories is developed with the aim to standardize NGS testing and accelerate implementation of NGS in clinical settings across China.

[A consensus on the standardization of the next generation sequencing process for the diagnosis of genetic diseases (3) - Data analysis].

Bioinformatic analysis and variant classification are the key components of high-throughput sequencing-based genetic diagnostic approach. This consensus is part of the effort to develop a standardized process for next generation sequencing (NGS)-based test for germline mutations underlying Mendelian disorders in China. The flow-chart, common software, key parameters of bioinformatics pipeline for data processing, annotation, storage and variant classification are reviewed, which is aimed to help improving and maintaining a high-quality process and obtaining consistent outcomes for NGS-based molecular diagnosis.

[A consensus on the standardization of the next generation sequencing process for the diagnosis of genetic diseases (4) - Report interpretation and genetic counseling].

Clinical genetic testing results are compiled into a standardized report by genetic specialists and provided to clinicians and patients (Should the patient be intellectually disabled or under 18, the report will be provided to his/her parents or legal guardians). The content of genetic testing report should conform to relevant guidelines, industry standards and consensus. The decisions of clinicians will be made based on the report and clinical indications. Genetic counselors should provide post-test counseling to clinicians and patients or their authorized family members. A mechanism of follow-up visit after the genetic testing should be established with informed consent. Data should be shared by clinical institutions and genome sequencing institutions. As findings upon follow-up visit can help with further evaluation of the results, genome sequencing institutions should regularly re-analyze historical and follow-up data, and the updated results should be shared with clinical institutions. All activities involving reporting, genetic counselling, follow-up visiting, and re-analyzing should follow the relevant guidelines and regulations.

Histone deacetylases HDA6 and HDA9 coordinately regulate valve cell elongation through affecting auxin signaling in Arabidopsis.

Both Histone Deacetylases HDA6 and HDA9 belong to class I subfamily of RPD3/HDA1 HDACs. Loss-of-function mutants of HDA9 form slightly blunt siliques. However, the involvement of HDA6 in regulating silique tip growth is unclear. In this study, we show that HDA6 acts redundantly with HDA9 in regulating the elongation of valve cells in the silique tip. Although the hda6 single mutant does not exhibit a detectable silique phenotype, the silique tip of hda6 hda9 double mutant displays a more severe bulge, a morphology we termed as "nock-shaped". The valve cells of the silique tip of hda9 are longer than wild-type, and loss of HDA6 in hda9 enhances the valve cell elongation phenotype. The transcript levels of auxin-signaling-related genes are mis-regulated in hda9 and hda6 hda9 siliques, and the GFP reporter driven by the auxin response promoter DR5 is weaker in hda9 or hda6 hda9 than wild-type or hda6. Thus, our findings reveal that HDA6 and HDA9 coordinately control the elongation of silique valve cells through regulating the expression of auxin-related genes in silique tips.

Genome-wide characterization of a SRO gene family involved in response to biotic and abiotic stresses in banana (Musa spp.).

BACKGROUND: Banana (Musa spp.) is one of the world's most important fruits and its production is largely limited by diverse stress conditions. SROs (SIMILAR TO RCD-ONE) have important functions in abiotic stress resistance and development of plants. They contain a catalytic core of the poly(ADP-ribose) polymerase (PARP) domain and a C-terminal RST (RCD-SRO-TAF4) domain. In addition, partial SROs also include an N-terminal WWE domain. Although a few of SROs have been characterized in some model plants, little is known about their functions in banana, especially in response to biotic stress. RESULTS: Six MaSRO genes in banana genome were identified using the PARP and RST models as a query. Phylogenetic analysis showed that 77 SROs from 15 species were divided into two structurally distinct groups. The SROs in the group I possessed three central regions of the WWE, PARP and RST domains. The WWE domain was lacking in the group II SROs. In the selected monocots, only MaSROs of banana were present in the group II. Most of MaSROs expressed in more than one banana tissue. The stress- and hormone-related cis-regulatory elements (CREs) in the promoter regions of MaSROs supported differential transcripts of MaSROs in banana roots treated with abiotic and biotic stresses. Moreover, expression profiles of MaSROs in the group I were clearly distinct with those observed in the group II after hormone treatment. Notably, the expression of MaSRO4 was significantly upregulated by the multiple stresses and hormones. The MaSRO4 protein could directly interact with MaNAC6 and MaMYB4, and the PARP domain was required for the protein-protein interaction. CONCLUSIONS: Six MaSROs in banana genome were divided into two main groups based on the characteristics of conserved domains. Comprehensive expression analysis indicated that MaSROs had positive responses to biotic and abiotic stresses via a complex interaction network with hormones. MaSRO4 could interact directly with MaNAC6 and MaMYB4 through the PARP domain to regulate downstream signaling pathway.

The LYSIN MOTIF-CONTAINING RECEPTOR-LIKE KINASE 1 protein of banana is required for perception of pathogenic and symbiotic signals.

How plants can distinguish pathogenic and symbiotic fungi remains largely unknown. Here, we characterized the role of MaLYK1, a lysin motif receptor kinase of banana. Live cell imaging techniques were used in localization studies. RNA interference (RNAi)-silenced transgenic banana plants were generated to analyze the biological role of MaLYK1. The MaLYK1 ectodomain, chitin beads, chitooligosaccharides (COs) and mycorrhizal lipochitooligosaccharides (Myc-LCOs) were used in pulldown assays. Ligand-induced MaLYK1 complex formation was tested in immunoprecipitation experiments. Chimeric receptors were expressed in Lotus japonicus to characterize the function of the MaLYK1 kinase domain. MaLYK1 was localized to the plasma membrane. MaLYK1 expression was induced by Foc4 (Fusarium oxysporum f. sp. cubense race 4) and diverse microbe-associated molecular patterns. MaLYK1-silenced banana lines showed reduced chitin-triggered defense responses, increased Foc4-induced disease symptoms and reduced mycorrhization. The MaLYK1 ectodomain was pulled down by chitin beads and LCOs or COs impaired this process. Ligand treatments induced MaLYK1 complex formation in planta. The kinase domain of MaLYK1 could functionally replace that of the chitin elicitor receptor kinase 1 (AtCERK1) in Arabidopsis thaliana and of a rhizobial LCO (Nod factor) receptor (LjNFR1) in L. japonicus. MaLYK1 represents a central molecular switch that controls defense- and symbiosis-related signaling.

The Arabidopsis SWI2/SNF2 chromatin Remodeler BRAHMA regulates polycomb function during vegetative development and directly activates the flowering repressor gene SVP.

The chromatin remodeler BRAHMA (BRM) is a Trithorax Group (TrxG) protein that antagonizes the functions of Polycomb Group (PcG) proteins in fly and mammals. Recent studies also implicate such a role for Arabidopsis (Arabidopsis thaliana) BRM but the molecular mechanisms underlying the antagonism are unclear. To understand the interplay between BRM and PcG during plant development, we performed a genome-wide analysis of trimethylated histone H3 lysine 27 (H3K27me3) in brm mutant seedlings by chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq). Increased H3K27me3 deposition at several hundred genes was observed in brm mutants and this increase was partially supressed by removal of the H3K27 methyltransferase CURLY LEAF (CLF) or SWINGER (SWN). ChIP experiments demonstrated that BRM directly binds to a subset of the genes and prevents the inappropriate association and/or activity of PcG proteins at these loci. Together, these results indicate a crucial role of BRM in restricting the inappropriate activity of PcG during plant development. The key flowering repressor gene SHORT VEGETATIVE PHASE (SVP) is such a BRM target. In brm mutants, elevated PcG occupancy at SVP accompanies a dramatic increase in H3K27me3 levels at this locus and a concomitant reduction of SVP expression. Further, our gain- and loss-of-function genetic evidence establishes that BRM controls flowering time by directly activating SVP expression. This work reveals a genome-wide functional interplay between BRM and PcG and provides new insights into the impacts of these proteins in plant growth and development.

[Discussion on the standard of clinical genetic testing report and the consensus of gene testing industry].

The widespread application of next generation sequencing (NGS) in clinical settings has enabled testing, diagnosis, treatment and prevention of genetic diseases. However, many issues have arisen in the meanwhile. One of the most pressing issues is the lack of standards for reporting genetic test results across different service providers. The First Forum on Standards and Specifications for Clinical Genetic Testing was held to address the issue in Shenzhen, China, on October 28, 2017. Participants, including geneticists, clinicians, and representatives of genetic testing service providers, discussed problems of clinical genetic testing services across in China and shared opinions on principles, challenges, and standards for reporting clinical genetic test results. Here we summarize expert opinions presented at the seminar and report the consensus, which will serve as a basis for the development of standards and guidelines for reporting of clinical genetic testing results, in order to promote the standardization and regulation of genetic testing services in China.

Ectopic expression of NnPER1, a Nelumbo nucifera 1-cysteine peroxiredoxin antioxidant, enhances seed longevity and stress tolerance in Arabidopsis.

Seed longevity, the maintenance of viability during storage, is a major factor for conservation of genetic resources and biodiversity. Seed longevity is an important trait of agriculture crop and is impaired by reactive oxygen species (ROS) during seed desiccation, storage and germination (C. R. Biol., 331, 2008 and 796). Seeds possess a wide range of systems (protection, detoxification, repair) allowing them to survive during storage and to preserve a high germination ability. In many plants, 1-cys peroxiredoxin (1-Cys Prx, also named PER1) is a seed-specific antioxidant which eliminates ROS with cysteine residues. Here we identified and characterized a seed-specific PER1 protein from seeds of sacred lotus (Nelumbo nucifera Gaertn.). Purified NnPER1 protein protects DNA against the cleavage by ROS in the mixed-function oxidation system. The transcription and protein accumulation of NnPER1 increased during seed desiccation and imbibition and under abiotic stress treatment. Ectopic expression of NnPER1 in Arabidopsis enhanced the seed germination ability after controlled deterioration treatment (CDT), indicating that NnPER1 improves seed longevity of transgenic plants. Consistent with the function of NnPER1 on detoxifying ROS, we found that the level of ROS release and lipid peroxidation was strikingly lower in transgenic seeds compared to wild-type with or without CDT. Furthermore, transgenic Arabidopsis seeds ectopic-expressing NnPER1 displayed enhanced tolerance to high temperature and abscisic acid (ABA), indicating that NnPER1 may participate in the thermotolerance and ABA signaling pathway.

Increase of a group of PTC(+) transcripts by curcumin through inhibition of the NMD pathway.

Nonsense-mediated mRNA decay (NMD), an mRNA surveillance mechanism, eliminates premature termination codon-containing (PTC⁺) transcripts. For instance, it maintains the homeostasis of splicing factors and degrades aberrant transcripts of human genetic disease genes. Here we examine the inhibitory effect on the NMD pathway and consequent increase of PTC+ transcripts by the dietary compound curcumin. We have found that several PTC⁺ transcripts including that of serine/arginine-rich splicing factor 1 (SRSF1) were specifically increased in cells by curcumin. We also observed a similar curcumin effect on the PTC⁺ mutant transcript from a Tay-Sachs-causing HEXA allele or from a beta-globin reporter gene. The curcumin effect was accompanied by significantly reduced expression of the NMD factors UPF1, 2, 3A and 3B. Consistently, in chromatin immunoprecipitation assays, curcumin specifically reduced the occupancy of acetyl-histone H3 and RNA polymerase II at the promoter region (-376 to -247nt) of human UPF1, in a time- and dosage-dependent way. Importantly, knocking down UPF1 abolished or substantially reduced the difference of PTC(+) transcript levels between control and curcumin-treated cells. The disrupted curcumin effect was efficiently rescued by expression of exogenous Myc-UPF1 in the knockdown cells. Together, our data demonstrate that a group of PTC⁺ transcripts are stabilized by a dietary compound curcumin through the inhibition of UPF factor expression and the NMD pathway.

HISTONE DEACETYLASE19 interacts with HSL1 and participates in the repression of seed maturation genes in Arabidopsis seedlings.

The seed maturation genes are specifically and highly expressed during late embryogenesis. In this work, yeast two-hybrid, bimolecular fluorescence complementation, and coimmunoprecipitation assays revealed that HISTONE DEACETYLASE19 (HDA19) interacted with the HIGH-LEVEL EXPRESSION OF SUGAR-INDUCIBLE GENE2-LIKE1 (HSL1), and the zinc-finger CW [conserved Cys (C) and Trp (W) residues] domain of HSL1 was responsible for the interaction. Furthermore, we found that mutations in HDA19 resulted in the ectopic expression of seed maturation genes in seedlings, which was associated with increased levels of gene activation marks, such as Histone H3 acetylation (H3ac), Histone H4 acetylation (H4ac), and Histone H3 Lys 4 tri-methylation (H3K4me3), but decreased levels of the gene repression mark Histone H3 Lys 27 tri-methylation (H3K27me3) in the promoter and/or coding regions. In addition, elevated transcription of certain seed maturation genes was also found in the hsl1 mutant seedlings, which was also accompanied by the enrichment of gene activation marks but decreased levels of the gene repression mark. Chromatin immunoprecipitation assays showed that HDA19 could directly bind to the chromatin of the seed maturation genes. These results suggest that HDA19 and HSL1 may act together to repress seed maturation gene expression during germination. Further genetic analyses revealed that the homozygous hsl1 hda19 double mutants are embryonic lethal, suggesting that HDA19 and HSL1 may play a vital role during embryogenesis.

[A novel indel NF1 mutation identified in a patient with neurofibromatosis type 1].

OBJECTIVE: To identify the genetic etiology in a Chinese patient with neurofibromatosis type 1 (NF-1). METHODS: All coding exons and the flanking sequences of neurofibromin 1 (NF1) gene from the patient were captured, individually barcoded and subjected to HiSeq2000 high-throughput sequencing. Suspected mutation was validated in the nuclear family members with Sanger sequencing. RESULTS: A novel indel mutation, c.789_790delAGinsT, was identified in the exon 8 of the NF1 gene in the patient but not in her asymptomatic parents. The mutation was predicted to have caused shifting of the reading frame and a premature downstream stop codon (p.K263Nfs*18). Two known polymorphisms, c.888+108 C>T (rs2953000) and c.888+118 G>T (rs2952999), was detected in the flanking of the indel mutation in the patient and her father. Sequencing chromatogram for the family indicates that above changes are located on the same chromosome. CONCLUSION: The c.789_790delAGinsT, as a de novo mutation occurring on the paternally derived chromosome, is most likely to be causative for the disease. Compared with Sanger sequencing, targeted next-generation sequencing is more efficient and can dramatically reduce the cost for the genetic testing of NF-1.

[Association study between candidate genes on transforming growth factor-ß signaling pathway and the risk of non-syndromic cleft lip with or without cleft palate in Chinese populations].

OBJECTIVE: To explore the association between 10 candidate genes on transforming growth factor-ß (TGFB) signaling pathway and non-syndromic cleft lip with or without cleft palate (NSCL/P) among Chinese populations, and to study the gene-environment interaction. METHODS: A total of 806 Chinese Han NSCL/P trios were ascertained from an international consortium, which conducted a genome-wide association study using a case-parent trio design to investigate the genes affecting risk to NSCL/P. The transmission disequilibrium test (TDT) was used to test for effects of 343 single nucleotide polymorphisms (SNPs) in 10 genes on TGFB signaling pathway including DCN, TGFB1, TGFB2, TGFB3, TGFBR1, TGFBR2, BAMBI, SMAD2, SMAD3 and SMAD4. The conditional regression models were used to test for gene-environment interaction. RESULTS: For TDT, although 19 SNPs showed nominal significant association with NSCL/P, no significant evidence of association was seen for all SNPs in 806 NSCL/P trios after Bonferroni correction. The interactions between genes and maternal smoking, environmental tobacco smoke, alcohol consumption and multi-vitamin supplementation during pregnancy did not attain statistical significance after correction for multiple comparisons. CONCLUSION: No evidence for SNP effect of genes on TGFB signaling pathway and significant gene-environment interaction was seen in our data.