Two rare autosomal recessive neurological disorders identified by combined genetic approaches in a single consanguineous family with multiple offspring

INTRODUCTION

Neurodevelopmental disorders (NDDs) refer to a broad range of diseases including developmental delay, intellectual disability, epilepsy, autism spectrum disorders, and attention-deficit/hyperactivity disorder caused by dysfunctions in tightly controlled brain development. The genetic backgrounds of NDDs are quite heterogeneous; to date, recessive or dominant variations in numerous genes have been implicated. Herein, we present a large consanguineous family from Turkiye, who has been suffering from NDDs with two distinct clinical presentations.

METHODS AND RESULTS

Combined in-depth genetic approaches led us to identify a homozygous frameshift variant in NALCN related to NDD and expansion of dodecamer repeat in CSTB related to Unverricht-Lundborg disease (ULD). Additionally, we sought to functionally analyze the NALCN variant in terms of mRNA expression level and current alteration. We have both detected a decrease in the level of premature stop codon-bearing mRNA possibly through nonsense-mediated mRNA decay mechanism and also an increased current in patch-clamp recordings for the expressed truncated protein.

CONCLUSION

In conclusion, increased consanguinity may lead to the revealing of distinct rare neurogenetic diseases in a single family. Exome sequencing is generally considered the first-tier diagnostic test in individuals with NDD. Yet we underline the fact that customized approaches other than exome sequencing may be used as in the case of ULD to aid diagnosis and better genetic counseling.

Did the human genome project affect research on Schizophrenia?

The Human Genome Project was undertaken primarily to discover genetic causes and better treatments for human diseases. Schizophrenia was targeted since three of the project`s principal architects had a personal interest and also because, based on family, adoption, and twin studies, schizophrenia was widely believed to be a genetic disorder. Extensive studies using linkage analysis, candidate genes, genome wide association studies [GWAS], copy number variants, exome sequencing and other approaches have failed to identify causal genes. Instead, they identified almost 300 single nucleotide polymorphisms [SNPs] associated with altered risks of developing schizophrenia as well as some rare variants associated with increased risk in a small number of individuals. Risk genes play a role in the clinical expression of most diseases but do not cause the disease in the absence of other factors. Increasingly, observers question whether schizophrenia is strictly a genetic disorder. Beginning in 1996 NIMH began shifting its research resources from clinical studies to basic research based on the promise of the Human Genome Project. Consequently, three decades later NIMH's genetics investment has yielded almost nothing clinically useful for individuals currently affected. It is time to review NIMH`s schizophrenia research portfolio.

Joint Linkage and Association Analysis Using GENEHUNTER-MODSCORE with an Application to Familial Pancreatic Cancer

INTRODUCTION

Joint linkage and association (JLA) analysis combines two disease gene mapping strategies: linkage information contained in families and association information contained in populations. Such a JLA analysis can increase mapping power, especially when the evidence for both linkage and association is low to moderate. Similarly, an association analysis based on haplotypes instead of single markers can increase mapping power when the association pattern is complex.

METHODS

In this paper, we present an extension to the GENEHUNTER-MODSCORE software package that enables a JLA analysis based on haplotypes and uses information from arbitrary pedigree types and unrelated individuals. Our new JLA method is an extension of the MOD score approach for linkage analysis, which allows the estimation of trait-model and linkage disequilibrium (LD) parameters, i.e., penetrance, disease-allele frequency, and haplotype frequencies. LD is modeled between alleles at a single diallelic disease locus and up to three diallelic test markers. Linkage information is contributed by additional multi-allelic flanking markers. We investigated the statistical properties of our JLA implementation using extensive simulations, and we compared our approach to another commonly used single-marker JLA test. To demonstrate the applicability of our new method in practice, we analyzed pedigree data from the German National Case Collection for Familial Pancreatic Cancer (FaPaCa).

RESULTS

Based on the simulated data, we demonstrated the validity of our JLA-MOD score analysis implementation and identified scenarios in which haplotype-based tests outperformed the single-marker test. The estimated trait-model and LD parameters were in good accordance with the simulated values. Our method outperformed another commonly used JLA single-marker test when the LD pattern was complex. The exploratory analysis of the FaPaCa families led to the identification of a promising genetic region on chromosome 22q13.33, which can serve as a starting point for future mutation analysis and molecular research in pancreatic cancer.

CONCLUSION

Our newly proposed JLA-MOD score method proves to be a valuable gene mapping and characterization tool, especially when either linkage or association information alone provide insufficient power to identify the disease-causing genetic variants.

Genetic Map Construction Using F2 and RIL/DH Mapping Population

Genetic mapping is the determination of the position and relative genetic distance between genes or molecular markers in the chromosomes of a particular species. The construction of genetic maps uses data from the genotyping of the mapping population. Among the different mapping populations used, two are relatively common: the F2 and recombinant inbred lines (RILs) obtained as a result of the controlled crossing of genetically diverse parental forms (e.g., inbred lines). Also, the dihaploid (DH) population is often used in plants, but obtaining DHs in different crops, including rye, is very difficult or even impossible. Any molecular marker system can be used for genotyping. Polymorphic markers are used for linkage analysis, differentiating parental forms with segregation in the mapping population, consistent with the appropriate single-gene model. A genetic map is a great source of information on a species and can be an exquisite tool for analyzing important quantitative traits (QT).This chapter presents the procedure of genetic map construction with two different algorithms using the JoinMap5.0 program. First, the Materials section briefly informs about the mapping program, showing how to obtain a mapping population and prepare data for mapping. Finally, the Methods section describes the protocol for the mapping procedure itself.

Improved haplotype resolution of highly duplicated MHC genes in a long-read genome assembly using MiSeq amplicons

Long-read sequencing offers a great improvement in the assembly of complex genomic regions, such as the major histocompatibility complex (MHC) region, which can contain both tandemly duplicated MHC genes (paralogs) and high repeat content. The MHC genes have expanded in passerine birds, resulting in numerous MHC paralogs, with relatively high sequence similarity, making the assembly of the MHC region challenging even with long-read sequencing. In addition, MHC genes show rather high sequence divergence between alleles, making diploid-aware assemblers incorrectly classify haplotypes from the same locus as sequences originating from different genomic regions. Consequently, the number of MHC paralogs can easily be over- or underestimated in long-read assemblies. We therefore set out to verify the MHC diversity in an original and a haplotype-purged long-read assembly of one great reed warbler Acrocephalus arundinaceus individual (the focal individual) by using Illumina MiSeq amplicon sequencing. Single exons, representing MHC class I (MHC-I) and class IIB (MHC-IIB) alleles, were sequenced in the focal individual and mapped to the annotated MHC alleles in the original long-read genome assembly. Eighty-four percent of the annotated MHC-I alleles in the original long-read genome assembly were detected using 55% of the amplicon alleles and likewise, 78% of the annotated MHC-IIB alleles were detected using 61% of the amplicon alleles, indicating an incomplete annotation of MHC genes. In the haploid genome assembly, each MHC-IIB gene should be represented by one allele. The parental origin of the MHC-IIB amplicon alleles in the focal individual was determined by sequencing MHC-IIB in its parents. Two of five larger scaffolds, containing 6-19 MHC-IIB paralogs, had a maternal and paternal origin, respectively, as well as a high nucleotide similarity, which suggests that these scaffolds had been incorrectly assigned as belonging to different loci in the genome rather than as alternate haplotypes of the same locus. Therefore, the number of MHC-IIB paralogs was overestimated in the haploid genome assembly. Based on our findings we propose amplicon sequencing as a suitable complement to long-read sequencing for independent validation of the number of paralogs in general and for haplotype inference in multigene families in particular.

Selection for, and characterization of, fluralaner resistance in the house fly, Musca domestica

House flies, Musca domestica (L), are the mechanical vector of >100 human and animal pathogens, including those that are antibiotic-resistant. Given that house flies are associated closely with human and livestock activity, they present medical and veterinary health concerns. Although there are numerous strategies for control of house fly populations, chemical control has been favored in many facilities. Products with pyrethroid active ingredients have been used predominantly for >35 years in space sprays. As a result, strong selection for pyrethroid resistance has led to reduced control of many populations. Reliance on a limited number of insecticides for decades has created fly control problems necessitating the discovery and formulation of new control chemistries. Fluralaner is a relatively new insecticide and acaricide (first reported in 2010), belonging to the isoxazoline class. These insecticides target the glutamate- and gamma-aminobutyric acid-gated (GABA) chloride channels, which is a different mode of action from other insecticides used against house flies. Although is it not currently registered for house fly control in the United States, previous work has shown that fluralaner is highly toxic to house flies and that there was limited cross-resistance found in laboratory strains having high levels of resistance to other insecticides. Herein, we characterized the time and age dependency of fluralaner toxicity, detected cross-resistance in populations from across the United States, and selected a highly resistant (>11,000-fold) house fly strain. We found that the fluralaner LD₅₀ of 18-24 h old flies was 2-fold higher than for 5-6 d old flies. This appears to be due to more rapid penetration of fluralaner into the 5-6 d old flies. Fluralaner resistance was inherited as an intermediate to incompletely dominant trait and was mapped to chromosomes 5 and 3. Resistance could be suppressed to 7-fold with piperonyl butoxide, suggesting that cytochrome P450 (CYP)-mediated detoxification was a major mechanism of resistance. Decreased penetration was also demonstrated as a mechanism of resistance. The utility of fluralaner for house fly control is discussed.

