QTL Intervals Research Articles

MQTL-seq delineates functionally relevant candidate gene harbouring a major QTL regulating pod number in chickpea.

The present study used a whole-genome, NGS resequencing-based mQTL-seq (multiple QTL-seq) strategy in two inter-specific mapping populations (Pusa 1103 × ILWC 46 and Pusa 256 × ILWC 46) to scan the major genomic region(s) underlying QTL(s) governing pod number trait in chickpea. Essentially, the whole-genome resequencing of low and high pod number-containing parental accessions and homozygous individuals (constituting bulks) from each of these two mapping populations discovered >8 million high-quality homozygous SNPs with respect to the reference kabuli chickpea. The functional significance of the physically mapped SNPs was apparent from the identified 2,264 non-synonymous and 23,550 regulatory SNPs, with 8–10% of these SNPs-carrying genes corresponding to transcription factors and disease resistance-related proteins. The utilization of these mined SNPs in Δ (SNP index)-led QTL-seq analysis and their correlation between two mapping populations based on mQTL-seq, narrowed down two (CaqaPN4.1: 867.8 kb and CaqaPN4.2: 1.8 Mb) major genomic regions harbouring robust pod number QTLs into the high-resolution short QTL intervals (CaqbPN4.1: 637.5 kb and CaqbPN4.2: 1.28 Mb) on chickpea chromosome 4. The integration of mQTL-seq-derived one novel robust QTL with QTL region-specific association analysis delineated the regulatory (C/T) and coding (C/A) SNPs-containing one pentatricopeptide repeat (PPR) gene at a major QTL region regulating pod number in chickpea. This target gene exhibited anther, mature pollen and pod-specific expression, including pronounced higher up-regulated (∼3.5-folds) transcript expression in high pod number-containing parental accessions and homozygous individuals of two mapping populations especially during pollen and pod development. The proposed mQTL-seq-driven combinatorial strategy has profound efficacy in rapid genome-wide scanning of potential candidate gene(s) underlying trait-associated high-resolution robust QTL(s), thereby expediting genomics-assisted breeding and genetic enhancement of crop plants, including chickpea.

Identification of the Bile Acid Transporter Slco1a6 as a Candidate Gene That Broadly Affects Gene Expression in Mouse Pancreatic Islets.

We surveyed gene expression in six tissues in an F2 intercross between mouse strains C57BL/6J (abbreviated B6) and BTBR T(+) tf/J (abbreviated BTBR) made genetically obese with the Leptin(ob) mutation. We identified a number of expression quantitative trait loci (eQTL) affecting the expression of numerous genes distal to the locus, called trans-eQTL hotspots. Some of these trans-eQTL hotspots showed effects in multiple tissues, whereas some were specific to a single tissue. An unusually large number of transcripts (∼8% of genes) mapped in trans to a hotspot on chromosome 6, specifically in pancreatic islets. By considering the first two principal components of the expression of genes mapping to this region, we were able to convert the multivariate phenotype into a simple Mendelian trait. Fine mapping the locus by traditional methods reduced the QTL interval to a 298-kb region containing only three genes, including Slco1a6, one member of a large family of organic anion transporters. Direct genomic sequencing of all Slco1a6 exons identified a nonsynonymous coding SNP that converts a highly conserved proline residue at amino acid position 564 to serine. Molecular modeling suggests that Pro564 faces an aqueous pore within this 12-transmembrane domain-spanning protein. When transiently overexpressed in HEK293 cells, BTBR organic anion transporting polypeptide (OATP)1A6-mediated cellular uptake of the bile acid taurocholic acid (TCA) was enhanced compared to B6 OATP1A6. Our results suggest that genetic variation in Slco1a6 leads to altered transport of TCA (and potentially other bile acids) by pancreatic islets, resulting in broad gene regulation.

Genome-wide insertion-deletion (InDel) marker discovery and genotyping for genomics-assisted breeding applications in chickpea.

