Automated Landmarking via Multiple Templates

Geometric morphometrics based on landmark data has been increasingly used in biomedical and biological researchers for quantifying complex phenotypes. However, manual landmarking can be laborious and subject to intra and interobserver errors. This has motivated the development of automated landmarking methods. We have recently introduced ALPACA (Automated Landmarking through Point cloud Alignment and Correspondence), a fast method to automatically annotate landmarks via use of a landmark template as part of the SlicerMorph toolkit. Yet, using a single template may not consistently perform well for large study samples, especially when the sample consists of specimen with highly variable morphology, as it is common evolutionary studies. In this study, we introduce a variation on our ALPACA pipeline that supports multiple specimen templates, which we call MALPACA. We show that MALPACA outperforms ALPACA consistently by testing on two different datasets. We also introduce a method of choosing the templates that can be used in conjunction with MALPACA, when no prior information is available. This K-means method uses an approximation of the total morphological variation in the dataset to suggest samples within the population to be used as landmark templates. While we advise investigators to pay careful attention to the template selection process in any of the template-based automated landmarking approaches, our analyses show that the introduced K-means based method of templates selection is better than randomly choosing the templates. In summary, MALPACA can accommodate larger morphological disparity commonly found in evolutionary studies with performance comparable to human observer.

Gene by Environment Interactions reveal new regulatory aspects of signaling network plasticity

Phenotypes can change during exposure to different environments through the regulation of signaling pathways that operate in integrated networks. How signaling networks produce different phenotypes in different settings is not fully understood. Here, Gene by Environment Interactions (GEIs) were used to explore the regulatory network that controls filamentous/invasive growth in the yeast Saccharomyces cerevisiae. GEI analysis revealed that the regulation of invasive growth is decentralized and varies extensively across environments. Different regulatory pathways were critical or dispensable depending on the environment, microenvironment, or time point tested, and the pathway that made the strongest contribution changed depending on the environment. Some regulators even showed conditional role reversals. Ranking pathways’ roles across environments revealed an under-appreciated pathway (OPI1) as the single strongest regulator among the major pathways tested (RAS, RIM101, and MAPK). One mechanism that may explain the high degree of regulatory plasticity observed was conditional pathway interactions, such as conditional redundancy and conditional cross-pathway regulation. Another mechanism was that different pathways conditionally and differentially regulated gene expression, such as target genes that control separate cell adhesion mechanisms (FLO11 and SFG1). An exception to decentralized regulation of invasive growth was that morphogenetic changes (cell elongation and budding pattern) were primarily regulated by one pathway (MAPK). GEI analysis also uncovered a round-cell invasion phenotype. Our work suggests that GEI analysis is a simple and powerful approach to define the regulatory basis of complex phenotypes and may be applicable to many systems.

NetREx – Network-based Rice Expression Analysis Server for Abiotic Stress Conditions

Recent focus on transcriptomic studies in food crops like rice, wheat and maize provide new opportunities to address issues related to agriculture and climate change. Re-analysis of such data available in public domain supplemented with annotations across molecular hierarchy can be of immense help to the plant research community, particularly co-expression networks representing transcriptionally coordinated genes that are often part of the same biological process. With this objective we have developed NetREx, a Network based Rice Expression Analysis Server, that hosts ranked co-expression networks of Oryza sativa using publicly available mRNA-seq data across uniform experimental conditions. It provides a range of interactable data viewers and modules for analysing user queried genes across different stress conditions (drought, flood, cold and osmosis) and hormonal treatments (abscisic and jasmonic acid) and tissues (root and shoot). Subnetworks of user-defined genes can be queried in preconstructed tissue-specific networks, allowing users to view the fold-change, module memberships, gene annotations and analysis of their neighborhood genes and associated pathways. The webserver also allows querying of orthologous genes from Arabidopsis, wheat, maize, barley, and sorghum. Here we demonstrate that NetREx can be used to identify novel candidate genes and tissue-specific interactions under stress conditions and can aid in the analysis and understanding of complex phenotypes linked to stress response in rice. Available at: https://bioinf.iiit.ac.in/netrex/index.html