Reanalysis of exome sequencing data reveals a treatable neurometabolic origin in two previously undiagnosed siblings with neurodevelopmental disorder

Neurodevelopmental disorders (NDDs) have broad heterogeneity both clinically and genetically. Inborn errors of metabolism can be one of the reasons of neurodevelopmental disruption causing specific NDDs. Although there is tremendous advance in molecular identification via next-generation sequencing (NGS), there are still many unsolved patients with NDD. Reanalysis of NGS data with different pipelines can at least partially accomplish this challenge. Herein, we report clinic and genetic components of an adult sib-pair with an undiagnosed NDD condition, which has been solved through reanalysis of whole-exome sequencing (WES). Parallel analysis of SNP-based genotyping and WES was performed to focus on variants only in loci with positive logarithm of the odds scores. WES data was analyzed through three different pipelines with two distinct bed files. Reanalysis of WES data led us to detect a homozygous FOLR1 variant (ENST00000393676.5:c.610C > T, p.(Arg204Ter), rs952165627) in the affected sib-pair. Surprisingly, the variant could not be detected in the first analysis as the variant region is not included in the first bed file which may frequently be used. Biochemical tests of CSF have confirmed the genetic analysis, CSF folic acid levels were detected low in sib-pair, and intravenous folinic acid treatment improved the disease course for the first 6 months of follow-up even at late diagnosis age. Although combined analysis of SNP-based genotyping and WES is a powerful tool to reveal the genetic components of heterogeneous diseases, reanalysis of genome data still should be considered in unsolved patients. Also, biochemical screening helps us to decipher undiagnosed NDD that may be a treatable neurometabolic condition.

Adaptive tracking algorithm for trajectory analysis of cells and layer-by-layer assessment of motility dynamics

Tracking biological objects such as cells or subcellular components imaged with time-lapse microscopy enables us to understand the molecular principles about the dynamics of cell behaviors. However, automatic object detection, segmentation and extracting trajectories remain as a rate-limiting step due to intrinsic challenges of video processing. This paper presents an adaptive tracking algorithm (Adtari) that automatically finds the optimum search radius and cell linkages to determine trajectories in consecutive frames. A critical assumption in most tracking studies is that displacement remains unchanged throughout the movie and cells in a few frames are usually analyzed to determine its magnitude. Tracking errors and inaccurate association of cells may occur if the user does not correctly evaluate the value or prior knowledge is not present on cell movement. The key novelty of our method is that minimum intercellular distance and maximum displacement of cells between frames are dynamically computed and used to determine the threshold distance. Since the space between cells is highly variable in a given frame, our software recursively alters the magnitude to determine all plausible matches in the trajectory analysis. Our method therefore eliminates a major preprocessing step where a constant distance was used to determine the neighbor cells in tracking methods. Cells having multiple overlaps and splitting events were further evaluated by using the shape attributes including perimeter, area, ellipticity and distance. The features were applied to determine the closest matches by minimizing the difference in their magnitudes. Finally, reporting section of our software were used to generate instant maps by overlaying cell features and trajectories. Adtari was validated by using videos with variable signal-to-noise, contrast ratio and cell density. We compared the adaptive tracking with constant distance and other methods to evaluate performance and its efficiency. Our algorithm yields reduced mismatch ratio, increased ratio of whole cell track, higher frame tracking efficiency and allows layer-by-layer assessment of motility to characterize single-cells. Adaptive tracking provides a reliable, accurate, time efficient and user-friendly open source software that is well suited for analysis of 2D fluorescence microscopy video datasets.

Genome-wide linkage search for cancer susceptibility loci in a cohort of non BRCA1/2 families in Sri Lanka

Objective

Although linkage studies have been utilized for the identification of variants associated with cancer in the world, little is known about their role in non BRCA1/2 individuals in the Sri Lankans. Hence we performed linkage analysis to identify susceptibility loci related to the inherited risk of cancer in a cohort of Sri Lankans affected with hereditary breast cancer. The Illumina global screening array having 654,027 single nucleotide polymorphism markers was performed in four families, in which at least three individuals within third degree relatives were affected by breast cancer. Two-point parametric linkage analysis was conducted assuming disease allele frequency of 1%. Penetrance was set at 90% for carriers with a 10% phenocopy rate.

Results

Thirty-one variants exhibited genome-wide suggestive HLODs. The top overall HLOD score was at rs1856277, an intronic variant in MYO16 on chromosome 13. The two most informative families also suggested several candidate linked loci in genes, including ERAP1, RPRM, WWOX, CDH1, EXOC1, HUS1B, STIM1 and TUSC1. This study provides the first step in identifying germline variants that may be involved in risk of cancer in cancer-aggregated non-BRCA1/2 families from the understudied Sri Lankan population. Several candidate linked regions showed suggestive evidence of linkage to cancer risk.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-06081-5.

High-resolution detection of quantitative trait loci for seven important yield-related traits in wheat (Triticum aestivum L.) using a high-density SLAF-seq genetic map

BACKGROUND

Yield-related traits including thousand grain weight (TGW), grain number per spike (GNS), grain width (GW), grain length (GL), plant height (PH), spike length (SL), and spikelet number per spike (SNS) are greatly associated with grain yield of wheat (Triticum aestivum L.). To detect quantitative trait loci (QTL) associated with them, 193 recombinant inbred lines derived from two elite winter wheat varieties Chuanmai42 and Chuanmai39 were employed to perform QTL mapping in six/eight environments.

RESULTS

A total of 30 QTLs on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 4A, 5A, 5B, 6A, 6D, 7A, 7B and 7D were identified. Among them, six major QTLs QTgw.cib-6A.1, QTgw.cib-6A.2, QGw.cib-6A, QGl.cib-3A, QGl.cib-6A, and QSl.cib-2D explaining 5.96-23.75% of the phenotypic variance were detected in multi-environments and showed strong and stable effects on corresponding traits. Three QTL clusters on chromosomes 2D and 6A containing 10 QTLs were also detected, which showed significant pleiotropic effects on multiple traits. Additionally, three Kompetitive Allele Specific PCR (KASP) markers linked with five of these major QTLs were developed. Candidate genes of QTgw.cib-6A.1/QGl.cib-6A and QGl.cib-3A were analyzed based on the spatiotemporal expression patterns, gene annotation, and orthologous search.

CONCLUSIONS

Six major QTLs for TGW, GL, GW and SL were detected. Three KASP markers linked with five of these major QTLs were developed. These QTLs and KASP markers will be useful for elucidating the genetic architecture of grain yield and developing new wheat varieties with high and stable yield in wheat.

QTL mapping for adult plant resistance to wheat stripe rust in M96-5 × Guixie 3 wheat population

Development of cultivars with multiple resistances has proven to be an effective way to prevent diseases in wheat breeding. The Guixie 3 variety (GX3) has shown excellent performance in resistance to stripe rust in field for many years. The purpose of this study was to detect quantitative trait loci (QTL) associated with resistance to stripe rust in the adult plant stage and determine closely linked molecular markers. A population of recombinant inbred lines (n = 228) was derived from a cross between the susceptible landrace Mian 96-5 (M96-5) and GX3 variety and evaluated in multiple field studies, and QTL analysis enabled to elucidate genetic architecture of wheat resistance to stripe rust. A total of 19 QTL for stripe rust resistance were mapped on 12 chromosomes using phenotypic data from multiple field tests over the course of 6 years. These chromosomes included 1B (2), 1D (2), 2A (2), 2B (2), 2D (1), 4B (2), 4D (1), 5A (3), 5B (1), 6A (1), 6B (1), and 7B (1). Two stable QTL on chromosomes 2AS (Qyr.gaas.2A) and 6AL (Qyr.gaas.6A) were detected in six and five different environments, respectively; in both QTL, positive allele was contributed by GX3 variety. Qyr.gaas.2A was found to be crucial for increasing adult plant resistance, which may explain the large phenotypic variation of 45.52%. Our results provide theoretical and molecular insight for wheat breeding and suggest the cloning of genes associated with the GX3 variety may be beneficial in future studies.

Gene Hunting Approaches through the Combination of Linkage Analysis with Whole-Exome Sequencing in Mendelian Diseases: From Darwin to the Present Day

