Functional Bph14 from Rathu Heenati promotes resistance to BPH at the early seedling stage of rice (Oryza sativa L.) as revealed by QTL-seq.

KEY MESSAGE: A QTL associated with BPH resistance at the early seedling stage was identified on chromosome 3. Functional Bph14 in Rathu Heenati was associated with BPH resistance at the early seedling stage. Brown planthopper (BPH; Nilaparvata lugens Stål) is considered the most important rice pest in many Asian countries. Several BPH resistance genes have previously been identified. However, there are few reports of genes specific for BPH resistance at the early seedling stage, a crucial stage for direct-seeding cultivation. In this study, we performed a QTL-seq analysis using two bulks (20 F2 lines in each bulk) of the F2 population (n = 300) derived from a cross of Rathu Heenati (RH) × HCS-1 to identify QTL/genes associated with BPH resistance at the early seedling stage. An important QTL was identified on chromosome 3 and Bph14 was identified as a potential candidate gene based on the differences in gene expression and sequence variation when compared with the two parents. All plants in the resistant bulks possessed the functional Bph14 from RH and all plants in the susceptible bulk and HCS-1 contained a large deletion (2703 bp) in Bph14. The functional Bph14 gene of RH appears to be important for BPH resistance at the early seedling stage of rice and could be used in conjunction with other BPH resistance genes in rice breeding programs that confer resistance to BPH at the early and later growth stages.

Reproductive phasiRNA loci and DICER-LIKE5, but not microRNA loci, diversified in monocotyledonous plants.

In monocots other than maize (Zea mays) and rice (Oryza sativa), the repertoire and diversity of microRNAs (miRNAs) and the populations of phased, secondary, small interfering RNAs (phasiRNAs) are poorly characterized. To remedy this, we sequenced small RNAs (sRNA) from vegetative and dissected inflorescence tissue in 28 phylogenetically diverse monocots and from several early-diverging angiosperm lineages, as well as publicly available data from 10 additional monocot species. We annotated miRNAs, small interfering RNAs (siRNAs) and phasiRNAs across the monocot phylogeny, identifying miRNAs apparently lost or gained in the grasses relative to other monocot families, as well as a number of transfer RNA fragments misannotated as miRNAs. Using our miRNA database cleaned of these misannotations, we identified conservation at the 8th, 9th, 19th, and 3'-end positions that we hypothesize are signatures of selection for processing, targeting, or Argonaute sorting. We show that 21-nucleotide (nt) reproductive phasiRNAs are far more numerous in grass genomes than other monocots. Based on sequenced monocot genomes and transcriptomes, DICER-LIKE5, important to 24-nt phasiRNA biogenesis, likely originated via gene duplication before the diversification of the grasses. This curated database of phylogenetically diverse monocot miRNAs, siRNAs, and phasiRNAs represents a large collection of data that should facilitate continued exploration of sRNA diversification in flowering plants.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

MS23, a master basic helix-loop-helix factor, regulates the specification and development of the tapetum in maize.

Successful male gametogenesis involves orchestration of sequential gene regulation for somatic differentiation in pre-meiotic anthers. We report here the cloning of Male Sterile23 (Ms23), encoding an anther-specific predicted basic helix-loop-helix (bHLH) transcription factor required for tapetal differentiation; transcripts localize initially to the precursor secondary parietal cells then predominantly to daughter tapetal cells. In knockout ms23-ref mutant anthers, five instead of the normal four wall layers are observed. Microarray transcript profiling demonstrates a more severe developmental disruption in ms23-ref than in ms32 anthers, which possess a different bHLH defect. RNA-seq and proteomics data together with yeast two-hybrid assays suggest that MS23 along with MS32, bHLH122 and bHLH51 act sequentially as either homo- or heterodimers to choreograph tapetal development. Among them, MS23 is the earliest-acting factor, upstream of bHLH51 and bHLH122, controlling tapetal specification and maturation. By contrast, MS32 is constitutive and independently regulated and is required later than MS23 in tapetal differentiation.

QTL-seq reveals genomic regions associated with spikelet fertility in response to a high temperature in rice (Oryza sativa L.).