The Complete Chloroplast Genome Sequences of Eight Fagopyrum Species: Insights Into Genome Evolution and Phylogenetic Relationships

Buckwheat (Fagopyrum genus, Polygonaceae), is an annual or perennial, herbaceous or semi-shrub dicotyledonous plant. There are mainly three cultivated buckwheat species, common buckwheat (Fagopyrum esculentum) is widely cultivated in Asia, Europe, and America, while Tartary buckwheat (F. tataricum) and F. cymosum (also known as F. dibotrys) are mainly cultivated in China. The genus Fagopyrum is taxonomically confusing due to the complex phenotypes of different Fagopyrum species. In this study, the chloroplast (cp) genomes of three Fagopyrum species, F. longistylum, F. leptopodum, F. urophyllum, were sequenced, and five published cp genomes of Fagopyrum were retrieved for comparative analyses. We determined the sequence differentiation, repeated sequences of the cp genomes, and the phylogeny of Fagopyrum species. The eight cp genomes ranged, gene number, gene order, and GC content were presented. Most of variations of Fagopyrum species cp genomes existed in the LSC and SSC regions. Among eight Fagopyrum chloroplast genomes, six variable regions (ndhF-rpl32, trnS-trnG, trnC, trnE-trnT, psbD, and trnV) were detected as promising DNA barcodes. In addition, a total of 66 different SSR (simple sequence repeats) types were found in the eight Fagopyrum species, ranging from 8 to 16 bp. Interestingly, many SSRs showed significant differences especially in some photosystem genes, which provided valuable information for understanding the differences in light adaptation among different Fagopyrum species. Genus Fagopyrum has shown a typical branch that is distinguished from the Rumex, Rheum, and Reynoutria, which supports the unique taxonomic status in Fagopyrum among the Polygonaceae. In addition, phylogenetic analysis based on the cp genomes strongly supported the division of eight Fagopyrum species into two independent evolutionary directions, suggesting that the separation of cymosum group and urophyllum group may be earlier than the flower type differentiation in Fagopyrum plants. The results of the chloroplast-based phylogenetic tree were further supported by the matK and Internal Transcribed Spacer (ITS) sequences of 17 Fagopyrum species, which may help to further anchor the taxonomic status of other members in the urophyllum group in Fagopyrum. This study provides valuable information and high-quality cp genomes for identifying species and evolutionary analysis for future Fagopyrum research.

A single nucleotide variant in the PPARγ-homolog Eip75B affects fecundity in Drosophila

ABSTRACTStill little is understood about the nucleotide changes that underlie variation in complex phenotypes. Variation in the PPARγ-homolog Eip75B has previously been suggested to be associated with longevity and life-history differences in the fruit fly Drosophila melanogaster. Using RNAi knockdown, we first demonstrate that reduced expression of Eip75B in adults affects lifespan, egg laying rate and egg volume. To then test the effect of a naturally occurring SNP variant within a cis-regulatory domain of Eip75B, we screened wildtype lines with alternative alleles and conducted precise genome editing using CRISPR/Cas9. These experiments revealed that this natural polymorphism has a significant effect on fecundity and egg-to-adult viability, but not on longevity or other life-history traits. These results provide a rare functional validation for the role of a natural allelic variant in adaptation of life-history traits directly linked to fitness at the single nucleotide level.

BXD Recombinant Inbred Mice as a Model to Study Neurotoxicity

BXD recombinant inbred (RI) lines represent a genetic reference population derived from a cross between C57BL/6J mice (B6) and DBA/2J mice (D2), which through meiotic recombination events possesses recombinant chromosomes containing B6 or D2 haplotype segments. The quantitative trait loci (QTLs) are the locations of segregating genetic polymorphisms and are fundamental to understanding genetic diversity in human disease susceptibility and severity. QTL mapping represents the typical approach for identifying naturally occurring polymorphisms that influence complex phenotypes. In this process, genotypic values at markers of known genomic locations are associated with phenotypic values measured in a segregating population. Indeed, BXD RI strains provide a powerful tool to study neurotoxicity induced by different substances. In this review, we describe the use of BXD RI lines to understand the underlying mechanisms of neurotoxicity in response to ethanol and cocaine, as well as metals and pesticide exposures.