BACKGROUND

In the context of medical genetics, gene hunting is the process of identifying and functionally characterizing genes or genetic variations that contribute to disease phenotypes. In this review, we would like to summarize gene hunting process in terms of historical aspects from Darwin to now. For this purpose, different approaches and recent developments will be detailed.

SUMMARY

Linkage analysis and association studies are the most common methods in use for explaining the genetic background of hereditary diseases and disorders. Although linkage analysis is a relatively old approach, it is still a powerful method to detect disease-causing rare variants using family-based data, particularly for consanguineous marriages. As is known that, consanguineous marriages or endogamy poses a social problem in developing countries, however, this same condition also provides a unique opportunity for scientists to identify and characterize pathogenic variants. The rapid advancements in sequencing technologies and their parallel implementation together with linkage analyses now allow us to identify the candidate variants related to diseases in a relatively short time. Furthermore, we can now go one step further and functionally characterize the causative variant through in vitro and in vivo studies and unveil the variant-phenotype relationships on a molecular level more robustly. Key Messages: Herein, we suggest that the combined analysis of linkage and exome analysis is a powerful and precise tool to diagnose clinically rare and recessively inherited conditions.

First report of a successful pregnancy by preimplantation genetic testing for Beckwith-Wiedemann syndrome

OBJECTIVE

Beckwith-Wiedemann syndrome (BWS) is a rare imprinting gene disorder. Maternal CDKN1C mutation comprises 5% of etiologies of BWS. There is no successful report of preventing BWS by preimplantation genetic testing for monogenic disease (PGT-M) in the literature. Is PGT-M applicable for preventing BWS ?

CASE REPORT

This 39-year-old woman conceived naturally and delivered a boy who was diagnosed of BWS. The genetic testing of her son revealed CDKN1C gene mutation, and of the mother showed a carrier of the same mutation. She underwent controlled ovarian stimulation, oocyte pickup, and intracytoplasmic sperm injection. Trophectoderm biopsies were performed and samples were checked for PGT. Two wild-type and euploid embryos were thawed and transferred. One intrauterine pregnancy was achieved. The patient delivered a healthy female baby at 37 weeks of gestation.

CONCLUSION

In this case, we first report a successful pregnancy with a wild-type CDKN1C gene baby achieved by PGT-M.

Mechanistic insights into the digestion of complex dietary fibre by the rumen microbiota using combinatorial high-resolution glycomics and transcriptomic analyses

There is a knowledge gap regarding the factors that impede the ruminal digestion of plant cell walls or if rumen microbiota possess the functional activities to overcome these constraints. Innovative experimental methods were adopted to provide a high-resolution understanding of plant cell wall chemistries, identify higher-order structures that resist microbial digestion, and determine how they interact with the functional activities of the rumen microbiota. We characterized the total tract indigestible residue (TTIR) from cattle fed a low-quality straw diet using two comparative glycomic approaches: ELISA-based glycome profiling and total cell wall glycosidic linkage analysis. We successfully detected numerous and diverse cell wall glycan epitopes in barley straw (BS) and TTIR and determined their relative abundance pre- and post-total tract digestion. Of these, xyloglucans and heteroxylans were of higher abundance in TTIR. To determine if the rumen microbiota can further saccharify the residual plant polysaccharides within TTIR, rumen microbiota from cattle fed a diet containing BS were incubated with BS and TTIR ex vivo in batch cultures. Transcripts coding for carbohydrate-active enzymes (CAZymes) were identified and characterized for their contribution to cell wall digestion based on glycomic analyses, comparative gene expression profiles, and associated CAZyme families. High-resolution phylogenetic fingerprinting of these sequences encoded CAZymes with activities predicted to cleave the primary linkages within heteroxylan and arabinan. This experimental platform provides unprecedented precision in the understanding of forage structure and digestibility, which can be extended to other feed-host systems and inform next-generation solutions to improve the performance of ruminants fed low-quality forages.

Secular change of interdependencies in the clusters of the Japanese motor vehicle industry: a case study in the Tokai region

The Tokai region, including the prefectures of Aichi, Shizuoka, Mie, Gifu, and Shiga, in the center of Japan’s largest island, Honshu, houses clusters of the motor vehicle industry (MVI) and accounts for more than 50% of the domestic value from the MVI’s manufactured goods shipments. The subcontracted factories in the MVI tend to locate their main hub close to the assembly factories, creating large supply-chain structures with wide transactions. Since the 1990s, new automotive assembly factories and subcontracted factories have been built outside Aichi, in the Kyushu and Tohoku regions. This may have had an influence on local industries and economy. This study applies the hypothetical extraction method to the Tokai region using the Chubu-region multi-regional input–output tables to understand the secular changes of interdependencies in the MVI as each prefecture in the Tokai region has automotive assembly factories. The study shows that the Aichi MVI has strengthened its ties with industries in other prefectures in the Tokai region and elsewhere because the rate of growth for Aichi dropped from 2005 to 2010. Still, other prefecture rates rose, thus expanding its supply-chain structure throughout Japan more than before. In addition, the study shows that the Shizuoka MVI has increased business relationships not only within its region but also with Aichi. These results highlight that the industrial linkage between the Aichi MVI and the Shizuoka MVI has deepened. However, the results show that the supply-chain structures of the Mie, Gifu, and Shiga MVI have remained inside of their respective regions.

Mapping-by-Sequencing of Point and Insertional Mutations with Easymap

Random mutagenesis followed by screening for phenotypes of interest is a widely used strategy for genetic dissection of biological pathways; however, identifying the causal gene traditionally required time-consuming mapping approaches based on iterative linkage analysis. Mapping-by-sequencing accelerates this process, efficiently linking the phenotype of a mutant to a narrow candidate genomic region, using next-generation sequencing (NGS) data from a mapping population segregating for the mutant phenotype. To enable researchers at any bioinformatics skill level to conduct mapping-by-sequencing, we developed the Easymap mapping software. In this protocol we break down the steps involved in mapping-by-sequencing. First, we describe different ways of obtaining a mapping population and the steps used to generate NGS data. Next, we show how to analyze the NGS data using Easymap and how to obtain a list of candidate mutations, along with comprehensive information for assessing the potential causality of each candidate. Thus, this protocol enables the user to conduct mapping-by-sequencing using Easymap, facilitating the identification of causal loci for a mutant phenotype of interest.

Expanding the phenotypic spectrum of mutations in LRP2: a novel candidate gene of non-syndromic familial comitant strabismus

Background

Comitant strabismus (CS) is a heterogeneous disorder that is a major contributing factor to unilateral childhood-onset visual impairment. Studies have confirmed that genetic factors play an important role in the development of CS. The aim of this study was to identify the genetic cause of non-syndromic familial CS.

Methods

Fourteen unrelated CS families were recruited for the study. Twelve affected and 2 unaffected individuals from a large four-generation family (CS08) were selected to perform whole genome-wide linkage analysis. Parallel whole-exome sequencing (WES) was conducted in the same family (9 patients and 1 unaffected member) and 31 additional CS cases from 13 other unrelated families. Sanger sequencing was used to determine whether any of the remaining variants co-segregated with the disease phenotype in the corresponding family.

Results

Based on linkage analysis, CS in family CS08 mapped to a novel region of 34.17 centimorgan (cM) on chromosome 2q22.3-2q32.1 between markers D2S151 and D2S364, with a maximum log odds (LOD) score of 3.54 (theta = 0) at D2S142. Parallel WES identified a heterozygous variant, LRP2 c.335 A > G (p.Q112R), located in such a linkage interval that completely co-segregated with the disease in the family. Furthermore, another novel heterozygous variant (c.7274A > G, p.D2425G) in LRP2 that co-segregated was detected in 2 additional affected individuals from another unrelated family by WES. Both variants are predicted to be damaging by PolyPhen-2, SIFT and MutationTaster, and were absent in 100 ethnically matched normal controls.

Conclusion

LRP2 is a novel candidate genetic cause of non-syndromic familial CS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03155-z.

Rapid identification of candidate genes controlling male-sterility in Foxtail millet (Setaria italica)

Photo-(thermo-) sensitive genic male-sterile line is the key component of two-line hybridization system in foxtail millet (Setaria italica), but the genetic basis of male sterility in most male-sterile lines is still unclear. In the present study, a large F2 population was developed derived from a cross between the photo-(thermo-) sensitive male-sterile line A2 and the fertile-line 1484-5. Thirty plants with extreme high and extreme low fertility were selected from the population to construct a sterile DNA pool and a fertile DNA pool, respectively. Sequencing both DNA pools and data analysis revealed that two QTLs conferred male-sterility, qSiMS6.1 with a major effect and qSiMS6.2 with a minor effect, on chromosome 6. Both QTLs exhibited complete dominance. The major QTL, qSiMS6.1, was delimited to a 186-kb interval between the markers SiM20 and SiM9 by the joint analysis of QTL-seq and QTL mapping with SSR and structure variation markers. Millet_GLEAN_10020454 in this region is the most likely candidate gene for qSiMS6.1 since it is predicted to encode a male-sterile 5 like protein. These results lay a solid foundation for qSiMS6.1 cloning and provided gene resources for breeding new male-sterile lines.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01269-2.

