Cloning and functional mechanism of the dwarf gene gba affecting stem elongation and cellulose biosynthesis in jute (Corchorus olitorius).

Plant height (PH) is an important factor affecting bast fiber yield in jute. Here, we report the mechanism of dwarfism in the 'Guangbaai' (gba) of jute. The mutant gba had shorter internode length and cell length compared to the standard cultivar 'TaiZi 4' (TZ4). Exogenous GA3 treatment indicated that gba is a GA-insensitive dwarf mutant. Quantitative trait locus (QTL) analysis of three PH-related traits via a high-density genetic linkage map according to re-seq showed that a total of 25 QTLs were identified, including 13 QTLs for PH, with phenotypic variation explained ranging from 2.42 to 74.16%. Notably, the functional mechanism of the candidate gene CoGID1a, the gibberellic acid receptor, of the major locus qPHIL5 was evaluated by transgenic analysis and virus-induced gene silencing. A dwarf phenotype-related single nucleotide mutation in CoGID1a was identified in gba, which was also unique to the dwarf phenotype of gba among 57 cultivars. Cogid1a was unable to interact with the growth-repressor DELLA even in the presence of highly accumulated gibberellins in gba. Differentially expressed genes between transcriptomes of gba and TZ4 after GA3 treatment indicated up-regulation of genes involved in gibberellin and cellulose synthesis in gba. Interestingly, it was found that up-regulation of CoMYB46, a key transcription factor in the secondary cell wall, by the highly accumulated gibberellins in gba promoted the expression of cellulose synthase genes CoCesA4 and CoCesA7. These findings provide valuable insights into fiber development affected by endogenous gibberellin accumulation in plants.

Genome-wide identification of the key Kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L).

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of Kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of Kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key Kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (Gossypium hirsutum L.). Results showed that 159 Kinesin genes, including 15 genes of the Kinesin-13 gene subfamily, were identified in upland cotton; of which 157 Kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 Kinesin genes in upland cotton, including 10 Kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven Kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 Kinesin genes were significantly associated with three fiber traits, among which a Kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one Kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the Kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited Kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

Genome-wide identification of the key kinesin genes during fiber and boll development in upland cotton (Gossypium hirsutum L.).

Kinesin is a kind of motor protein, which interacts with microtubule filaments and regulates cellulose synthesis. Cotton fiber is a natural model for studying the cellular development and cellulose synthesis. Therefore, a systematic research of kinesin gene family in cotton (Gossypium spp.) will be beneficial for both understanding the function of kinesin protein and assisting the fiber improvement. Here, we aimed to identify the key kinesin genes present in cotton by combining genome-wide expression profile data, association mapping, and public quantitative trait loci (QTLs) in upland cotton (G. hirsutum L.). Results showed that 159 kinesin genes, including 15 genes of the kinesin-13 gene subfamily, were identified in upland cotton; of which 157 kinesin genes can be traced back to the diploid ancestors, G. raimondii and G. arboreum. Using a combined analysis of public QTLs and genome-wide expression profile information, there were 29 QTLs co-localized together with 28 kinesin genes in upland cotton, including 10 kinesin-13 subfamily genes. Genome-wide expression profile data indicated that, among the 28 co-localized genes, seven kinesin genes were predominantly expressed in fibers or ovules. By association mapping analysis, 30 kinesin genes were significantly associated with three fiber traits, among which a kinesin-13 gene, Ghir_A11G028430, was found to be associated with both cotton boll length and lint weight, and one kinesin-7 gene, Ghir_D04G017880 (Gh_Kinesin7), was significantly associated with fiber strength. In addition, two missense mutations were identified in the motor domain of the Gh_Kinesin7 protein. Overall, the kinesin gene family seemingly plays an important role in cotton fiber and boll development. The exploited kinesin genes will be beneficial for the genetic improvement of fiber quality and yield.

Multiseriate cortical sclerenchyma enhance root penetration in compacted soils.

