Search

Background: Plant long non-coding RNAs (lncRNAs) have important regulatory roles in responses to various biotic and abiotic stresses, including light quality. However, no lncRNAs have been specifically linked to the Shade Avoidance Response (SAS). Results: To better understand the involvement of lncRNAs in shade avoidance, we examined RNA-seq libraries for lncRNAs with the potential to function in the neighbor proximity phenomenon in Arabidopsis thaliana (A. thaliana). Using transcriptomes generated from seedlings exposed to high and low red/far-red (R/FR) light conditions, we identified 13 lncRNA genes differentially expressed in cotyledons and 138 in hypocotyls. To infer possible functions for these lncRNAs, we used a 'guilt-by-association' approach to identify genes co-expressed with lncRNAs in a weighted gene co-expression network. Of 34 co-expression modules, 10 showed biological functions related to differential growth. We identified three potential lncRNAs co-regulated with genes related to SAS. T-DNA insertions in two of these lncRNAs were correlated with morphological differences in seedling responses to increased FR light, supporting our strategy for computational identification of lncRNAs involved in SAS. Conclusions: Using a computational approach, we identified multiple lncRNAs in Arabidopsis involved in SAS. T-DNA insertions caused altered phenotypes under low R/FR light, suggesting functional roles in shade avoidance. Further experiments are needed to determine the specific mechanisms of these lncRNAs in SAS. [ABSTRACT FROM AUTHOR]

As a Brassica crop, Brassica napus typically has single flowers that contain four petals. The double-flower phenotype of rapeseed has been a desirable trait in China because of its potential commercial value in ornamental tourism. However, few double-flowered germplasms have been documented in B. napus, and knowledge of the underlying genes is limited. Here, B. napus D376 was characterized as a double-flowered strain that presented an average of 10.92 ± 1.40 petals and other normal floral organs. F1, F2 and BC1 populations were constructed by crossing D376 with a single-flowered line reciprocally. Genetic analysis revealed that the double-flower trait was a recessive trait controlled by multiple genes. To identify the key genes controlling the double-flower trait, bulk segregant analysis sequencing (BSA-seq) and RNA-seq analyses were conducted on F2 individual bulks with opposite extreme phenotypes. Through BSA-seq, one candidate interval was mapped at the region of chromosome C05: 14.56–16.17 Mb. GO and KEGG enrichment analyses revealed that the DEGs were significantly enriched in carbohydrate metabolic processes, notably starch and sucrose metabolism. Interestingly, five and thirty-six DEGs associated with floral development were significantly up- and down-regulated, respectively, in the double-flowered plants. A combined analysis of BSA-seq and RNA-seq data revealed that five genes were candidates associated with the double flower trait, and BnaC05.ERS2 was the most promising gene. These findings provide novel insights into the breeding of double-flowered varieties and lay a theoretical foundation for unveiling the molecular mechanisms of floral development in B. napus. [ABSTRACT FROM AUTHOR]

Background: The Lin−Sca1+c-Kit+ (LSK) fraction of the bone marrow (BM) comprises multipotent hematopoietic stem cells (HSCs), which are vital to tissue homeostasis and vascular repair. While diabetes affects HSC homeostasis overall, the molecular signature of mRNA and miRNA transcriptomic under the conditions of long-standing type 2 diabetes (T2D;>6 months) remains unexplored. Methods: In this study, we assessed the transcriptomic signature of HSCs in db/db mice, a well-known and widely used model for T2D. LSK cells of db/db mice enriched using a cell sorter were subjected to paired-end mRNA and single-end miRNA seq library and sequenced on Illumina NovaSeq 6000. The mRNA sequence reads were mapped using STAR (Spliced Transcripts Alignment to a Reference), and the miRNA sequence reads were mapped to the designated reference genome using the Qiagen GeneGlobe RNA-seq Analysis Portal with default parameters for miRNA. Results: We uncovered 2076 out of 13,708 mRNAs and 35 out of 191 miRNAs that were expressed significantly in db/db animals; strikingly, previously unreported miRNAs (miR-3968 and miR-1971) were found to be downregulated in db/db mice. Furthermore, we observed a molecular shift in the transcriptome of HSCs of diabetes with an increase in pro-inflammatory cytokines (Il4, Tlr4, and Tnf11α) and a decrease in anti-inflammatory cytokine IL10. Pathway mapping demonstrated inflammation mediated by chemokine, cytokine, and angiogenesis as one of the top pathways with a significantly higher number of transcripts in db/db mice. These molecular changes were reflected in an overt defect in LSK mobility in the bone marrow. miRNA downstream target analysis unveils several mRNAs targeting leukocyte migration, microglia activation, phagosome formation, and macrophage activation signaling as their primary pathways, suggesting a shift to an inflammatory phenotype. Conclusion: Our findings highlight that chronic diabetes adversely alters HSCs' homeostasis at the transcriptional level, thus potentially contributing to the inflammatory phenotype of HSCs under long-term diabetes. We also believe that identifying HSCs-based biomarkers in miRNAs or mRNAs could serve as diagnostic markers and potential therapeutic targets for diabetes and associated vascular complications. [ABSTRACT FROM AUTHOR]

Background: In the face of contemporary climatic vulnerabilities and escalating global temperatures, the prevalence of maydis leaf blight (MLB) poses a potential threat to maize production. This study endeavours to discern marker-trait associations and elucidate the candidate genes that underlie resistance to MLB in maize by employing a diverse panel comprising 336 lines. The panel was screening for MLB across four environments, employing standard artificial inoculation techniques. Genome-wide association studies (GWAS) and haplotype analysis were conducted utilizing a total of 128,490 SNPs obtained from genotyping-by-sequencing (GBS). Results: GWAS identified 26 highly significant SNPs associated with MLB resistance, among the markers examined. Seven of these SNPs, reported in novel chromosomal bins (9.06, 5.01, 9.01, 7.04, 4.06, 1.04, and 6.05) were associated with genes: bzip23, NAGS1, CDPK7, aspartic proteinase NEP-2, VQ4, and Wun1, which were characterized for their roles in diminishing fungal activity, fortifying defence mechanisms against necrotrophic pathogens, modulating phyto-hormone signalling, and orchestrating oxidative burst responses. Gene mining approach identified 22 potential candidate genes associated with SNPs due to their functional relevance to resistance against necrotrophic pathogens. Notably, bin 8.06, which hosts five SNPs, showed a connection to defense-regulating genes against MLB, indicating the potential formation of a functional gene cluster that triggers a cascade of reactions against MLB. In silico studies revealed gene expression levels exceeding ten fragments per kilobase million (FPKM) for most genes and demonstrated coexpression among all candidate genes in the coexpression network. Haplotype regression analysis revealed the association of 13 common significant haplotypes at Bonferroni ≤ 0.05. The phenotypic variance explained by these significant haplotypes ranged from low to moderate, suggesting a breeding strategy that combines multiple resistance alleles to enhance resistance to MLB. Additionally, one particular haplotype block (Hap_8.3) was found to consist of two SNPs (S8_152715134, S8_152460815) identified in GWAS with 9.45% variation explained (PVE). Conclusion: The identified SNPs/ haplotypes associated with the trait of interest contribute to the enrichment of allelic diversity and hold direct applicability in Genomics Assisted Breeding for enhancing MLB resistance in maize. [ABSTRACT FROM AUTHOR]

Background: Acinetobacter baumannii (A. baumannii) is a common opportunistic pathogen in hospitals that causes nosocomial infection. In order to understand the phenotypic and genotypic characteristics of A. baumannii isolates, we sequenced and analyzed 62 A. baumannii isolates from a hospital in Gansu province. Results: Non-repeated 62 A. baumannii isolates were collected from August 2015 to November 2021. Most isolates (56/62) were resistant to multiple drugs. All the 62 A. baumannii isolates were resistant to aztreonam and contained blaADC-25 gene which exists only on chromosome contigs. The 62 isolates in this study were not clustered in a single clade, but were dispersed among multiple clades in the common genome. Seven sequence types were identified by Multilocus sequence type (MLST) analysis and most isolates (52/62) belonged to ST2. The plasmids were grouped into 11 clusters by MOB-suite. Conclusions: This study furthers the understanding of A. baumannii antimicrobial-resistant genotypes, and may aid in prevention and control nosocomial infection caused by drug-resistant A. baumannii. [ABSTRACT FROM AUTHOR]

Background: Heritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium. Results: Here, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively). Conclusions: Heritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used. [ABSTRACT FROM AUTHOR]

