Signatures of copy number alterations in human cancer.

Gains and losses of DNA are prevalent in cancer and emerge as a consequence of inter-related processes of replication stress, mitotic errors, spindle multipolarity and breakage-fusion-bridge cycles, among others, which may lead to chromosomal instability and aneuploidy1,2. These copy number alterations contribute to cancer initiation, progression and therapeutic resistance3-5. Here we present a conceptual framework to examine the patterns of copy number alterations in human cancer that is widely applicable to diverse data types, including whole-genome sequencing, whole-exome sequencing, reduced representation bisulfite sequencing, single-cell DNA sequencing and SNP6 microarray data. Deploying this framework to 9,873 cancers representing 33 human cancer types from The Cancer Genome Atlas6 revealed a set of 21 copy number signatures that explain the copy number patterns of 97% of samples. Seventeen copy number signatures were attributed to biological phenomena of whole-genome doubling, aneuploidy, loss of heterozygosity, homologous recombination deficiency, chromothripsis and haploidization. The aetiologies of four copy number signatures remain unexplained. Some cancer types harbour amplicon signatures associated with extrachromosomal DNA, disease-specific survival and proto-oncogene gains such as MDM2. In contrast to base-scale mutational signatures, no copy number signature was associated with many known exogenous cancer risk factors. Our results synthesize the global landscape of copy number alterations in human cancer by revealing a diversity of mutational processes that give rise to these alterations.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Locomotion of bovine spermatozoa during the transition from individual cells to bundles.

Various locomotion strategies employed by microorganisms are observed in complex biological environments. Spermatozoa assemble into bundles to improve their swimming efficiency compared to individual cells. However, the dynamic mechanisms for the formation of sperm bundles have not been fully characterized. In this study, we numerically and experimentally investigate the locomotion of spermatozoa during the transition from individual cells to bundles of two cells. Three consecutive dynamic behaviors are found across the course of the transition: hydrodynamic attraction/repulsion, alignment, and synchronization. The hydrodynamic attraction/repulsion depends on the relative orientation and distance between spermatozoa as well as their flagellar wave patterns and phase shift. Once the heads are attached, we find a stable equilibrium of the rotational hydrodynamics resulting in the alignment of the heads. The synchronization results from the combined influence of hydrodynamic and mechanical cell-to-cell interactions. Additionally, we find that the flagellar beat is regulated by the interactions during the bundle formation, whereby spermatozoa can synchronize their beats to enhance their swimming velocity.

Assigning mutational signatures to individual samples and individual somatic mutations with SigProfilerAssignment.

MOTIVATION: Analysis of mutational signatures is a powerful approach for understanding the mutagenic processes that have shaped the evolution of a cancer genome. To evaluate the mutational signatures operative in a cancer genome, one first needs to quantify their activities by estimating the number of mutations imprinted by each signature. RESULTS: Here we present SigProfilerAssignment, a desktop and an online computational framework for assigning all types of mutational signatures to individual samples. SigProfilerAssignment is the first tool that allows both analysis of copy-number signatures and probabilistic assignment of signatures to individual somatic mutations. As its computational engine, the tool uses a custom implementation of the forward stagewise algorithm for sparse regression and nonnegative least squares for numerical optimization. Analysis of 2700 synthetic cancer genomes with and without noise demonstrates that SigProfilerAssignment outperforms four commonly used approaches for assigning mutational signatures. AVAILABILITY AND IMPLEMENTATION: SigProfilerAssignment is available under the BSD 2-clause license at https://github.com/AlexandrovLab/SigProfilerAssignment with a web implementation at https://cancer.sanger.ac.uk/signatures/assignment/.

Insect meal in aquafeeds: A sustainable path to enhanced mucosal immunity in fish.

The mucosal surfaces of fish, including their intestines, gills, and skin, are constantly exposed to various environmental threats, such as water quality fluctuations, pollutants, and pathogens. However, various cells and microbiota closely associated with these surfaces work in tandem to create a functional protective barrier against these conditions. Recent research has shown that incorporating specific feed ingredients into fish diets can significantly boost their mucosal and general immune response. Among the various ingredients being investigated, insect meal has emerged as one of the most promising options, owing to its high protein content and immunomodulatory properties. By positively influencing the structure and function of mucosal surfaces, insect meal (IM) has the potential to enhance the overall immune status of fish. This review provides a comprehensive overview of the potential benefits of incorporating IM into aquafeed as a feed ingredient for augmenting the mucosal immune response of fish.

