Genome-wide association analysis and QTL mapping reveal the genetic control of cadmium accumulation in maize leaf

Background

Accumulation of cadmium (Cd) in maize (Zea mays L.) poses a significant risk to human health as it is ingested via the food chain. A genome-wide association study (GWAS) was conducted in a population of 269 maize accessions with 43,737 single nucleotide polymorphisms (SNPs) to identify candidate genes and favorable alleles for controlling Cd accumulation in maize.

Results

When grown in contaminated soil, accessions varied significantly in leaf Cd concentration at both the seeding and maturing stages with phenotypic variation and the coefficient of variation all above 48%. The co-localized region between SYN27837 (147,034,650 bp) and SYN36598 (168,551,327 bp) on chromosome 2 was associated with leaf Cd under three soil conditions varying in Cd content in 2015 and 2016. The significant SNP (SYN25051) at position 161,275,547 could explained 27.1% of the phenotype variation. Through QTL mapping using the IBMSyn10 double haploid (DH) population, we validated the existence of a major QTL identified by GWAS; qLCd2 could explain the 39.8% average phenotype variation across the experiments. Expression of GRMZM2G175576 encoding a cadmium/zinc-transporting ATPase underlying the QTL was significantly increased in roots, stems and leaves of B73, a low Cd accumulation line in response to Cd stress.

Conclusions

Our findings provide new insights into the genetic control of Cd accumulation and could aid rapid development of maize genotypes with low-Cd accumulation by manipulation of the favorable alleles.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4395-x) contains supplementary material, which is available to authorized users.

Genetic Dissection of Maize Embryonic Callus Regenerative Capacity Using Multi-Locus Genome-Wide Association Studies

The regenerative capacity of the embryonic callus, a complex quantitative trait, is one of the main limiting factors for maize transformation. This trait was decomposed into five traits, namely, green callus rate (GCR), callus differentiating rate (CDR), callus plantlet number (CPN), callus rooting rate (CRR), and callus browning rate (CBR). To dissect the genetic foundation of maize transformation, in this study multi-locus genome-wide association studies (GWAS) for the five traits were performed in a population of 144 inbred lines genotyped with 43,427 SNPs. Using the phenotypic values in three environments and best linear unbiased prediction (BLUP) values, as a result, a total of 127, 56, 160, and 130 significant quantitative trait nucleotides (QTNs) were identified by mrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB, respectively. Of these QTNs, 63 QTNs were commonly detected, including 15 across multiple environments and 58 across multiple methods. Allele distribution analysis showed that the proportion of superior alleles for 36 QTNs was <50% in 31 elite inbred lines. Meanwhile, these superior alleles had obviously additive effect on the regenerative capacity. This indicates that the regenerative capacity-related traits can be improved by proper integration of the superior alleles using marker-assisted selection. Moreover, a total of 40 candidate genes were found based on these common QTNs. Some annotated genes were previously reported to relate with auxin transport, cell fate, seed germination, or embryo development, especially, GRMZM2G108933 (WOX2) was found to promote maize transgenic embryonic callus regeneration. These identified candidate genes will contribute to a further understanding of the genetic foundation of maize embryonic callus regeneration.

Melatonin at pharmacological concentrations suppresses osteoclastogenesis via the attenuation of intracellular ROS

Osteoporosis is linked to age-related decline of melatonin production; however, the direct effects of melatonin on osteoclastogenesis remain unknown. Our study demonstrates that melatonin at pharmacological concentrations, rather than at physiological concentrations, significantly inhibits osteoclastogenesis. Melatonin-mediated anti-osteoclastogenesis involves a reactive oxygen species (ROS)-mediated but not a silent information regulator type 1 (SIRT1)-independent pathway.

INTRODUCTION

Osteoporosis is a bone disorder linked to impaired bone formation and excessive bone resorption. Melatonin has been suggested to treat osteoporosis due to its beneficial actions on osteoblast differentiation. However, the direct effects of melatonin on osteoclastogenesis in bone marrow monocytes (BMMs) remain unknown. This study was to investigate whether melatonin at either physiological or pharmacological concentrations could affect osteoclast differentiation.