Background: The number of spikelets per spike is a key trait that affects the yield of bread wheat (Triticum aestivum L.). Identification of the QTL for spikelets per spike and its genetic effects that could be used in molecular assistant breeding in the future. Results: In this study, four recombinant inbred line (RIL) populations were generated and used, having YuPi branching wheat (YP), with Supernumerary Spikelets (SS) phenotype, as a common parent. QTL (QSS.sicau-2 A and QSS.sicau-2D) related to SS trait were mapped on chromosomes 2 A and 2D through bulked segregant exome sequencing (BSE-Seq). Fourteen molecular markers were further developed within the localization interval, and QSS.sicau-2 A was narrowed to 3.0 cM covering 7.6 Mb physical region of the reference genome, explaining 13.7 − 15.9% the phenotypic variance. Similarly, the QSS.sicau-2D was narrowed to 1.8 cM covering 2.4 Mb physical region of the reference genome, and it explained 27.4 − 32.9% the phenotypic variance. These two QTL were validated in three different genetic backgrounds using the linked markers. QSS.sicau-2 A was identified as WFZP-A, and QSS.sicau-2D was identified a novel locus, different to the previously identified WFZP-D. Based on the gene expression patterns, gene annotation and sequence analysis, TraesCS2D03G0260700 was predicted to be a potential candidate gene for QSS.sicau-2D. Conclusion: Two significant QTL for SS, namely QSS.sicau-2 A and QSS.sicau-2D were identified in multiple environments were identified and their effect in diverse genetic populations was assessed. QSS.sicau-2D is a novel QTL associated with the SS trait, with TraesCS2D03G0260700 predicted as its candidate gene. [ABSTRACT FROM AUTHOR]

Background: Breeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals. Methods: The present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats. Results: The analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats. Conclusion: We have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits. [ABSTRACT FROM AUTHOR]

Background: The breeding of layers emphasizes the continual selection of egg-related traits, such as egg production, egg quality and eggshell, which enhance their productivity and meet the demand of market. As the breeding process continued, the genomic homozygosity of layers gradually increased, resulting in the emergence of runs of homozygosity (ROH). Therefore, ROH analysis can be used in conjunction with other methods to detect selection signatures and identify candidate genes associated with various important traits in layer breeding. Results: In this study, we generated whole-genome sequencing data from 686 hens in a Rhode Island Red population that had undergone fifteen consecutive generations of intensive artificial selection. We performed a genome-wide ROH analysis and utilized multiple methods to detect signatures of selection. A total of 141,720 ROH segments were discovered in whole population, and most of them (97.35%) were less than 3 Mb in length. Twenty-three ROH islands were identified, and they overlapped with some regions bearing selection signatures, which were detected by the De-correlated composite of multiple signals methods (DCMS). Sixty genes were discovered and functional annotation analysis revealed the possible roles of them in growth, development, immunity and signaling in layers. Additionally, two-tailed analyses including DCMS and ROH for 44 phenotypes of layers were conducted to find out the genomic differences between subgroups of top and bottom 10% phenotype of individuals. Combining the results of GWAS, we observed that regions significantly associated with traits also exhibited selection signatures between the high and low subgroups. We identified a region significantly associated with egg weight near the 25 Mb region of GGA 1, which exhibited selection signatures and has higher genomic homozygosity in the low egg weight subpopulation. This suggests that the region may be play a role in the decline in egg weight. Conclusions: In summary, through the combined analysis of ROH, selection signatures, and GWAS, we identified several genomic regions that associated with the production traits of layers, providing reference for the study of layer genome. [ABSTRACT FROM AUTHOR]

We investigated whole blood and hepatic mRNA expressions of immune genes and rumen microbiome of crossbred beef steers with divergent residual feed intake phenotype to identify relevant biological processes underpinning feed efficiency in beef cattle. Low-RFI beef steers (n = 20; RFI = − 1.83 kg/d) and high-RFI beef steers (n = 20; RFI = + 2.12 kg/d) were identified from a group of 108 growing crossbred beef steers (average BW = 282 ± 30.4 kg) fed a high-forage total mixed ration after a 70-d performance testing period. At the end of the 70-d testing period, liver biopsies and blood samples were collected for total RNA extraction and cDNA synthesis. Rumen fluid samples were also collected for analysis of the rumen microbial community. The mRNA expression of 84 genes related to innate and adaptive immunity was analyzed using pathway-focused PCR-based arrays. Differentially expressed genes were determined using P-value ≤ 0.05 and fold change (FC) ≥ 1.5 (in whole blood) or ≥ 2.0 (in the liver). Gene ontology analysis of the differentially expressed genes revealed that pathways related to pattern recognition receptor activity, positive regulation of phagocytosis, positive regulation of vitamin metabolic process, vascular endothelial growth factor production, positive regulation of epithelial tube formation and T-helper cell differentiation were significantly enriched (FDR < 0.05) in low-RFI steers. In the rumen, the relative abundance of PeH15, Arthrobacter, Moryella, Weissella, and Muribaculaceae was enriched in low-RFI steers, while Methanobrevibacter, Bacteroidales_BS11_gut_group, Bacteroides and Clostridium_sensu_stricto_1 were reduced. In conclusion, our study found that low-RFI beef steers exhibit increased mRNA expression of genes related to immune cell functions in whole blood and liver tissues, specifically those involved in pathogen recognition and phagocytosis regulation. Additionally, these low-RFI steers showed differences in the relative abundance of some microbial taxa which may partially account for their improved feed efficiency compared to high-RFI steers. [ABSTRACT FROM AUTHOR]

Background: Alpaca is a domestic South American camelid probably arising from the domestication of two wild camelids, the vicugna and the guanaco. Two phenotypes are described for alpaca, known as huacaya and suri. Huacaya fleece is characterized by compact, soft, and highly crimped fibers, while suri fleece is longer, straight, less crimped, and lustrous. The gene variants determining these phenotypes are still unknown, although previous studies suggested a dominant inheritance of the suri. Based on that, the aim of this study was the identification of the gene variants determining alpaca coat phenotypes through whole genome sequencing (WGS) analysis. Results: The sample used includes two test-cross alpaca families, suri × huacaya, which produced two offspring, one with the suri phenotype and one with the huacaya phenotype. The analyzed sample was expanded through the addition of WGS data from six vicugnas and six guanacos; this because we assumed the absence of the gene variants linked to the suri phenotype in these wild species. The analysis of gene variant segregation with the suri phenotype, coupled with the filtering of gene variants present in the wild species, disclosed the presence in all the suri samples of a premature termination codon (PTC) in TRPV3 (transient receptor potential cation channel subfamily V member 3), a gene known to be involved in hair growth and cycling, thermal sensation, cold tolerance and adaptation in several species. Mutations in TRPV3 were previously associated with the alteration of hair structure leading to an impaired formation of the hair canal and the hair shaft in mouse. This PTC in TRPV3, due to a G > T substitution (p.Glu475*), results in a loss of 290 amino acids from the canonical translated protein, plausibly leading to a physiological dysfunction. Conclusion: The present results suggest that the suri phenotype may arise from a TRPV3 gene variant which may explain some of the suri features such as its longer hair fibre with lower number of cuticular scales compared to huacaya. [ABSTRACT FROM AUTHOR]

Background: Understanding the molecular underpinnings of phenotypic variations is critical for enhancing poultry breeding programs. The Brazilian broiler (TT) and laying hen (CC) lines exhibit striking differences in body weight, growth potential, and muscle mass. Our work aimed to compare the global transcriptome of wing and pectoral tissues during the early development (days 2.5 to 3.5) of these chicken lines, unveiling disparities in gene expression and regulation. Results: Different and bona-fide transcriptomic profiles were identified for the compared lines. A similar number of up- and downregulated differentially expressed genes (DEGs) were identified, considering the broiler line as a reference. Upregulated DEGs displayed an enrichment of protease-encoding genes, whereas downregulated DEGs exhibited a prevalence of receptors and ligands. Gene Ontology analysis revealed that upregulated DEGs were mainly associated with hormone response, mitotic cell cycle, and different metabolic and biosynthetic processes. In contrast, downregulated DEGs were primarily linked to communication, signal transduction, cell differentiation, and nervous system development. Regulatory networks were constructed for the mitotic cell cycle and cell differentiation biological processes, as their contrasting roles may impact the development of distinct postnatal traits. Within the mitotic cell cycle network, key upregulated DEGs included CCND1 and HSP90, with central regulators being NF-κB subunits (RELA and REL) and NFATC2. The cell differentiation network comprises numerous DEGs encoding transcription factors (e.g., HOX genes), receptors, ligands, and histones, while the main regulatory hubs are CREB, AR and epigenetic modifiers. Clustering analyses highlighted PIK3CD as a central player within the differentiation network. Conclusions: Our study revealed distinct developmental transcriptomes between Brazilian broiler and layer lines. The gene expression profile of broiler embryos seems to favour increased cell proliferation and delayed differentiation, which may contribute to the subsequent enlargement of pectoral tissues during foetal and postnatal development. Our findings pave the way for future functional studies and improvement of targeted traits of economic interest in poultry. [ABSTRACT FROM AUTHOR]

Genetic pleiotropy refers to the simultaneous association of a gene with multiple phenotypes. It is widely distributed in the whole genome and can help to understand the common genetic mechanism of diseases or traits. In this study, a multivariate response best-subset selection (MRBSS) model based pleiotropic association analysis method is proposed. Different from the traditional genetic association model, the high-dimensional genotypic data are viewed as response variables while the multiple phenotypic data as predictor variables. Moreover, the response best-subset selection procedure is converted into an 0-1 integer optimization problem by introducing a separation parameter and a tuning parameter. Furthermore, the model parameters are estimated by using the curve search under the modified Bayesian information criterion. Simulation experiments show that the proposed method MRBSS remarkably reduces the computational time, obtains higher statistical power under most of the considered scenarios, and controls the type I error rate at a low level. The application studies in the datasets of maize yield traits and pig lipid traits further verifies the effectiveness. [ABSTRACT FROM AUTHOR]