We developed 21,499 genome-wide insertion–deletion (InDel) markers (2- to 54-bp in silico fragment length polymorphism) by comparing the genomic sequences of four (desi, kabuli and wild C. reticulatum) chickpea [Cicer arietinum (L.)] accessions. InDel markers showing 2- to 6-bp fragment length polymorphism among accessions were abundant (76.8%) in the chickpea genome. The physically mapped 7,643 and 13,856 markers on eight chromosomes and unanchored scaffolds, respectively, were structurally and functionally annotated. The 4,506 coding (23% large-effect frameshift mutations) and regulatory InDel markers were identified from 3,228 genes (representing 11.7% of total 27,571 desi genes), suggesting their functional relevance for trait association/genetic mapping. High amplification (97%) and intra-specific polymorphic (60–83%) potential and wider genetic diversity (15–89%) were detected by genome-wide 6,254 InDel markers among desi, kabuli and wild accessions using even a simpler cost-effective agarose gel-based assay. This signifies added advantages of this user-friendly genetic marker system for manifold large-scale genotyping applications in laboratories with limited infrastructure and resources. Utilizing 6,254 InDel markers-based high-density (inter-marker distance: 0.212 cM) inter-specific genetic linkage map (ICC 4958 × ICC 17160) of chickpea as a reference, three major genomic regions harboring six flowering and maturity time robust QTLs (16.4–27.5% phenotypic variation explained, 8.1–11.5 logarithm of odds) were identified. Integration of genetic and physical maps at these target QTL intervals mapped on three chromosomes delineated five InDel markers-containing candidate genes tightly linked to the QTLs governing flowering and maturity time in chickpea. Taken together, our study demonstrated the practical utility of developing and high-throughput genotyping of such beneficial InDel markers at a genome-wide scale to expedite genomics-assisted breeding applications in chickpea.

Mapping of seedling resistance in barley to Puccinia striiformis f. sp. pseudohordei.

The barley grass stripe rust (BGYR) pathogen Puccinia striiformis f. sp. pseudohordei was first detected in Australia in 1997. While studies have established that it is virulent on wild barley grass, and can infect several barley cultivars, the basis of genetic resistance to this pathogen in barley is largely unknown. Understanding the genetic basis of host resistance and ensuring the selection of germplasm with multiple resistance genes are important to mitigate the potential impact of BGYR in barley production. Genetic analysis of seedling resistance to BGYR in two barley doubled haploid populations, Amaji Nijo/WI2585 (AN/WI) and Galleon/Haruna Nijo (GL/HN), indicated that resistance is governed by several genes. Marker regression analysis of the seedling resistance data from the AN/WI population detected a major QTL, BGYR_WI1 (resistance contributed by WI2585 with the closest marker explaining 52 % of the total phenotypic effect) on chromosome 1HS, flanked by the loci Xabg59 and Xabc310b at map positions 0.0 and 6.9 cM, respectively. Similarly, a major QTL, BGYR_HN1, (resistance contributed by Haruna Nijo with the closest marker explaining 70 % of the total phenotypic effect) was detected in the GL/HN population and was mapped to 1HS, flanked by the loci Xbcd135 and XHOR1 at map positions 12.8 and 24.5 cM, respectively. In addition, several minor loci that provided resistance against BGYR were detected in both populations. While defined QTL intervals were large, the analysis nonetheless provides new information on sources of major QTL controlling resistance to BGYR.

Targeted Recombinant Progeny: a design for ultra-high resolution mapping of Quantitative Trait Loci in crosses between inbred or pure lines.