Mechanical impedance limits soil exploration and resource capture by plant roots. We examine the role of root anatomy in regulating plant adaptation to mechanical impedance and identify a root anatomical phene in maize (Zea mays) and wheat (Triticum aestivum) associated with penetration of hard soil: Multiseriate cortical sclerenchyma (MCS). We characterize this trait and evaluate the utility of MCS for root penetration in compacted soils. Roots with MCS had a greater cell wall-to-lumen ratio and a distinct UV emission spectrum in outer cortical cells. Genome-wide association mapping revealed that MCS is heritable and genetically controlled. We identified a candidate gene associated with MCS. Across all root classes and nodal positions, maize genotypes with MCS had 13% greater root lignin concentration compared to genotypes without MCS. Genotypes without MCS formed MCS upon exogenous ethylene exposure. Genotypes with MCS had greater lignin concentration and bending strength at the root tip. In controlled environments, MCS in maize and wheat was associated improved root tensile strength and increased penetration ability in compacted soils. Maize genotypes with MCS had root systems with 22% greater depth and 49% greater shoot biomass in compacted soils in the field compared to lines without MCS. Of the lines we assessed, MCS was present in 30 to 50% of modern maize, wheat, and barley cultivars but was absent in teosinte and wild and landrace accessions of wheat and barley. MCS merits investigation as a trait for improving plant performance in maize, wheat, and other grasses under edaphic stress.

Xylem systems genetics analysis reveals a key regulator of lignin biosynthesis in Populus deltoides.

Despite the growing resources and tools for high-throughput characterization and analysis of genomic information, the discovery of the genetic elements that regulate complex traits remains a challenge. Systems genetics is an emerging field that aims to understand the flow of biological information that underlies complex traits from genotype to phenotype. In this study, we used a systems genetics approach to identify and evaluate regulators of the lignin biosynthesis pathway in Populus deltoides by combining genome, transcriptome, and phenotype data from a population of 268 unrelated individuals of P. deltoides The discovery of lignin regulators began with the quantitative genetic analysis of the xylem transcriptome and resulted in the detection of 6706 and 4628 significant local- and distant-eQTL associations, respectively. Among the locally regulated genes, we identified the R2R3-MYB transcription factor MYB125 (Potri.003G114100) as a putative trans-regulator of the majority of genes in the lignin biosynthesis pathway. The expression of MYB125 in a diverse population positively correlated with lignin content. Furthermore, overexpression of MYB125 in transgenic poplar resulted in increased lignin content, as well as altered expression of genes in the lignin biosynthesis pathway. Altogether, our findings indicate that MYB125 is involved in the control of a transcriptional coexpression network of lignin biosynthesis genes during secondary cell wall formation in P. deltoides.

Weak genetic signal for phenotypic integration implicates developmental processes as major regulators of trait covariation.

Phenotypic integration is an important metric that describes the degree of covariation among traits in a population, and is hypothesized to arise due to selection for shared functional processes. Our ability to identify the genetic and/or developmental underpinnings of integration is marred by temporally overlapping cell-, tissue- and structure-level processes that serve to continually 'overwrite' the structure of covariation among traits through ontogeny. Here, we examine whether traits that are integrated at the phenotypic level also exhibit a shared genetic basis (e.g. pleiotropy). We micro-CT scanned two hard tissue traits, and two soft tissue traits (mandible, pectoral girdle, atrium and ventricle, respectively) from an F5 hybrid population of Lake Malawi cichlids, and used geometric morphometrics to extract 3D shape information from each trait. Given the large degree of asymmetric variation that may reflect developmental instability, we separated symmetric from asymmetric components of shape variation. We then performed quantitative trait loci (QTL) analysis to determine the degree of genetic overlap between shapes. While we found ubiquitous associations among traits at the phenotypic level, except for a handful of notable exceptions, our QTL analysis revealed few overlapping genetic regions. Taken together, this indicates developmental interactions can play a large role in determining the degree of phenotypic integration among traits, and likely obfuscate the genotype to phenotype map, limiting our ability to gain a comprehensive picture of the genetic contributors responsible for phenotypic divergence.