Background: Single synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed. Results: Synonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo. Conclusion: In summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations. [ABSTRACT FROM AUTHOR]

Background: The prevalence of Acinetobacter baumannii in nosocomial infections and its remarkable ability to develop antimicrobial resistance have been a critical issue in hospital settings. Here, we examined the genomic features related to resistance phenotype displayed by carbapenem-resistant A. baumannii (CRAB) MTC1106 (ST2) and MTC0619 (ST25). Results: Resistome analysis of both strains revealed that MTC1106 possessed higher numbers of antimicrobial resistance genes compared to MTC0619. Some of those genetic determinants were present in accordance with the susceptibility profile of the isolates. The predicted ISAba1 region upstream of blaOXA-23 gene was related to carbapenem resistance since this IS element was well-characterized to mediate overexpression of carbapenemase genes and eventually provided capability to confer resistance. Unlike MTC0619 strain, which only carried class B and D β-lactamase genes, MTC1106 strain also possessed blaTEM-1D, a class A β-lactamase. Regarding to aminoglycosides resistance, MTC0619 contained 5 related genes in which all of them belonged to three groups of aminoglycosides modifying enzyme (AME), namely, N-acetyltransferase (AAC), O-nucleotidyltransferase (ANT), and O-phosphotransferase (APH). On the other hand, MTC1106 lacked only the AAC of which found in MTC0619, yet it also carried an armA gene encoding for 16S rRNA methyltransferase. Two macrolides resistance genes, mph(E) and msr(E), were identified next to the armA gene of MTC1106 isolate in which they encoded for macrolide 2'-phosphotransferase and ABC-type efflux pump, respectively. Besides acquired resistance genes, some chromosomal genes and SNPs associated with resistance to fluoroquinolones (i.e. gyrA and parC) and colistin (i.e. pmrCAB, eptA, and emrAB) were observed. However, gene expression analysis suggested that the genetic determinants significantly contributing to low-level colistin resistance remained unclear. In addition, similar number of efflux pumps genes were identified in both lineages with only the absence of adeC, a part of adeABC RND-type multidrug efflux pump in MTC0619 strain. Conclusions: We found that MTC1106 strain harbored more antimicrobial resistance genes and showed higher resistance to antibiotics than MTC0619 strain. Regarding genomic characterization, this study was likely the first genome comparative analysis of CARB that specifically included isolates belonging to ST2 and ST25 which were widely spread in Thailand. Taken altogether, this study suggests the importance to monitor the resistance status of circulating A. baumannii clones and identify genes that may contribute to shifting the resistance trend among isolates. [ABSTRACT FROM AUTHOR]

Background: Selective constraint, the depletion of variation due to negative selection, provides insights into the functional impact of variants and disease mechanisms. However, its characterization in mice, the most commonly used mammalian model, remains limited. This study aims to quantify mouse gene constraint using a new metric called the nonsynonymous observed expected ratio (NOER) and investigate its relationship with gene function. Results: NOER was calculated using whole-genome sequencing data from wild mouse populations (Mus musculus sp and Mus spretus). Positive correlations were observed between mouse gene constraint and the number of associated knockout phenotypes, indicating stronger constraint on pleiotropic genes. Furthermore, mouse gene constraint showed a positive correlation with the number of pathogenic variant sites in their human orthologues, supporting the relevance of mouse models in studying human disease variants. Conclusions: NOER provides a resource for assessing the fitness consequences of genetic variants in mouse genes and understanding the relationship between gene constraint and function. The study's findings highlight the importance of pleiotropy in selective constraint and support the utility of mouse models in investigating human disease variants. Further research with larger sample sizes can refine constraint estimates in mice and enable more comprehensive comparisons of constraint between mouse and human orthologues. [ABSTRACT FROM AUTHOR]

The H9c2 myoblast cell line, isolated from the left ventricular tissue of rat, is currently used in vitro as a mimetic for skeletal and cardiac muscle due to its biochemical, morphological, and electrical/hormonal signaling properties. During culture, H9c2 cells acquire a myotube phenotype, where a critical component is the inclusion of retinoic acid (RA). The results from some authors on H9c2 suggested that thousands of genes respond to RA stimuli, while others report hundreds of genes responding to RA over different cell types. In this article, using a more appropriate experimental design, we first confirm the H9c2 cardiac phenotype with and without RA and report transcriptomic and physiological changes regarding calcium handling, bioenergetics, and other biological concepts. Interestingly, of the 2360 genes showing a transcriptional change, 622 genes were statistically associated with the RA response. Of these genes, only 305 were RA-specific, and the rest also showed a culture-time component. Thus, the major expression changes (from 74 to 87%) were indeed due to culture conditions over time. Unexpectedly, only a few components of the retinol pathway in KEGG responded to RA. Our results show the role of RA in the H9c2 cultures impacting the interpretation using H9c2 as an in vitro model. [ABSTRACT FROM AUTHOR]

Background: Genome-wide association studies (GWASs) have identified many single-nucleotide polymorphisms (SNPs) associated with complex phenotypes in the European (EUR) population; however, the extent to which EUR-associated SNPs can be generalized to other populations such as East Asian (EAS) is not clear. Results: By leveraging summary statistics of 31 phenotypes in the EUR and EAS populations, we first evaluated the difference in heritability between the two populations and calculated the trans-ethnic genetic correlation. We observed the heritability estimates of some phenotypes varied substantially across populations and 53.3% of trans-ethnic genetic correlations were significantly smaller than one. Next, we examined whether EUR-associated SNPs of these phenotypes could be identified in EAS using the trans-ethnic false discovery rate method while accounting for winner's curse for SNP effect in EUR and difference of sample sizes in EAS. We found on average 54.5% of EUR-associated SNPs were also significant in EAS. Furthermore, we discovered non-significant SNPs had higher effect heterogeneity, and significant SNPs showed more consistent linkage disequilibrium and allele frequency patterns between the two populations. We also demonstrated non-significant SNPs were more likely to undergo natural selection. Conclusions: Our study revealed the extent to which EUR-associated SNPs could be significant in the EAS population and offered deep insights into the similarity and diversity of genetic architectures underlying phenotypes in distinct ancestral groups. [ABSTRACT FROM AUTHOR]

Background: In light of previous studies that profiled breed-specific traits or used genome-wide association studies to refine loci associated with characteristic morphological features in dogs, the field has gained tremendous genetic insights for known dog traits observed among breeds. Here we aim to address the question from a reserve perspective: whether there are breed-specific genotypes that may underlie currently unknown phenotypes. This study provides a complete set of breed-specific genetic signatures (BSGS). Several novel BSGS with significant protein-altering effects were highlighted and validated. Results: Using the next generation whole-genome sequencing technology coupled with unsupervised machine learning for pattern recognitions, we constructed and analyzed a high-resolution sequence map for 76 breeds of 412 dogs. Genomic structures including novel single nucleotide polymorphisms (SNPs), SNP clusters, insertions, deletions (INDELs) and short tandem repeats (STRs) were uncovered mutually exclusively among breeds. We also partially validated some novel nonsense variants by Sanger sequencing with additional dogs. Four novel nonsense BSGS were found in the Bernese Mountain Dog, Samoyed, Bull Terrier, and Basset Hound, respectively. Four INDELs resulting in either frame-shift or codon disruptions were found in the Norwich Terrier, Airedale Terrier, Chow Chow and Bernese Mountain Dog, respectively. A total of 15 genomic regions containing three types of BSGS (SNP-clusters, INDELs and STRs) were identified in the Akita, Alaskan Malamute, Chow Chow, Field Spaniel, Keeshond, Shetland Sheepdog and Sussex Spaniel, in which Keeshond and Sussex Spaniel each carried one amino-acid changing BSGS in such regions. Conclusion: Given the strong relationship between human and dog breed-specific traits, this study might be of considerable interest to researchers and all. Novel genetic signatures that can differentiate dog breeds were uncovered. Several functional genetic signatures might indicate potentially breed-specific unknown phenotypic traits or disease predispositions. These results open the door for further investigations. Importantly, the computational tools we developed can be applied to any dog breeds as well as other species. This study will stimulate new thinking, as the results of breed-specific genetic signatures may offer an overarching relevance of the animal models to human health and disease. [ABSTRACT FROM AUTHOR]