METHODS

Primary BMMs were isolated from the femurs and tibias of C57BL/6 mice and were induced toward multinucleated osteoclasts, in the presence of melatonin at either physiological (0.01 to 10 nM) or pharmacological (1 to 100 μM) concentrations. Tartrate-resistant acid phosphatase (TRAP) staining was used to label multinucleated osteoclasts and the levels of osteoclast-specific genes were evaluated. To further explore the underlying mechanisms, the roles of silent information regulator type 1 (SIRT1) and reactive oxygen species (ROS) were evaluated.

RESULTS

We found that melatonin at pharmacological concentrations, rather than at physiological concentrations, significantly inhibited osteoclast formation in a dose-dependent manner. The number of TRAP-positive cells and the gene expression of osteoclast-specific markers were significantly downregulated in melatonin-treated BMMs. The melatonin-mediated repression of osteoclast differentiation involved the inhibition of the nuclear factor κ-light-chain-enhancer of activated B cells (NF-κB) signaling pathway. The treatment with SIRT1 inhibitors did not affect osteoclast differentiation but, when supplemented with exogenous hydrogen peroxide, a partial rescue of melatonin-suppressed osteoclastogenesis was observed.

CONCLUSION

Melatonin at pharmacological doses directly inhibited osteoclastogenesis of BMMs by a ROS-mediated but not a SIRT1-independent pathway.

[Analysis of gene mutation of early onset epileptic spasm with unknown reason]

Objective: To summarize the gene mutation of early onset epileptic spasm with unknown reason. Method: In this prospective study, data of patients with early onset epileptic spasm with unknown reason were collected from neurological department of Children's Hospital of Fudan University between March 2016 and December 2016. Patients with known disorders such as infection, metabolic, structural, immunological problems and known genetic mutations were excluded. Patients with genetic disease that can be diagnosed by clinical manifestations and phenotypic characteristics were also excluded. Genetic research methods included nervous system panel containing 1 427 epilepsy genes, whole exome sequencing (WES), analysis of copy number variation (CNV) and karyotype analysis of chromosome. The basic information, phenotypes, genetic results and the antiepileptic treatment of patients were analyzed. Result: Nine of the 17 cases with early onset epileptic spasm were boys and eight were girls. Patients' age at first seizure onset ranged from 1 day after birth to 8 months (median age of 3 months). The first hospital visit age ranged from 1 month to 2 years (median age of 4.5 months). The time of following-up ranged from 8 months to 3 years and 10 months. All the 17 patients had early onset epileptic spasm. Video electroencephalogram was used to monitor the spasm seizure. Five patients had Ohtahara syndrome, 10 had West syndrome, two had unclear classification. In 17 cases, 10 of them had detected pathogenic genes. Nine cases had point mutations, involving SCN2A, ARX, UNC80, KCNQ2, and GABRB3. Except one case of mutations in GABRB3 gene have been reported, all the other cases had new mutations. One patient had deletion mutation in CDKL5 gene. One CNV case had 6q 22.31 5.5MB repeats. Ten cases out of 17 were using 2-3 antiepileptic drugs (AEDs) and the drugs had no effect. Seven cases used adrenocorticotropic hormone (ACTH) and prednisone besides AEDs (a total course for 8 weeks). Among them, five cases had no effect and two cases were seizure free recently. A case with GABRB3 (C.905A>G) had seizure controlled for 3 mouths. A case with ARX (C.700G>A) had seizure controlled for 6 mouths. Conclusion: The early onset epileptic spasm with unknown reason is highly related to genetic disorders. A variety of genetic mutations, especially new mutations were found. Genetic heterogeneity of epileptic spasm is obvious.

Endogenous small interfering RNAs associated with maize embryonic callus formation