Integration of genome-wide association study and expression quantitative trait loci data identifies AIM2 as a risk gene of periodontitis.

AIM: To identify risk variants associated with gene expression in peripheral blood and to identify genes whose expression change may contribute to the susceptibility to periodontitis. MATERIAL AND METHODS: We systematically integrated the genetic associations from a recent large-scale periodontitis GWAS and blood expression quantitative trait loci (eQTL) data using Sherlock, a Bayesian statistical framework. We then validated the potential causal genes in independent gene expression data sets. Gene co-expression analysis was used to explore the functional relationship for the identified causal genes. RESULTS: Sherlock analysis identified 10 genes (rs7403881 for MT1L, rs12459542 for SIGLEC5, rs12459542 for SIGLEC14, rs6680386 for S100A12, rs10489524 for TRIM33, rs11962642 for HIST1H3E, rs2814770 for AIM2, rs7593959 for FASTKD2, rs10416904 for PKN1, and rs10508204 for WDR37) whose expression may influence periodontitis. Among these genes, AIM2 was consistent significantly upregulated in periodontium of periodontitis patients across four data sets. The cis-eQTL (rs2814770, ~16 kb upstream of AIM2) showed significant association with AIM2 (p = 6.63 × 10-6 ) and suggestive association with periodontitis (p = 7.52 × 10-4 ). We also validated the significant association between rs2814770 and AIM2 expression in independent expression data set. Pathway analysis revealed that genes co-expressed with AIM2 were significantly enriched in immune- and inflammation-related pathways. CONCLUSION: Our findings implicate that AIM2 is a susceptibility gene, expression of which in gingiva may influence periodontitis risk. Further functional investigation of AIM2 may provide new insight for periodontitis pathogenesis.

Genetic architecture underlying the lignin biosynthesis pathway involves noncoding RNAs and transcription factors for growth and wood properties in Populus.

Lignin provides structural support in perennial woody plants and is a complex phenolic polymer derived from phenylpropanoid pathway. Lignin biosynthesis is regulated by coordinated networks involving transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). However, the genetic networks underlying the lignin biosynthesis pathway for tree growth and wood properties remain unknown. Here, we used association genetics (additive, dominant and epistasis) and expression quantitative trait nucleotide (eQTN) mapping to decipher the genetic networks for tree growth and wood properties in 435 unrelated individuals of Populus tomentosa. We detected 124 significant associations (P ≤ 6.89E-05) for 10 growth and wood property traits using 30 265 single nucleotide polymorphisms from 203 lignin biosynthetic genes, 81 TF genes, 36 miRNA genes and 71 lncRNA loci, implying their common roles in wood formation. Epistasis analysis uncovered 745 significant pairwise interactions, which helped to construct proposed genetic networks of lignin biosynthesis pathway and found that these regulators might affect phenotypes by linking two lignin biosynthetic genes. eQTNs were used to interpret how causal genes contributed to phenotypes. Lastly, we investigated the possible functions of the genes encoding 4-coumarate: CoA ligase and cinnamate-4-hydroxylase in wood traits using epistasis, eQTN mapping and enzymatic activity assays. Our study provides new insights into the lignin biosynthesis pathway in poplar and enables the novel genetic factors as biomarkers for facilitating genetic improvement of trees.

In silico pathway analysis and tissue specific cis-eQTL for colorectal cancer GWAS risk variants.