Understanding the socio-demographic and programmatic factors associated with adolescent motherhood and its association with child undernutrition in Bangladesh.

BACKGROUND: Worldwide, a significant number of girls become mothers during adolescence. In Bangladesh, adolescent childbirth is highly prevalent and has adverse effects on children's health and undernutrition. We aimed to identify the relationship between the undernutrition of children and adolescent motherhood, the factors associated with adolescent mothers' age at first birth, and to examine the programmatic factors and gaps influencing children's undernutrition in Bangladesh. METHODS: We analysed the 'Bangladesh Demographic and Health Survey' BDHS-17-18 data and desk review. To examine the factors associated with adolescent motherhood and its impact on child undernutrition, data from 7,643 mother-child pairs were selected. Child stunting, wasting, and underweight were measured according to the World Health Organisation (WHO) median growth guidelines based on z-scores - 2. Univariate, bivariate, simple, and multiple logistic regressions were used for analyse. We followed the systematic procedures for the literature review. RESULTS: Approximately, 89% of adolescents aged ≤ 19 years were married and 71% of them gave their first childbirth. Children of adolescent mothers (≤ 19 years) were significantly 1.68 times more wasted (aOR: 1.68; 95% CI: 1.08 to 2.64), 1.37 times more underweight (aOR: 1.37; 95% CI: 1.01 to 1.86) and either form 1.32 times more stunting, wasting or underweight (aOR:1.32; 95% Cl: 1.05 to 1.66) compared to the children of adult mothers (> 19 years) after adjusting potential confounders. The factors associated with mothers' first childbirth during adolescence were the age gap between husband and wife 5-10 years (aOR: 1.81; 95% Cl: 1.57-2.10) and age gap > 10 years (aOR: 2.41; 95% Cl: 1.96-2.97) compared with the age group < 5 years, and husbands' education (aOR: 1.29; 95% Cl: 1.04-1.61) compared with the uneducated husbands. In the literature review, we found potential gaps in focusing on the Adolescent Sexual and Reproductive Health (ASRH) program in Bangladesh, from thirty-two programmes only half of them focused on adolescents aged 10-19 years, and eleven programmes focused only on girls. CONCLUSION: Children of adolescent mothers are at risk of wasting, underweight, and any form of undernutrition. For effective policies and interventions in Bangladesh, it is important to emphasise delaying adolescent pregnancy and prioritising child undernutrition.

Elevated temperature impairs gonadal development by suppressing the expression of the genes for kisspeptin, GnRH1 and GTH subunits in Nile tilapia Oreochromis niloticus.

Temperature is a preeminent factor in the regulation of fish reproduction and hinders gonadal development beyond a specific threshold. To comprehend the molecular mechanism responsible for reproductive suppression at different temperature, expression of the genes encoding kisspeptin (kiss2), gonadotropin-releasing hormone (gnrh1) and their receptors (gpr54, gnrh1r) in the brain, and the gonadotropin (GTH) subunits (fshb and lhb) in the pituitary were studied in juvenile Nile tilapia (Oreochromis niloticus) along with gonadal histology. Fish were acclimatized to three distinct temperatures, including 31 °C, 34 °C and 37 °C for 14 days. The mRNA levels of kiss2, gpr54, gnrh1, and gnrh1r were significantly decreased at 37 °C compared to 31 °C and 34 °C in the both sexes. In parallel, the expression level of fshb in the both sexes and lhb in the female were significantly lower at 37 °C in the pituitary. Histologically, the gonads of both sexes had normal growth of gametes at control temperature (31 °C), whereas the spermatogenesis and oocyte maturation were slowed down and atretic oocytes were found in the ovary at 37 °C acclimation temperature. Taken together, the results imply that elevated temperature beyond the specific threshold may have a negative impact on reproduction by suppressing the gene expressions of kisspeptin/GnRH1/GTH system and eventually restrains normal growth and maturation of gametes in the both sexes of Nile tilapia.

Hypoxia alters the upper thermal limits and blood physiology in zebrafish, Danio rerio.