BackgroundHigh-resolution mapping of the loci (QTN) responsible for genetic variation in quantitative traits is essential for positional cloning of candidate genes, and for effective marker assisted selection. The confidence interval (QTL) flanking the point estimate of QTN-location is proportional to the number of individuals in the mapping population carrying chromosomes recombinant in the given interval. Consequently, many designs for high resolution QTN mapping are based on increasing the proportion of recombinants in the mapping population.The "Targeted Recombinant Progeny" (TRP) design is a new design for high resolution mapping of a target QTN in crosses between pure, or inbred lines. It is a three-generation procedure generating a large number of recombinant individuals within a QTL previously shown to contain a QTN. This is achieved by having individuals that carry chromosomes recombinant across the target QTL interval as parents of a large mapping population; most of whom will therefore carry recombinant chromosomes targeted to the given QTL. The TRP design is particularly useful for high resolution mapping of QTN that differentiate inbred or pure lines, and hence are not amenable to high resolution mapping by genome-wide association tests.ResultsIn the absence of residual polygenic variation, population sizes required for achieving given mapping resolution by the TRP-F2 design relative to a standard F2 design ranged from 0.289 for a QTN with standardized allele substitution effect = 0.2, mapped to an initial QTL of 0.2 Morgan to 0.041 for equivalent QTN mapped to an initial QTL of 0.02 M. In the presence of residual polygenic variation, the relative effectiveness of the TRP design ranges from 1.068 to 0.151 for the same initial QTL intervals and QTN effect. Thus even in the presence of polygenic variation, the TRP can still provide major savings. Simulation showed that mapping by TRP should be based on 30-50 markers spanning the initial interval; and on at least 50 or more G2 families representing this number of recombination points,.ConclusionsThe TRP design can be an effective procedure for achieving high and ultra-high mapping resolution of a target QTN previously mapped to a known confidence interval (QTL).

OsGRAS23, a rice GRAS transcription factor gene, is involved in drought stress response through regulating expression of stress-responsive genes.

BackgroundDrought is a major abiotic stress factors that reduces agricultural productivity. GRAS transcription factors are plant-specific proteins that play diverse roles in plant development. However, the functions of a number of GRAS genes identified in rice are unknown, especially the GRAS genes related to rice drought resistance have not been characterized.ResultsIn this study, a novel GRAS transcription factor gene named OsGRAS23, which is located in a drought-resistant QTL interval on chromosome 4 of rice, was isolated. The expression of OsGRAS23 was induced by drought, NaCl, and jasmonic acid treatments. The OsGRAS23-GFP fused protein was localized in the nucleus of tobacco epidermal cells. A trans-activation assay in yeast cells demonstrated that the OsGRAS23 protein possessed a strong transcriptional activation activity. OsGRAS23-overexpressing rice plants showed improved drought resistance and oxidative stress tolerance as well as less H2O2 accumulation compared with the wild-type plants. Furthermore, microarray analysis showed that several anti-oxidation related genes were up-regulated in the OsGRAS23-overexpressing rice plants. The yeast one hybrid test indicated that OsGRAS23 could bind to the promoters of its potential target genes.ConclusionsOur results demonstrate that OsGRAS23 encodes a stress-responsive GRAS transcription factor and positively modulates rice drought tolerance via the induction of a number of stress-responsive genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0532-3) contains supplementary material, which is available to authorized users.

Genetic analysis of inflorescence and plant height components in sorghum (Panicoidae) and comparative genetics with rice (Oryzoidae).

BackgroundDomestication has played an important role in shaping characteristics of the inflorescence and plant height in cultivated cereals. Taking advantage of meta-analysis of QTLs, phylogenetic analyses in 502 diverse sorghum accessions, GWAS in a sorghum association panel (n = 354) and comparative data, we provide insight into the genetic basis of the domestication traits in sorghum and rice.ResultsWe performed genome-wide association studies (GWAS) on 6 traits related to inflorescence morphology and 6 traits related to plant height in sorghum, comparing the genomic regions implicated in these traits by GWAS and QTL mapping, respectively. In a search for signatures of selection, we identify genomic regions that may contribute to sorghum domestication regarding plant height, flowering time and pericarp color. Comparative studies across taxa show functionally conserved ‘hotspots’ in sorghum and rice for awn presence and pericarp color that do not appear to reflect corresponding single genes but may indicate co-regulated clusters of genes. We also reveal homoeologous regions retaining similar functions for plant height and flowering time since genome duplication an estimated 70 million years ago or more in a common ancestor of cereals. In most such homoeologous QTL pairs, only one QTL interval exhibits strong selection signals in modern sorghum.ConclusionsIntersections among QTL, GWAS and comparative data advance knowledge of genetic determinants of inflorescence and plant height components in sorghum, and add new dimensions to comparisons between sorghum and rice.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0477-6) contains supplementary material, which is available to authorized users.