BACKGROUND: Genome-wide association studies have identified 55 genetic variants associated with colorectal cancer risk to date. However, potential causal genes and pathways regulated by these risk variants remain to be characterized. Therefore, we performed gene ontology enrichment and pathway analyses to determine if there was an enrichment of genes in proximity to the colorectal cancer risk variants that could further elucidate the probable causal genes and pathways involved in colorectal cancer biology. RESULTS: For the 65 unique genes that either contained, or were immediately neighboring up- and downstream, of these variants there was a significant enrichment for the KEGG pathway, Pathways in Cancer (p-value = 2.67 × 10-5) and an enrichment for multiple biological processes (FDR < 0.05), such as cell junction organization, tissue morphogenesis, regulation of SMAD protein phosphorylation, and odontogenesis identified through Gene Ontology analysis. To identify potential causal genes, we conducted a cis-expression quantitative trait loci (cis-eQTL) analysis using gene expression and genotype data from the Genotype-Tissue Expression (GTEx) Project portal in normal sigmoid (n = 124) and transverse (n = 169) colon tissue. In addition, we also did a cis-eQTL analysis on colorectal tumor tissue (n = 147) from The Cancer Genome Atlas (TCGA). We identified two risk alleles that were significant cis-eQTLs for FADS2 (rs1535) and COLCA1 and 2 (rs3802842) genes in the normal transverse colon tissue and two risk alleles that were significant cis-eQTLs for the CABLES2 (rs2427308) and LIPG (rs7229639) genes in the normal sigmoid colon tissue, but not tumor tissue. CONCLUSIONS: Our data reaffirm the potential to identify an enrichment for biological processes and candidate causal genes based on expression profiles correlated with genetic risk alleles of colorectal cancer, however, the identification of these significant cis-eQTLs is context and tissue specific.

Genetic complexity of miscanthus cell wall composition and biomass quality for biofuels.

BACKGROUND: Miscanthus sinensis is a high yielding perennial grass species with great potential as a bioenergy feedstock. One of the challenges that currently impedes commercial cellulosic biofuel production is the technical difficulty to efficiently convert lignocellulosic biomass into biofuel. The development of feedstocks with better biomass quality will improve conversion efficiency and the sustainability of the value-chain. Progress in the genetic improvement of biomass quality may be substantially expedited by the development of genetic markers associated to quality traits, which can be used in a marker-assisted selection program. RESULTS: To this end, a mapping population was developed by crossing two parents of contrasting cell wall composition. The performance of 182 F1 offspring individuals along with the parents was evaluated in a field trial with a randomized block design with three replicates. Plants were phenotyped for cell wall composition and conversion efficiency characters in the second and third growth season after establishment. A new SNP-based genetic map for M. sinensis was built using a genotyping-by-sequencing (GBS) approach, which resulted in 464 short-sequence uniparental markers that formed 16 linkage groups in the male map and 17 linkage groups in the female map. A total of 86 QTLs for a variety of biomass quality characteristics were identified, 20 of which were detected in both growth seasons. Twenty QTLs were directly associated to different conversion efficiency characters. Marker sequences were aligned to the sorghum reference genome to facilitate cross-species comparisons. Analyses revealed that for some traits previously identified QTLs in sorghum occurred in homologous regions on the same chromosome. CONCLUSION: In this work we report for the first time the genetic mapping of cell wall composition and bioconversion traits in the bioenergy crop miscanthus. These results are a first step towards the development of marker-assisted selection programs in miscanthus to improve biomass quality and facilitate its use as feedstock for biofuel production.

Hedgehog signaling mediates adaptive variation in a dynamic functional system in the cichlid feeding apparatus.

Adaptive variation in the craniofacial skeleton is a key component of resource specialization and habitat divergence in vertebrates, but the proximate genetic mechanisms that underlie complex patterns of craniofacial variation are largely unknown. Here we demonstrate that the Hedgehog (Hh) signaling pathway mediates widespread variation across a complex functional system that affects the kinematics of lower jaw depression--the opercular four-bar linkage apparatus--among Lake Malawi cichlids. By using a combined quantitative trait locus mapping and population genetics approach, we show that allelic variation in the Hh receptor, ptch1, affects the development of distinct bony elements in the head that represent two of three movable links in this functional system. The evolutionarily derived allele is found in species that feed from the water column, and is associated with shifts in anatomy that translate to a four-bar system capable of faster jaw rotation. Alternatively, the ancestral allele is found in species that feed on attached algae, and is associated with the development of a four-bar system that predicts slower jaw movement. Experimental manipulation of the Hh pathway during cichlid development recapitulates functionally salient natural variation in craniofacial geometry. In all, these results significantly extend our understanding of the mechanisms that fine-tune the craniofacial skeletal complex during adaptation to new foraging niches.