Background: Considering relatives' health history in logistic regression for case–control genome-wide association studies (CC-GWAS) may provide new information that increases accuracy and power to detect disease associated genetic variants. We conducted simulations and analyzed type 2 diabetes (T2D) data from the Framingham Heart Study (FHS) to compare two methods, liability threshold model conditional on both case–control status and family history (LT-FH) and Fam-meta, which incorporate family history into CC-GWAS. Results: In our simulation scenario of trait with modest T2D heritability (h2 = 0.28), variant minor allele frequency ranging from 1% to 50%, and 1% of phenotype variance explained by the genetic variants, Fam-meta had the highest overall power, while both methods incorporating family history were more powerful than CC-GWAS. All three methods had controlled type I error rates, while LT-FH was the most conservative with a lower-than-expected error rate. In addition, we observed a substantial increase in power of the two familial history methods compared to CC-GWAS when the prevalence of the phenotype increased with age. Furthermore, we showed that, when only the phenotypes of more distant relatives were available, Fam-meta still remained more powerful than CC-GWAS, confirming that leveraging disease history of both close and distant relatives can increase power of association analyses. Using FHS data, we confirmed the well-known association of TCF7L2 region with T2D at the genome-wide threshold of P-value < 5 × 10–8, and both familial history methods increased the significance of the region compared to CC-GWAS. We identified two loci at 5q35 (ADAMTS2) and 5q23 (PRR16), not previously reported for T2D using CC-GWAS and Fam-meta; both genes play a role in cardiovascular diseases. Additionally, CC-GWAS detected one more significant locus at 13q31 (GPC6) reported associated with T2D-related traits. Conclusions: Overall, LT-FH and Fam-meta had higher power than CC-GWAS in simulations, especially using phenotypes that were more prevalent in older age groups, and both methods detected known genetic variants with lower P-values in real data application, highlighting the benefits of including family history in genetic association studies. [ABSTRACT FROM AUTHOR]

Background: Functional genomics uses unbiased systematic genome-wide gene disruption or analyzes natural variations such as gene expression profiles of different tissues from multicellular organisms to link gene functions to particular phenotypes. Functional genomics approaches are of particular importance to identify large sets of genes that are specifically important for a particular biological process beyond known candidate genes, or when the process has not been studied with genetic methods before. Results: Here, we present a large set of genes whose disruption interferes with the function of the odoriferous defensive stink glands of the red flour beetle Tribolium castaneum. This gene set is the result of a large-scale systematic phenotypic screen using RNA interference applied in a genome-wide forward genetics manner. In this first-pass screen, 130 genes were identified, of which 69 genes could be confirmed to cause phenotypic changes in the glands upon knock-down, which vary from necrotic tissue and irregular reservoir size to irregular color or separation of the secreted gland compounds. Gene ontology analysis revealed that many of those genes are encoding enzymes (peptidases and cytochromes P450) as well as proteins involved in membrane trafficking with an enrichment in lysosome and mineral absorption pathways. The knock-down of 13 genes caused specifically a strong reduction of para-benzoquinones in the gland reservoirs, suggesting a specific function in the synthesis of these toxic compounds. Only 14 of the 69 confirmed gland genes are differentially overexpressed in stink gland tissue and thus could have been detected in a transcriptome-based analysis. However, only one out of eight genes identified by a transcriptomics approach known to cause phenotypic changes of the glands upon knock-down was recognized by this phenotypic screen, indicating the limitation of such a non-redundant first-pass screen. Conclusion: Our results indicate the importance of combining diverse and independent methodologies to identify genes necessary for the function of a certain biological tissue, as the different approaches do not deliver redundant results but rather complement each other. The presented phenotypic screen together with a transcriptomics approach are now providing a set of close to hundred genes important for odoriferous defensive stink gland physiology in beetles. [ABSTRACT FROM AUTHOR]

Background: Morinda officinalis How (MO) is a perennial herb distributed in tropical and subtropical regions, which known as one of the "Four Southern Herbal Medicines". The extent of genetic variability and the population structure of MO are presently little understood. Here, nine morphological traits, six chemical components and Single nucleotide polymorphism (SNP) markers were used in integrative research of MO germplasm variation among 88 individuals collected from ten populations across four geographical provinces of China. Results: Both phenotype and chemical composition have significant genetic variation, and there is a certain correlation between them such as root diameter and the nystose content, as well as geographical distribution. The principal component analysis (PCA) showed the leaf length, leaf width, nystose, 1F-furanosaccharide nystose, and the section color were the major contributors to diversity. The cluster analysis based on phenotypic and oligosaccharide data distinguished three significant groups, which was consistent with the result of a corresponding analysis with 228,615 SNP markers, and importantly, they all showed a significant correlation with geographical origin. However, there was little similarity between two cluster results. The Shannon's information index (I) varied from 0.17 to 0.53 with a mean of 0.37, suggesting a high level of genetic diversity in MO populations, which mainly existed among individuals within populations, accounting for 99.66% of the total according to the analysis of molecular variance (AMOVA) results. Each population also maintains the connection because of certain gene communication, so that the genetic differentiation between populations was not very significant. The STRUCTURE software was used to analyse the population structure and the result showed that 88 accessions were clustered into three groups, and 67% of them were pure type, which was also confirmed through PCA. Conclusions: The comprehensive study of phenotypic, chemical and molecular markers will provide valuable information for future breeding plans and understanding the phylogenetic relationship of MO population. [ABSTRACT FROM AUTHOR]

Background: Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes. Results: Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation. Conclusions: The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes. [ABSTRACT FROM AUTHOR]

Population-scale expression profiling studies can provide valuable insights into biological and disease-underlying mechanisms. The availability of phenotypic traits is essential for studying clinical effects. Therefore, missing, incomplete, or inaccurate phenotypic information can make analyses challenging and prevent RNA-seq or other omics data to be reused. A possible solution are predictors that infer clinical or behavioral phenotypic traits from molecular data. While such predictors have been developed based on different omics data types and are being applied in various studies, metabolomics-based surrogates are less commonly used than predictors based on DNA methylation profiles.In this study, we inferred 17 traits, including diabetes status and exposure to lipid medication, using previously trained metabolomic predictors. We evaluated whether these metabolomic surrogates can be used as an alternative to reported information for studying the respective phenotypes using expression profiling data of four population cohorts. For the majority of the 17 traits, the metabolomic surrogates performed similarly to the reported phenotypes in terms of effect sizes, number of significant associations, replication rates, and significantly enriched pathways.The application of metabolomics-derived surrogate outcomes opens new possibilities for reuse of multi-omics data sets. In studies where availability of clinical metadata is limited, missing or incomplete information can be complemented by these surrogates, thereby increasing the size of available data sets. Additionally, the availability of such surrogates could be used to correct for potential biological confounding. In the future, it would be interesting to further investigate the use of molecular predictors across different omics types and cohorts. [ABSTRACT FROM AUTHOR]

Background: Multi-omics represent a promising link between phenotypes and genome variation. Few studies yet address their integration to understand genetic architecture and improve predictability. Results: Our study used 241 poplar genotypes, phenotyped in two common gardens, with xylem and cambium RNA sequenced at one site, yielding large phenotypic, genomic (SNP), and transcriptomic datasets. Prediction models for each trait were built separately for SNPs and transcripts, and compared to a third model integrated by concatenation of both omics. The advantage of integration varied across traits and, to understand such differences, an eQTL analysis was performed to characterize the interplay between the genome and transcriptome and classify the predicting features into cis or trans relationships. A strong, significant negative correlation was found between the change in predictability and the change in predictor ranking for trans eQTLs for traits evaluated in the site of transcriptomic sampling. Conclusions: Consequently, beneficial integration happens when the redundancy of predictors is decreased, likely leaving the stage to other less prominent but complementary predictors. An additional gene ontology (GO) enrichment analysis appeared to corroborate such statistical output. To our knowledge, this is a novel finding delineating a promising method to explore data integration. [ABSTRACT FROM AUTHOR]

Background: The pigmentation phenotype diversity is rich in domestic goats, and identification of the genetic loci affecting coat color in goats has long been of interest. Via the detections of selection signatures, a duplication upstream ASIP was previously reported to be a variant affecting the Swiss markings depigmentation phenotype in goats. Results: We conducted a genome-wide association study using whole-genome sequencing (WGS) data to identify the genetic loci and causal variants affecting the pigmentation phenotype in 65 Jintang black (JT) goats (i.e., 48 solid black vs. 17 non-classic Swiss markings). Although a single association peak harboring the ASIP gene at 52,619,845–72,176,538 bp on chromosome 13 was obtained using a linear mixed model approach, all the SNPs and indels in this region were excluded as causal variants for the pigmentation phenotype. We then found that all 17 individuals with non-classic Swiss markings carried a 13,420-bp duplication (CHI13:63,129,198–63,142,617 bp) nearly 101 kb upstream of ASIP, and this variant was strongly associated (P = 1.48 × 10− 12) with the coat color in the 65 JT goats. The copy numbers obtained from the WGS data also showed that the duplication was present in all 53 goats from three European breeds with Swiss markings and absent in 45 of 51 non-Swiss markings goats from four other breeds and 21 Bezoars, which was further validated in 314 samples from seven populations based on PCR amplification. The copy numbers of the duplication vary in different goat breeds with Swiss markings, indicating a threshold effect instead of a dose-response effect at the molecular level. Furthermore, breakpoint flanking repeat analysis revealed that the duplication was likely to be a result of the Bov-B-mediated nonallelic homologous recombination. Conclusion: We confirmed that a genomic region harboring the ASIP gene is a major locus affecting the coat color phenotype of Swiss markings in goats. Although the molecular genetic mechanisms remain unsolved, the 13,420-bp duplication upstream of ASIP is a necessary but not sufficient condition for this phenotype in goats. Moreover, the variations in the copy number of the duplication across different goat breeds do not lead to phenotypic heterogeneity. [ABSTRACT FROM AUTHOR]