Identification of novel QTL governing root architectural traits in an interspecific soybean population.

Cultivated soybean (Glycine max L.) cv. Dunbar (PI 552538) and wild G. soja (PI 326582A) exhibited significant differences in root architecture and root-related traits. In this study, phenotypic variability for root traits among 251 BC2F5 backcross inbred lines (BILs) developed from the cross Dunbar/PI 326582A were identified. The root systems of the parents and BILs were evaluated in controlled environmental conditions using a cone system at seedling stage. The G. max parent Dunbar contributed phenotypically favorable alleles at a major quantitative trait locus on chromosome 8 (Satt315-I locus) that governed root traits (tap root length and lateral root number) and shoot length. This QTL accounted for >10% of the phenotypic variation of both tap root and shoot length. This QTL region was found to control various shoot- and root-related traits across soybean genetic backgrounds. Within the confidence interval of this region, eleven transcription factors (TFs) were identified. Based on RNA sequencing and Affymetrix expression data, key TFs including MYB, AP2-EREBP and bZIP TFs were identified in this QTL interval with high expression in roots and nodules. The backcross inbred lines with different parental allelic combination showed different expression pattern for six transcription factors selected based on their expression pattern in root tissues. It appears that the marker interval Satt315–I locus on chromosome 8 contain an essential QTL contributing to early root and shoot growth in soybean.

High-resolution genetic mapping of allelic variants associated with cell wall chemistry in Populus.

BackgroundQTL cloning for the discovery of genes underlying polygenic traits has historically been cumbersome in long-lived perennial plants like Populus. Linkage disequilibrium-based association mapping has been proposed as a cloning tool, and recent advances in high-throughput genotyping and whole-genome resequencing enable marker saturation to levels sufficient for association mapping with no a priori candidate gene selection. Here, multiyear and multienvironment evaluation of cell wall phenotypes was conducted in an interspecific P. trichocarpa x P. deltoides pseudo-backcross mapping pedigree and two partially overlapping populations of unrelated P. trichocarpa genotypes using pyrolysis molecular beam mass spectrometry, saccharification, and/ or traditional wet chemistry. QTL mapping was conducted using a high-density genetic map with 3,568 SNP markers. As a fine-mapping approach, chromosome-wide association mapping targeting a QTL hot-spot on linkage group XIV was performed in the two P. trichocarpa populations. Both populations were genotyped using the 34 K Populus Infinium SNP array and whole-genome resequencing of one of the populations facilitated marker-saturation of candidate intervals for gene identification.ResultsFive QTLs ranging in size from 0.6 to 1.8 Mb were mapped on linkage group XIV for lignin content, syringyl to guaiacyl (S/G) ratio, 5- and 6-carbon sugars using the mapping pedigree. Six candidate loci exhibiting significant associations with phenotypes were identified within QTL intervals. These associations were reproducible across multiple environments, two independent genotyping platforms, and different plant growth stages. cDNA sequencing for allelic variants of three of the six loci identified polymorphisms leading to variable length poly glutamine (PolyQ) stretch in a transcription factor annotated as an ANGUSTIFOLIA C-terminus Binding Protein (CtBP) and premature stop codons in a KANADI transcription factor as well as a protein kinase. Results from protoplast transient expression assays suggested that each of the polymorphisms conferred allelic differences in the activation of cellulose, hemicelluloses, and lignin pathway marker genes.ConclusionThis study illustrates the utility of complementary QTL and association mapping as tools for gene discovery with no a priori candidate gene selection. This proof of concept in a perennial organism opens up opportunities for discovery of novel genetic determinants of economically important but complex traits in plants.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1215-z) contains supplementary material, which is available to authorized users.