Identification of quantitative trait loci controlling fibre length and lignin content in Arabidopsis thaliana stems.

Fibre properties and the biochemical composition of cell walls are important traits in many applications. For example, the lengths of fibres define the strength and quality of paper, and lignin content is a critical parameter for the use of biomass in biofuel production. Identifying genes controlling these traits is comparatively difficult in woody species, because of long generation times and limited amenability to high-resolution genetic mapping. To address this problem, this study mapped quantitative trait loci (QTLs) defining fibre length and lignin content in the Arabidopsis recombinant inbred line population Col-4 × Ler-0. Adapting high-throughput phenotyping techniques for both traits for measurements in Arabidopsis inflorescence stems identified significant QTLs for fibre length on chromosomes 2 and 5, as well as one significant QTL affecting lignin content on chromosome 2. For fibre length, total variation within the population was 208% higher than between parental lines and the identified QTLs explained 50.58% of the observed variation. For lignin content, the values were 261 and 26.51%, respectively. Bioinformatics analysis of the associated intervals identified a number of candidate genes for fibre length and lignin content. This study demonstrates that molecular mapping of QTLs pertaining to wood and fibre properties is possible in Arabidopsis, which substantially broadens the use of Arabidopsis as a model species for the functional characterization of plant genes.

Genetic mapping of craniofacial traits in the Mexican tetra reveals loci associated with bite differences between cave and surface fish.

BACKGROUND: The Mexican tetra, Astyanax mexicanus, includes interfertile surface-dwelling and cave-dwelling morphs, enabling powerful studies aimed at uncovering genes involved in the evolution of cave-associated traits. Compared to surface fish, cavefish harbor several extreme traits within their skull, such as a protruding lower jaw, a wider gape, and an increase in tooth number. These features are highly variable between individual cavefish and even across different cavefish populations. RESULTS: To investigate these traits, we created a novel feeding behavior assay wherein bite impressions could be obtained. We determined that fish with an underbite leave larger bite impressions with an increase in the number of tooth marks. Capitalizing on the ability to produce hybrids from surface and cavefish crosses, we investigated genes underlying these segregating orofacial traits by performing Quantitative Trait Loci (QTL) analysis with F2 hybrids. We discovered significant QTL for bite (underbite vs. overbite) that mapped to a single region of the Astyanax genome. Within this genomic region, multiple genes exhibit coding region mutations, some with known roles in bone development. Further, we determined that there is evidence that this genomic region is under natural selection. CONCLUSIONS: This work highlights cavefish as a valuable genetic model for orofacial patterning and will provide insight into the genetic regulators of jaw and tooth development.

Genome-Wide Analysis of Dental Caries Variability Reveals Genotype-by-Environment Interactions.

Genotype-by-environment interactions (GEI) may influence dental caries, although their effects are difficult to detect. Variance quantitative trait loci (vQTL) may serve as an indicator of underlying GEI effects. The aim of this study was to investigate GEI effects on dental caries by prioritizing variants from genome-wide vQTL analysis. First, we identified vQTLs from ~4.3 M genome-wide variants in three cohorts of white children aged 3-5 (n = 396, n = 328, n = 773) using Levene's test. A total of 39 independent vQTLs with p < 1 × 10-6 were identified, some of which were located in or near genes with plausible biological roles in dental caries (IGFBP7, SLC5A8, and SHH involved in tooth development and enamel mineralization). Next, we used linear regression to test GEI effects on dental caries with the 39 prioritized variants and self-reported environmental factors (demographic, socioeconomic, behavioral, and dietary factors) in the three cohorts separately. We identified eight significant GEIs indicating that children with vQTL risk genotypes had higher caries experience if they had less educated parents, lower household/parental income, brushed their teeth less frequently, consumed sugar-sweetened beverages more frequently, were not breastfed, and were female. We reported the first genome-wide vQTL analysis of dental caries in children nominating several novel genes and GEI for further investigations.