Hypoxic aquatic environments occur more frequently as a result of climate change, thereby exerting challenges on the physiological and metabolic functions of aquatic animals. In this study, a model fish, zebrafish (Danio rerio) was used to observe the climate-induced hypoxic effect on the upper thermal limit (critical thermal maximum; CTmax), hemoglobin, and blood glucose levels, and abnormalities of erythrocytes at cellular and nuclear level. The value of CTmax decreased significantly under hypoxia (39.10 ± 0.96 °C) compared to normoxia (43.70 ± 0.91 °C). At CTmax, hemoglobin levels were much lower (9.33 ± 0.60 g/dL) and blood glucose levels were significantly higher (194.20 ± 11.33 mg/L) under hypoxia than they were under normoxia and at the beginning of the experiment. Increased frequencies of abnormalities in the erythrocytes at both cellular (fusion, twin, elongated, spindle and tear drop shaped) and nuclear (micronucleus, karyopyknosis, binuclei, nuclear degeneration and notched nuclei) levels were also found under hypoxia compared to normoxia. These results suggest that hypoxic conditions significantly alter the temperature tolerance and subsequent physiology in zebrafish. Our findings will aid in the development of effective management techniques for aquatic environments with minimum oxygen availability.

Phenotypic and microarray analysis reveals salinity stress-induced oxidative tolerance in transgenic rice expressing a DEAD-box RNA helicase, OsDB10.

Helicases are the motor proteins not only involved in the process of mRNA metabolism but also played a significant role in providing abiotic stresses tolerance. In this study, a DEAD-box RNA helicase OsDB10 was cloned and functionally characterized. The transcript levels of OsDB10 were increased both in shoot and root upon salt, heat, cold, and ABA application and was more prominent in shoot compared to root. Genomic integration of OsDB10 in transgenic rice was confirmed by PCR, Southern blot and qRT-PCR analysis. The transgenic plants showed quicker seed germination, reduced necrosis, higher chlorophyll, more survival rate, better seedling growth, and produced more grain yield under salinity stress. Furthermore, transgenic lines also accumulated less Na+ and high K+ ions and salinity tolerance of the transgenic were also assayed by measuring different bio-physiological indices. Moreover, the OsDB10 transgenic plants showed enhanced tolerance to salinity-induced oxidative stress by scavenging ROS and increased activity of antioxidants enzymes. Microarray analysis showed upregulation of transcriptional regulations and metabolic reprogramming as OsDB10 overexpression modulates the expression of many other genes. Altogether, our results confirmed that OsDB10 is a functional DEAD-box RNA helicase and played vital roles in plant defence response against salinity stress.

Characterization of immune responses associated with ERAP-1 expression in HSV-induced Behçet's disease mouse model.

Behçet's disease (BD) is a chronic multisystem inflammatory disorder. Endoplasmic reticulum aminopeptidase 1 (ERAP1) polymorphism has been reported as a risk factor for BD. However, the immunological role of ERAP1 in BD remains unclear. Therefore, the purpose of this study was to investigate the immunological role of ERAP1 in BD using a mouse model. ERAP1 incomplete expressing mice (ERAP1 hetero, +/-) were generated and inoculated with herpes simplex virus 1 to produce a BD mouse model. In these mice, dendritic cell activation markers and other immune response-related markers were analyzed. Among them, the factor showing a significant difference between ERAP1+/- BD mice and WT BD mice was IL-17. In ERAP1+/-, BD had significantly different expression levels of CD80, CD11b, Ly6G, RORγt, IFNγ, and IL-17 compared to asymptomatic controls. This study demonstrates ERAP1 defective expressions play an important role in BD development through inappropriate regulation of Th17.

SARS-CoV-2 infection alters the gut microbiome in diabetes patients: A cross-sectional study from Bangladesh.