Genetic Analysis of Immunity to Soybean Cyst Nematode Race 3 in Elite Line Zhongpin 03-5373

Soybean cyst nematode(SCN) race 3, one of the eight races identified in China, is widely distributed and severely reduced soybean yield. Zhongpin 03-5373(ZP03-5373) is an elite line immune to SCN race 3. In this study, a recombinant inbred line(RIL) population was developed from a cross between ZP03-5373 and Zhonghuang 13(ZH13). A genetic linkage map was constructed using a total of 506 molecular markers, including SSRs, EST-SSRs, In Del, and SNP. The total length of the genetic map was 2651.90 c M with an average marker spacing of 5.24 c M. Based on the phenotyping data, we detected three QTL intervals to dominate SCN3, including SCN3-7(Gm07), SCN3-11(Gm11), and SCN3-18(Gm11). Main effect QTL, SCN3-18 could explain 29.5% of resistant variation. Two minor effect QTLs, SCN3-7 and SCN3-11, explained 6.2% and 5.5% of resistant variation, respectively. And it further showed there was significant epistatic interaction between SCN3-7 and SCN3-18 for resistance to SCN3. Both SCN3-7 and SCN3-11 were confirmed to be resistant to SCN3 by tracking flanking marker in the ancestors of ZP03-5373. These markers will be helpful for developing SCN resistant cultivars and cloning resistant genes by marker assisted selection.

Phenomics: A Challenge for Crop Improvement in Genomic Era

The last two decades have observed tremendous progress in the genomics for plant breeding research, especially the availability of a large number of high-throughput cost-effective molecular markers and genotyping platforms, advances in sequencing technologies leading the reduction in sequencing costs and making available genome sequences of major crop genomes, etc. But the advance in phenomics is lagging far behind our present day capacity to generate high-throughput molecular genotyping data, thus creating phenotypic bottleneck”. The accurate and precise phenotypic data is an essential component in the discovery of genes/QTLs of important agronomic traits using modern genomic approaches; otherwise it leads to false positives and false negatives. Therefore, it necessitates the development of high-throughput phenotyping facilities. Although, several phenotyping facilities have been developed around the world that can scan and record precise and accurate data for thousands of plants quickly by making use of non-invasive imaging, spectroscopy, image analysis, robotics and high-performance computing facilities, but more efforts and funds are required to be allocated in this field to achieve fruitful results from genomics/molecular breeding approaches like QTL interval mapping, association mapping, genome-wide association studies (GWAS), QTL cloning, QTL validation, marker-assisted selection (MAS), marker-assisted recurrent selection (MARS), TILLING (Targeting Induced Local Lesions in Genomes) and genomic selection (GS) or genome-wide selection (GWS).

Resistances to anthracnose (Colletotrichum acutatum) of Capsicum mature green and ripe fruit are controlled by a major dominant cluster of QTLs on chromosome P5

Anthracnose (Colletotrichum spp.) is a serious disease worldwide in pepper (Capsicum) production. Inheritance of resistance to Colletotrichum acutatum from a Capsicum chinense accession (PBC932) was studied in a BC1 population derived from a hybrid with Capsicum annuum line 77013 (susceptible) using a QTL analysis method. Resistance test was performed on detached mature green and mature red fruit under laboratory conditions by evaluated in disease incidence, true lesion diameter and overall lesion diameter. Based on a linkage map with 14 linkage groups, 385 markers (SSR, InDel and CAPS), 1310.2cM in length, inclusive Composite Interval Mapping (ICIM) revealed main effect QTLs located in a close marker interval on P5 chromosome for all fruit stages and resistance criteria, and four minor-effect QTLs only at green mature stage. Identification of recombinant individuals suggested that resistance in green versus red fruit may be controlled by distinct genes within the QTL interval on P5.

