ALBERT

Description: Background Historically, geneticists have relied on genotyping arrays and imputation to study human genetic variation. However, an underrepresentation of diverse populations has resulted in arrays that poorly capture global genetic variation, and a lack of reference panels. This has contributed to deepening global health disparities. Whole genome sequencing (WGS) better captures genetic variation but remains prohibitively expensive. Thus, we explored WGS at “mid-pass” 1-7x coverage. Results Here, we developed and benchmarked methods for mid-pass sequencing. When applied to a population without an existing genomic reference panel, 4x mid-pass performed consistently well across ethnicities, with high recall (98%) and precision (97.5%). Conclusion Compared to array data imputed into 1000 Genomes, mid-pass performed better across all metrics and identified novel population-specific variants with potential disease relevance. We hope our work will reduce financial barriers for geneticists from underrepresented populations to characterize their genomes prior to biomedical genetic applications.

Description: Background Kuruma shrimp, a major commercial shrimp species in the world, has two cryptic or sibling species, Marsupenaeus japonicus and Marsupenaeus pulchricaudatus. Codon usage analysis would contribute to our understanding of the genetic and evolutionary characteristics of the two Marsupenaeus species. In this study, we analyzed codon usage and related indices using coding sequences (CDSs) from RNA-seq data. Results Using CodonW 1.4.2 software, we performed the codon bias analysis of transcriptomes obtained from hepatopancreas tissues, which indicated weak codon bias. Almost all parameters had similar correlations for both species. The gene expression level (FPKM) was negatively correlated with A/T3s. We determined 12 and 14 optimal codons for M. japonicus and M. pulchricaudatus, respectively, and all optimal codons have a C/G-ending. The two Marsupenaeus species had different usage frequencies of codon pairs, which contributed to further analysis of transcriptional differences between them. Orthologous genes that underwent positive selection (ω 〉 1) had a higher correlation coefficient than that of experienced purifying selection (ω 

Description: Background The molecular complexity of neural retina development remains poorly studied. Knowledge of retinal neurogenesis regulation sheds light on retinal degeneration therapy exploration. Therefore, we integrated the time-series circRNA, lncRNA, miRNA, and mRNA expression profiles of the developing retina through whole-transcriptome sequencing. The key functional ncRNAs and the ceRNA network regulating retinal neurogenesis were identified. Results Transcriptomic analysis identified circRNA as the most variable ncRNA subtype. We screened a series of neurogenesis-related circRNAs, lncRNAs, and miRNAs using different strategies based on their diversified molecular functions. The expression of circCDYL, circATXN1, circDYM, circPRGRIP, lncRNA Meg3, and lncRNA Vax2os was validated by quantitative real-time PCR. These circRNAs and lncRNAs participate in neurotransmitter transport and multicellular organism growth through the intricate circRNA/lncRNA-miRNA-mRNA network. Conclusion Whole-transcriptome sequencing and bioinformatics analysis systematically screened key ncRNAs in retinal neurogenesis. The validated ncRNAs and their circRNA/lncRNA-miRNA-mRNA network involve neurotransmitter transport and multicellular organism growth during retinal development.

Description: Background The evaluation of alternative splicing, including differential isoform expression and differential exon usage, can provide some insights on the transcriptional changes that occur in response to environmental perturbations. Maternal nutrition is considered a major intrauterine regulator of fetal developmental programming. The objective of this study was to assess potential changes in splicing events in the longissimus dorsi muscle of beef calves gestated under control or methionine-rich diets. RNA sequencing and whole-genome bisulfite sequencing were used to evaluate muscle transcriptome and methylome, respectively. Results Alternative splicing patterns were significantly altered by maternal methionine supplementation. Most of the altered genes were directly implicated in muscle development, muscle physiology, ATP activities, RNA splicing and DNA methylation, among other functions. Interestingly, there was a significant association between DNA methylation and differential exon usage. Indeed, among the set of genes that showed differential exon usage, significant differences in methylation level were detected between significant and non-significant exons, and between contiguous and non-contiguous introns to significant exons. Conclusions Overall, our findings provide evidence that a prenatal diet rich in methyl donors can significantly alter the offspring transcriptome, including changes in isoform expression and exon usage, and some of these changes are mediated by changes in DNA methylation.

Description: Background The gut microbiome is a diverse network of bacteria which inhabit our digestive tract and is crucial for efficient cellular metabolism, nutrient absorption, and immune system development. Spinal cord injury (SCI) disrupts autonomic function below the level of injury and can alter the composition of the gut microbiome. Studies in rodent models have shown that SCI-induced bacterial imbalances in the gut can exacerbate the spinal cord damage and impair recovery. In this study we, for the first time, characterized the composition of the gut microbiome in a Yucatan minipig SCI model. We compared the relative abundance of the most dominant bacterial phyla in control samples to those collected from animals who underwent a contusion-compression SCI at the 2nd or 10th Thoracic level. Results We identify specific bacterial fluctuations that are unique to SCI animals, which were not found in uninjured animals given the same dietary regimen or antibiotic administration. Further, we identified a specific time-frame, “SCI-acute stage”, during which many of these bacterial fluctuations occur before returning to “baseline” levels. Conclusion This work presents a dynamic view of the microbiome changes that accompany SCI, establishes a resource for future studies and to understand the changes that occur to gut microbiota after spinal cord injury and may point to a potential therapeutic target for future treatment.

Description: Background Grapevine (Vitis vinifera) productivity has been severely affected by various bacterial, viral and fungal diseases worldwide. When a plant is infected with the pathogen, various defense mechanisms are subsequently activated in plants at various molecular levels. Thus, for substantiating the disease control in an eco-friendly way, it is essential to understand the molecular mechanisms governing pathogen resistance in grapes. Results In our study, we performed genome-wide identification of various defensive genes expressed during powdery mildew (PM) and downy mildew (DM) infections in grapevine. Consequently, we identified 6, 21, 2, 5, 3 and 48 genes of Enhanced Disease Susceptibility 1 (EDS1), Non-Race-specific Disease Resistance (NDR1), Phytoalexin deficient 4 (PAD4), Nonexpressor of PR Gene (NPR), Required for Mla-specified resistance (RAR) and Pathogenesis Related (PR), respectively, in the grapevine genome. The phylogenetic study revealed that V. vinifera defensive genes are evolutionarily related to Arabidopsis thaliana. Differential expression analysis resulted in identification of 2, 4, 7, 2, 4, 1 and 7 differentially expressed Nucleotide-binding leucine rich repeat receptor (NLR), EDS1, NDR1, PAD4, NPR, RAR1 and PR respectively against PM infections and 28, 2, 5, 4, 1 and 19 differentially expressed NLR, EDS1, NDR1, NPR, RAR1 and PR respectively against DM infections in V. vinifera. The co-expression study showed the occurrence of closely correlated defensive genes that were expressed during PM and DM stress conditions. Conclusion The PM and DM responsive defensive genes found in this study can be characterized in future for impelling studies relaying fungal and oomycete resistance in plants, and the functionally validated genes would then be available for conducting in-planta transgenic gene expression studies for grapes.

Description: Background Members of the basic helix-loop-helix (bHLH) transcription factor family perform indispensable functions in various biological processes, such as plant growth, seed maturation, and abiotic stress responses. However, the bHLH family in foxtail millet (Setaria italica), an important food and feed crop, has not been thoroughly studied. Results In this study, 187 bHLH genes of foxtail millet (SibHLHs) were identified and renamed according to the chromosomal distribution of the SibHLH genes. Based on the number of conserved domains and gene structure, the SibHLH genes were divided into 21 subfamilies and two orphan genes via phylogenetic tree analysis. According to the phylogenetic tree, the subfamilies 15 and 18 may have experienced stronger expansion in the process of evolution. Then, the motif compositions, gene structures, chromosomal spread, and gene duplication events were discussed in detail. A total of sixteen tandem repeat events and thirty-eight pairs of segment duplications were identified in bHLH family of foxtail millet. To further investigate the evolutionary relationship in the SibHLH family, we constructed the comparative syntenic maps of foxtail millet associated with representative monocotyledons and dicotyledons species. Finally, the gene expression response characteristics of 15 typical SibHLH genes in different tissues and fruit development stages, and eight different abiotic stresses were analysed. The results showed that there were significant differences in the transcription levels of some SibHLH members in different tissues and fruit development stages, and different abiotic stresses, implying that SibHLH members might have different physiological functions. Conclusions In this study, we identified 187 SibHLH genes in foxtail millet and further analysed the evolution and expression patterns of the encoded proteins. The findings provide a comprehensive understanding of the bHLH family in foxtail millet, which will inform further studies on the functional characteristics of SibHLH genes.

Description: Background Faced with the ongoing global pandemic of coronavirus disease, the ‘National Reference Centre for Whole Genome Sequencing of microbial pathogens: database and bioinformatic analysis’ (GENPAT) formally established at the ‘Istituto Zooprofilattico Sperimentale dell’Abruzzo e del Molise’ (IZSAM) in Teramo (Italy) is in charge of the SARS-CoV-2 surveillance at the genomic scale. In a context of SARS-CoV-2 surveillance requiring correct and fast assessment of epidemiological clusters from substantial amount of samples, the present study proposes an analytical workflow for identifying accurately the PANGO lineages of SARS-CoV-2 samples and building of discriminant minimum spanning trees (MST) bypassing the usual time consuming phylogenomic inferences based on multiple sequence alignment (MSA) and substitution model. Results GENPAT constituted two collections of SARS-CoV-2 samples. The first collection consisted of SARS-CoV-2 positive swabs collected by IZSAM from the Abruzzo region (Italy), then sequenced by next generation sequencing (NGS) and analyzed in GENPAT (n = 1592), while the second collection included samples from several Italian provinces and retrieved from the reference Global Initiative on Sharing All Influenza Data (GISAID) (n = 17,201). The main results of the present work showed that (i) GENPAT and GISAID detected the same PANGO lineages, (ii) the PANGO lineages B.1.177 (i.e. historical in Italy) and B.1.1.7 (i.e. ‘UK variant’) are major concerns today in several Italian provinces, and the new MST-based method (iii) clusters most of the PANGO lineages together, (iv) with a higher dicriminatory power than PANGO lineages, (v) and faster that the usual phylogenomic methods based on MSA and substitution model. Conclusions The genome sequencing efforts of Italian provinces, combined with a structured national system of NGS data management, provided support for surveillance SARS-CoV-2 in Italy. We propose to build phylogenomic trees of SARS-CoV-2 variants through an accurate, discriminant and fast MST-based method avoiding the typical time consuming steps related to MSA and substitution model-based phylogenomic inference.

Description: Background Viral vectors, including adenovirus (Ad) and modified vaccinia Ankara (MVA), have gained increasing attention as vaccine platforms in recent years due to their capacity to express antigens from a wide array of pathogens, their rapid induction of humoral and cellular protective immune responses, and their relatively low production costs. In particular, the chimpanzee Ad vector, ChAdOx1, has taken centre stage as a leading COVID-19 vaccine candidate. However, despite mounting data, both clinical and pre-clinical, demonstrating effective induction of adaptive immune responses, the innate immune signals that precede the protective responses that make these vectors attractive vaccine platforms remain poorly understood. Results In this study, a mouse immunisation model was used to evaluate whole blood gene expression changes 24 h after either a single dose or heterologous prime-boost regimen of an Ad and/or MVA vaccine. We demonstrate through comparative analysis of Ad vectors encoding different antigens that a transgene product-specific gene signature can be discerned from the vector-induced transcriptional response. Expression of genes involved in TLR2 stimulation and γδ T cell and natural killer cell activation were induced after a single dose of Ad, while MVA led to greater expression of type I interferon genes. The order of prime-boost combinations was found to influence the magnitude of the gene expression changes, with MVA/Ad eliciting greater transcriptional perturbation than Ad/MVA. Contrasting the two regimens revealed significant enrichment of epigenetic regulation pathways and augmented expression of MHC class I and II molecules associated with MVA/Ad. Conclusion These data demonstrate that the order in which vaccines from heterologous prime-boost regimens are administered leads to distinct transcriptional responses and may shape the immune response induced by such combinations. The characterisation of early vaccine-induce responses strengthens our understanding of viral vector vaccine mechanisms of action ahead of their characterisation in human clinical trials and are a valuable resource to inform the pre-clinical design of appropriate vaccine constructs for emerging infectious diseases.

Description: Background High-density SNP arrays are now available for a wide range of crop species. Despite the development of many tools for generating genetic maps, the genome position of many SNPs from these arrays is unknown. Here we propose a linkage disequilibrium (LD)-based algorithm to allocate unassigned SNPs to chromosome regions from sparse genetic maps. This algorithm was tested on sugarcane, wheat, and barley data sets. We calculated the algorithm’s efficiency by masking SNPs with known locations, then assigning their position to the map with the algorithm, and finally comparing the assigned and true positions. Results In the 20-fold cross-validation, the mean proportion of masked mapped SNPs that were placed by the algorithm to a chromosome was 89.53, 94.25, and 97.23% for sugarcane, wheat, and barley, respectively. Of the markers that were placed in the genome, 98.73, 96.45 and 98.53% of the SNPs were positioned on the correct chromosome. The mean correlations between known and new estimated SNP positions were 0.97, 0.98, and 0.97 for sugarcane, wheat, and barley. The LD-based algorithm was used to assign 5920 out of 21,251 unpositioned markers to the current Q208 sugarcane genetic map, representing the highest density genetic map for this species to date. Conclusions Our LD-based approach can be used to accurately assign unpositioned SNPs to existing genetic maps, improving genome-wide association studies and genomic prediction in crop species with fragmented and incomplete genome assemblies. This approach will facilitate genomic-assisted breeding for many orphan crops that lack genetic and genomic resources.