The induction efficiency of maize embryonic callus is highly dependent on the genotype, and only a few lines possess a high capacity for callus formation. Although certain genes and pathways have been reported to contribute to the regulation of callus induction, to the best of our knowledge, the functions of the small interfering RNAs (siRNAs) involved in this process remain unknown. In this study, we identified 861 differentially expressed siRNAs and 576 target genes in the callus induction process. These target genes were classified into 3 clusters, and their functions involve controlling metalloexopeptidase activity, catalase activity, transcription regulation, and O-methyltransferase activity. In addition, certain genes related to auxin transport and stem cell or meristem development (e.g., PLT5-like, ARF15, SAUR-like, FAS1-like, Fea3, SCL5, and Zmwox2A) were regulated by the differentially expressed siRNAs. Moreover, zma-siR004119-2 directly cleaves the 5' UTR of Homeobox-transcription factor 25, which further leads to the down-regulation of its expression. Twelve 24-nt-siRNAs led to the hyper-methylation of GRMZM2G013465, which further decreases its expression. These results suggest that differentially expressed siRNAs regulate callus formation by controlling the expression of their target genes.

Integrative analysis of DNA methylation, mRNAs, and small RNAs during maize embryo dedifferentiation

BACKGROUND

Maize (Zea mays) is an important model crop for transgenic studies. However, genetic transformation of maize requires embryonic calli derived from immature embryo, and the impact of utilizing tissue culture methods on the maize epigenome is poorly understood. Here, we generated whole-genome MeDIP-seq data examining DNA methylation in dedifferentiated and normal immature maize embryos.

RESULTS

We observed that most of the dedifferentiated embryos exhibited a methylation increase compared to normal embryos. Increased methylation at promoters was associated with down-regulated protein-coding gene expression; however, the correlation was not strong. Analysis of the callus and immature embryos indicated that the methylation increase was induced during induction of embryonic callus, suggesting phenotypic consequences may be caused by perturbations in genomic DNA methylation levels. The correlation between the 21-24nt small RNAs and DNA methylation regions were investigated but only a statistically significant correlation for 24nt small RNAs was observed.

CONCLUSIONS

These data extend the significance of epigenetic changes during maize embryo callus formation, and the methylation changes might explain some of the observed embryonic callus variation in callus formation.

Metabolomic and Proteomic Analysis of Maize Embryonic Callus induced from immature embryo

The low ratio of embryonic callus (EC) induction has inhibited the rapid development of maize genetic engineering. Still, little is known to explain the genotype-dependence of EC induction. Here, we performed a large-scale, quantitative analysis of the maize EC metabolome and proteome at three typical induction stages in two inbred lines with a range of EC induction capabilities. Comparison of the metabolomes and proteomes suggests that the differential molecular responses begin at an early stage of development and continue throughout the process of EC formation. The two inbred lines show different responses under various conditions, such as metal ion binding, cell enlargement, stem cell formation, meristematic activity maintenance, somatic embryogenesis, cell wall synthesis, and hormone signal transduction. Furthermore, the differences in hormone (auxin, cytokinin, gibberellin, salicylic acid, jasmonic acid, brassinosteroid and ethylene) synthesis and transduction ability could partially explain the higher EC induction ratio in the inbred line 18-599R. During EC formation, repression of the "histone deacetylase 2 and ERF transcription factors" complex in 18-599R activated the expression of downstream genes, which further promoted EC induction. Together, our data provide new insights into the molecular regulatory mechanism responsible for efficient EC induction in maize.

Genome-wide analysis of transcription factors involved in maize embryonic callus formation

In this study, a maize inbred line with a strong capacity to induce embryonic callus, 18-599R, was used to analyze the transcription factors expressed during embryonic callus formation. A total of 1180 transcription factors were found to be expressed during three key stages of callus induction. Of these, compared with control, 361, 346 and 328 transcription factors were significantly downregulated during stages I, II and III, respectively. In contrast, 355, 372 and 401 transcription factors (TFs) were upregulated during the respective stages. We constructed a transcription factor-mediated regulatory network and found that plant hormone signal transduction was the pathway most significantly enriched among TFs. This pathway includes 48 TFs regulating cell enlargement, cell differentiation, cell division and cell dedifferentiation via the response to plant hormones. Through real-time polymerase chain reaction (PCR) and degradome sequencing, we identified 23 transcription factors that are regulated by miRNA. Through further analysis, ZmMYB138, a member of the MYB transcription factor family localized in the nucleus, was verified to promote embryonic callus formation in the maize embryo through GA signal transduction.