Genome-wide association mapping of quantitative resistance to sudden death syndrome in soybean.

BackgroundSudden death syndrome (SDS) is a serious threat to soybean production that can be managed with host plant resistance. To dissect the genetic architecture of quantitative resistance to the disease in soybean, two independent association panels of elite soybean cultivars, consisting of 392 and 300 unique accessions, respectively, were evaluated for SDS resistance in multiple environments and years. The two association panels were genotyped with 52,041 and 5,361 single nucleotide polymorphisms (SNPs), respectively. Genome-wide association mapping was carried out using a mixed linear model that accounted for population structure and cryptic relatedness.ResultA total of 20 loci underlying SDS resistance were identified in the two independent studies, including 7 loci localized in previously mapped QTL intervals and 13 novel loci. One strong peak of association on chromosome 18, associated with all disease assessment criteria across the two panels, spanned a physical region of 1.2 Mb around a previously cloned SDS resistance gene (GmRLK18-1) in locus Rfs2. An additional variant independently associated with SDS resistance was also found in this genomic region. Other peaks were within, or close to, sequences annotated as homologous to genes previously shown to be involved in plant disease resistance. The identified loci explained an average of 54.5% of the phenotypic variance measured by different disease assessment criteria.ConclusionsThis study identified multiple novel loci and refined the map locations of known loci related to SDS resistance. These insights into the genetic basis of SDS resistance can now be used to further enhance durable resistance to SDS in soybean. Additionally, the associations identified here provide a basis for further efforts to pinpoint causal variants and to clarify how the implicated genes affect SDS resistance in soybean.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-809) contains supplementary material, which is available to authorized users.

Microarray expression analysis of the main inflorescence in Brassica napus.

The effect of the number of pods on the main inflorescence (NPMI) on seed yield in Brassica napus plants grown at high density is a topic of great economic and scientific interest. Here, we sought to identify patterns of gene expression that determine the NPMI during inflorescence differentiation. We monitored gene expression profiles in the main inflorescence of two B. napus F6 RIL pools, each composed of nine lines with a low or high NPMI, and their parental lines, Zhongshuang 11 (ZS11) and 73290, using a Brassica 90K elements oligonucleotide array. We identified 4,805 genes that were differentially expressed (≥1.5 fold-change) between the low- and high-NPMI samples. Of these, 82.8% had been annotated and 17.2% shared no significant homology with any known genes. About 31 enriched GO clusters were identified amongst the differentially expressed genes (DEGs), including those involved in hormone responses, development regulation, carbohydrate metabolism, signal transduction, and transcription regulation. Furthermore, 92.8% of the DEGs mapped to chromosomes that originated from B. rapa and B. oleracea, and 1.6% of the DEGs co-localized with two QTL intervals (PMI10 and PMI11) known to be associated with the NPMI. Overexpression of BnTPI, which co-localized with PMI10, in Arabidopsis suggested that this gene increases the NPMI. This study provides insight into the molecular factors underlying inflorescence architecture, NPMI determination and, consequently, seed yield in B. napus.

NEW INSIGHT INTO THE GENETICS OF COLOR IN GRAPE