Testing the accuracy of 3D automatic landmarking via genome-wide association studies.

Various advances in 3D automatic phenotyping and landmark-based geometric morphometric methods have been made. While it is generally accepted that automatic landmarking compromises the capture of the biological variation, no studies have directly tested the actual impact of such landmarking approaches in analyses requiring a large number of specimens and for which the precision of phenotyping is crucial to extract an actual biological signal adequately. Here, we use a recently developed 3D atlas-based automatic landmarking method to test its accuracy in detecting QTLs associated with craniofacial development of the house mouse skull and lower jaws for a large number of specimens (circa 700) that were previously phenotyped via a semiautomatic landmarking method complemented with manual adjustment. We compare both landmarking methods with univariate and multivariate mapping of the skull and the lower jaws. We find that most significant SNPs and QTLs are not recovered based on the data derived from the automatic landmarking method. Our results thus confirm the notion that information is lost in the automated landmarking procedure although somewhat dependent on the analyzed structure. The automatic method seems to capture certain types of structures slightly better, such as lower jaws whose shape is almost entirely summarized by its outline and could be assimilated as a 2D flat object. By contrast, the more apparent 3D features exhibited by a structure such as the skull are not adequately captured by the automatic method. We conclude that using 3D atlas-based automatic landmarking methods requires careful consideration of the experimental question.

Fine mapping of dental fluorosis quantitative trait loci in mice.

Genetic factors underlie the susceptibility and the resistance to dental fluorosis (DF). The A/J (DF susceptible) and 129P3/J (DF resistant) mouse strains have previously been used to detect quantitative trait loci (QTLs) associated with DF on chromosome (Chr) 2 and Chr 11. In the present study, increased marker density genotyping followed by interval mapping was performed to narrow the QTL intervals and improve the logarithm of the odds (to the base 10) (LOD) scores. Narrower intervals were obtained on Chr 2 where LOD ≥ 6.0 (57-84 cM or ≈ 51 Mb), LOD ≥ 7.0 (62-79 cM or ≈ 32 Mb), and LOD ≥ 8.0 (65-74 cM or ≈ 17 Mb); and on Chr 11 where LOD ≥ 6.0 (18-51 cM or ≈ 53 Mb), LOD ≥ 7.0 (28-48 cM or ≈ 34 Mb), and LOD ≥ 8.0 (31-45 cM or ≈ 22 Mb). Haplotype analysis between A/J and 129P3/J mice further reduced the QTL intervals. Accn1 was selected as a candidate gene based upon its location near the peak LOD score on Chr 11 and distant homology with the Caenorhabditis elegans fluoride-resistance gene, flr1. The severity of DF between Accn1(-/-) and wild-type mice was not significantly different. Hence, the loss of ACCN1 function does not modify DF severity in mice. Narrowing the DF QTL intervals will facilitate additional candidate gene selections and interrogation.

The cichlid oral and pharyngeal jaws are evolutionarily and genetically coupled.

Evolutionary constraints may significantly bias phenotypic change, while "breaking" from such constraints can lead to expanded ecological opportunity. Ray-finned fishes have broken functional constraints by developing two jaws (oral-pharyngeal), decoupling prey capture (oral jaw) from processing (pharyngeal jaw). It is hypothesized that the oral and pharyngeal jaws represent independent evolutionary modules and this facilitated diversification in feeding architectures. Here we test this hypothesis in African cichlids. Contrary to our expectation, we find integration between jaws at multiple evolutionary levels. Next, we document integration at the genetic level, and identify a candidate gene, smad7, within a pleiotropic locus for oral and pharyngeal jaw shape that exhibits correlated expression between the two tissues. Collectively, our data show that African cichlid evolutionary success has occurred within the context of a coupled jaw system, an attribute that may be driving adaptive evolution in this iconic group by facilitating rapid shifts between foraging habitats, providing an advantage in a stochastic environment such as the East African Rift-Valley.