Characterizing the population structure and genetic diversity of maize breeding germplasm in Southwest China using genome-wide SNP markers

Background

Maize breeding germplasm used in Southwest China has high complexity because of the diverse ecological features of this area. In this study, the population structure, genetic diversity, and linkage disequilibrium decay distance of 362 important inbred lines collected from the breeding program of Southwest China were characterized using the MaizeSNP50 BeadChip with 56,110 single nucleotide polymorphisms (SNPs).

Results

With respect to population structure, two (Tropical and Temperate), three (Tropical, Stiff Stalk and non-Stiff Stalk), four [Tropical, group A germplasm derived from modern U.S. hybrids (PA), group B germplasm derived from modern U.S. hybrids (PB) and Reid] and six (Tropical, PB, Reid, Iowa Stiff Stalk Synthetic, PA and North) subgroups were identified. With increasing K value, the Temperate group showed pronounced hierarchical structure with division into further subgroups. The Genetic Diversity of each group was also estimated, and the Tropical group was more diverse than the Temperate group. Seven low-genetic-diversity and one high-genetic-diversity regions were collectively identified in the Temperate, Tropical groups, and the entire panel. SNPs with significant variation in allele frequency between the Tropical and Temperate groups were also evaluated. Among them, a region located at 130 Mb on Chromosome 2 showed the highest genetic diversity, including both number of SNPs with significant variation and the ratio of significant SNPs to total SNPs. Linkage disequilibrium decay distance in the Temperate group was greater (2.5–3 Mb) than that in the entire panel (0.5–0.75 Mb) and the Tropical group (0.25–0.5 Mb). A large region at 30–120 Mb of Chromosome 7 was concluded to be a region conserved during the breeding process by comparison between S37, which was considered a representative tropical line in Southwest China, and its 30 most similar derived lines.

Conclusions

For the panel covered most of widely used inbred lines in Southwest China, this work representatively not only illustrates the foundation and evolution trend of maize breeding resource as a theoretical reference for the improvement of heterosis, but also provides plenty of information for genetic researches such as genome-wide association study and marker-assisted selection in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3041-3) contains supplementary material, which is available to authorized users.

Haploid Strategies for Functional Validation of Plant Genes

Increasing knowledge of plant genome sequences requires the development of more reliable and efficient genetic approaches for genotype-phenotype validation. Functional identification of plant genes is generally achieved by a combination of creating genetic modifications and observing the according phenotype, which begins with forward-genetic methods represented by random physical and chemical mutagenesis and move towards reverse-genetic tools as targeted genome editing. A major bottleneck is time need to produce modified homozygous genotypes that can actually be used for phenotypic validation. Herein, we comprehensively address and compare available experimental approaches for functional validation of plant genes, and propose haploid strategies to reduce the time needed and cost consumed for establishing gene function.

Comparative RNA-Seq Analysis Reveals That Regulatory Network of Maize Root Development Controls the Expression of Genes in Response to N Stress

Nitrogen (N) is an essential nutrient for plants, and it directly affects grain yield and protein content in cereal crops. Plant root systems are not only critical for anchorage in the soil, but also for N acquisition. Therefore, genes controlling root development might also affect N uptake by plants. In this study, the responses of nitrogen on root architecture of mutant rtcs and wild-type of maize were investigated by morphological and physiological analysis. Subsequently, we performed a comparative RNA-Seq analysis to compare gene expression profiles between mutant rtcs roots and wild-type roots under different N conditions. We identified 786 co-modulated differentially expressed genes (DEGs) related to root development. These genes participated in various metabolic processes. A co-expression cluster analysis and a cis-regulatory motifs analysis revealed the importance of the AP2-EREBP transcription factor family in the rtcs-dependent regulatory network. Some genotype-specific DEGs contained at least one LBD motif in their promoter region. Further analyses of the differences in gene transcript levels between rtcs and wild-type under different N conditions revealed 403 co-modulated DEGs with distinct functions. A comparative analysis revealed that the regulatory network controlling root development also controlled gene expression in response to N-deficiency. Several AP2-EREBP family members involved in multiple hormone signaling pathways were among the DEGs. These transcription factors might play important roles in the rtcs-dependent regulatory network related to root development and the N-deficiency response. Genes encoding the nitrate transporters NRT2-1, NAR2.1, NAR2.2, and NAR2.3 showed much higher transcript levels in rtcs than in wild-type under normal-N conditions. This result indicated that the LBD gene family mainly functions as transcriptional repressors, as noted in other studies. In summary, using a comparative RNA-Seq-based approach, we identified DEGs related to root development that also participated in the N-deficiency response in maize. These findings will increase our understanding of the molecular regulatory networks controlling root development and N-stress responses.