A single synonymous nucleotide change impacts the male-killing phenotype of prophage WO gene wmk

Wolbachia are the most widespread bacterial endosymbionts in animals. Within arthropods, these maternally-transmitted bacteria can selfishly hijack host reproductive processes to increase the relative fitness of their transmitting females. One such form of reproductive parasitism called male killing, or the selective killing of infected males, is recapitulated to degrees by transgenic expression of the WO-mediated killing (wmk) gene. Here, we characterize the genotype-phenotype landscape of wmk-induced male killing in D. melanogaster using transgenic expression. While phylogenetically distant wmk homologs induce no sex-ratio bias, closely-related homologs exhibit complex phenotypes spanning no death, male death, or death of all hosts. We demonstrate that alternative start codons, synonymous codons, and notably a single synonymous nucleotide in wmk can ablate killing. These findings reveal previously unrecognized features of transgenic wmk-induced killing and establish new hypotheses for the impacts of post-transcriptional processes in male killing variation. We conclude that synonymous sequence changes are not necessarily silent in nested endosymbiotic interactions with life-or-death consequences.

Intergenerational hormesis is regulated by heritable 18S rRNA methylation

SummaryHeritable non-genetic information can regulate a variety of complex phenotypes. However, what specific non-genetic cues are transmitted from parents to their descendants are poorly understood. Here, we perform metabolic methyl-labelling experiments to track the heritable transmission of methylation from ancestors to their descendants in the nematode Caenorhabditis elegans. We find that methylation is transmitted to descendants in proteins, RNA, DNA and lipids. We further find that in response to parental starvation, fed naïve progeny display reduced fertility, increased heat stress resistance, and extended longevity. This intergenerational hormesis is accompanied by a heritable increase in N6’-dimethyl adenosine (m6,2A) on the 18S ribosomal RNA at adenosines 1735 and 1736. We identified the conserved DIMT-1 as the m6,2A methyltransferase in C. elegans and find that dimt-1 is required for the intergenerational hormesis phenotypes. This study provides the first labeling and tracking of heritable non-genetic material across generations and demonstrates the importance of rRNA methylation for regulating the heritable response to starvation.

Neurodevelopmental Disorders in Patients With Complex Phenotypes and Potential Complex Genetic Basis Involving Non-Coding Genes, and Double CNVs

Neurodevelopmental disorders (NDDs) are a heterogeneous class of brain diseases, with a complex genetic basis estimated to account for up to 50% of cases. Nevertheless, genetic diagnostic yield is about 20%. Array-comparative genomic hybridization (array-CGH) is an established first-level diagnostic test able to detect pathogenic copy number variants (CNVs), however, most identified variants remain of uncertain significance (VUS). Failure of interpretation of VUSs may depend on various factors, including complexity of clinical phenotypes and inconsistency of genotype-phenotype correlations. Indeed, although most NDD-associated CNVs are de novo, transmission from unaffected parents to affected children of CNVs with high risk for NDDs has been observed. Moreover, variability of genetic components overlapped by CNVs, such as long non-coding genes, genomic regions with long-range effects, and additive effects of multiple CNVs can make CNV interpretation challenging. We report on 12 patients with complex phenotypes possibly explained by complex genetic mechanisms, including involvement of antisense genes and boundaries of topologically associating domains. Eight among the 12 patients carried two CNVs, either de novo or inherited, respectively, by each of their healthy parents, that could additively contribute to the patients’ phenotype. CNVs overlapped either known NDD-associated or novel candidate genes (PTPRD, BUD13, GLRA3, MIR4465, ABHD4, and WSCD2). Bioinformatic enrichment analyses showed that genes overlapped by the co-occurring CNVs have synergistic roles in biological processes fundamental in neurodevelopment. Double CNVs could concur in producing deleterious effects, according to a two-hit model, thus explaining the patients’ phenotypes and the incomplete penetrance, and variable expressivity, associated with the single variants. Overall, our findings could contribute to the knowledge on clinical and genetic diagnosis of complex forms of NDD.