Clinics and genetic background of hereditary gingival fibromatosis

BACKGROUND

Hereditary gingival fibromatosis (HGF) is a rare condition characterized by slowly progressive overgrowth of the gingiva. The severity of overgrowth may differ from mild causing phonetic and masticatory issues, to severe resulting in diastemas or malposition of teeth. Both, autosomal-dominant and autosomal-recessive forms of HGF are described. The aim of this review is a clinical overview, as well as a summary and discussion of the involvement of candidate chromosomal regions, pathogenic variants of genes, and candidate genes in the pathogenesis of HGF. The loci related to non-syndromic HGF have been identified on chromosome 2 (GINGF, GINGF3), chromosome 5 (GINGF2), chromosome 11 (GINGF4), and 4 (GINGF5). Of these loci, pathogenic variants of the SOS-1 and REST genes inducing HGF have been identified in the GINGF and the GINGF5, respectively. Furthermore, among the top 10 clusters of genes ranked by enrichment score, ATP binding, and fibronectin encoding genes were proposed as related to HGF.

CONCLUSION

The analysis of clinical reports as well as translational genetic studies published since the late'90s indicate the clinical and genetic heterogeneity of non-syndromic HGF and point out the importance of genetic studies and bioinformatics of more numerous unrelated families to identify novel pathogenic variants potentially inducing HGF. This strategy will help to unravel the molecular  mechanisms as well as uncover specific targets for novel and less invasive therapies of this rare, orphan condition.

In vitro immunomodulatory activity of water-soluble glucans from fresh and dried Longan (Dimocarpus longan Lour.)

Longan (Dimocarpus longan Lour.) is a seasonal tropical fruit used by Chinese medicine in both fresh and dried pulp forms. Their polysaccharides have been reported to have biological activity. However, their composition and immune activity have not yet been disclosed. To fulfil this aim, hot water-soluble polysaccharides of fresh and dried longan pulp were fractionated according to their molecular weight by ultrafiltration (10, 50, 100 kDa cut off). The main polysaccharides recovered were 1,6-linked glucans branched at O-3 (4-8%), O-2 (1%), O-2,4 (0.1%), and O-3,4 (0.1%). The drying process promotes the solubility of the polysaccharides. These glucans from fresh and dried longan pulp have immunomodulatory activity, shown by in vitro phagocytosis, NO, TNF-α, and IL-6 macrophages production assays. They showed also to inhibit the inflammatory response induced by LPS. The immunological activity of these glucans seems to have different responses dependent on their molecular weight, related to the immune regulatory pathways.

Linkage Analysis of Oligosaccharides and Polysaccharides: A Tutorial

Polysaccharides and oligosaccharides are a diverse group of natural polymers with important biological functions. The diversity of carbohydrate polymers is vast, ranging from small oligosaccharides of defined composition decorating proteins, to large, complex heteropolymers comprising integral cell wall components of plants, fungi and bacteria. An important step in the elucidation of unknown carbohydrate structures in a sample is the assessment of the various linkages present. This is accomplished by performing linkage analysis of the sample. The analysis proceeds as a successive series of chemical steps in which unlinked carbohydrate hydroxyls are marked with methyl groups, the sample is hydrolyzed into monosaccharides and reduced to alditols, and finally free hydroxyls are acetylated. Gas chromatography-mass spectrometry (GC-MS) analysis is employed to analyze the resultant partially methylated alditol acetates (PMAAs). The following paper reviews the major literature pertaining to the specific protocol for linkage analysis of carbohydrates outlined herein. The review details additional steps necessary for the completion of uronic acid linkage analysis, as well as analysis of chitin containing polymers. It also gives chromatographic examples of common erroneous results which the first time practitioner will want to be aware of. Our hope is that this protocol will serve as a definitive guide, allowing novice researchers to perform linkage analysis of carbohydrates in their own lab.

An intronic variant in the CELF4 gene is associated with risk for colorectal cancer

BACKGROUND

Germline predisposition variants associated with colorectal cancer (CRC) have been identified but all are not yet identified. We sought to identify the responsible predisposition germline variant in an extended high-risk CRC pedigree that exhibited evidence of linkage to the 18q12.2 region (TLOD = +2.81).

METHODS

DNA from two distantly related carriers of the hypothesized predisposition haplotype on 18q12.2 was sequenced to identify candidate variants. The candidate rare variants shared by the related sequenced subjects were screened in 3,094 CRC cases and 5x population-matched controls from UKBiobank to test for association. Further segregation of the variant was tested via Taqman assay in other sampled individuals in the pedigree.

RESULTS

Analysis of whole genome sequence data for the two related hypothesized predisposition haplotype carriers, restricted to the shared haplotype boundaries, identified multiple (n = 6) rare candidate non-coding variants that were tested for association with CRC risk in UKBiobank. A rare intronic variant ofCELF4 gene, rs568643870, was significantly associated with CRC (p = 0.004, OR = 5.0), and segregated with CRC in other members of the linked pedigree.

CONCLUSION

Evidence of segregation in a high-risk pedigree, case-control association in an external dataset, and identification of additional CRC-affected carriers in the linked pedigree support a role for a rareCELF4 intronic variant in CRC risk.

Structural elucidation and immuno-stimulatory activity of a novel polysaccharide containing glucuronic acid from the fungus Echinodontium tinctorium

An immuno-stimulatory polysaccharide (EtISPFa) was purified from water extract of the fungus Echinodontium tinctorium. EtISPFa has an estimated weight average molecular weight (M_w) of 1354 kDa and is composed of glucose (66.2 %), glucuronic acid (10.1 %), mannose (6.7 %), galactose (6.4 %), xylose (5.6 %), rhamnose (3.1 %), fucose (1.8 %), and arabinose (0.2 %). It has multiple glycosidic linkages, with 3-Glcp (19.8 %), 4-GlcpA (10.8 %), 6-Glcp (10.7 %), and 3,6-Glcp (8.7 %) being the most prominent. NMR analysis showed that EtISPFa has a backbone containing mostly of 3-substituted β-glucopyranose with 4-substituted glucopyranosyluronic acid. Short side chains consisting of an average of two β-glycopyranose residues, connected through 1→6 linkages, are attached to the 6-position of about every 4th or 5th backbone glucose residue. EtISPFa is a novel glucuronic acid-containing β-glucan capable of significantly inducing the production of cytokines IL-17, IL-16, MIP-2, G-CSF,GM-CSF, LIF, MIP-1α, MIP-1β, and RANTES in vitro. EtISPFa should be further explored for its immuno-stimulatory activity in vivo.

Evidence of genetic epistasis in autoimmune diabetes susceptibility revealed by mouse congenic sublines

Susceptibility to autoimmune diabetes is a complex genetic trait. Linkage analyses exploiting the NOD mouse, which spontaneously develops autoimmune diabetes, have proved to be a useful tool for the characterization of some of these traits. In a linkage analysis using 3A9 TCR transgenic mice on both B10.BR and NOD.H2^k backgrounds, we previously determined that both the Idd2 and Idd13 loci were linked to the proportion of immunoregulatory CD4^-CD8^- double negative (DN) T cells. In addition to Idd2 and Idd13, five other loci showed weak linkage to the proportion of DN T cells. Of interest, in an interim analysis, a locus on chromosome 12 is linked to DN T cell proportion in both the spleen and the lymph nodes. To determine the impact of this locus on DN T cells, we generated two congenic sublines, which we named Chr12P and Chr12D for proximal and distal, respectively. While 3A9 TCR:insHEL NOD.H2^k-Chr12D mice were protected from diabetes, 3A9 TCR:insHEL NOD.H2^k-Chr12P showed an increase in diabetes incidence. Yet, the proportion of DN T cells was similar to the parental 3A9 TCR NOD.H2^k strain for both of these congenic sublines. A genome-wide two dimensional LOD score analysis reveals genetic epistasis between chromosome 12 and the Idd13 locus. Altogether, this study identified further complex genetic interactions in defining the proportion of DN T cells, along with evidence of genetic epistasis within a locus on chromosome 12 influencing autoimmune susceptibility.

Genome-wide association and whole exome sequencing studies reveal a novel candidate locus for restless legs syndrome

INTRODUCTION

The restless legs syndrome (RLS) is a common heritable neurologic disorder which is characterized by an irresistible desire to move and unpleasant sensations in the legs.

METHODS