KEY MESSAGE: The QTL-seq approach was used to identify QTLs for spikelet fertility under heat stress in rice. QTLs were detected on chromosomes 1, 2 and 3. Rice is a staple food of more than half of the global population. Rice production is increasingly affected by extreme environmental fluctuations caused by climate change. Increasing temperatures that exceed the optimum temperature adversely affect rice growth and development, especially during reproductive stages. Heat stress during the reproductive stages has a large effect on spikelet fertility; hence, the yield decreases. To sustain rice yields under increasing temperatures, the development of rice varieties for heat tolerance is necessary. In this study, we applied the QTL-seq approach to rapidly identify QTLs for spikelet fertility under heat stress (air temperature of 40-45 °C) based on two DNA pools, each consisting of 25 individual plants that exhibited a heat-tolerant or heat-sensitive phenotype from an F2 population of a cross between M9962 (heat tolerant) and Sinlek (heat sensitive). Three QTLs, qSF1, qSF2 and qSF3, were detected on chromosomes 1, 2 and 3, respectively, according to the highest contrasting SNP index between the two bulks. The QTLs identified in this study were found to overlap or were linked to QTLs previously identified in other crosses using conventional QTL mapping. A few highly abundant and anther-specific genes that contain nonsynonymous variants were identified within the QTLs and were proposed to be potential candidate genes. These genes could be targets in rice breeding programs for heat tolerance.

Efficiency and precision of microRNA biogenesis modes in plants.

Many evolutionarily conserved microRNAs (miRNAs) in plants regulate transcription factors with key functions in development. Hence, mutations in the core components of the miRNA biogenesis machinery cause strong growth defects. An essential aspect of miRNA biogenesis is the precise excision of the small RNA from its precursor. In plants, miRNA precursors are largely variable in size and shape and can be processed by different modes. Here, we optimized an approach to detect processing intermediates during miRNA biogenesis. We characterized a miRNA whose processing is triggered by a terminal branched loop. Plant miRNA processing can be initiated by internal bubbles, small terminal loops or branched loops followed by dsRNA segments of 15-17 bp. Interestingly, precision and efficiency vary with the processing modes. Despite the various potential structural determinants present in a single a miRNA precursor, DCL1 is mostly guided by a predominant structural region in each precursor in wild-type plants. However, our studies in fiery1, hyl1 and se mutants revealed the existence of cleavage signatures consistent with the recognition of alternative processing determinants. The results provide a general view of the mechanisms underlying the specificity of miRNA biogenesis in plants.

Spatiotemporally dynamic, cell-type-dependent premeiotic and meiotic phasiRNAs in maize anthers.

Maize anthers, the male reproductive floral organs, express two classes of phased small-interfering RNAs (phasiRNAs). PhasiRNA precursors are transcribed by RNA polymerase II and map to low-copy, intergenic regions similar to PIWI-interacting RNAs (piRNAs) in mammalian testis. From 10 sequential cohorts of staged maize anthers plus mature pollen we find that 21-nt phased siRNAs from 463 loci appear abruptly after germinal and initial somatic cell fate specification and then diminish, whereas 24-nt phasiRNAs from 176 loci coordinately accumulate during meiosis and persist as anther somatic cells mature and haploid gametophytes differentiate into pollen. Male-sterile ocl4 anthers defective in epidermal signaling lack 21-nt phasiRNAs. Male-sterile mutants with subepidermal defects--mac1 (excess meiocytes), ms23 (defective pretapetal cells), and msca1 (no normal soma or meiocytes)--lack 24-nt phasiRNAs. ameiotic1 mutants (normal soma, no meiosis) accumulate both 21-nt and 24-nt phasiRNAs, ruling out meiotic cells as a source or regulator of phasiRNA biogenesis. By in situ hybridization, miR2118 triggers of 21-nt phasiRNA biogenesis localize to epidermis; however, 21-PHAS precursors and 21-nt phasiRNAs are abundant subepidermally. The miR2275 trigger, 24-PHAS precursors, and 24-nt phasiRNAs all accumulate preferentially in tapetum and meiocytes. Therefore, each phasiRNA type exhibits independent spatiotemporal regulation with 21-nt premeiotic phasiRNAs dependent on epidermal and 24-nt meiotic phasiRNAs dependent on tapetal cell differentiation. Maize phasiRNAs and mammalian piRNAs illustrate putative convergent evolution of small RNAs in male reproduction.