Background: As genomic sequencing moves closer to clinical implementation, there has been an increasing acceptance of returning incidental findings to research participants and patients for mutations in highly penetrant, medically actionable genes. A curated list of genes has been recommended by the American College of Medical Genetics and Genomics (ACMG) for return of incidental findings. However, the pleiotropic effects of these genes are not fully known. Such effects could complicate genetic counseling when returning incidental findings. In particular, there has been no systematic evaluation of psychiatric manifestations associated with rare variation in these genes. Results: Here, we leveraged a targeted sequence panel and real-world electronic health records from the eMERGE network to assess the burden of rare variation in the ACMG-56 genes and two psychiatric-associated genes (CACNA1C and TCF4) across common mental health conditions in 15,181 individuals of European descent. As a positive control, we showed that this approach replicated the established association between rare mutations in LDLR and hypercholesterolemia with no visible inflation from population stratification. However, we did not identify any genes significantly enriched with rare deleterious variants that confer risk for common psychiatric disorders after correction for multiple testing. Suggestive associations were observed between depression and rare coding variation in PTEN (P = 1.5 × 10–4), LDLR (P = 3.6 × 10–4), and CACNA1S (P = 5.8 × 10–4). We also observed nominal associations between rare variants in KCNQ1 and substance use disorders (P = 2.4 × 10–4), and APOB and tobacco use disorder (P = 1.1 × 10–3). Conclusions: Our results do not support an association between psychiatric disorders and incidental findings in medically actionable gene mutations, but power was limited with the available sample sizes. Given the phenotypic and genetic complexity of psychiatric phenotypes, future work will require a much larger sequencing dataset to determine whether incidental findings in these genes have implications for risk of psychopathology. [ABSTRACT FROM AUTHOR]

Identifying the genes underlying fitness-related traits such as body size and male ornamentation can provide tools for conservation and management and are often subject to various selective pressures. Here we performed high-depth whole genome re-sequencing of pools of individuals representing the phenotypic extremes for antler and body size in white-tailed deer (Odocoileus virginianus). Samples were selected from a tissue repository containing phenotypic data for 4,466 male white-tailed deer from Anticosti Island, Quebec, with four pools representing the extreme phenotypes for antler and body size after controlling for age. Our results revealed a largely homogenous population but detected highly divergent windows between pools for both traits, with the mean allele frequency difference of 14% for and 13% for antler and body SNPs in outlier windows, respectively. Genes in outlier antler windows were enriched for pathways associated with cell death and protein metabolism and some of the most differentiated windows included genes associated with oncogenic pathways and reproduction, processes consistent with antler evolution and growth. Genes associated with body size were more nuanced, suggestive of a highly complex trait. Overall, this study revealed the complex genomic make-up of both antler morphology and body size in free-ranging white-tailed deer and identified target loci for additional analyses. [ABSTRACT FROM AUTHOR]

Background: Gamma-irradiated mutants of Triticum aestivum L., hexaploid wheat, provide novel and agriculturally important traits and are used as breeding materials. However, the identification of causative genomic regions of mutant phenotypes is challenging because of the large and complicated genome of hexaploid wheat. Recently, the combined use of high-quality reference genome sequences of common wheat and cost-effective resequencing technologies has made it possible to evaluate genome-wide polymorphisms, even in complex genomes. Results: To investigate whether the genome sequencing approach can effectively detect structural variations, such as deletions, frequently caused by gamma irradiation, we selected a grain-hardness mutant from the gamma-irradiated population of Japanese elite wheat cultivar "Kitahonami." The Hardness (Ha) locus, including the puroindoline protein-encoding genes Pina-D1 and Pinb-D1 on the short arm of chromosome 5D, primarily regulates the grain hardness variation in common wheat. We performed short-read genome sequencing of wild-type and grain-hardness mutant plants, and subsequently aligned their short reads to the reference genome of the wheat cultivar "Chinese Spring." Genome-wide comparisons of depth-of-coverage between wild-type and mutant strains detected ~ 130 Mbp deletion on the short arm of chromosome 5D in the mutant genome. Molecular markers for this deletion were applied to the progeny populations generated by a cross between the wild-type and the mutant. A large deletion in the region including the Ha locus was associated with the mutant phenotype, indicating that the genome sequencing is a powerful and efficient approach for detecting a deletion marker of a gamma-irradiated mutant phenotype. In addition, we investigated a pre-harvest sprouting tolerance mutant and identified a 67.8 Mbp deletion on chromosome 3B where Viviparous-B1 and GRAS family transcription factors are located. Co-dominant markers designed to detect the deletion-polymorphism confirmed the association with low germination rate, leading to pre-harvest sprouting tolerance. Conclusions: Short read-based genome sequencing of gamma-irradiated mutants facilitates the identification of large deletions linked to mutant phenotypes when combined with segregation analyses in progeny populations. This method allows effective application of mutants with agriculturally important traits in breeding using marker-assisted selection. [ABSTRACT FROM AUTHOR]

Background: The characterization of genetic diversity and population differentiation for maize inbred lines from breeding programs is of great value in assisting breeders in maintaining and potentially increasing the rate of genetic gain. In our study, we characterized a set of 187 tropical maize inbred lines from the public breeding program of the Universidade Federal de Viçosa (UFV) in Brazil based on 18 agronomic traits and 3,083 single nucleotide polymorphisms (SNP) markers to evaluate whether this set of inbred lines represents a panel of tropical maize inbred lines for association mapping analysis and investigate the population structure and patterns of relationships among the inbred lines from UFV for better exploitation in our maize breeding program. Results: Our results showed that there was large phenotypic and genotypic variation in the set of tropical maize inbred lines from the UFV maize breeding program. We also found high genetic diversity (GD = 0.34) and low pairwise kinship coefficients among the maize inbred lines (only approximately 4.00 % of the pairwise relative kinship was above 0.50) in the set of inbred lines. The LD decay distance over all ten chromosomes in the entire set of maize lines with r2 = 0.1 was 276,237 kb. Concerning the population structure, our results from the model-based STRUCTURE and principal component analysis methods distinguished the inbred lines into three subpopulations, with high consistency maintained between both results. Additionally, the clustering analysis based on phenotypic and molecular data grouped the inbred lines into 14 and 22 genetic divergence clusters, respectively. Conclusions: Our results indicate that the set of tropical maize inbred lines from UFV maize breeding programs can comprise a panel of tropical maize inbred lines suitable for a genome-wide association study to dissect the variation of complex quantitative traits in maize, mainly in tropical environments. In addition, our results will be very useful for assisting us in the assignment of heterotic groups and the selection of the best parental combinations for new breeding crosses, mapping populations, mapping synthetic populations, guiding crosses that target highly heterotic and yielding hybrids, and predicting untested hybrids in the public breeding program UFV. [ABSTRACT FROM AUTHOR]

Abstract Background Efforts to elucidate the function of enhancers in vivo are underway but their vast numbers alongside differing enhancer architectures make it difficult to determine their impact on gene activity. By systematically annotating multiple mouse tissues with super- and typical-enhancers, we have explored their relationship with gene function and phenotype. Results Though super-enhancers drive high total- and tissue-specific expression of their associated genes, we find that typical-enhancers also contribute heavily to the tissue-specific expression landscape on account of their large numbers in the genome. Unexpectedly, we demonstrate that both enhancer types are preferentially associated with relevant ‘tissue-type’ phenotypes and exhibit no difference in phenotype effect size or pleiotropy. Modelling regulatory data alongside molecular data, we built a predictive model to infer gene-phenotype associations and use this model to predict potentially novel disease-associated genes. Conclusion Overall our findings reveal that differing enhancer architectures have a similar impact on mammalian phenotypes whilst harbouring differing cellular and expression effects. Together, our results systematically characterise enhancers with predicted phenotypic traits endorsing the role for both types of enhancers in human disease and disorders.

Intellectual disability (ID) can be caused by non-genetic and genetic factors, the latter being responsible for more than 1700 ID-related disorders. The broad ID phenotypic and genetic heterogeneity, as well as the difficulty in the establishment of the inheritance pattern, often result in a delay in the diagnosis. It has become apparent that massive parallel sequencing can overcome these difficulties. In this review we address: (i) ID genetic aetiology, (ii) clinical/medical settings testing, (iii) massive parallel sequencing, (iv) variant filtering and prioritization, (v) variant classification guidelines and functional studies, and (vi) ID diagnostic yield. Furthermore, the need for a constant update of the methodologies and functional tests, is essential. Thus, international collaborations, to gather expertise, data and resources through multidisciplinary contributions, are fundamental to keep track of the fast progress in ID gene discovery. [ABSTRACT FROM AUTHOR]