We aim to identify new variants associated with RLS by performing genome-wide linkage and subsequent association analysis of forty member's family with history of RLS.

RESULTS

We found evidence of linkage for three loci 7q21.11 (HLOD = 3.02), 7q21.13-7q21.3 (HLOD = 3.02) and 7q22.3 (HLOD = 3.09). Fine-mapping of those regions in association study using exome sequencing identified SEMA3A (p-value = 8.5·10^-4), PPP1R9A (p-value = 7.2·10^-4), PUS7 (p-value = 8.7·10^-4), CDHR3 (p-value = 7.2·10^-4), HBP1 (p-value = 1.5·10^-4) and COG5 (p-value = 1.5·10^-4) genes with p-values below significance threshold.

CONCLUSION

Linkage analysis with subsequent association study of exome variants identified six new genes associated with RLS mapped on 7q21 and q22.

Linkage analysis, GWAS, transcriptome analysis to identify candidate genes for rice seedlings in response to high temperature stress

BACKGROUND

Rice plants suffer from the rising temperature which is becoming more and more prominent. Mining heat-resistant genes and applying them to rice breeding is a feasible and effective way to solve the problem.

RESULT

Three main biomass traits, including shoot length, dry weight, and fresh weight, changed after abnormally high-temperature treatment in the rice seedling stage of a recombinant inbred lines and the natural indica germplasm population. Based on a comparison of the results of linkage analysis and genome-wide association analysis, two loci with lengths of 57 kb and 69 kb in qDW7 and qFW6, respectively, were associated with the rice response to abnormally high temperatures at the seedling stage. Meanwhile, based on integrated transcriptome analysis, some genes are considered as important candidate genes. Combining with known genes and analysis of homologous genes, it was found that there are eight genes in candidate intervals that need to be focused on in subsequent research.

CONCLUSIONS

The results indicated several relevant loci, which would help researchers to further discover beneficial heat-resistant genes that can be applied to rice heat-resistant breeding.

Linkage analysis, GWAS, transcriptome analysis to identify candidate genes for rice seedlings in response to high temperature stress

BACKGROUND

Rice plants suffer from the rising temperature which is becoming more and more prominent. Mining heat-resistant genes and applying them to rice breeding is a feasible and effective way to solve the problem.

RESULT

Three main biomass traits, including shoot length, dry weight, and fresh weight, changed after abnormally high-temperature treatment in the rice seedling stage of a recombinant inbred lines and the natural indica germplasm population. Based on a comparison of the results of linkage analysis and genome-wide association analysis, two loci with lengths of 57 kb and 69 kb in qDW7 and qFW6, respectively, were associated with the rice response to abnormally high temperatures at the seedling stage. Meanwhile, based on integrated transcriptome analysis, some genes are considered as important candidate genes. Combining with known genes and analysis of homologous genes, it was found that there are eight genes in candidate intervals that need to be focused on in subsequent research.

CONCLUSIONS

The results indicated several relevant loci, which would help researchers to further discover beneficial heat-resistant genes that can be applied to rice heat-resistant breeding.

Multiparental Population in Crops: Methods of Development and Dissection of Genetic Traits

Multiparental populations are located midway between association mapping that relies on germplasm collections and classic linkage analysis, based upon biparental populations. They provide several key advantages such as the possibility to include a higher number of alleles and increased level of recombination with respect to biparental populations, and more equilibrated allelic frequencies than association mapping panels. Moreover, in these populations new allele's combinations arise from recombination that may reveal transgressive phenotypes and make them a useful pre-breeding material. Here we describe the strategies for working with multiparental populations, focusing on nested association mapping populations (NAM) and multiparent advanced generation intercross populations (MAGIC). We provide details from the selection of founders, population development, and characterization to the statistical methods for genetic mapping and quantitative trait detection.

Genomic regions associated with important seed quality traits in food-grade soybeans

BACKGROUND

The production of soy-based food products requires specific physical and chemical characteristics of the soybean seed. Identification of quantitative trait loci (QTL) associated with value-added traits, such as seed weight, seed protein and sucrose concentration, could accelerate the development of competitive high-protein soybean cultivars for the food-grade market through marker-assisted selection (MAS). The objectives of this study were to identify and validate QTL associated with these value-added traits in two high-protein recombinant inbred line (RIL) populations.

RESULTS

The RIL populations were derived from the high-protein cultivar 'AC X790P' (49% protein, dry weight basis), and two high-yielding commercial cultivars, 'S18-R6' (41% protein) and 'S23-T5' (42% protein). Fourteen large-effect QTL (R2 > 10%) were identified associated with seed protein concentration. Of these QTL, seven QTL were detected in both populations, and eight of them were co-localized with QTL associated with either seed sucrose concentration or seed weight. None of the protein-related QTL was found to be associated with seed yield in either population. Sixteen candidate genes with putative roles in protein metabolism were identified within seven of these protein-related regions: qPro_Gm02-3, qPro_Gm04-4, qPro_Gm06-1, qPro_Gm06-3, qPro_Gm06-6, qPro_Gm13-4 and qPro-Gm15-3.

CONCLUSION

The use of RIL populations derived from high-protein parents created an opportunity to identify four novel QTL that may have been masked by large-effect QTL segregating in populations developed from diverse parental cultivars. In total, we have identified nine protein QTL that were detected either in both populations in the current study or reported in other studies. These QTL may be useful in the curated selection of new soybean cultivars for optimized soy-based food products.

Biallelic variants in four genes underlying recessive osteogenesis imperfecta

Osteogenesis imperfecta (OI) is an inherited heterogeneous rare skeletal disorder characterized by increased bone fragility and low bone mass. The disorder mostly segregates in an autosomal dominant manner. However, several rare autosomal recessive and X-linked forms, caused by mutations in 18 different genes, have also been described in the literature. Here, we present five consanguineous families segregating OI in an autosomal recessive pattern. Affected individuals in the five families presented severe forms of skeletal deformities. It included frequent bone fractures with abnormal healing, short stature, facial dysmorphism, osteopenia, joint laxity, and severe scoliosis. In order to search for the causative variants, DNA of at least one affected individual in three families (A-C) were subjected to whole exome sequencing (WES). In two other families (D-E), linkage analysis using highly polymorphic microsatellite markers was followed by Sanger sequencing. Sequence analysis revealed two novels and three previously reported disease-causing variants. The two novel homozygous variants including [c.824G > A; p.(Cys275Tyr)] in the SP7 gene and [c.397C > T, p.(Gln133*)] in the SERPINF1 gene were identified in families A and B, respectively. The three previously reported homozygous variants including [c.497G > A; p.(Arg166His)] in the SPARC gene, (c.359-3C > G; intron 2) and [c.677C > T; p.(Ser226Leu)] in the WNT1 gene were identified in family C, D, and E. In conclusion, our findings provided additional evidence of involvement of homozygous sequence variants in the SP7, SERPINF1, SPARC and WNT1 genes causing severe OI. It also highlights the importance of extensive genetic investigations to search for the culprit gene in each case of skeletal deformity.

Genetic analysis of loose smut (Ustilago tritici) resistance in Sonop spring wheat

BACKGROUND

The genetics of resistance to loose smut of wheat (Triticum aestivum L.) caused by the fungus Ustilago tritici (Pers.) Rostr. is not well understood. This study examines loose smut resistance in Sonop (TD-14), a South African spring wheat variety. A doubled haploid (DH) population of 163 lines derived from the cross Diamont/TD-14 was studied. The parents and progenies were inoculated with U. tritici races T2, T9, and T39 individually in growth facilities at Morden and Swift Current, Canada. Loose smut incidence (LSI) and partial loose smut resistance (PLSR) were assessed.

RESULTS

A whole genome linkage map was developed consisting of 11,519 SNP loci found on 31 linkage groups spanning 2845 cM. A new major resistance gene Ut11 was located to the distal end of chromosome arm 7BS. Ut11 conferred resistance to U. tritici race T2, but not races T9 and T39. Quantitative trait locus (QTL) mapping identified four QTL controlling LSI in the Diamont/TD-14 DH population on chromosomes 3B, 4B, 5B, and 7B (at Ut11) with TD-14 contributing the resistance alleles at three of these loci. The major QTL QUt.mrc-5B was effective against all three races and explained up to 81% of the phenotypic variation. The only QTL identified for PLSR coincided with the LSI QTL QUt.mrc-5B indicating that this locus affected both loose smut incidence and partial smutting of spikes.

CONCLUSIONS

A race-specific resistance gene Ut11 and a broadly effective resistance QTL QUt.mrc-5B were the main loci controlling loose smut resistance in the differential line TD-14 (cultivar Sonop). This study provides insight into the genetics of loose smut resistance in spring wheat Sonop and the single nucleotide polymorphism (SNP) markers linked to the resistance gene Ut11 and QTL QUt.mrc-5B will be useful for selecting loose smut resistance in breeding programs.