An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes.

Small RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNA-dependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

A deletion of the gene encoding amino aldehyde dehydrogenase enhances the "pandan-like" aroma of winter melon (Benincasa hispida) and is a functional marker for the development of the aroma.

KEY MESSAGE: The gene conferring a "pandan-like" aroma of winter melon was identified. The sequence variation (804-bp deletion) found in the gene was used as the target for functional marker development. Winter melon (Benincasa hispida), a member of the Cucurbitaceae family, is a commonly consumed vegetable in Asian countries that is popular for its nutritional and medicinal value. A "pandan-like" aroma, which is economically important in crops including rice and soybean, is rarely found in most commercial varieties of winter melon, but is present in some landraces. This aroma is a value-added potential trait in breeding winter melon with a higher economic value. In this study, we confirmed that the aroma of winter melon is due to the potent volatile compound 2-acetyl-1-pyrroline (2AP) as previously identified in other plants. Based on an analysis of public transcriptome data, BhAMADH encoding an aminoaldehyde dehydrogenase (AMADH) was identified as a candidate gene conferring aroma of winter melon. A sequence comparison of BhAMADH between the aromatic and non-aromatic accessions revealed an 804-bp deletion encompassing exons 11-13 in the aromatic accession. The deletion caused several premature stop codons and could result in a truncated protein with a length of only 208 amino acids compared with 503 amino acids in the normal protein. A functional marker was successfully developed based on the 804-bp deletion and validated in 237 F2 progenies. A perfect association of the marker genotypes and aroma phenotypes indicates that BhAMADH is the major gene conferring the aroma. The recently developed functional marker could be efficiently used in breeding programs for the aroma trait in winter melon.

A new approach for annotation of transposable elements using small RNA mapping.

Transposable elements (TEs) are mobile genomic DNA sequences found in most organisms. They so densely populate the genomes of many eukaryotic species that they are often the major constituents. With the rapid generation of many plant genome sequencing projects over the past few decades, there is an urgent need for improved TE annotation as a prerequisite for genome-wide studies. Analogous to the use of RNA-seq for gene annotation, we propose a new method for de novo TE annotation that uses as a guide 24 nt-siRNAs that are a part of TE silencing pathways. We use this new approach, called TASR (for Transposon Annotation using Small RNAs), for de novo annotation of TEs in Arabidopsis, rice and soybean and demonstrate that this strategy can be successfully applied for de novo TE annotation in plants.Executable PERL is available for download from: http://tasr-pipeline.sourceforge.net/.

High-resolution identification and abundance profiling of cassava (Manihot esculenta Crantz) microRNAs.

BACKGROUND: Small RNAs (sRNAs) are endogenous sRNAs that play regulatory roles in plant growth, development, and biotic and abiotic stress responses. In plants, one subset of sRNAs, microRNAs (miRNAs) exhibit tissue-differential expression and regulate gene expression mainly through direct cleavage of mRNA or indirectly via production of secondary phased siRNAs (phasiRNAs) that silence cognate target transcripts in trans. RESULTS: Here, we have identified cassava (Manihot esculenta Crantz) miRNAs using high resolution sequencing of sRNA libraries from leaf, stem, callus, male and female flower tissues. To analyze the data, we built a cassava genome database and, via sequence analysis and secondary structure prediction, 38 miRNAs not previously reported in cassava were identified. These new cassava miRNAs included two miRNAs not previously been reported in any plant species. The miRNAs exhibited tissue-differential accumulation as confirmed by quantitative RT-PCR and Northern blot analysis, largely reflecting levels observed in sequencing data. Some of the miRNAs identified were predicted to trigger production of secondary phased siRNAs (phasiRNAs) from 80 PHAS loci. CONCLUSIONS: Cassava is a woody perennial shrub, grown principally for its starch-rich storage roots, which are rich in calories. In this study, new miRNAs were identified and their expression was validated using qRT-PCR of RNA from five different tissues. The data obtained expand the list of annotated miRNAs and provide additional new resources for cassava improvement research.

Extensive Families of miRNAs and PHAS Loci in Norway Spruce Demonstrate the Origins of Complex phasiRNA Networks in Seed Plants.