Populations of different South Asian nations including Bangladesh reportedly have a high risk of developing diabetes in recent years. This study aimed to investigate the differences in the gut microbiome of COVID-19-positive participants with or without type 2 diabetes mellitus (T2DM) compared with healthy control subjects. Microbiome data of 30 participants with T2DM were compared with 22 age-, sex-, and body mass index (BMI)-matched individuals. Clinical features were recorded while fecal samples were collected aseptically from the participants. Amplicon-based (16S rRNA) metagenome analyses were employed to explore the dysbiosis of gut microbiota and its correlation with genomic and functional features in COVID-19 patients with or without T2DM. Comparing the detected bacterial genera across the sample groups, 98 unique genera were identified, of which 9 genera had unique association with COVID-19 T2DM patients. Among different bacterial groups, Shigella (25%), Bacteroides (23.45%), and Megamonas (15.90%) had higher mean relative abundances in COVID-19 patients with T2DM. An elevated gut microbiota dysbiosis in T2DM patients with COVID-19 was observed while some metabolic functional changes correlated with bidirectional microbiome dysbiosis between diabetes and non-diabetes humans gut were also found. These results further highlight the possible association of COVID-19 infection that might be linked with alteration of gut microbiome among T2DM patients.

Microbiome pattern and diversity of an anadromous fish, hilsa shad (Tenualosa ilisha).

BACKGROUND: The host-microbe interactions are complex, dynamic and context-dependent. In this regard, migratory fish species like hilsa shad (Tenualosa ilisha), which migrates from seawater to freshwater for spawning, provides a unique system for investigating the microbiome under an additional change in fish's habitat. This work was undertaken to detect taxonomic variation of microbiome and their function in the migration of hilsa. METHODS AND RESULTS: The study employed 16S rRNA amplicon-based metagenomic analysis to scrutinize bacterial diversity in hilsa gut, skin mucus and water. Thus, a total of 284 operational taxonomic units (OTUs), 9 phyla, 35 orders and 121 genera were identified in all samples. More than 60% of the identified bacteria were Proteobacteria with modest abundance (> 5%) of Firmicutes, Bacteroidetes and Actinobacteria. Leucobacter in gut and Serratia in skin mucus were the core bacterial genera, while Acinetobacter, Pseudomonas and Psychrobacter exhibited differential compositions in gut, skin mucus and water. CONCLUSIONS: Representative fresh-, brackish- and seawater samples of hilsa habitats were primarily composed of Vibrio, Serratia and Psychrobacter, and their diversity in seawater was significantly higher (P < 0.05) than freshwater. Overall, salinity and water microbiota had an influence on the microbial composition of hilsa shad, contributing to host metabolism and adaptation processes. This pioneer exploration of hilsa gut and skin mucus bacteria across habitats will advance our insights into microbiome assembly in migratory fish populations.

An implicit electrolyte model for plane wave density functional theory exhibiting nonlinear response and a nonlocal cavity definition.

We have developed and implemented an implicit electrolyte model in the Vienna Ab initio Simulation Package (VASP) that includes nonlinear dielectric and ionic responses as well as a nonlocal definition of the cavities defining the spatial regions where these responses can occur. The implementation into the existing VASPsol code is numerically efficient and exhibits robust convergence, requiring computational effort only slightly higher than the original linear polarizable continuum model. The nonlinear + nonlocal model is able to reproduce the characteristic "double hump" shape observed experimentally for the differential capacitance of an electrified metal interface while preventing "leakage" of the electrolyte into regions of space too small to contain a single water molecule or solvated ion. The model also gives a reasonable prediction of molecular solvation free energies as well as the self-ionization free energy of water and the absolute electron chemical potential of the standard hydrogen electrode. All of this, combined with the additional ability to run constant potential density functional theory calculations, should enable the routine computation of activation barriers for electrocatalytic processes.

Open-label, randomised controlled trial found that a green banana mixed rice suji diet was most effective for persistent diarrhoea in children in Bangladesh.

AIM: Persistent diarrhoea continues for at least 14 days and kills more children than acute diarrhoea. We assessed whether rice suji, green banana mixed rice suji or 75% rice suji improved persistent diarrhoea in young children. METHODS: This open-labelled randomised controlled trial was carried out between December 2017 and August 2019 at the Dhaka Hospital of icddr,b, Bangladesh, with 135 children aged 6-35 months with persistent diarrhoea. The children were randomly assigned to green banana mixed rice suji, rice suji or 75% rice suji, with 45 in each group. The primary outcome was the percentage who recovered from diarrhoea by day 5 using an intention-to-treat analysis. RESULTS: The children's median age was 8 months (interquartile range: 7-10 months). By day 5, the recovery rate was 58%, 31% and 58% for children in the green banana mixed rice suji, rice suji and 75% rice suji groups, respectively. The green banana mixed rice suji group had fewer relapses (7%) than the 75% rice suji group (24%). Enteroaggregative Escherichia coli, rotavirus, norovirus, Enteropathogenic Escherichia coli, astrovirus and Campylobacter were the major pathogens for persistent diarrhoea. CONCLUSION: Green banana mixed rice suji was the most effective option for managing persistent diarrhoea in young children.