Identification of QTLs associated with curd architecture in cauliflower

BACKGROUND

Curd architecture is one of the most important characters determining the curd morphology of cauliflower. However, the genetic mechanism dissection of this complex trait at molecular level is lacking. Genes/QTLs responsible for the morphological differences between present-day loose-curd and compact-curd cauliflower haven't been well revealed.

RESULTS

Herein, by using a common compact-curd parent and two loose-curd parents, we developed two double haploid (DH) populations including 122 and 79 lines, respectively. For each population, we decomposed the curd architecture concept into four parameters (basal diameter, stalk length, stalk angle and curd solidity), and collected corresponding phenotypic data for each parameter across two environments. The Kosambi function and composite interval mapping algorithm were conducted to construct the linkage map and analyze the QTLs associated with curd architecture parameters. A total of 20 QTLs were detected with the minimum likelihood of odd (LOD) values ranging from 2.61 to 8.38 and the percentage of the phenotypic variance explained by each QTL (PVE) varying between 7.69 and 25.10%. Of these, two QTLs controlling stalk length (qSL.C6-1, qSL.C6-2) and two QTLs controlling curd solidity (qCS.C6-1 and qCS.C6-2) were steadily expressed in both environments. Further, qSL.C6-1, qSL.C6-2, qCS.C6-1 and qCS.C6-4 fell into the same chromosomal region of the reference genome, indicating that these loci are involved in pleiotropic effects or are tightly linked.

CONCLUSION

The current study identified a series of QTLs associated with curd architecture parameters, which might contribute essentially to the formation of present-day loose-curd cauliflower that is widely cultivated in China. These results may pave the way for intensive deciphering the molecular mechanisms of curd development and for marker-assisted selection of curd morphology in cauliflower breeding.

Cellulose, pectin and water in cell walls determine apple flesh viscoelastic mechanical properties

The viscoelastic mechanical properties are important quality traits for fleshy fruit uses. The contribution of cell wall polysaccharides chemistry and organization on their variability was studied in six varieties of apple. Correlation between damping and storage modulus of plasmolyzed tissue distinguished better apple varieties on their viscoelasticity than fresh samples. Galactose, arabinose and uronic acids correlated positively with the storage modulus of fresh apple samples (E'_f). These corresponded to 4-linked galactan but no specific arabinose linkage. Galacturonic acid branched on O-3 and terminal rhamnose correlated negatively with E'_f. These correlations formed two groups of fruit except for branched methyl-esterified galacturonic. Solid-state ¹³C NMR spectroscopy analyses showed that E'_f correlated negatively with cellulose C4 T_1ρ^H relaxation and positively with pectin methyl esters T_HH proton diffusion. The results point to the key roles of pectin structure and hydration and cellulose microfibrils distribution on apple mechanical properties.

Genome-wide linkage analysis of carotid artery traits in exceptionally long-lived families

BACKGROUND AND AIMS

Atherosclerosis develops with age and is partially controlled by genetics. Research to date has identified common variants with small effects on atherosclerosis related traits. We aimed to use family-based genome-wide linkage analysis to identify chromosomal regions potentially harboring rare variants with larger effects for atherosclerosis related traits.

METHODS

Participants included 2205 individuals from the Long Life Family Study (LLFS), which recruited families with exceptional longevity from Boston, New York, Pittsburgh, and Denmark. Participants underwent B-mode ultrasonography of the carotid arteries to measure intima-media thickness (IMT), inter-adventitial diameter (IAD), and plaque presence and severity. We conducted residual heritability and genome-wide linkage analyses adjusted for age, age², sex, and field center using pedigree-based maximum-likelihood methods in SOLAR.

RESULTS

All carotid traits were significantly heritable with a range of 0.68 for IAD to 0.38 for IMT. We identified three chromosomal regions with linkage to IAD (3q13; max LOD 5.3), plaque severity (17q22-q23, max LOD 3.2), and plaque presence (17q24, max LOD 3.1). No common allelic variants within these linkage peaks were associated with the carotid artery traits.

CONCLUSIONS

We identified three chromosomal regions with evidence of linkage to carotid artery diameter and atherosclerotic plaque in exceptionally long-lived families. Since common allelic variants within our linkage peaks did not account for our findings, future follow-up resequencing of these regions in LLFS families should help advance our understanding of atherosclerosis, CVD, and healthy vascular aging.

Enhancing Upland cotton for drought resilience, productivity, and fiber quality: comparative evaluation and genetic dissection

To provision the world sustainably, modern society must increase overall crop production, while conserving and preserving natural resources. Producing more with diminishing water resources is an especially daunting endeavor. Toward the goal of genetically improving drought resilience of cultivated Upland cotton (Gossypium hirsutum L.), this study addresses the genetics of differential yield components referred to as productivity and fiber quality traits under regular-water versus low-water (LW) field conditions. We used ten traits to assess water stress deficit, which included six productivity and four fiber quality traits on two recombinant inbred line (RIL) populations from reciprocally crossed cultivars, Phytogen 72 and Stoneville 474. To facilitate genetic inferences, we genotyped RILs with the CottonSNP63K array, assembled high-density linkage maps of over 7000 SNPs and then analyzed quantitative trait variations. Analysis of variance revealed significant differences for all traits (p < 0.05) in these RIL populations. Although the LW irrigation regime significantly reduced all traits, except lint percent, the RILs exhibited a broad phenotypic spectrum of heritable differences across the water regimes. Transgressive segregation occurred among the RILs, suggesting the possibility of genetic gain through phenotypic selection for drought resilience and perhaps through marker-based selection. Analyses revealed more than 150 quantitative trait loci (QTLs) associated with productivity and fiber quality traits (p < 0.005) on different genomic regions of the cotton genome. The multiple-QTL models analysis with LOD > 3.0 detected 21 QTLs associated with productivity and 22 QTLs associated with fiber quality. For fiber traits, strong clustering and QTL associations occurred in c08 and its homolog c24 as well as c10, c14, and c21. Using contemporary genome sequence assemblies and bioinformatically related information, the identification of genomic regions associated with responses to plant stress/drought elevates the possibility of using marker-assisted and omics-based selection to enhance breeding for drought resilient cultivars and identifying candidate genes and networks. RILs with different responses to drought indicated that it is possible to maintain high fiber quality under LW conditions or reduce the of LW impact on quality. The heritable variation among elite bi-parental RILs for productivity and quality under field drought conditions, and their association of QTLs, and thus specific genomic regions, indicate opportunities for breeding-based gains in water resource conservation, i.e., enhancing cotton's agricultural sustainability.

Genetic Analysis of a Large Family with Migraine, Vertigo, and Motion Sickness

BACKGROUND

Migraine is a common disorder most typically presenting as headache and often associated with vertigo and motion sickness. It is a genetically complex condition with multiple genes ultimately contributing to the predisposition and development of this episodic neurological disorder. We identified a large American family of 29 individuals of which 17 members suffered from at least one of these disorders, migraine, vertigo, or motion sickness. Many of these individuals suffered from several simultaneously. We hypothesized that vertigo and motion sickness may involve genes that are independent to those directly contributing to migraine susceptibility.

METHODS

Genome-wide linkage analysis performed using 400 microsatellite repeat markers spaced at 10 cM throughout the genome. The members of this family were phenotyped for each condition, migraine, vertigo, and motion sickness and analyzed separately. Statistical analysis was performed using two-point and multipoint linkage analysis employing a number of models including autosomal recessive or dominant patterns of inheritance with high and low genetic penetrance.

RESULTS

We identified a novel locus for migraine, 9q13-q22 (maximum two-point logarithm of odds [LOD] score-2.51). In addition, there are suggestive LOD scores that localize to different chromosomes for each phenotype; vertigo (chromosome 18, LOD score of 1.82) and motion sickness (chromosome 4, LOD score of 2.09).

CONCLUSIONS

Our analysis supports our hypothesis that the migraine-associated vertigo and motion sickness may involve distinct susceptibility genes.

Linkage and sequence analysis of neutral oligosaccharides by negative-ion MALDI tandem mass spectrometry with laser-induced dissociation

Mass spectrometry (MS) has become the primary method for high-sensitivity structural determination of oligosaccharides. Fragmentation in the negative-ion MS can provide a wealth of structural information and these can be used for sequence determination. However, although negative-ion MS of neutral oligosaccharide using the deprotonated molecule [M-H]^- as the precursor has been very successful for electrospray ionization (ESI), it has only limited success for matrix-assisted laser desorption/ionization (MALDI). In the present study, the features of negative-ion MALDI primary spectra were investigated in detail and the product-ion spectra using [M-H]^- and [M+Cl]^- as the precursors were carefully compared. The formation of [M-H]^- was the main difficulty for MALDI while [M+Cl]^- was proved to be useful as alternative precursor anion for MALDI-MS/MS to produce similar fragmentation for sequencing of neutral oligosaccharides. N-(1-naphthyl)ethylenediamine dihydrochloride was then used as both the matrix and the Cl^- dopant to evaluate the extent of structural information that can be obtained by negative-ion fragmentation from [M+Cl]^- using laser-induced dissociation (LID)-MS/MS for linkage assignment of gluco-oligosaccharides and for typing of blood-group ABO(H) and Lewis antigens on either type 1 or type 2 backbone-chains.

Statistical Methods and Software for Substance Use and Dependence Genetic Research

BACKGROUND

Substantial substance use disorders and related health conditions emerged dur-ing the mid-20th century and continue to represent a remarkable 21st century global burden of disease. This burden is largely driven by the substance-dependence process, which is a complex process and is influenced by both genetic and environmental factors. During the past few decades, a great deal of pro-gress has been made in identifying genetic variants associated with Substance Use and Dependence (SUD) through linkage, candidate gene association, genome-wide association and sequencing studies.

METHODS

Various statistical methods and software have been employed in different types of SUD ge-netic studies, facilitating the identification of new SUD-related variants.

CONCLUSION

In this article, we review statistical methods and software that are currently available for SUD genetic studies, and discuss their strengths and limitations.

HCV genotype 1-6 NS3 residue 80 substitutions impact protease inhibitor activity and promote viral escape

BACKGROUND & AIMS

Protease inhibitors (PIs) are of central importance in the treatment of patients with chronic hepatitis C virus (HCV) infection. HCV NS3 protease (NS3P) position 80 displays polymorphisms associated with resistance to the PI simeprevir for HCV genotype 1a. We investigated the effects of position-80-substitutions on fitness and PI-resistance for HCV genotypes 1-6, and analyzed evolutionary mechanisms underlying viral escape mediated by pre-existing Q80K.

METHODS

The fitness of infectious NS3P recombinants of HCV genotypes 1-6, with engineered position-80-substitutions, was studied by comparison of viral spread kinetics in Huh-7.5 cells in culture. Median effective concentration (EC50) and fold resistance for PIs simeprevir, asunaprevir, paritaprevir, grazoprevir, glecaprevir and voxilaprevir were determined in short-term treatment assays. Viral escape was studied by long-term treatment of genotype 1a recombinants with simeprevir, grazoprevir, glecaprevir and voxilaprevir and of genotype 3a recombinants with glecaprevir and voxilaprevir, next generation sequencing, NS3P substitution linkage and haplotype analysis.

RESULTS

Among tested PIs, only glecaprevir and voxilaprevir showed pan-genotypic activity against the original genotype 1-6 culture viruses. Variants with position-80-substitutions were all viable, but fitness depended on the specific substitution and the HCV isolate. Q80K conferred resistance to simeprevir across genotypes but had only minor effects on the activity of the remaining PIs. For genotype 1a, pre-existing Q80K mediated accelerated escape from simeprevir, grazoprevir and to a lesser extent glecaprevir, but not voxilaprevir. For genotype 3a, Q80K mediated accelerated escape from glecaprevir and voxilaprevir. Escape was mediated by rapid and genotype-, PI- and PI-concentration-dependent co-selection of clinically relevant resistance associated substitutions.

CONCLUSIONS

Position-80-substitutions had relatively low fitness cost and the potential to promote HCV escape from clinically relevant PIs in vitro, despite having a minor impact on results in classical short-term resistance assays.

LAY SUMMARY

Among all clinically relevant hepatitis C virus protease inhibitors, voxilaprevir and glecaprevir showed the highest and most uniform activity against cell culture infectious hepatitis C virus with genotype 1-6 proteases. Naturally occurring amino acid changes at protease position 80 had low fitness cost and influenced sensitivity to simeprevir, but not to other protease inhibitors in short-term treatment assays. Nevertheless, the pre-existing change Q80K had the potential to promote viral escape from protease inhibitors during long-term treatment by rapid co-selection of additional resistance changes, detected by next generation sequencing.

Structural elucidation of fucoidan from Cladosiphon okamuranus (Okinawa mozuku)

Fucoidan is a sulphated polysaccharide, made up mainly of l-fucose, which is found in brown seaweeds. Its chemical structure is diverse and depends on maturity, species and geographical location. The objective of this study was to elucidate the chemical structure of fucoidan from Cladosiphon okamuranus harvested in Japan. The fucoidan was subject to purification prior to monosaccharide profiling, sulphate content determination, and linkage analysis. Our results showed that Japanese Cladosiphon okamuranus fucoidan contained 70.13 ± 0.22 wt% fucose and 15.16 ± 1.17 wt% sulphate. Other minor monosaccharides found were d-xylose, d-galactose, d-mannose, d-glucose, d-arabinose, d-rhamnose and d-glucuronic acid. Linkage analysis revealed that fucopyranoside units along the backbone are linked, through α-1,3-glycosidic bonds, with fucose branching at C-2, and one sulphate group at C-4 per every three fucose units, i.e. the structure of fucoidan from Japanese Cladosiphon okamuranus is [→3)-α-fuc(1→]0.52[→3)-α-fuc-4-OSO3-(1→]0.33[→2)-α-fuc]0.14.