In eudicot plants, the miR482/miR2118 superfamily regulates and instigates the production of phased secondary small interfering RNAs (siRNAs) from NB-LRR (nucleotide binding leucine-rich repeat) genes that encode disease resistance proteins. In grasses, this miRNA family triggers siRNA production specifically in reproductive tissues from long noncoding RNAs. To understand this functional divergence, we examined the small RNA population in the ancient gymnosperm Norway spruce (Picea abies). As many as 41 miRNA families in spruce were found to trigger phasiRNA (phased, secondary siRNAs) production from diverse PHAS loci, with a remarkable 19 miRNA families capable of targeting over 750 NB-LRR genes to generate phasiRNAs. miR482/miR2118, encoded in spruce by at least 24 precursor loci, targets not only NB-LRR genes to trigger phasiRNA production (as in eudicots) but also noncoding PHAS loci, generating phasiRNAs preferentially in male or female cones, reminiscent of its role in the grasses. These data suggest a dual function of miR482/miR2118 present in gymnosperms that was selectively yet divergently retained in flowering plants. A few MIR482/MIR2118 precursors possess an extremely long stem-loop structure, one arm of which shows significant sequence similarity to spruce NB-LRR genes, suggestive of an evolutionary origin from NB-LRR genes through gene duplication. We also characterized an expanded miR390-TAS3 (TRANS-ACTING SIRNA GENE 3)-ARF (AUXIN RESPONSIVE FACTOR) pathway, comprising 18 TAS3 genes of diverse features. Finally, we annotated spruce miRNAs and their targets. Taken together, these data expand our understanding of phasiRNA network in plants and the evolution of plant miRNAs, particularly miR482/miR2118 and its functional diversification.

Plant microRNAs display differential 3' truncation and tailing modifications that are ARGONAUTE1 dependent and conserved across species.

Plant small RNAs are 3' methylated by the methyltransferase HUA1 ENHANCER1 (HEN1). In plant hen1 mutants, 3' modifications of small RNAs, including oligo-uridylation (tailing), are associated with accelerated degradation of microRNAs (miRNAs). By sequencing small RNAs of the wild type and hen1 mutants from Arabidopsis thaliana, rice (Oryza sativa), and maize (Zea mays), we found 3' truncation prior to tailing is widespread in these mutants. Moreover, the patterns of miRNA truncation and tailing differ substantially among miRNA families but are conserved across species. The same patterns are also observable in wild-type libraries from a broad range of species, only at lower abundances. ARGONAUTE (AGO1), even with defective slicer activity, can bind these truncated and tailed variants of miRNAs. An ago1 mutation in hen1 suppressed such 3' modifications, indicating that they occur while miRNAs are in association with AGO1, either during or after RNA-induced silencing complex assembly. Our results showed AGO1-bound miRNAs are actively 3' truncated and tailed, possibly reflecting the activity of cofactors acting in conserved patterns in miRNA degradation.

Roles of small RNAs in soybean defense against Phytophthora sojae infection.

The genus Phytophthora consists of many notorious pathogens of crops and forestry trees. At present, battling Phytophthora diseases is challenging due to a lack of understanding of their pathogenesis. We investigated the role of small RNAs in regulating soybean defense in response to infection by Phytophthora sojae, the second most destructive pathogen of soybean. Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are universal regulators that repress target gene expression in eukaryotes. We identified known and novel small RNAs that differentially accumulated during P. sojae infection in soybean roots. Among them, miR393 and miR166 were induced by heat-inactivated P. sojae hyphae, indicating that they may be involved in soybean basal defense. Indeed, knocking down the level of mature miR393 led to enhanced susceptibility of soybean to P. sojae; furthermore, the expression of isoflavonoid biosynthetic genes was drastically reduced in miR393 knockdown roots. These data suggest that miR393 promotes soybean defense against P. sojae. In addition to miRNAs, P. sojae infection also resulted in increased accumulation of phased siRNAs (phasiRNAs) that are predominantly generated from canonical resistance genes encoding nucleotide binding-leucine rich repeat proteins and genes encoding pentatricopeptide repeat-containing proteins. This work identifies specific miRNAs and phasiRNAs that regulate defense-associated genes in soybean during Phytophthora infection.