Genome-wide association studies provide insights into the genetic determination of fruit traits of pear.

Pear is a major fruit tree crop distributed worldwide, yet its breeding is a very time-consuming process. To facilitate molecular breeding and gene identification, here we have performed genome-wide association studies (GWAS) on eleven fruit traits. We identify 37 loci associated with eight fruit quality traits and five loci associated with three fruit phenological traits. Scans for selective sweeps indicate that traits including fruit stone cell content, organic acid and sugar contents might have been under continuous selection during breeding improvement. One candidate gene, PbrSTONE, identified in GWAS, has been functionally verified to be involved in the regulation of stone cell formation, one of the most important fruit quality traits in pear. Our study provides insights into the complex fruit related biology and identifies genes controlling important traits in pear through GWAS, which extends the genetic resources and basis for facilitating molecular breeding in perennial trees.

Translation of mouse model to human gives insights into periodontitis etiology.

To suggest candidate genes involved in periodontitis, we combined gene expression data of periodontal biopsies from Collaborative Cross (CC) mouse lines, with previous reported quantitative trait loci (QTL) in mouse and with human genome-wide association studies (GWAS) associated with periodontitis. Periodontal samples from two susceptible, two resistant and two lines that showed bone formation after periodontal infection were collected during infection and naïve status. Differential expressed genes (DEGs) were analyzed in a case-control and case-only design. After infection, eleven protein-coding genes were significantly stronger expressed in resistant CC lines compared to susceptible ones. Of these, the most upregulated genes were MMP20 (P = 0.001), RSPO4 (P = 0.032), CALB1 (P = 1.06×10-4), and AMTN (P = 0.05). In addition, human orthologous of candidate genes were tested for their association in a case-controls samples of aggressive (AgP) and chronic (CP) periodontitis (5,095 cases, 9,908 controls). In this analysis, variants at two loci, TTLL11/PTGS1 (rs9695213, P = 5.77×10-5) and RNASE2 (rs2771342, P = 2.84×10-5) suggested association with both AgP and CP. In the association analysis with AgP only, the most significant associations were located at the HLA loci HLA-DQH1 (rs9271850, P = 2.52×10-14) and HLA-DPA1 (rs17214512, P = 5.14×10-5). This study demonstrates the utility of the CC RIL populations as a suitable model to investigate the mechanism of periodontal disease.

Detection of dental fluorosis-associated quantitative trait Loci on mouse chromosomes 2 and 11.

Systemic exposure to greater than optimal fluoride (F) can lead to dental fluorosis (DF). Parental A/J (DF-susceptible) and 129P3/J (DF-resistant) inbred mice were used for histological studies and to generate F2 progeny. Mice were treated with 0 or 50 ppm F in their drinking water for 60 days. A clinical criterion (modified Thylstrup and Fejerskov categorical scale) was used to assess the severity of DF for each individual F2 animal. Parental strains were subjected to histological examination of maturing enamel. F treatment resulted in accumulation of amelogenins in the maturing enamel of A/J mice. Quantitative trait loci (QTL) detection was performed using phenotypic extreme F2 animals genotyped for 354 single nucleotide polymorphism-based markers distributed throughout the mouse genome followed by chi(2) analysis. Significant evidence of association was observed on chromosomes 2 and 11 for a series of consecutive markers (p < 0.0001). Further analyses were performed to examine whether the phenotypic effects were found in both male and female F2 mice or whether there was evidence for gender-specific effects. Analyses performed using the markers on chromosomes 2 and 11 which were significant in the mixed-gender mice were also significant when analyses were limited to only the male or female mice. The QTL detected on chromosomes 2 and 11 which influence the variation in response to fluorosis have their effect in mice of both genders. Finally, the QTL in both chromosomes appear to have an additive effect.