Background: Pyramiding different resistance genes into one plant genotype confers enhanced resistance at the phenotypic level, but the molecular mechanisms underlying this effect are not well-understood. In soybean, aphid resistance is conferred by Rag genes. We compared the transcriptional response of four soybean genotypes to aphid feeding to assess how the combination of Rag genes enhanced the soybean resistance to aphid infestation. Results: A strong synergistic interaction between Rag1 and Rag2, defined as genes differentially expressed only in the pyramid genotype, was identified. This synergistic effect in the Rag1/2 phenotype was very evident early (6 h after infestation) and involved unique biological processes. However, the response of susceptible and resistant genotypes had a large overlap 12 h after aphid infestation. Transcription factor (TF) analyses identified a network of interacting TF that potentially integrates signaling from Rag1 and Rag2 to produce the unique Rag1/2 response. Pyramiding resulted in rapid induction of phytochemicals production and deposition of lignin to strengthen the secondary cell wall, while repressing photosynthesis. We also identified Glyma.07G063700 as a novel, strong candidate for the Rag1 gene. Conclusions: The synergistic interaction between Rag1 and Rag2 in the Rag1/2 genotype can explain its enhanced resistance phenotype. Understanding molecular mechanisms that support enhanced resistance in pyramid genotypes could facilitate more directed approaches for crop improvement. [ABSTRACT FROM AUTHOR]