PRAME promotes in vitro leukemia cells death by regulating S100A4/p53 signaling

OBJECTIVE

PRAME (Preferentially Expressed Antigen in Melanoma) is a tumor-associated antigen recognized by immunocytes, and it induces cytotoxic T cell-mediated responses in melanoma. PRAME is expressed in a wide variety of tumors, but in contrast with most other tumor-associated antigens, it is also expressed in leukemias. The physiologic role of PRAME remains elusive. Recently, it has found PRAME could be involved in the regulation of cell death in leukemias, but the mechanism of the function is unclear. Here, we confirm that PRAME induces leukemias cell death by regulation of S100A4/p53 signaling.

MATERIALS AND METHODS

The pCDNA3-PRAME plasmid and its control were transfected with the KG-1 cells. The pCDNA3-PRAME transfected KG-1 cells were then transiently transfected with S100A4 cDNA or wt-p53 siRNA. The PRAME siRNA and its control were transfected with the K562 cells. The PRAME siRNA transfected K562 cells were then transiently transfected with S100A4 siRNA or pGMp53-Lu. PRAME, S100A4 and P53 were detected by Western blot assay in different time point. Annexin V/propidium iodide and MTT methods were used to detect apoptosis and cell survival rate.

RESULTS

KG-1 cells overexpressing the PRAME gene significantly induces apoptosis and decreases proliferation in vitro, followed by down-regulation of S100A4 and up-regulation of p53. Up-regulation of S100A4 by S100A4 transfection inhibits PRAME-induced p53 up-regulation. Furthermore, up-regulation of S100A4 by S100A4 transfection or down-regulation of p53 by p53 siRNA transfection reduces apoptosis and increases proliferation in vitro. Knockdown of PRAME in K562 cells significantly increases proliferation in vitro, followed by up-regulation of S100A4 and down-regulation of p53. The downregulation of S100A4 by S100A4 siRNA transfection increased p53 expression. Furthermore, downregulation of S100A4 by S100A4 siRNA transfection or up-regulation of p53 by p53 transfection decreases proliferation in vitro.

CONCLUSIONS

Our results suggest that the leukemias expressing high levels of PRAME has a favorable prognosis. PRAME promotes in vitro leukemia cells death by regulating S100A4/p53 signaling.

An ultra-high-density map as a community resource for discerning the genetic basis of quantitative traits in maize

BACKGROUND

To safeguard the food supply for the growing human population, it is important to understand and exploit the genetic basis of quantitative traits. Next-generation sequencing technology performs advantageously and effectively in genetic mapping and genome analysis of diverse genetic resources. Hence, we combined re-sequencing technology and a bin map strategy to construct an ultra-high-density bin map with thousands of bin markers to precisely map a quantitative trait locus.

RESULTS