Grapes may either be white or colored, ranging from the lightest pink to the darkest purple tone depending on the quantity and composition of anthocyanins accumulated in the berry skin. QTL mapping on cross-derived progeny and associa-tion mapping on natural diversity were combined to decipher the genetic architecture of anthocyanin content and composition in berry skin, a crucial trait for both wine quality and human nutrition. For the total anthocyanin content in the skin of mature berries on a segregating population of 191 progenies, we identified 2 QTLs accounting for 62% and 7.1% of the total variation. Further studying candidate genes from these QTL intervals on a core collection of natural resources (141 individuals), we identified Myb genes highly associated with total content of anthocyanins. Using a multivariate regression method, we demonstrated that four polymorphisms accounted for 86% of the observed variation, with all these polymorphisms leading to structural changes in the MYB proteins. Skin color tones are not only due to the overall pigment concentration but also to variation in pigment composition, mostly due to hydroxylation, methylation and acylation differences. For all these variations, the major QTL on LG2 played a key role. In addition, specific QTLs were found for each of the composition variable: on LG6, explaining 13% of the hydroxylation variation, on LG1, explaining up to 17% of the methylation variation and on LG7 and LG13, explaining 11 and 7% of the variation in p-coumaric acylation. For the hydroxylation level, the specific QTL on LG6 co-localized with candidate VvF3?5?H genes and two genic SNPs within two of those genes were identified as associated to the variation. For the methylation variation, the specific QTL on LG1 co-localized with a cluster of O-methyltransferase genes among which VvAOMT3.2 showed a strong association.

An Expression QTL of Closely Linked Candidate Genes Affects pH of Meat in Chickens

Genetical genomics has been suggested as a powerful approach to study the genotype-phenotype gap. However, the relatively low power of these experiments (usually related to the high cost) has hindered fulfillment of its promise, especially for loci (QTL) of moderate effects.One strategy with which to overcome the issue is to use a targeted approach. It has two clear advantages: (i) it reduces the problem to a simple comparison between different genotypic groups at the QTL and (ii) it is a good starting point from which to investigate downstream effects of the QTL. In this study, from 698 F2 birds used for QTL mapping, gene expression profiles of 24 birds with divergent homozygous QTL genotypes were investigated. The targeted QTL was on chromosome 1 and affected initial pH of breast muscle. The biological mechanisms controlling this trait can be similar to those affecting malignant hyperthermia or muscle fatigue in humans. The gene expression study identified 10 strong local signals that were markedly more significant compared to any genes on the rest of the genome. The differentially expressed genes all mapped to a region <1 Mb, suggesting a remarkable reduction of the QTL interval. These results, combined with analysis of downstream effect of the QTL using gene network analysis, suggest that the QTL is controlling pH by governing oxidative stress. The results were reproducible with use of as few as four microarrays on pooled samples (with lower significance level). The results demonstrate that this cost-effective approach is promising for characterization of QTL.

Expression patterns of microRNAs in different organs and developmental stages of a superhybrid rice LYP9 and its parental lines

Heterosis is an important phenomenon, and the molecular mechanisms underlying heterosis are still enigmatic. microRNAs (miRNAs) play vital roles in many aspects of plant development. A set of miRNAs was selected to investigate the roles of miRNAs in heterosis displayed in a superhybrid rice. We analysed the expression patterns of miRNAs in different organs and developmental stages of the superhybrid rice and its parental lines. All possible modes of miRNA action were observed, including additive, high- and low-parent value, above high- and below low-parent value. Different organs and developmental stages exhibited different modes of miRNA expression. Overall, the non-additive mode is the predominant expression pattern of miRNAs observed in this superhybrid. Many heterotic QTL intervals harbour miRNAs, whose expression patterns reveal their specific roles in different organs and developmental stages. miRNAs regulate the expression levels of target genes that have important functions in plant development. The predominant non-additive mode of miRNA expression pattern in the hybrid suggests that miRNAs play critical roles in hybrid development, in particular, those miRNAs located in the heterotic QTL intervals may have important roles in heterosis. Our research sheds new light on understanding of the molecular mechanisms of heterosis.

Fox domestication: Molecular mechanisms involved in selection for behavior