Identification of microRNAs and their mRNA targets during soybean nodule development: functional analysis of the role of miR393j-3p in soybean nodulation.

Plant microRNAs (miRNAs) play important regulatory roles in a number of developmental processes. The present work investigated the roles of miRNAs during nodule development in the crop legume soybean (Glycine max). Fifteen soybean small RNA libraries were sequenced from different stages of nodule development, including young nodules, mature nodules and senescent nodules. In order to identify the regulatory targets of the miRNAs, five parallel analysis of RNA ends (PARE) libraries were also sequenced from the same stages of nodule development. Sequencing identified 284 miRNAs, including 178 novel soybean miRNAs. Analysis of miRNA abundance identified 139 miRNAs whose expression was significantly regulated during nodule development, including 12 miRNAs whose expression changed > 10-fold. Analysis of the PARE libraries identified 533 miRNA targets, including three nodulation-related genes and eight nodule-specific genes. miR393j-3p was selected for detailed analysis as its expression was significantly regulated during nodule formation, and it targeted a nodulin gene, Early Nodulin 93 (ENOD93). Strong, ectopic expression of miR393j-3p, as well as RNAi silencing of ENOD93 expression, significantly reduced nodule formation. The data indicate that miR393j-3p regulation of ENOD93 mRNA abundance is a key control point for soybean nodule formation.

Rapid construction of parallel analysis of RNA end (PARE) libraries for Illumina sequencing.

MicroRNAs (miRNAs) are ∼21nt small RNAs that pair to their target mRNAs and in many cases trigger cleavage, particularly in plants. Although many computational tools can predict miRNA:mRNA interactions, it remains critical to validate cleavage events, due to miRNA function in translational repression or due to high rates of false positives (over 90%) for unvalidated target predictions. A few years ago, three laboratories described similar methods to validate cleavage of miRNA targets by the cloning en masse of 5' ends of cleaved or uncapped mRNAs. To take advantage of the recent progress in high-throughput sequencing technology, we have devised an updated protocol to (1) enable much faster library preparation, and (2) reduce the cost by pooling indexed samples together for sequencing. Here we provide a step-by-step protocol for PARE library construction, starting from total RNA. This protocol has been successfully used in our laboratory to validate miRNA targets in a variety of plant species. We also provide advice for troubleshooting on some common issues.

Rice mutants, selected under severe drought stress, show reduced stomatal density and improved water use efficiency under restricted water conditions.

Introduction: Rice is among the least water-use-efficient crops, and rice plants utilise most of their water uptake for transpirational cooling via stomata. To improve water-use efficiency (WUE) in rice, reducing stomatal density and size could help optimise transpiration and photosynthesis. Methodology: In this study, we compared two series of purple rice stomata mutants: the Stomatal Model Mutant (SMM) identified by microscopic observation of flag-leaf stomata, and the Drought-selected Model Mutant (DMM) generated through screening under severe water stress. After undergoing two rounds of severe water stress between -60 to -80 Ym, right before the R1-2 reproductive stage, three DMMs were selected based on their rapid recovery rate and % filled-grain percentage. Result: The three DMMs displayed 618-697 stomatal units per mm2, similar to the SMMs low-density stomata mutant (JHN 8756 (LD)). Furthermore, the four SMMs, three DMMs and the Jao Hom Nin wild type (JHN WT) were treated with two restricted water condition schemes from seedlings to harvest. The total amount of irrigation and precipitation during the experiment was 78.1 L/plant (69.1 mm/plant) for the less restricted water condition (LR) and 47.5 L/plant (42 mm/plant) for the more restricted water condition (MR). Water condition treatments had no effects on stomatal density and stomatal index. In contrast, genotypes and restricted water condition schemes affected plant height, tillers/plant, % filled grains and shoot dry weight (SDW). The three DMMs and the JHN 8756 (LD), the SMM's low-density stomata mutant, displayed greater resilience towards more restricted water conditions than the SMMs and the JHN wild type. Particularly, DMMs were tolerant to more restricted water condition treatments, showing no SDW penalties. Together, the DMMs and the JHN 8756 (LD) displayed higher WUE under these conditions of more restricted water conditions. Conclusion: A rigorous screening process to distinguish tolerant mutants with a rapid drought recovery rate from severe water stress could pave the way to isolate more mutants with better stomatal functionality and resilience in preparation for imminent climate changes.