Background: Recently, ProFAT and BATLAS studies identified brown and white adipocytes marker genes based on analysis of large databases. They offered scores to determine the thermogenic status of adipocytes using the gene-expression data of these markers. In this work, we investigated the functional context of these genes. Results: Gene Set Enrichment Analyses (KEGG, Reactome) of the BATLAS and ProFAT marker-genes identified pathways deterministic in the formation of brown and white adipocytes. The collection of the annotated proteins of the defined pathways resulted in expanded white and brown characteristic protein-sets, which theoretically contain all functional proteins that could be involved in the formation of adipocytes. Based on our previously obtained RNA-seq data, we visualized the expression profile of these proteins coding genes and found patterns consistent with the two adipocyte phenotypes. The trajectory of the regulatory processes could be outlined by the transcriptional profile of progenitor and differentiated adipocytes, highlighting the importance of suppression processes in browning. Protein interaction network-based functional genomics by STRING, Cytoscape and R-Igraph platforms revealed that different biological processes shape the brown and white adipocytes and highlighted key regulatory elements and modules including GAPDH-CS, DECR1, SOD2, IL6, HRAS, MTOR, INS-AKT, ERBB2 and 4-NFKB, and SLIT-ROBO-MAPK. To assess the potential role of a particular protein in shaping adipocytes, we assigned interaction network location-based scores (betweenness centrality, number of bridges) to them and created a freely accessible platform, the AdipoNET (https//adiponet.com), to conveniently use these data. The Eukaryote Promoter Database predicted the response elements in the UCP1 promoter for the identified, potentially important transcription factors (HIF1A, MYC, REL, PPARG, TP53, AR, RUNX, and FoxO1). Conclusion: Our integrative approach-based results allowed us to investigate potential regulatory elements of thermogenesis in adipose tissue. The analyses revealed that some unique biological processes form the brown and white adipocyte phenotypes, which presumes the existence of the transitional states. The data also suggests that the two phenotypes are not mutually exclusive, and differentiation of thermogenic adipocyte requires induction of browning as well as repressions of whitening. The recognition of these simultaneous actions and the identified regulatory modules can open new direction in obesity research. [ABSTRACT FROM AUTHOR]

Background: Advancements in statistical methods and sequencing technology have led to numerous novel discoveries in human genetics in the past two decades. Among phenotypes of interest, most attention has been given to studying genetic associations with continuous or binary traits. Efficient statistical methods have been proposed and are available for both types of traits under different study designs. However, for multinomial categorical traits in related samples, there is a lack of efficient statistical methods and software. Results: We propose an efficient score test to analyze a multinomial trait in family samples, in the context of genome-wide association/sequencing studies. An alternative Wald statistic is also proposed. We also extend the methodology to be applicable to ordinal traits. We performed extensive simulation studies to evaluate the type-I error of the score test, Wald test compared to the multinomial logistic regression for unrelated samples, under different allele frequency and study designs. We also evaluate the power of these methods. Results show that both the score and Wald tests have a well-controlled type-I error rate, but the multinomial logistic regression has an inflated type-I error rate when applied to family samples. We illustrated the application of the score test with an application to the Framingham Heart Study to uncover genetic variants associated with diabesity, a multi-category phenotype. Conclusion: Both proposed tests have correct type-I error rate and similar power. However, because the Wald statistics rely on computer-intensive estimation, it is less efficient than the score test in terms of applications to large-scale genetic association studies. We provide computer implementation for both multinomial and ordinal traits. [ABSTRACT FROM AUTHOR]

Background: With the emphasis on analysing genotype-by-environment interactions within the framework of genomic selection and genome-wide association analysis, there is an increasing demand for reliable tools that can be used to simulate large-scale genomic data in order to assess related approaches. Results: We proposed a theory to simulate large-scale genomic data on genotype-by-environment interactions and added this new function to our developed tool GPOPSIM. Additionally, a simulated threshold trait with large-scale genomic data was also added. The validation of the simulated data indicated that GPOSPIM2.0 is an efficient tool for mimicking the phenotypic data of quantitative traits, threshold traits, and genetically correlated traits with large-scale genomic data while taking genotype-by-environment interactions into account. Conclusions: This tool is useful for assessing genotype-by-environment interactions and threshold traits methods. [ABSTRACT FROM AUTHOR]

Background: The genus Xanthomonas has long been considered to consist predominantly of plant pathogens, but over the last decade there has been an increasing number of reports on non-pathogenic and endophytic members. As Xanthomonas species are prevalent pathogens on a wide variety of important crops around the world, there is a need to distinguish between these plant-associated phenotypes. To date a large number of Xanthomonas genomes have been sequenced, which enables the application of machine learning (ML) approaches on the genome content to predict this phenotype. Until now such approaches to the pathogenomics of Xanthomonas strains have been hampered by the fragmentation of information regarding pathogenicity of individual strains over many studies. Unification of this information into a single resource was therefore considered to be an essential step. Results: Mining of 39 papers considering both plant-associated phenotypes, allowed for a phenotypic classification of 578 Xanthomonas strains. For 65 plant-pathogenic and 53 non-pathogenic strains the corresponding genomes were available and de novo annotated for the presence of Pfam protein domains used as features to train and compare three ML classification algorithms; CART, Lasso and Random Forest. Conclusion: The literature resource in combination with recursive feature extraction used in the ML classification algorithms provided further insights into the virulence enabling factors, but also highlighted domains linked to traits not present in pathogenic strains. [ABSTRACT FROM AUTHOR]

Background: The salmon louse (Lepeophtheirus salmonis) is an obligate ectoparasitic copepod living on Atlantic salmon and other salmonids in the marine environment. Salmon lice cause a number of environmental problems and lead to large economical losses in aquaculture every year. In order to develop novel parasite control strategies, a better understanding of the mechanisms of moulting and development of the salmon louse at the transcriptional level is required. Methods: Three weighted gene co-expression networks were constructed based on the pairwise correlations of salmon louse gene expression profiles at different life stages. Network-based approaches and gene annotation information were applied to identify genes that might be important for the moulting and development of the salmon louse. RNA interference was performed for validation. Regulatory impact factors were calculated for all the transcription factor genes by examining the changes in co-expression patterns between transcription factor genes and deferentially expressed genes in middle stages and moulting stages. Results: Eight gene modules were predicted as important, and 10 genes from six of the eight modules have been found to show observable phenotypes in RNA interference experiments. We knocked down five hub genes from three modules and observed phenotypic consequences in all experiments. In the infection trial, no copepodids with a RAB1A-like gene knocked down were found on fish, while control samples developed to chalimus-1 larvae. Also, a FOXO-like transcription factor obtained highest scores in the regulatory impact factor calculation. Conclusions: We propose a gene co-expression network-based approach to identify genes playing an important role in the moulting and development of salmon louse. The RNA interference experiments confirm the effectiveness of our approach and demonstrated the indispensable role of a RAB1A-like gene in the development of the salmon louse. We propose that our approach could be generalized to identify important genes associated with a phenotype of interest in other organisms. [ABSTRACT FROM AUTHOR]

Background: Rainbow trout (Oncorhynchus mykiss) is a salmonid species with a complex life-history. Wild populations are naturally divided into freshwater residents and sea-run migrants. Migrants undergo an energy-demanding adaptation for life in seawater, known as smoltification, while freshwater residents display these changes in an attenuated magnitude and rate. Despite this, in seawater rainbow trout farming all fish are transferred to seawater. Under these circumstances, weeks after seawater transfer, a significant portion of the fish die (around 10%) or experience growth stunting (GS; around 10%), which represents an important profitability and welfare issue. The underlying causes leading to GS in seawater-transferred rainbow trout remain unknown. In this study, we aimed at characterising the GS phenotype in seawater-transferred rainbow trout using untargeted and targeted approaches. To this end, the liver proteome (LC-MS/MS) and lipidome (LC-MS) of GS and fast-growing phenotypes were profiled to identify molecules and processes that are characteristic of the GS phenotype. Moreover, the transcription, abundance or activity of key proteins and hormones related to osmoregulation (Gill Na+, K + –ATPase activity), growth (plasma IGF-I, and liver igf1, igfbp1b, ghr1 and ctsl) and stress (plasma cortisol) were measured using targeted approaches. Results: No differences in Gill Na+, K + –ATPase activity and plasma cortisol were detected between the two groups. However, a significant downregulation in plasma IGF-I and liver igf1 transcription pointed at this growth factor as an important pathomechanism for GS. Changes in the liver proteome revealed reactive-oxygen-species-mediated endoplasmic reticulum stress as a key mechanism underlying the GS phenotype. From the lipidomic analysis, key observations include a reduction in triacylglycerols and elevated amounts of cardiolipins, a characteristic lipid class associated with oxidative stress, in GS phenotype. Conclusion: While the triggers to the activation of endoplasmic reticulum stress are still unknown, data from this study point towards a nutritional deficiency as an underlying driver of this phenotype. [ABSTRACT FROM AUTHOR]

Background: Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds. Results: Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes. Conclusions: Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds. [ABSTRACT FROM AUTHOR]

Background: The aerial organs of most terrestrial plants are covered by cuticular waxes, which impart plants a glaucous appearance and play important roles in protecting against various biotic and abiotic stresses. Despite many glossy green (wax-defective) mutants being well characterized in model plants, little is known about the genetic basis of glossy green mutant in broccoli. Results: B156 is a spontaneous broccoli mutant showing a glossy green phenotype. Detection by scanning electron microscopy (SEM) and chromatography-mass spectrometry (GC-MS) revealed that B156 is a cuticular wax-defective mutant, lacking waxes mostly longer than C28. Inheritance analysis revealed that this trait was controlled by a single recessive gene, BoGL5. Whole-genome InDel markers were developed, and a segregating F2 population was constructed to map BoGL5. Ultimately, BoGL5 was mapped to a 94.1 kb interval on C01. The BoCER2 gene, which is homologous to the Arabidopsis CER2 gene, was identified as a candidate of BoGL5 from the target interval. Sequence analyses revealed that Bocer2 in B156 harbored a G-to-T SNP mutation at the 485th nucleotide of the CDS, resulting in a W-to-L transition at the 162nd amino acid, a conserved site adjacent to an HXXXD motif of the deduced protein sequence. Expression analysis revealed that BoCER2 was significantly down-regulated in the leaves, stems, and siliques of B156 mutant than that of B3. Last, ectopic expression of BoCER2 in A. thaliana could, whereas Bocer2 could not, rescue the phenotype of cer2 mutant. Conclusions: Overall, this study mapped the locus determining glossy phenotype of B156 and proved BoCER2 is functional gene involved in cuticular wax biosynthesis which would promotes the utilization of BoCER2 to enhance plant resistance to biotic and abiotic stresses, and breeding of B. oleracea cultivars with glossy traits. [ABSTRACT FROM AUTHOR]

Background: There is a long-term interest in investigating the genetic basis of the horned/polled phenotype in domestic goats. Here, we report a genome-wide association study (GWAS) to detect the genetic loci affecting the polled phenotype in goats. Results: We obtained a total of 13,980,209 biallelic SNPs, using the genotyping-by-sequencing data from 45 Jintang Black (JT) goats, which included 32 female and nine male goats, and four individuals with the polled intersex syndrome (PIS). Using a mixed-model based GWAS, we identified two association signals, which were located at 150,334,857–150,817,260 bp (P = 5.15 × 10− 119) and 128,286,704–131,306,537 bp (P = 2.74 × 10− 15) on chromosome 1. The genotype distributions of the 14 most significantly associated SNPs were completely correlated with horn status in goats, based on the whole-genome sequencing (WGS) data from JT and two other Chinese horned breeds. However, variant annotation suggested that none of the detected SNPs within the associated regions were plausible causal mutations. Via additional read-depth analyses and visual inspections of WGS data, we found a 10.1-kb deletion (CHI1:g. 129424781_129434939del) and a 480-kb duplication (CHI1:150,334,286–150,818,098 bp) encompassing two genes KCNJ15 and ERG in the associated regions of polled and PIS-affected goats. Notably, the 10.1-kb deletion also served as the insertion site for the 480-kb duplication, as validated by PCR and Sanger sequencing. Our WGS genotyping showed that all horned goats were homozygous for the reference alleles without either the structural variants (SVs), whereas the PIS-affected goats were homozygous for both the SVs. We also demonstrated that horned, polled, and PIS-affected individuals among 333 goats from JT and three other Chinese horned breeds can be accurately classified via PCR amplification and agarose gel electrophoresis of two fragments in both SVs. Conclusion: Our results revealed that two genomic regions on chromosome 1 are major loci affecting the polled phenotypes in goats. We provided a diagnostic PCR to accurately classify horned, polled, and PIS-affected goats, which will enable a reliable genetic test for the early-in-life prediction of horn status in goats. [ABSTRACT FROM AUTHOR]

Background: Plants grown under shade are exposed to low red/far-red ratio, thereby triggering an array of altered phenotypes called shade avoidance syndrome (SAS). Shade negatively influences plant growth, leading to a reduction in agricultural productivity. Understanding of SAS is crucial for sustainable agricultural practices, especially for high-density indoor farming. Brassicaceae vegetables are widely consumed around the world and are commonly cultivated in indoor farms. However, our understanding of SAS in Brassicaceae vegetables and their genome-wide transcriptional regulatory networks are still largely unexplored. Results: Shade induced common signs of SAS, including hypocotyl elongation and reduced carotenoids/anthocyanins biosynthesis, in two different Brassicaceae species: Brassica rapa (Choy Sum and Pak Choy) and Brassica oleracea (Kai Lan). Phenotype-assisted transcriptome analysis identified a set of genes induced by shade in these species, many of which were related to auxin biosynthesis and signaling [e.g. YUCCA8 (YUC8), YUC9, and INDOLE-3-ACETIC ACID INDUCIBLE (IAAs)] and other phytohormones signaling pathways including brassinosteroids and ethylene. The genes functioning in plant defense (e.g. MYB29 and JASMONATE-ZIM-DOMAIN PROTEIN 9) as well as in biosynthesis of anthocyanins and glucosinolates were repressed upon shade. Besides, each species also exhibited distinct SAS phenotypes. Shade strongly reduced primary roots and elongated petioles of B. oleracea, Kai Lan. However, these SAS phenotypes were not clearly recognized in B. rapa, Choy Sum and Pak Choy. Some auxin signaling genes (e.g. AUXIN RESPONSE FACTOR 19, IAA10, and IAA20) were specifically induced in B. oleracea, while homologs in B. rapa were not up-regulated under shade. Contrastingly, shade-exposed B. rapa vegetables triggered the ethylene signaling pathway earlier than B. oleracea, Kai Lan. Interestingly, shade induced the transcript levels of LONG HYPOCOTYL IN FAR-RED 1 (HFR1) homolog in only Pak Choy as B. rapa. As HFR1 is a key negative regulator of SAS in Arabidopsis, our finding suggests that Pak Choy HFR1 homolog may also function in conferring higher shade tolerance in this variety. Conclusions: Our study shows that two Brassicaceae species not only share a conserved SAS mechanism but also exhibit distinct responses to shade, which will provide comprehensive information to develop new shade-tolerant cultivars that are suitable for high-density indoor farms. [ABSTRACT FROM AUTHOR]

Background: The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on identifying genomic regions affected by different selection pressures, mainly highlighting the recent positive selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the Mangalarga Marchador genome. Besides, we seek to direct the discussion about genes and traits of importance in this breed, especially traits related to the type and quality of gait, temperament, conformation, and locomotor system. Results: Three different methods were used to search for signals of selection: Tajima's D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n = 62) and females (n = 130) that were initially chosen considering well-defined phenotypes for gait: picada (n = 86) and batida (n = 106). All horses were genotyped using a 670 k Axiom® Equine Genotyping Array​ (Axiom MNEC670). In total, 27, 104 (chosen), and 38 candidate genes were observed within the signatures of selection identified in TD, iHS, and ROH analyses, respectively. The genes are acting in essential biological processes. The enrichment analysis highlighted the following functions: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system. Conclusions: Our findings revealed evidence of signatures of selection in the MM breed that encompass genes acting on athletic performance, limb development, and energy to muscle activity, with the particular involvement of the HOX family genes. The genome of MM is marked by recent positive selection. However, Tajima's D and iHS results point also to the presence of balancing selection in specific regions of the genome. [ABSTRACT FROM AUTHOR]

Background: Despite increasing interest in γδ T cells and their non-classical behaviour, most studies focus on animals with low numbers of circulating γδ T cells, such as mice and humans. Arguably, γδ T cell functions might be more prominent in chickens where these cells form a higher proportion of the circulatory T cell compartment. The TCR repertoire defines different subsets of γδ T cells, and such analysis is facilitated by well-annotated TCR loci. γδ T cells are considered at the cusp of innate and adaptive immunity but most functions have been identified in γδ low species. A deeper understanding of TCR repertoire biology in γδ high and γδ low animals is critical for defining the evolution of the function of γδ T cells. Repertoire dynamics will reveal populations that can be classified as innate-like or adaptive-like as well as those that straddle this definition. Results: Here, a recent discrepancy in the structure of the chicken TCR gamma locus is resolved, demonstrating that tandem duplication events have shaped the evolution of this locus. Importantly, repertoire sequencing revealed large differences in the usage of individual TRGV genes, a pattern conserved across multiple tissues, including thymus, spleen and the gut. A single TRGV gene, TRGV3.3, with a highly diverse private CDR3 repertoire dominated every tissue in all birds. TRGV usage patterns were partly explained by the TRGV-associated recombination signal sequences. Public CDR3 clonotypes represented varying proportions of the repertoire of TCRs utilising different TRGVs, with one TRGV dominated by super-public clones present in all birds. Conclusions: The application of repertoire analysis enabled functional annotation of the TCRG locus in a species with a high circulating γδ phenotype. This revealed variable usage of TCRGV genes across multiple tissues, a pattern quite different to that found in γδ low species (human and mouse). Defining the repertoire biology of avian γδ T cells will be key to understanding the evolution and functional diversity of these enigmatic lymphocytes in an animal that is numerically more reliant on them. Practically, this will reveal novel ways in which these cells can be exploited to improve health in medical and veterinary contexts. [ABSTRACT FROM AUTHOR]

Background: The flesh pigmentation of farmed Atlantic salmon is formed by accumulation of carotenoids derived from commercial diets. In the salmon gastrointestinal system, the hindgut is considered critical in the processes of carotenoids uptake and metabolism. In Tasmania, flesh color depletion can noticeably affect farmed Atlantic salmon at different levels of severity following extremely hot summers. In this study, RNA sequencing (RNA-Seq) was performed to investigate the reduction in flesh pigmentation. Library preparation is a key step that significantly impacts the effectiveness of RNA sequencing (RNA-Seq) experiments. Besides the commonly used whole transcript RNA-Seq method, the 3' mRNA-Seq method is being applied widely, owing to its reduced cost, enabling more repeats to be sequenced at the expense of lower resolution. Therefore, the output of the Illumina TruSeq kit (whole transcript RNA-Seq) and the Lexogen QuantSeq kit (3' mRNA-Seq) was analyzed to identify genes in the Atlantic salmon hindgut that are differentially expressed (DEGs) between two flesh color phenotypes. Results: In both methods, DEGs between the two color phenotypes were associated with metal ion transport, oxidation-reduction processes, and immune responses. We also found DEGs related to lipid metabolism in the QuantSeq method. In the TruSeq method, a missense mutation was detected in DEGs in different flesh color traits. The number of DEGs found in the TruSeq libraries was much higher than the QuantSeq; however, the trend of DEGs in both library methods was similar and validated by qPCR. Conclusions: Flesh coloration in Atlantic salmon is related to lipid metabolism in which apolipoproteins, serum albumin and fatty acid-binding protein genes are hypothesized to be linked to the absorption, transport and deposition of carotenoids. Our findings suggest that Grp could inhibit the feeding behavior of low color-banded fish, resulting in the dietary carotenoid shortage. Several SNPs in genes involving in carotenoid-binding cholesterol and oxidative stress were detected in both flesh color phenotypes. Regarding the choice of the library preparation method, the selection criteria depend on the research design and purpose. The 3' mRNA-Seq method is ideal for targeted identification of highly expressed genes, while the whole RNA-Seq method is recommended for identification of unknown genes, enabling the identification of splice variants and trait-associated SNPs, as we have found for duox2 and duoxa1. [ABSTRACT FROM AUTHOR]