Although animal domestication has enjoyed the attention of geneticists and evolutionary biologists since Darwin’s times, the focal question about the genetic basis of this process has not been duly considered. We present some results of experimental modeling of historical domestication in silver foxes (Vulpes vulpes), an object of commercial breeding. Attention is focused on the role of artificial selection in the transformation of the behavior of these animals, which has brought about unique populations of tame and aggressive foxes. Additional resources were developed to analyze the molecular nature of differences in the behavior of these unique foxes: the meiotic map of the fox genome was constructed, informative segregating subpopulations were obtained by crosses of tame and aggressive animals, and the formerly applied method of behavioral phenotype quantification was improved. The integrated behavioral phenotypes (principal components PC1 and PC2) used in the study were obtained by analysis of ethological parameters recorded with a camcorder. The most important result is that the region most closely associated with tame behavior was identified on fox chromosome 12 by QTL interval mapping. The principle significance of this result is strengthened by the fact that the region is similar to the region on dog chromosome 5, presumed to be responsible for the early domestication of wolves and their evolutionary transformation to primitive dogs.

Identification and characterization of a repertoire of genes differentially expressed in developing top ear shoots between a superior hybrid and its parental inbreds in Zea mays L.

Heterosis has been widely used in crop breeding and production; however, little is known about the genes controlling trait heterosis. The shortage of genes known to function in heterosis significantly limits our understanding of the molecular basis underlying heterosis. Here, we report 748 genes differentially expressed (DG) in the developing top ear shoots between a maize heterotic F1 hybrid (Mo17×B73) and its parental inbreds identified using maize microarrays containing 28,608 unigene features. Of the 748 DG, over 600 were new for the inbred and hybrid combination. The DG were enriched for 35 of the total 213 maize gene ontology (GO) terms, including those describing photosynthesis, respiration, DNA replication, metabolism, and hormone biosynthesis. From the DG, we identified six genes involved in glycolysis, three genes in the citrate cycle, and four genes in the C4-dicarboxylic acid cycle. We mapped 533 of the 748 DG to the maize B73 genome, 298 (55.9%) of which mapped to the QTL intervals of 11 maize ear traits. Moreover, we compared the repertoire of the DG with that of 14-day seedlings of the same inbred and hybrid combination. Only approximately 5% of the DG was shared between the two organs and developmental stages. Furthermore, we mapped 417 (55.7%) of the 748 maize DG to the QTL intervals of 26 rice yield-related traits. Therefore, this study provides a repertoire of genes useful for identification of genes involved in maize ear trait heterosis and information for a better understanding of the molecular basis underlying heterosis in maize.

Genomic Correlates of Relationship QTL Involved in Fore- versus Hind Limb Divergence in Mice

Divergence of serially homologous elements of organisms is a common evolutionary pattern contributing to increased phenotypic complexity. Here, we study the genomic intervals affecting the variational independence of fore- and hind limb traits within an experimental mouse population. We use an advanced intercross of inbred mouse strains to map the loci associated with the degree of autonomy between fore- and hind limb long bone lengths (loci affecting the relationship between traits, relationship quantitative trait loci [rQTL]). These loci have been proposed to interact locally with the products of pleiotropic genes, thereby freeing the local trait from the variational constraint due to pleiotropic mutations. Using the known polymorphisms (single nucleotide polymorphisms [SNPs]) between the parental strains, we characterized and compared the genomic regions in which the rQTL, as well as their interaction partners (intQTL), reside. We find that these two classes of QTL intervals harbor different kinds of molecular variation. SNPs in rQTL intervals more frequently reside in limb-specific cis-regulatory regions than SNPs in intQTL intervals. The intQTL loci modified by the rQTL, in contrast, show the signature of protein-coding variation. This result is consistent with the widely accepted view that protein-coding mutations have broader pleiotropic effects than cis-regulatory polymorphisms. For both types of QTL intervals, the underlying candidate genes are enriched for genes involved in protein binding. This finding suggests that rQTL effects are caused by local interactions among the products of the causal genes harbored in rQTL and intQTL intervals. This is the first study to systematically document the population-level molecular variation underlying the evolution of character individuation.