Family-Based Association Study of Pulmonary Function in a Population in Northeast Asia

Large-scale genome sequencing redefines the genetic footprints of high-altitude adaptation in Tibetans

BACKGROUND

Tibetans are genetically adapted to high-altitude environments. Though many studies have been conducted, the genetic basis of the adaptation remains elusive due to the poor reproducibility for detecting selective signatures in the Tibetan genomes.

RESULTS

Here, we present whole-genome sequencing (WGS) data of 1001 indigenous Tibetans, covering the major populated areas of the Qinghai-Tibetan Plateau in China. We identify 35 million variants, and more than one-third of them are novel variants. Utilizing the large-scale WGS data, we construct a comprehensive map of allele frequency and linkage disequilibrium and provide a population-specific genome reference panel, referred to as 1KTGP. Moreover, with the use of a combined approach, we redefine the signatures of Darwinian-positive selection in the Tibetan genomes, and we characterize a high-confidence list of 4320 variants and 192 genes that have undergone selection in Tibetans. In particular, we discover four new genes, TMEM132C, ATP13A3, SANBR, and KHDRBS2, with strong signals of selection, and they may account for the adaptation of cardio-pulmonary functions in Tibetans. Functional annotation and enrichment analysis indicate that the 192 genes with selective signatures are likely involved in multiple organs and physiological systems, suggesting polygenic and pleiotropic effects.

CONCLUSIONS

Overall, the large-scale Tibetan WGS data and the identified adaptive variants/genes can serve as a valuable resource for future genetic and medical studies of high-altitude populations.

Differential expression of the Tmem132 family genes in the developing mouse nervous system

The family of novel transmembrane proteins (TMEM) 132 have been associated with multiple neurological disorders and cancers in humans, but have hardly been studied in vivo. Here we report the expression patterns of the five Tmem132 genes (a, b, c, d and e) in developing mouse nervous system with RNA in situ hybridization in wholemount embryos and tissue sections. Our results reveal differential and partially overlapping expression of multiple Tmem132 family members in both the central and peripheral nervous system, suggesting potential partial redundancy among them.

Elastic Correlation Adjusted Regression (ECAR) scores for high dimensional variable importance measuring

Investigation of the genetic basis of traits or clinical outcomes heavily relies on identifying relevant variables in molecular data. However, characteristics such as high dimensionality and complex correlation structures of these data hinder the development of related methods, resulting in the inclusion of false positives and negatives. We developed a variable importance measure method, termed the ECAR scores, that evaluates the importance of variables in the dataset. Based on this score, ranking and selection of variables can be achieved simultaneously. Unlike most current approaches, the ECAR scores aim to rank the influential variables as high as possible while maintaining the grouping property, instead of selecting the ones that are merely predictive. The ECAR scores' performance is tested and compared to other methods on simulated, semi-synthetic, and real datasets. Results showed that the ECAR scores improve the CAR scores in terms of accuracy of variable selection and high-rank variables' predictive power. It also outperforms other classic methods such as lasso and stability selection when there is a high degree of correlation among influential variables. As an application, we used the ECAR scores to analyze genes associated with forced expiratory volume in the first second in patients with lung cancer and reported six associated genes.

Efficient multiplexed genome engineering with a polycistronic tRNA and CRISPR guide-RNA reveals an important role of detonator in reproduction of Drosophila melanogaster

Genome association studies in human and genetic studies in mouse implicated members of the transmembrane protein 132 (TMEM132) family in multiple conditions including panic disorder, hearing loss, limb and kidney malformation. However, the presence of five TMEM132 paralogs in mammalian genomes makes it extremely challenging to reveal the full requirement for these proteins in vivo. In contrast, there is only one TMEM132 homolog, detonator (dtn), in the genome of fruit fly Drosophila melanogaster, enabling straightforward research into its in vivo function. In the current study, we generate multiple loss-of-function dtn mutant fly strains through a polycistronic tRNA-gRNA approach, and show that most embryos lacking both maternal and paternal dtn fail to hatch into larvae, indicating an essential role of dtn in Drosophila reproduction.

CG4928 Is Vital for Renal Function in Fruit Flies and Membrane Potential in Cells: A First In-Depth Characterization of the Putative Solute Carrier UNC93A

The number of transporter proteins that are not fully characterized is immense. Here, we used Drosophila melanogaster and human cell lines to perform a first in-depth characterization of CG4928, an ortholog to the human UNC93A, of which little is known. Solute carriers regulate and maintain biochemical pathways important for the body, and malfunctioning transport is associated with multiple diseases. Based on phylogenetic analysis, CG4928 is closely related to human UNC93A and has a secondary and a tertiary protein structure and folding similar to major facilitator superfamily transporters. Ubiquitous knockdown of CG4928 causes flies to have a reduced secretion rate from the Malpighian tubules; altering potassium content in the body and in the Malpighian tubules, homologous to the renal system; and results in the development of edema. The edema could be rescued by using amiloride, a common diuretic, and by maintaining the flies on ion-free diets. CG4928-overexpressing cells did not facilitate the transport of sugars and amino acids; however, proximity ligation assay revealed that CG4928 co-localized with TASK¹ channels. Overexpression of CG4928 resulted in induced apoptosis and cytotoxicity, which could be restored when cells were kept in high-sodium media. Furthermore, the basal membrane potential was observed to be disrupted. Taken together, the results indicate that CG4928 is of importance for generating the cellular membrane potential by an unknown manner. However, we speculate that it most likely acts as a regulator or transporter of potassium flows over the membrane.

The Neuronal and Peripheral Expressed Membrane-Bound UNC93A Respond to Nutrient Availability in Mice

Many transporters such as the solute carriers belonging to the Major facilitator superfamily Pfam clan are orphans in that their tissue and cellular localization as well as substrate profile and function are still unknown. Here we have characterized the putative solute carrier UNC93A. We aimed to investigate the expression profile on both protein and mRNA level of UNC93A in mouse since it has not been clarified. UNC93A staining was found in cortex, hippocampus and cerebellum. It was found to be expressed in many neurons, but not all, with staining located in close proximity to the plasma membrane. Furthermore, we aimed to extend the starvation data available for Unc93a in hypothalamic cell cultures from mouse. We investigated the Unc93a alterations with focus on amino acid deprivation in embryonic cortex cells from mice as well as 24 h starvation in adult male mice and compared it to recently studied putative and known solute carriers. Unc93a expression was found both in the brain and peripheral organs, in low to moderate levels in the adult mice and was affected by amino acid deprivation in embryonic cortex cultures and starvation in in vivo samples. In conclusion, the membrane-bound UNC93A is expressed in both the brain and peripheral tissues and responds to nutrient availability in mice.