Background: Isatidis Radix, the root of Isatis indigotica Fort. (Chinese woad) can produce a variety of efficacious compound with medicinal properties. The tetraploid I. indigotica plants exhibit superior phenotypic traits, such as greater yield, higher bioactive compounds accumulation and enhanced stress tolerance. In this study, a comparative transcriptomic and metabolomic study on Isatidis Radix autotetraploid and its progenitor was performed. Results: Through the targeted metabolic profiling, 283 metabolites were identified in Isatidis Radix, and 70 polyploidization-altered metabolites were obtained. Moreover, the production of lignans was significantly increased post polyploidization, which implied that polyploidization-modulated changes in lignan biosynthesis. Regarding the transcriptomic shift, 2065 differentially expressed genes (DEGs) were identified as being polyploidy-responsive genes, and the polyploidization-altered DEGs were enriched in phenylpropanoid biosynthesis and plant hormone signal transduction. The further integrative analysis of polyploidy-responsive metabolome and transcriptome showed that 1584 DEGs were highly correlated with the 70 polyploidization-altered metabolites, and the transcriptional factors TFs-lignans network highlighted 10 polyploidy-altered TFs and 17 fluctuated phenylpropanoid pathway compounds. Conclusions: These results collectively indicated that polyploidization contributed to the high content of active compounds in autotetraploid roots, and the gene–lignan pathway network analysis highlighted polyploidy–responsive key functional genes and regulators. [ABSTRACT FROM AUTHOR]

Background: The human microbiome is inherently dynamic and its dynamic nature plays a critical role in maintaining health and driving disease. With an increasing number of longitudinal microbiome studies, scientists are eager to learn the comprehensive characterization of microbial dynamics and their implications to the health and disease-related phenotypes. However, due to the challenging structure of longitudinal microbiome data, few analytic methods are available to characterize the microbial dynamics over time. Results: We propose a microbial trend analysis (MTA) framework for the high-dimensional and phylogenetically-based longitudinal microbiome data. In particular, MTA can perform three tasks: 1) capture the common microbial dynamic trends for a group of subjects at the community level and identify the dominant taxa; 2) examine whether or not the microbial overall dynamic trends are significantly different between groups; 3) classify an individual subject based on its longitudinal microbial profiling. Our extensive simulations demonstrate that the proposed MTA framework is robust and powerful in hypothesis testing, taxon identification, and subject classification. Our real data analyses further illustrate the utility of MTA through a longitudinal study in mice. Conclusions: The proposed MTA framework is an attractive and effective tool in investigating dynamic microbial pattern from longitudinal microbiome studies. [ABSTRACT FROM AUTHOR]

Background: RNA polymerase II plays critical roles in transcription in eukaryotic organisms. C-terminal Domain Phosphatase-like 1 (CPL1) regulates the phosphorylation state of the C-terminal domain of RNA polymerase II subunit B1, which is critical in determining RNA polymerase II activity. CPL1 plays an important role in miRNA biogenesis, plant growth and stress responses. Although cpl1 mutant showes delayed-flowering phenotype, the molecular mechanism behind CPL1's role in floral transition is still unknown. Results: To study the role of CPL1 during the floral transition, we first tested phenotypes of cpl1-3 mutant, which harbors a point-mutation. The cpl1-3 mutant contains a G-to-A transition in the second exon, which results in an amino acid substitution from Glu to Lys (E116K). Further analyses found that the mutated amino acid (Glu) was conserved in these species. As a result, we found that the cpl1-3 mutant experienced delayed flowering under both long- and short-day conditions, and CPL1 is involved in the vernalization pathway. Transcriptome analysis identified 109 genes differentially expressed in the cpl1 mutant, with 2 being involved in floral transition. Differential expression of the two flowering-related DEGs was further validated by qRT-PCR. Conclusions: Flowering genetic pathways analysis coupled with transciptomic analysis provides potential genes related to floral transition in the cpl1-3 mutant, and a framework for future studies of the molecular mechanisms behind CPL1's role in floral transition. [ABSTRACT FROM AUTHOR]

Background: Association mapping studies of quantitative trait loci (QTL) for canine hip dysplasia (CHD) can contribute to the understanding of the genetic background of this common and debilitating disease and might contribute to its genetic improvement. The power of association studies for CHD is limited by relatively small sample numbers for CHD records within countries, suggesting potential benefits of joining data across countries. However, this is complicated due to the use of different scoring systems across countries. In this study, we incorporated routinely assessed CHD records and genotype data of German Shepherd dogs from two countries (UK and Sweden) to perform genome-wide association studies (GWAS) within populations using different variations of CHD phenotypes. As phenotypes, dogs were either classified into cases and controls based on the Fédération Cynologique Internationale (FCI) five-level grading of the worst hip or the FCI grade was treated as an ordinal trait. In a subsequent meta-analysis, we added publicly available data from a Finnish population and performed the GWAS across all populations. Genetic associations for the CHD phenotypes were evaluated in a linear mixed model using 62,089 SNPs. Results: Multiple SNPs with genome-wide significant and suggestive associations were detected in single-population GWAS and the meta-analysis. Few of these SNPs overlapped between populations or between single-population GWAS and the meta-analysis, suggesting that many CHD-related QTL are population-specific. More significant or suggestive SNPs were identified when FCI grades were used as phenotypes in comparison to the case-control approach. MED13 (Chr 9) and PLEKHA7 (Chr 21) emerged as novel positional candidate genes associated with hip dysplasia. Conclusions: Our findings confirm the complex genetic nature of hip dysplasia in dogs, with multiple loci associated with the trait, most of which are population-specific. Routinely assessed CHD information collected across countries provide an opportunity to increase sample sizes and statistical power for association studies. While the lack of standardisation of CHD assessment schemes across countries poses a challenge, we showed that conversion of traits can be utilised to overcome this obstacle. [ABSTRACT FROM AUTHOR]