Accelerating haploid induction rate and haploid validation through marker-assisted selection for qhir1 and qhir8 in maize.

Doubled haploid (DH) technology becomes more routinely applied in maize hybrid breeding. However, some issues in haploid induction and identification persist, requiring resolution to optimize DH production. Our objective was to implement simultaneous marker-assisted selection (MAS) for qhir1 (MTL/ZmPLA1/NLD) and qhir8 (ZmDMP) using TaqMan assay in F2 generation of four BHI306-derived tropical × temperate inducer families. We also aimed to assess their haploid induction rate (HIR) in the F3 generation as a phenotypic response to MAS. We highlighted remarkable increases in HIR of each inducer family. Genotypes carrying qhir1 and qhir8 exhibited 1 - 3-fold higher haploid frequency than those carrying only qhir1. Additionally, the qhir1 marker was employed for verifying putative haploid seedlings at 7 days after planting. Flow cytometric analysis served as the gold standard test to assess the accuracy of the R1-nj and the qhir1 marker. The qhir1 marker showed high accuracy and may be integrated in multiple haploid identifications at early seedling stage succeeding pre-haploid sorting via R1-nj marker.

Candidate genes affecting stomatal density in rice (Oryza sativa L.) identified by genome-wide association.

Stomata regulate photosynthesis and water loss. They have been an active subject of research for centuries, but our knowledge of the genetic components that regulate stomatal development in crops remains very limited in comparison to the model plant Arabidopsis thaliana. Leaf stomatal density was found to vary by over 2.5-fold across a panel of 235 rice accessions. Using GWAS, we successfully identified five different QTLs associated with stomatal density on chromosomes 2, 3, 9, and 12. Forty-two genes were identified within the haplotype blocks corresponding to these QTLs. Of these, nine genes contained haplotypes that were associated with different stomatal densities. These include a gene encoding a trehalose-6-phosphate synthase, an enzyme that has previously been associated with altered stomatal density in Arabidopsis, and genes encoding a B-BOX zinc finger family protein, a leucine-rich repeat family protein, and the 40 S ribosomal protein S3a, none of which have previously been linked to stomatal traits. We investigated further and show that a closely related B-BOX protein regulates stomatal development in Arabidopsis. The results of this study provide information on genetic associations with stomatal density in rice. The QTLs and candidate genes may be useful in future breeding programs for low or high stomatal density and, consequently, improved photosynthetic capacity, water use efficiency, or drought tolerance.

Genetic Divergence of Thai Indigenous Pigs from Three Distinct Geographic Regions Revealed by Microsatellite Marker Analysis.

Thai indigenous pigs (TIPs) are important genetic resources. Crosses with exotic pig breeds and wild boars may cause genetic losses. To date, the physical characteristics of TIPs have been inconsistent. The classification of TIPs by genetic information is needed to pursue an appropriate conservation program. In this study, the genetic diversity, cluster analysis, and phylogenetic relationship of TIPs were investigated using twenty-nine pig microsatellite markers. Blood samples were collected from TIPs from three regions of Thailand: north (NT, n = 118), northeast (NE, n = 61), and south (ST, n = 75). The mean total number of distinct alleles and the effective number of alleles per locus were 11.851 and 5.497, respectively. The mean observed heterozygosity (Ho) and mean expected heterozygosity (He) were 0.562 and 0.837, respectively. The F values of the microsatellite loci were positive under Hardy-Weinberg Equilibrium at p < 0.001, with overall mean values of Fis, Fit, and Fst of 0.247, 0.281, and 0.046, respectively. A total of 5, 5, and 17 private alleles were found at frequencies greater than 0.050 in the NT, NE, and ST pigs, respectively. Three optimal clusters (K = 3) were proposed within the TIP populations. Pigs from the NT and NE regions were mixed in two clusters, while members of the ST region were clearly separated. The phylogenetic tree confirmed that the pigs from NT and NE were each divided into two subgroups, while the pigs from ST were clustered into one group. A microsatellite analysis revealed the high genetic diversity of the TIP populations and confirmed the genetic divergence of the TIPs from the European and Chinese breeds. A genetic admixture of the TIP with the local wild boars was detected.