A comparison of cosegregation analysis methods for the clinical setting

Quantitative cosegregation analysis can help evaluate the pathogenicity of genetic variants. However, genetics professionals without statistical training often use simple methods, reporting only qualitative findings. We evaluate the potential utility of quantitative cosegregation in the clinical setting by comparing three methods. One thousand pedigrees each were simulated for benign and pathogenic variants in BRCA1 and MLH1 using United States historical demographic data to produce pedigrees similar to those seen in the clinic. These pedigrees were analyzed using two robust methods, full likelihood Bayes factors (FLB) and cosegregation likelihood ratios (CSLR), and a simpler method, counting meioses. Both FLB and CSLR outperform counting meioses when dealing with pathogenic variants, though counting meioses is not far behind. For benign variants, FLB and CSLR greatly outperform as counting meioses is unable to generate evidence for benign variants. Comparing FLB and CSLR, we find that the two methods perform similarly, indicating that quantitative results from either of these methods could be combined in multifactorial calculations. Combining quantitative information will be important as isolated use of cosegregation in single families will yield classification for less than 1% of variants. To encourage wider use of robust cosegregation analysis, we present a website ( http://www.analyze.myvariant.org ) which implements the CSLR, FLB, and Counting Meioses methods for ATM, BRCA1, BRCA2, CHEK2, MEN1, MLH1, MSH2, MSH6, and PMS2. We also present an R package, CoSeg, which performs the CSLR analysis on any gene with user supplied parameters. Future variant classification guidelines should allow nuanced inclusion of cosegregation evidence against pathogenicity.

Multiple genetic loci define Ca++ utilization by bloodstream malaria parasites

Background

Bloodstream malaria parasites require Ca⁺⁺ for their development, but the sites and mechanisms of Ca⁺⁺ utilization are not well understood. We hypothesized that there may be differences in Ca⁺⁺ uptake or utilization by genetically distinct lines of P. falciparum. These differences, if identified, may provide insights into molecular mechanisms.

Results

Dose response studies with the Ca⁺⁺ chelator EGTA (ethylene glycol-bis(β-aminoethyl ether)-N,N,N′,N′-tetraacetic acid) revealed stable differences in Ca⁺⁺ requirement for six geographically divergent parasite lines used in previous genetic crosses, with the largest difference seen between the parents of the HB3 x Dd2 cross. Genetic mapping of Ca⁺⁺ requirement yielded complex inheritance in 34 progeny clones with a single significant locus on chromosome 7 and possible contributions from other loci. Although encoded by a gene in the significant locus and a proposed Ca⁺⁺ target, PfCRT (P. falciparum chloroquine resistance transporter), the primary determinant of clinical resistance to the antimalarial drug chloroquine, does not appear to contribute to this quantitative trait. Stage-specific application of extracellular EGTA also excluded determinants associated with merozoite egress and erythrocyte reinvasion.

Conclusions

We have identified differences in Ca⁺⁺ utilization amongst P. falciparum lines. These differences are under genetic regulation, segregating as a complex trait in genetic cross progeny. Ca⁺⁺ uptake and utilization throughout the bloodstream asexual cycle of malaria parasites represents an unexplored target for therapeutic intervention.

Clinical and genetic characterization of an Italian family with slow-channel syndrome

INTRODUCTION

The slow-channel congenital myasthenic syndrome (SCCMS) is a postsynaptic form of congenital myasthenic syndromes (CMSs), a clinically heterogeneous group of disorders caused by genetic defects leading to an abnormal signal transmission at the endplate.

METHODS

We report clinical and molecular data of a multigenerational family in which the presentation of a progressive proximal-distal weakness with ocular involvement led to a number of different clinical diagnoses.

RESULTS

A comprehensive genetic study which included whole-genome linkage analysis and whole-exome sequencing identified a heterozygous missense substitution (c.721C>T, p.L241F) in the ε subunit of the acetylcholine receptor (CHRNE) that was consistent with clinical weakness in all patients.

DISCUSSION

SCCMS is characterized by a broad and heterogeneous clinical phenotype in which disease onset, symptoms, severity, and progression can be highly variable even between family members. The identification of a CHRNE mutation allowed to make the definitive diagnosis of CMS in this family and contributed to define the clinical spectrum of this disease.

Construction of a high-density genetic map: genotyping by sequencing (GBS) to map purple seed coat color (Psc) in hulless barley

Background

Colored hulless barley are more suitable in food processing compared to normal (yellow) varieties because it is rich in bioactive compounds and produces higher extraction pearling fractions. Therefore, seed coat color is an important agronomic trait for the breeding and study of hulless barley.

Results

Genotyping-by-sequencing single-nucleotide polymorphism (GBS-SNP) analysis of a doubled haploid (DH) mapping population (Nierumuzha × Kunlun10) was conducted to map the purple seed coat color genes (Psc). A high-density genetic map of hulless barley was constructed, which contains 3662 efficient SNP markers with 1129 bin markers. Seven linkage groups were resolved, which had a total length of 645.56 cM. Chromosome length ranged from 60.21 cM to 127.21 cM, with average marker density of 0.57 cM. A total of five loci accounting for 3.79% to 23.86% of the observed phenotypic variation for Psc were detected using this high-density map. Five structural candidate genes (F3’M, HID, UF3GT, UFGT and 5MAT) and one regulatory factor (Ant1) related to flavonoid or anthocyanin biosynthesis were identified..

Conclusions

Five structural candidate genes and one regulatory factor related to flavonoid or anthocyanin biosynthesis have been identified using a high-density genetic map of hulless barley. This study lays the foundation for map-based cloning of Psc but provides a valuable tool for studying marker-trait associations and its application to marker-assisted breeding of hulless barley.

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0072-6) contains supplementary material, which is available to authorized users.