CD83 Regulates the Immune Responses in Inflammatory Disorders.

Activating the immune system plays an important role in maintaining physiological homeostasis and defending the body against harmful infections. However, abnormalities in the immune response can lead to various immunopathological responses and severe inflammation. The activation of dendritic cells (DCs) can influence immunological responses by promoting the differentiation of T cells into various functional subtypes crucial for the eradication of pathogens. CD83 is a molecule known to be expressed on mature DCs, activated B cells, and T cells. Two isotypes of CD83, a membrane-bound form and a soluble form, are subjects of extensive scientific research. It has been suggested that CD83 is not only a ubiquitous co-stimulatory molecule but also a crucial player in monitoring and resolving inflammatory reactions. Although CD83 has been involved in immunological responses, its functions in autoimmune diseases and effects on pathogen immune evasion remain unclear. Herein, we outline current immunological findings and the proposed function of CD83 in inflammatory disorders.

Mitigating greenhouse gas emissions from irrigated rice cultivation through improved fertilizer and water management.

Greenhouse gas (GHG) emissions from agriculture sector play an important role for global warming and climate change. Thus, it is necessary to find out GHG emissions mitigation strategies from rice cultivation. The efficient management of nitrogen fertilizer using urea deep placement (UDP) and the use of the water-saving alternate wetting and drying (AWD) irrigation could mitigate greenhouse gas (GHG) emissions and reduce environmental pollution. However, there is a dearth of studies on the impacts of UDP and the integrated plant nutrient system (IPNS) which combines poultry manure and prilled urea (PU) with different irrigation regimes on GHG emissions, nitrogen use efficiency (NUE) and rice yields. We conducted field experiments during the dry seasons of 2018, 2019, and 2020 to compare the effects of four fertilizer treatments including control (no N), PU, UDP, and IPNS in combination with two irrigation systems- (AWD and continuous flooding, CF) on GHG emissions, NUE and rice yield. Fertilizer treatments had significant (p < 0.05) interaction effects with irrigation regimes on methane (CH4) and nitrous oxide (N2O) emissions. PU reduced CH4 and N2O emissions by 6% and 20% compared to IPNS treatment, respectively under AWD irrigation, but produced similar emissions under CF irrigation. Similarly, UDP reduced cumulative CH4 emissions by 9% and 15% under AWD irrigation, and 9% and 11% under CF condition compared to PU and IPNS treatments, respectively. Across the year and fertilizer treatments, AWD irrigation significantly (p < 0.05) reduced cumulative CH4 emissions and GHG intensity by 28%, and 26%, respectively without significant yield loss compared to CF condition. Although AWD irrigation increased cumulative N2O emissions by 73%, it reduced the total global warming potential by 27% compared to CF irrigation. The CH4 emission factor for AWD was lower (1.67 kg ha-1 day-1) compared to CF (2.33 kg ha-1 day-1). Across the irrigation regimes, UDP increased rice yield by 21% and N recovery efficiency by 58% compared to PU. These results suggest that both UDP and AWD irrigation might be considered as a carbon-friendly technology.

Hypoxia-inducible factor-1 drives divergent immunomodulatory functions in the pathogenesis of autoimmune diseases.