In this study, we generated a linkage map containing 1,151,856 high quality SNPs between Mo17 and B73, which were verified in the maize intermated B73 × Mo17 (IBM) Syn10 population. This resource is an excellent complement to existing maize genetic maps available in an online database  (iPlant, http://data.maizecode.org/maize/qtl/syn10/ ). Moreover, in this population combined with the IBM Syn4 RIL population, we detected 135 QTLs for flowering time and plant height traits across the two populations. Eighteen known functional genes and twenty-five candidate genes for flowering time and plant height trait were fine-mapped into a 2.21-4.96 Mb interval. Map expansion and segregation distortion were also analyzed, and evidence for inadvertent selection of early flowering time in the process of mapping population development was observed. Furthermore, an updated integrated map with 1,151,856 high-quality SNPs, 2,916 traditional markers and 6,618 bin markers was constructed. The data were deposited into the iPlant Discovery Environment (DE), which provides a fundamental resource of genetic data for the maize genetic research community.

CONCLUSIONS

Our findings provide basic essential genetic data for the maize genetic research community. An updated IBM Syn10 population and a reliable, verified high-quality SNP set between Mo17 and B73 will aid in future molecular breeding efforts.

Characterization of Sugarcane Mosaic Virus Scmv1 and Scmv2 Resistance Regions by Regional Association Analysis in Maize

Sugarcane Mosaic Virus (SCMV) causes one of the most severe virus diseases in maize worldwide, resulting in reduced grain and forage yield in susceptible cultivars. In this study, two association panels consisting of 94 inbred lines each, from China and the U.S., were characterized for resistance to two isolates: SCMV-Seehausen and SCMV-BJ. The population structure of both association panels was analyzed using 3072 single nucleotide polymorphism (SNP) markers. The Chinese and the U.S. panel were both subdivided into two sub-populations, the latter comprised of Stiff Stalk Synthetic (SS) lines and Non Stiff Stalk Synthetic (NSS). The relative kinships were calculated using informative 2947 SNPs with minor allele frequency ≥ 5% and missing data ≤ 20% for the Chinese panel and 2841 SNPs with the same characteristics were used for the U.S. panel. The Scmv1 region was genotyped using 7 single sequence repeat (SSR) and sequence-tagged site (STS) markers, and 12 SSR markers were used for the Scmv2 region in the U.S. panel, while 5 of them were used for the Chinese panel. For all traits, a MLM (Mix Linear Model) controlling both population structure and relative kinship (Q + K) was used for association analysis. Three markers Trx-1, STS-11, and STS-12 located in the Scmv1 region were strongly associated (P = 0.001) with SCMV resistance, and explained more than 16.0%, 10.6%, and 19.7% of phenotypic variation, respectively. 207FG003 located in the Scmv2 region was significantly associated (P = 0.001) with SCMV resistance, and explained around 18.5% of phenotypic variation.

Reduced Penetrance of PRRT2 Mutation in a Chinese Family With Infantile Convulsion and Choreoathetosis Syndrome

Paroxysmal kinesigenic dyskinesia is a rare episodic movement disorder that can be isolated or associated with benign infantile seizures as part of choreoathetosis syndrome. Mutations in the PRRT2 gene have been recently identified as a cause of paroxysmal kinesigenic dyskinesia and infantile convulsion and choreoathetosis (ICCA). We reported a PRRT2 heterozygous mutation (c.604-607delTCAC, p.S202Hfs*25) in a 3-generation Chinese family with infantile convulsion and choreoathetosis and paroxysmal kinesigenic dyskinesia. The mutation was present in 5 family members, of which 4 were clinically affected and 1 was an obligate carrier with reduced penetrance of PRRT2. The affected carriers of this mutation presented with a similar type of infantile convulsion during early childhood and developed additional paroxysmal kinesigenic dyskinesia symptoms later in life. In addition, they all had a dramatic clinical response to oxcarbazepine/phenytoin therapy. Reduced penetrance of the PRRT2 mutation in this family could warrant genetic counseling.

Transcriptomic changes during maize roots development responsive to Cadmium (Cd) pollution using comparative RNAseq-based approach

The heavy metal cadmium (Cd), acts as a widespread environmental contaminant, which has shown to adversely affect human health, food safety and ecosystem safety in recent years. However, research on how plant respond to various kinds of heavy metal stress is scarcely reported, especially for understanding of complex molecular regulatory mechanisms and elucidating the gene networks of plant respond to Cd stress. Here, transcriptomic changes during Mo17 and B73 seedlings development responsive to Cd pollution were investigated and comparative RNAseq-based approach in both genotypes were performed. 115 differential expression genes (DEGs) with significant alteration in expression were found co-modulated in both genotypes during the maize seedling development; of those, most of DGEs were found comprised of stress and defense responses proteins, transporters, as well as transcription factors, such as thaumatin-like protein, ZmOPR2 and ZmOPR5. More interestingly, genotype-specific transcriptional factors changes induced by Cd stress were found contributed to the regulatory mechanism of Cd sensitivity in both different genotypes. Moreover, 12 co-expression modules associated with specific biological processes or pathways (M1 to M12) were identified by consensus co-expression network. These results will expand our understanding of complex molecular mechanism of response and defense to Cd exposure in maize seedling roots.

The Impact of Genetic Relationship and Linkage Disequilibrium on Genomic Selection

Genomic selection is a promising research area due to its practical application in breeding. In this study, impact of realized genetic relationship and linkage disequilibrium (LD) on marker density and training population size required was investigated and their impact on practical application was further discussed. This study is based on experimental data of two populations derived from the same two founder lines (B73, Mo17). Two populations were genotyped with different marker sets at different density: IBM Syn4 and IBM Syn10. A high-density marker set in Syn10 was imputed into the Syn4 population with low marker density. Seven different prediction scenarios were carried out with a random regression best linear unbiased prediction (RR-BLUP) model. The result showed that the closer the real genetic relationship between training and validation population, the fewer markers were required to reach a good prediction accuracy. Taken the short-term cost for consideration, relationship information is more valuable than LD information. Meanwhile, the result indicated that accuracies based on high LD between QTL and markers were more stable over generations, thus LD information would provide more robust prediction capacity in practical applications.

Identification and characterisation of tobacco microRNA transcriptome using high-throughput sequencing

MicroRNAs (miRNAs) are post-transcriptional regulators that are involved in numerous biological processes in plants. In this study, we investigate miRNAs in Honghua Dajinyuan, an agronomically important species of tobacco in China. Here, we report a comprehensive analysis of miRNA expression profiles in the leaf, stem and root using a high-throughput sequencing approach. A total of 165 miRNAs, representing 55 conserved families, and 50 novel miRNAs, representing 19 families, were identified in three libraries. In addition, 12 miRNAs were randomly selected from a differentially expressed conserved miRNA family in three libraries with expression alterations and subjected to qRT-PCR validation. Of these, the expression level of nta-miR167d is highly enriched in the leaf tissue. In addition, the expression level of nta-miR319a is prominently enriched in the stem, while nta-miR160c is highly enriched in the root. Moreover, the target prediction showed that most of the targets coded for transcription factors that are involved in cellular and metabolic processes. GO analysis showed that most of the targets were involved in organelle function, served binding functions, and take part in cellular and metabolic processes. This study helps shed new light on understanding the role of miRNAs in different parts of the tobacco plant and adds a significant number of novel miRNAs to the tobacco miRNA transcriptome.

Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development

The present study profiled and analyzed gene expression of the maize ear at four key developmental stages. Based on genome-wide profile analysis, we detected differential mRNA of maize genes. Some of the differentially expressed genes (DEGs) were predicted to be potential candidates of maize ear development. Several well-known genes were found with reported mutant analyses, such as, compact plant2 (ct2), zea AGAMOUS homolog1 (zag1), bearded ear (bde), and silky1 (si1). MicroRNAs such as microRNA156 were predicted to target genes involved in maize ear development. Antisense transcripts were widespread throughout all the four stages, and are suspected to play important roles in maize ear development. Thus, identification and characterization of important genes and regulators at all the four developmental stages will contribute to an improved understanding of the molecular mechanisms responsible for maize ear development.

Genome-wide identification and analysis of drought-responsive genes and microRNAs in tobacco

Drought stress response is a complex trait regulated at transcriptional and post-transcriptional levels in tobacco. Since the 1990s, many studies have shown that miRNAs act in many ways to regulate target expression in plant growth, development and stress response. The recent draft genome sequence of Nicotiana benthamiana has provided a framework for Digital Gene Expression (DGE) and small RNA sequencing to understand patterns of transcription in the context of plant response to environmental stress. We sequenced and analyzed three Digital Gene Expression (DGE) libraries from roots of normal and drought-stressed tobacco plants, and four small RNA populations from roots, stems and leaves of control or drought-treated tobacco plants, respectively. We identified 276 candidate drought responsive genes (DRGs) with sequence similarities to 64 known DRGs from other model plant crops, 82 were transcription factors (TFs) including WRKY, NAC, ERF and bZIP families. Of these tobacco DRGs, 54 differentially expressed DRGs included 21 TFs, which belonged to 4 TF families such as NAC (6), MYB (4), ERF (10), and bZIP (1). Additionally, we confirmed expression of 39 known miRNA families (122 members) and five conserved miRNA families, which showed differential regulation under drought stress. Targets of miRNAs were further surveyed based on a recently published study, of which ten targets were DRGs. An integrated gene regulatory network is proposed for the molecular mechanisms of tobacco root response to drought stress using differentially expressed DRGs, the changed expression profiles of miRNAs and their target transcripts. This network analysis serves as a reference for future studies on tobacco response stresses such as drought, cold and heavy metals.

Prevalence and awareness of diabetes mellitus among a rural population in China: results from Liaoning Province

AIM

To clarify the diabetes prevalence trends among the rural population in northern China.

METHODS

All eligible permanent residents aged ≥ 35 years in selected rural villages of Liaoning province were invited to participate in the study. A total of 11 600 people completed all questionnaires and were included in the study. The response rate was 85.3%. Fasting plasma glucose levels were measured after at least 12 h of fasting and diabetes was diagnosed according to WHO criteria, i.e. fasting plasma glucose ≥ 7 mmol/l and/or being on treatment for diabetes. Impaired fasting glucose was defined according to the 1997 and the 2010 American Diabetes Association (ADA) criteria (6.1-6.9 and 5.6-6.9 mmol/l, respectively). Previous diagnoses of diabetes were assessed on the basis of self-reports.

RESULTS

The prevalence of diabetes among adults in the rural population was 10.6% (10.0% in men and 11.1% in women). The prevalence of impaired fasting glucose was 13.0 and 36.1% according to the 1997 and the 2010 ADA criteria, respectively. The prevalence of previously diagnosed diabetes was 4.3% among the whole population (3.3% in men and 5.1% in women). The prevalence of previously diagnosed diabetes was 34.8% in men and 50.2% in women. Only 29.6% of men and 42% of women with diabetes had taken oral hypoglycaemic agents or insulin to lower their blood sugar. In multivariate analysis, age, drinking habits, BMI, dyslipidaemia and family history of diabetes were identified as independent risk factors for diabetes, and occupational physical activity, smoking and lean meat intake were identified as independent protective factors for diabetes.

CONCLUSIONS

The prevalence of diabetes and impaired fasting glucose were found to be high in this rural population. Although the rate of treatment of people with diabetes has increased, the glycaemic control rate was still at a low level.

The development dynamics of the maize root transcriptome responsive to heavy metal Pb pollution

Lead (Pb), as a heavy metal element, has become the most important metal pollutant of the environment. With allocating a relatively higher proportion of its biomass in roots, maize could be a potential important model to study the phytoremediation of Pb-contaminated soil. Here we analyzed the maize root transcriptome of inbred lines 9782 under heavy metal lead (Pb) pollution, which was identified as a non-hyperaccumulator for Pb in roots. In the present study, more than 98 millions reads were mapped to define gene structure and detect polymorphism, thereby to qualify transcript abundance along roots development under Pb treatment. A total of 17,707, 17,440, 16,998 and 16,586 genes were identified in maize roots at four developmental stages (0, 12 h, 24 h and 48 h) respectively and 2,825, 2,626, 2161 and 2260 stage-specifically expressed genes were also identified respectively. In addition, based on our RNA-Seq data, transcriptomic changes during maize root development responsive to Pb were investigated. A total of 384 differentially expressed genes (DEGs) (log2Ratio ≥ 1, FDR ≤ 0.001) were identified, of which, 36 genes with significant alteration in expression were detected in four developmental stages; 12 DEGs were randomly selected and successful validated by qRT-PCR. Additionally, many transcription factor families might act as the important regulators at different developmental stages, such as bZIP, ERF and GARP et al. These results will expand our understanding of the complex molecular and cellular events in maize root development and provide a foundation for future study on root development in maize under heavy metal pollution and other cereal crops.