Contribution of Inbred Singletons to Variance Component Estimation of Heritability and Linkage

OBJECTIVES

An interesting consequence of consanguinity is that the inbred singleton becomes informative for genetic variance. We determine the contribution of an inbred singleton to variance component analysis of heritability and linkage.

METHODS

Statistical theory for the power of variance component analysis of quantitative traits is used to determine the expected contribution of an inbred singleton to likelihood-ratio tests of heritability and linkage.

RESULTS

In variance component models, an inbred singleton contributes relatively little to a test of heritability but can contribute substantively to a test of linkage. For small-to-moderate quantitative trait locus (QTL) effects and a level of inbreeding comparable to matings between first cousins (the preferred form of union in many human populations), an inbred singleton can carry nearly 25% of the information of a non-inbred sib pair. In more highly inbred contexts available with experimental animal populations, nonhuman primate colonies, and some human subpopulations, the contribution of an inbred singleton relative to a sib pair can exceed 50%.

CONCLUSIONS

Inbred individuals, even in isolation from other members of a sample, can contribute to variance component estimation and tests of heritability and linkage. Under certain conditions, the informativeness of the inbred singleton can approach that of a non-inbred sib pair.

Genome-wide scan for circulating vascular adhesion protein-1 levels: MACROD2 as a potential transcriptional regulator of adipogenesis

AIMS/INTRODUCTION

Vascular adhesion protein-1 (VAP-1) is a membrane-bound amine oxidase highly expressed in mature adipocytes and released into the circulation. VAP-1 has been strongly implicated in several pathological processes, including diabetes, inflammation, hypertension, hepatic steatosis and renal diseases, and is an important disease marker and therapeutic target. Here, we aimed to identify the genetic loci for circulating VAP-1 levels.

MATERIALS AND METHODS

We carried out a genomic-wide linkage scan for the quantitative trait locus of circulating VAP-1 levels in 1,100 Han Chinese individuals from 398 families in the Stanford Asian Pacific Program for Hypertension and Insulin Resistance study. Regional association fine mapping was carried out using additional single-nucleotide polymorphisms.

RESULTS

The estimated heritability of circulating VAP-1 levels is high (h2 = 69%). The most significant quantitative trait locus for circulating VAP-1 was located at 38 cM on chromosome 20, with a maximum empirical logarithm of odds score of 4.11 (P = 6.86 × 10-6 ) in females. Regional single-nucleotide polymorphism fine mapping within a 1-unit support region showed the strongest association signals in the MACRO domain containing 2 (MACROD2) gene in females (P = 5.38 × 10-6 ). Knockdown of MACROD2 significantly suppressed VAP-1 expression in human adipocytes, as well as the expression of key adipogenic genes. Furthermore, MACROD2 expression was found to be positively associated with VAP-1 in human visceral adipose tissue.

CONCLUSION

MACROD2 is a potential genetic determinant of serum VAP-1 levels, probably through transcriptional regulation of adipogenesis.

Linkage analysis of multiplex Caribbean Hispanic families loaded for unexplained early-onset cases identifies novel Alzheimer's disease loci

Introduction

Less than 10% of early-onset Alzheimer's disease (EOAD) is explained by known mutations.

Methods

We conducted genetic linkage analysis of 68 well-phenotyped Caribbean Hispanic families without clear inheritance patterns or mutations in APP, PSEN1, and PSEN2 and with two or more individuals with EOAD.

Results

We identified 16 (logarithm of odds > 3.6) linked regions, including eight novel loci for EOAD (2p15, 5q14.1, 11p15.1, 13q21.22, 13q33.1, 16p12.1, 20p12.1, and 20q11.21) and eight regions previously associated with late-onset Alzheimer's disease. The strongest signal was observed at 16p12.1 (25 cM, 33 Mb; heterogeneity logarithm of odds = 5.3), ∼3 Mb upstream of the ceroid lipofuscinosis 3 (CLN3) gene associated with juvenile neuronal ceroid lipofuscinosis (JNCL), which functions in retromer trafficking and has been reported to alter intracellular processing of the amyloid precursor protein.

Discussion

This study supports the notion that the genetic architectures of unexplained EOAD and late-onset AD overlap partially, but not fully.

Identifying Specific and Differentially Linked Glycosyl Residues in Mammalian Glycans by Targeted LC-MS Analysis

Glycans, which are widespread in nature, consist of a large number of monosaccharides linked via glycosidic bonds. Due to the complex nature of glycan structures on glycoproteins, assessing the configuration and positions of the glycosidic linkages of a glycan is a subject of considerable interest. In this study, a method for accomplishing this using partially O-methylated alditols (PMAs) from glycans combined with LC-MS analysis is reported. N-Glycans were first per-methylated with methyl iodide, and the levels of methylation were further confirmed by MALDI-TOF. PMAs were then produced via complete hydrolysis and reduction. PMAs derived from Fetuin N-glycan and Lewis^a antigen carbohydrates were successfully detected by LC-MS analysis. This analysis can be performed without the need for an additional derivatization step for GC analysis, and should be suitable for use in conjunction with a LC-MS-based analysis platform.

Enhancement of Glen Moy x Latham raspberry linkage map using GbS to further understand control of developmental processes leading to fruit ripening

BACKGROUND

The changing climate is altering timing of key fruit ripening processes and increasing the occurrence of fruit defects. To improve our understanding of the genetic control of raspberry fruit development an enhanced genetic linkage map was developed and used to examine ripening phenotypic data.

RESULTS

In this study we developed an enhanced genetic linkage map for the raspberry cvs. Glen Moy x Latham reference mapping population using genotyping by sequencing (GbS). Alignment to a newly sequenced draft reference genome of red raspberry, cultivar (cv.) Glen Moy, identified 8019 single nucleotide polymorphisms (SNPs). After stringent filtering to take account of read coverage over all the progeny individuals, association with a single chromosome, heterozygosity and marker regression mapping, 2348 high confidence SNPs were retained and integrated with an existing raspberry genetic map. The linkage map contained many more SNPs segregating in Latham than in Glen Moy. This caused difficulties in quantitative trait loci (QTL) mapping with standard software and a novel analysis based on a hidden Markov model was used to improve the mapping. QTL mapping using the newly generated dense genetic map not only corroborated previously identified genetic locations but also provided additional genetic elements controlling fruit ripening in raspberry.

CONCLUSION

The high-density GbS map located the QTL peaks more precisely than in earlier studies, aligned the QTLs with Glen Moy genome scaffolds, narrowed the range of potential candidate genes to these regions that can be utilised in other populations or in gene expression studies to confirm their role and increased the repertoire of markers available to understand the genetic control of fruit ripening traits.