Hypoxia-inducible factor-1 (HIF-1) is a heterodimeric (HIF-1α/ HIF-1ß) transcription factor in which the oxygen-sensitive HIF-1α subunit regulates gene transcription to mediate adaptive tissue responses to hypoxia. HIF-1 is a key mediator in both regulatory and pathogenic immune responses, because ongoing inflammation in localized tissues causes increased oxygen consumption and consequent hypoxia within the inflammatory lesions. In autoimmune diseases, HIF-1 plays complex and divergent roles within localized inflammatory lesions by orchestrating a critical immune interplay sponsoring the pathogenesis of the disease. In this review, we have summarized the role of HIF-1 in lymphoid and myeloid immunomodulation in autoimmune diseases. HIF-1 drives inflammation by controlling the Th17/Treg /Tr1 balance through the tipping of the differentiation of CD4+ T cells in favour of pro-inflammatory Th17 cells while suppressing the development of anti-inflammatory Treg /Tr1 cells. On the other hand, HIF-1 plays a protective role by facilitating the expression of anti-inflammatory cytokine IL-10 in and expansion of CD1dhi CD5+ B cells, known as regulatory B cells or B10 cells. Apart from lymphoid cells, HIF-1 also controls the activation of macrophages, neutrophils and dendritic cells, thus eventually further influences the activation and development of effector/regulatory T cells by facilitating the creation of a pro/anti-inflammatory microenvironment within the autoinflammatory lesions. Based on the critical immunomodulatory roles that HIF-1 plays, this master transcription factor seems to be a potent druggable target for the treatment of autoimmune diseases.

Vascular and immunopathological role of Asymmetric Dimethylarginine (ADMA) in Experimental Autoimmune Encephalomyelitis.

Asymmetric dimethylarginine (ADMA) is an endogenous nitric oxide synthase (NOS) inhibitor/uncoupler inducing vascular pathology. Vascular pathology is an important factor for the development and progression of CNS pathology of MS, yet the role of ADMA in MS remains elusive. Patients with multiple sclerosis (MS) are reported to have elevated blood levels of ADMA, and mice with experimental autoimmune encephalomyelitis (EAE, an animal model of MS) generated by auto-immunization of myelin oligodendrocyte glycoprotein (MOG) and blood-brain barrier (BBB) disruption by pertussis toxin also had increased blood ADMA levels in parallel with induction of clinical disease. To explore the role of ADMA in EAE pathogenesis, EAE mice were treated with a daily dose of ADMA. It is of special interest that ADMA treatment enhanced the BBB disruption in EAE mice and exacerbated the clinical and CNS disease of EAE. ADMA treatment also induced the BBB disruption and EAE disease in MOG-immunized mice even without pertussis toxin treatment, suggesting the role of ADMA in BBB dysfunction in EAE. T-cell polarization studies also documented that ADMA treatment promotes TH 1- and TH 17-mediated immune responses but without affecting Treg-mediated immune response in EAE mice as well as in in vitro T-cell culture. Taken together, these data, for the first time, document the vascular and immunopathogenic roles of ADMA in EAE, thus pointing to the potential of ADMA-mediated mechanism as a new target of potential therapy for MS.

Robust regression based genome-wide multi-trait QTL analysis.

In genome-wide quantitative trait locus (QTL) mapping studies, multiple quantitative traits are often measured along with the marker genotypes. Multi-trait QTL (MtQTL) analysis, which includes multiple quantitative traits together in a single model, is an efficient technique to increase the power of QTL identification. The two most widely used classical approaches for MtQTL mapping are Gaussian Mixture Model-based MtQTL (GMM-MtQTL) and Linear Regression Model-based MtQTL (LRM-MtQTL) analyses. There are two types of LRM-MtQTL approach known as least squares-based LRM-MtQTL (LS-LRM-MtQTL) and maximum likelihood-based LRM-MtQTL (ML-LRM-MtQTL). These three classical approaches are equivalent alternatives for QTL detection, but ML-LRM-MtQTL is computationally faster than GMM-MtQTL and LS-LRM-MtQTL. However, one major limitation common to all the above classical approaches is that they are very sensitive to outliers, which leads to misleading results. Therefore, in this study, we developed an LRM-based robust MtQTL approach, called LRM-RobMtQTL, for the backcross population based on the robust estimation of regression parameters by maximizing the ß-likelihood function induced from the ß-divergence with multivariate normal distribution. When ß = 0, the proposed LRM-RobMtQTL method reduces to the classical ML-LRM-MtQTL approach. Simulation studies showed that both ML-LRM-MtQTL and LRM-RobMtQTL methods identified the same QTL positions in the absence of outliers. However, in the presence of outliers, only the proposed method was able to identify all the true QTL positions. Real data analysis results revealed that in the presence of outliers only our LRM-RobMtQTL approach can identify all the QTL positions as those identified in the absence of outliers by both methods. We conclude that our proposed LRM-RobMtQTL analysis approach outperforms the classical MtQTL analysis methods.