From K[ZnBP2O8] to K[(Zn0.5Al0.5)2P2O8]: Substitution-Induced Microporous Zincoaluminophosphate with a Congruent-Melting Behavior.

An open-frame aluminophosphate, K[(Zn0.5Al0.5)2P2O8] (KZAPO), was rationally designed by a substitution design strategy and synthesized by a high-temperature molten salt method. Compared with the parent crystal of K[ZnBP2O8], KZAPO was characterized by similar 4 × 8 × 8 networks, a comparable short-wave ultraviolet transparency and a more regular tetrahedral frame with the mixing of (ZnO4)6- and (AlO4)5- anionic groups, highlighting the multifunctional roles that anionic group mixing played in structural and property modulations. In particular, KZAPO was characterized by a high thermal stability (over 850 °C) and a congruent-melting behavior, being conducive to practical applications.

Multi-omics study in monozygotic twins confirm the contribution of de novo mutation to psoriasis.

BACKGROUND: Genome-wide association studies have identified over 120 risk loci for psoriasis. However, most of the variations are located in non-coding region with high frequency and small effect size. Pathogenetic variants are rarely reported except HLA-C*0602 with the odds ratio being approximately 4.0 in Chinese population. Although rare variations still account for a small proportion of phenotypic variances in complex diseases, their effect on phenotypes is large. Recently, more and more studies focus on the low-frequency functional variants and have achieved a certain amount of success. METHOD: Whole genome sequencing and sanger sequencing was performed on 8 MZ twin pairs discordant for psoriasis to scan and verified the de novo mutations (DNMs). Additionally, 665 individuals with about 20 years' medical history versus 2054 healthy controls and two published large population studies which had about 8 years' medical history (including 10,727 cases versus 10,582 controls) were applied to validate the enrichment of rare damaging mutations in two DNMs genes. Besides, to verify the pathogenicity of candidate DNM in C3, RNA-sequencing for CD4+, CD8+ T cells of twins and lesion, non-lesion skin of psoriasis patients were carried out. Meanwhile, the enzyme-linked immunosorbent assay kit was used to detect the level of C3, C3b in the supernatant of peripheral blood. RESULT: A total of 27 DNMs between co-twins were identified. We found six of eight twins carry HLA-C∗0602 allele which have large effects on psoriasis. And it is interesting that a missense mutation in SPRED1 and a splice region mutation in C3 are found in the psoriasis individuals in the other two MZ twin pairs without carrying HLA-C*0602 allele. In the replication stage, we found 2 loss-of-function (LOF) variants of C3 only in 665 cases with about 20 years' medical history and gene-wise analysis in 665 cases and 2054 controls showed that the rare missense mutations in C3 were enriched in cases (OR = 1.91, P = 0.0028). We further scanned the LOF mutations of C3 in two published studies (about 8 years' medical history), and found one LOF mutation in the case without carrying HLA-C*0602. In the individual with DNM in C3, RNA sequencing showed the expression level of C3 in skin was significant higher than healthy samples in public database (TPM fold change = 1.40, P = 0.000181) and ELISA showed protein C3 in peripheral blood was higher (~2.2-fold difference) than the other samples of twins without DNM in C3. CONCLUSION: To the best of our knowledge, this is the first report that DNM in C3 is the likely pathological mutations, and it provided a better understanding of the genetic etiology of psoriasis and additional treatments for this disease.

Whole-Genome Sequencing Reveals Diverse Models of Structural Variations in Esophageal Squamous Cell Carcinoma.

Comprehensive identification of somatic structural variations (SVs) and understanding their mutational mechanisms in cancer might contribute to understanding biological differences and help to identify new therapeutic targets. Unfortunately, characterization of complex SVs across the whole genome and the mutational mechanisms underlying esophageal squamous cell carcinoma (ESCC) is largely unclear. To define a comprehensive catalog of somatic SVs, affected target genes, and their underlying mechanisms in ESCC, we re-analyzed whole-genome sequencing (WGS) data from 31 ESCCs using Meerkat algorithm to predict somatic SVs and Patchwork to determine copy-number changes. We found deletions and translocations with NHEJ and alt-EJ signature as the dominant SV types, and 16% of deletions were complex deletions. SVs frequently led to disruption of cancer-associated genes (e.g., CDKN2A and NOTCH1) with different mutational mechanisms. Moreover, chromothripsis, kataegis, and breakage-fusion-bridge (BFB) were identified as contributing to locally mis-arranged chromosomes that occurred in 55% of ESCCs. These genomic catastrophes led to amplification of oncogene through chromothripsis-derived double-minute chromosome formation (e.g., FGFR1 and LETM2) or BFB-affected chromosomes (e.g., CCND1, EGFR, ERBB2, MMPs, and MYC), with approximately 30% of ESCCs harboring BFB-derived CCND1 amplification. Furthermore, analyses of copy-number alterations reveal high frequency of whole-genome duplication (WGD) and recurrent focal amplification of CDCA7 that might act as a potential oncogene in ESCC. Our findings reveal molecular defects such as chromothripsis and BFB in malignant transformation of ESCCs and demonstrate diverse models of SVs-derived target genes in ESCCs. These genome-wide SV profiles and their underlying mechanisms provide preventive, diagnostic, and therapeutic implications for ESCCs.

Genomic analyses reveal mutational signatures and frequently altered genes in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is one of the most common cancers worldwide and the fourth most lethal cancer in China. However, although genomic studies have identified some mutations associated with ESCC, we know little of the mutational processes responsible. To identify genome-wide mutational signatures, we performed either whole-genome sequencing (WGS) or whole-exome sequencing (WES) on 104 ESCC individuals and combined our data with those of 88 previously reported samples. An APOBEC-mediated mutational signature in 47% of 192 tumors suggests that APOBEC-catalyzed deamination provides a source of DNA damage in ESCC. Moreover, PIK3CA hotspot mutations (c.1624G>A [p.Glu542Lys] and c.1633G>A [p.Glu545Lys]) were enriched in APOBEC-signature tumors, and no smoking-associated signature was observed in ESCC. In the samples analyzed by WGS, we identified focal (<100 kb) amplifications of CBX4 and CBX8. In our combined cohort, we identified frequent inactivating mutations in AJUBA, ZNF750, and PTCH1 and the chromatin-remodeling genes CREBBP and BAP1, in addition to known mutations. Functional analyses suggest roles for several genes (CBX4, CBX8, AJUBA, and ZNF750) in ESCC. Notably, high activity of hedgehog signaling and the PI3K pathway in approximately 60% of 104 ESCC tumors indicates that therapies targeting these pathways might be particularly promising strategies for ESCC. Collectively, our data provide comprehensive insights into the mutational signatures of ESCC and identify markers for early diagnosis and potential therapeutic targets.

Multiomic analysis of mice epilepsy models suggest that miR-21a expression modulates mRNA and protein levels related to seizure deterioration.

Epilepsy is now recognized as the second most common neurological disease in China. To determine the genetic cause of epileptic encephalopathy, we performed a multiomics study using mouse models of controls, anticonvulsant mice treated with five drugs and epileptic mice. Based on genome-wide profiling analysis, we discovered four genes in the epileptic mouse group with differentially-expressed mRNA. After isobaric tags for relative and absolute quantification (iTRAQ) validation, only one gene, SNCA, remained, which was associated with apoptotic response of neuronal cells, and regulation of dopamine release and transport. We also identified three miRNAs targeting SNCA, out of which mmu-miR-21a-3p demonstrated a seven-fold change in expression between control and epileptic mice.

A rare novel mutation in TECTA causes autosomal dominant nonsyndromic hearing loss in a Mongolian family.

BACKGROUND: The genetic basis of autosomal dominant nonsyndromic hearing loss is complex. Genetic factors are responsible for approximately 50% of cases with congenital hearing loss. However, no previous studies have documented the clinical phenotype and genetic basis of autosomal dominant nonsyndromic hearing loss in Mongolians. METHODS: In this study, we performed exon capture sequencing of a Mongolian family with hereditary hearing loss and identified a novel mutation in TECTA gene, which encodes α -tectorin, a major component of the inner ear extracellular matrix that contacts the specialized sensory hair cells. RESULTS: The novel G → T missense mutation at nucleotide 6016 results in a substitution of amino acid aspartate at 2006 with tyrosine (Asp2006Tyr) in a highly conserved zona pellucida (ZP) domain of α-tectorin. The mutation is not found in control subjects from the same family with normal hearing and a genotype-phenotype correlation is observed. CONCLUSION: A novel missense mutation c.6016 G > T (p.Asp2006Tyr) of TECTA gene is a characteristic TECTA-related mutation which causes autosomal dominant nonsyndromic hearing loss. Our result indicated that mutation in TECTA gene is responsible for the hearing loss in this Mongolian family.

Chromosome-scale genome assembly of Astragalus membranaceus using PacBio and Hi-C technologies.

Astragalus membranaceus (Fisch.) Bge (AM) is a medicinal herb plant belonging to the Leguminosae family. In this study, we present a chromosome-scale genome assembly of AM, aiming to enhance the molecular biology and functional studies of Astragali Radix. The genome size of AM is about 1.43 Gb, with a contig N50 value of 1.67 Mb. A total of 98.16% of the assembly anchored to 9 pseudochromosomes using Hi-C technology. The assembly completeness was estimated to be 97.27% using BUSCO with the long terminal repeat assembly index (LAI) of 16.22 and quality value (QV) of 48.58. Additionally, the genome contained 67.98% repetitive sequences. Genome annotation predicted 29,914 protein-coding genes, including 73 genes involved in the flavonoid biosynthetic pathway and 2,048 transcription factors. The high-quality genome assembly and gene annotation resources will greatly facilitate future functional genomic studies in Leguminosae species.

Telomere-to-telomere genome assembly of sorghum.

"Cuohu Bazi" (CHBZ) is an ancient sorghum variety collected from the fields of China, known for its agronomic traits like dwarf stature, early maturation. In this study, we present the first telomere-to-telomere (T2T) and gap-free genome assembly of CHBZ using PacBio HiFi reads, Oxford Nanopore Technologies, and Hi-C data. The assembled genome comprises 724.85 Mb, effectively resolving all 3,913 gaps that were present in the previous sorghum BTx623 reference genome. Notably, the T2T assembly captures 10 centromeres and all 20 telomeres, providing strong support for their integrity. This assembly is of high quality in terms of contiguity (contig N50: 71.1 Mb), completeness (BUSCO score: 99.01%, k-mer completeness: 98.88%), and correctness (QV: 61.60). Repetitive sequences accounted for 70.41% of the genome and a total of 32,855 protein-coding genes have been annotated. Furthermore, 161 CHBZ-specific presence/absence variants genes have been identified when comparing to BTx623 genome. This study provides valuable insights for future research on sorghum genetics, genomics, and evolutionary history.

A complete reference genome of broomcorn millet.

Broomcorn millet (Panicum miliaceum L.), known for its traits of drought resistance, adaptability to poor soil, short growth period, and high photosynthetic efficiency as a C4 plant, represents one of the earliest domesticated crops globally. This study reports the telomere-to-telomere (T2T) gap-free reference genome for broomcorn millet (AJ8) using PacBio high-fidelity (HiFi) long reads, Oxford Nanopore long-read technologies and high-throughput chromosome conformation capture (Hi-C) sequencing data. The size of AJ8 genome was approximately 834.7 Mb, anchored onto 18 pseudo-chromosomes. Notably, 18 centromeres and 36 telomeres were obtained. The assembled genome showed high quality in terms of completeness (BUSCO score: 99.6%, QV: 61.7, LAI value: 20.4). In addition, 63,678 protein-coding genes and 433.8 Mb (~52.0%) repetitive sequences were identified. The complete reference genome for broomcorn millet provides a valuable resource for genetic studies and breeding of this important cereal crop.

A near complete genome assembly of the East Friesian sheep genome.

Advancements in sequencing have enabled the assembly of numerous sheep genomes, significantly advancing our understanding of the link between genetic variation and phenotypic traits. However, the genome of East Friesian sheep (Ostfriesisches Milchschaf), a key high-yield milk breed, remains to be fully assembled. Here, we constructed a near-complete and gap-free East Friesian genome assembly using PacBio HiFi, ultra-long ONT and Hi-C sequencing. The resulting genome assembly spans approximately 2.96 Gb, with a contig N50 length of 104.1 Mb and only 164 unplaced sequences. Remarkably, our assembly has captured 41 telomeres and 24 centromeres. The assembled sequence is of high quality on completeness (BUSCO score: 97.1%) and correctness (QV: 69.1). In addition, a total of 24,580 protein-coding genes were predicted, of which 97.2% (23,891) carried at least one conserved functional domain. Collectively, this assembly provides not only a near T2T gap-free genome, but also provides a valuable genetic resource for comparative genome studies of sheep and will serve as an important tool for the sheep research community.

Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca.

BACKGROUND: Despite their regional economic importance and being increasingly reared globally, the origins and evolution of the llama and alpaca remain poorly understood. Here we report reference genomes for the llama, and for the guanaco and vicuña (their putative wild progenitors), compare these with the published alpaca genome, and resequence seven individuals of all four species to better understand domestication and introgression between the llama and alpaca. RESULTS: Phylogenomic analysis confirms that the llama was domesticated from the guanaco and the alpaca from the vicuña. Introgression was much higher in the alpaca genome (36%) than the llama (5%) and could be dated close to the time of the Spanish conquest, approximately 500 years ago. Introgression patterns are at their most variable on the X-chromosome of the alpaca, featuring 53 genes known to have deleterious X-linked phenotypes in humans. Strong genome-wide introgression signatures include olfactory receptor complexes into both species, hypertension resistance into alpaca, and fleece/fiber traits into llama. Genomic signatures of domestication in the llama include male reproductive traits, while in alpaca feature fleece characteristics, olfaction-related and hypoxia adaptation traits. Expression analysis of the introgressed region that is syntenic to human HSA4q21, a gene cluster previously associated with hypertension in humans under hypoxic conditions, shows a previously undocumented role for PRDM8 downregulation as a potential transcriptional regulation mechanism, analogous to that previously reported at high altitude for hypoxia-inducible factor 1α. CONCLUSIONS: The unprecedented introgression signatures within both domestic camelid genomes may reflect post-conquest changes in agriculture and the breakdown of traditional management practices.

Comparative genomics and transcriptomics of Chrysolophus provide insights into the evolution of complex plumage coloration.

Background: As one of the most recognizable characteristics in birds, plumage color has a high impact on understanding the evolution and mechanisms of coloration. Feather and skin are ideal tissues to explore the genomics and complexity of color patterns in vertebrates. Two species of the genus Chrysolophus, golden pheasant (Chrysolophus pictus) and Lady Amherst's pheasant (Chrysolophus amherstiae), exhibit brilliant colors in their plumage, but with extreme phenotypic differences, making these two species great models to investigate plumage coloration mechanisms in birds. Results: We sequenced and assembled a genome of golden pheasant with high coverage and annotated 15,552 protein-coding genes. The genome of Lady Amherst's pheasant is sequenced with low coverage. Based on the feather pigment identification, a series of genomic and transcriptomic comparisons were conducted to investigate the complex features of plumage coloration. By identifying the lineage-specific sequence variations in Chrysolophus and golden pheasant against different backgrounds, we found that four melanogenesis biosynthesis genes and some lipid-related genes might be candidate genomic factors for the evolution of melanin and carotenoid pigmentation, respectively. In addition, a study among 47 birds showed some candidate genes related to carotenoid coloration in a broad range of birds. The transcriptome data further reveal important regulators of the two colorations, particularly one splicing transcript of the microphthalmia-associated transcription factor gene for pheomelanin synthesis. Conclusions: Analysis of the golden pheasant and its sister pheasant genomes, as well as comparison with other avian genomes, are helpful to reveal the underlying regulation of their plumage coloration. The present study provides important genomic information and insights for further studies of avian plumage evolution and diversity.

Whole-genome sequencing of 175 Mongolians uncovers population-specific genetic architecture and gene flow throughout North and East Asia.

The genetic variation in Northern Asian populations is currently undersampled. To address this, we generated a new genetic variation reference panel by whole-genome sequencing of 175 ethnic Mongolians, representing six tribes. The cataloged variation in the panel shows strong population stratification among these tribes, which correlates with the diverse demographic histories in the region. Incorporating our results with the 1000 Genomes Project panel identifies derived alleles shared between Finns and Mongolians/Siberians, suggesting that substantial gene flow between northern Eurasian populations has occurred in the past. Furthermore, we highlight that North, East, and Southeast Asian populations are more aligned with each other than these groups are with South Asian and Oceanian populations.

High-throughput deep screening and identification of four peripheral leucocyte microRNAs as novel potential combination biomarkers for preeclampsia.

OBJECTIVE: To identify the specific microRNA (miRNA) biomarkers of preeclampsia (PE), the miRNA profiles analysis were performed. STUDY DESIGN: The blood samples were obtained from five PE patients and five normal healthy pregnant women. The small RNA profiles were analyzed to identify miRNA expression levels and find out miRNAs that may associate with PE. The quantitative reverse transcriptase-PCR (qRT-PCR) assay was used to validate differentially expressed peripheral leucocyte miRNAs in a new cohort. RESULT: The data analysis showed that 10 peripheral leucocyte miRNAs were significantly differently expressed in severe PE patients. Four differently expressed miRNAs were successfully validated using qRT-PCR method. CONCLUSION: We successfully constructed a model with high accuracy to predict PE. A combination of four peripheral leucocyte miRNAs has great potential to serve as diagnostic biomarkers of PE.

MicroRNAome in decidua: a new approach to assess the maintenance of pregnancy.

OBJECTIVE: To comparatively analyze the human microRNAomes between normal pregnant and miscarriage deciduas by an in-depth sequencing of microRNA (miRNA); and to specifically examine miRNA-199b-5p and serum/glucocorticoid regulated kinase 1 (SGK1) in vivo and in vitro for their possible roles in pregnancy maintenance. DESIGN: Samples of deciduas from 6-8-week spontaneous miscarriages and normal pregnant women were irrespectively collected and comparatively analyzed by miRNA sequencing. The miR-199b-5p and SGK1 expressions were validated in vivo and in vitro. SETTING: University research and clinical institutes. PATIENT(S): In this experimental study, samples of deciduas were obtained from October 2011 to April 2012 from 29 women with spontaneous miscarriages and 35 normal pregnant women (control group) who underwent pregnancy termination at 6-8 weeks at our university gynecology unit. INTERVENTION(S): Endometrial biopsies, cell transfection, and production of an miR-199b-5p transgenic mouse model. MAIN OUTCOME MEASURE(S): In-depth sequencing of the miRNAome on human deciduas was performed for statistically significant differences in miRNA expression. Expression levels of SGK1 were detected by quantitative polymerase chain reaction and immunoblotting (Western blot) in vitro while miR-199b-5p is overexpressed or knockdown in miR-199b-5p transgenic mice. RESULT(S): Expression of the 1,921 known miRNAs was analyzed in the study. In aborted deciduas, 0.57% of the miRNAs were expressed abundantly (>10,000 transcripts per million) and represented 86.38% of all the miRNA reads. Six miRNAs were down-regulated (let-7a-5p, let-7f-5p, let-7g-5p, let-7e-5p, let-7d-5p, and miR-98), whereas miR-199b-5p was significantly up-regulated. Overexpression or knockdown of miR-199b-5p in HEK293T and Ishikawa cells decreased or increased SGK1 expression. Furthermore, overexpression of miR-199b-5p in human endometrial stromal cells or in transgenic mouse decreased SGK1 expression at the mRNA and protein levels, respectively. CONCLUSION(S): Among the miRNAomes, the abundant expression of the let-7 members was decreased in aborted samples, whereas miR-199b-5p expression was consistently increased. A significant inverse correlation was found between miR-199b-5p and SGK1 in vivo and in vitro.

Camelid genomes reveal evolution and adaptation to desert environments.

Bactrian camel (Camelus bactrianus), dromedary (Camelus dromedarius) and alpaca (Vicugna pacos) are economically important livestock. Although the Bactrian camel and dromedary are large, typically arid-desert-adapted mammals, alpacas are adapted to plateaus. Here we present high-quality genome sequences of these three species. Our analysis reveals the demographic history of these species since the Tortonian Stage of the Miocene and uncovers a striking correlation between large fluctuations in population size and geological time boundaries. Comparative genomic analysis reveals complex features related to desert adaptations, including fat and water metabolism, stress responses to heat, aridity, intense ultraviolet radiation and choking dust. Transcriptomic analysis of Bactrian camels further reveals unique osmoregulation, osmoprotection and compensatory mechanisms for water reservation underpinned by high blood glucose levels. We hypothesize that these physiological mechanisms represent kidney evolutionary adaptations to the desert environment. This study advances our understanding of camelid evolution and the adaptation of camels to arid-desert environments.

The genome of a Mongolian individual reveals the genetic imprints of Mongolians on modern human populations.

Mongolians have played a significant role in modern human evolution, especially after the rise of Genghis Khan (1162[?]-1227). Although the social cultural impacts of Genghis Khan and the Mongolian population have been well documented, explorations of their genome structure and genetic imprints on other human populations have been lacking. We here present the genome of a Mongolian male individual. The genome was de novo assembled using a total of 130.8-fold genomic data produced from massively parallel whole-genome sequencing. We identified high-confidence variation sets, including 3.7 million single nucleotide polymorphisms (SNPs) and 756,234 short insertions and deletions. Functional SNP analysis predicted that the individual has a pathogenic risk for carnitine deficiency. We located the patrilineal inheritance of the Mongolian genome to the lineage D3a through Y haplogroup analysis and inferred that the individual has a common patrilineal ancestor with Tibeto-Burman populations and is likely to be the progeny of the earliest settlers in East Asia. We finally investigated the genetic imprints of Mongolians on other human populations using different approaches. We found varying degrees of gene flows between Mongolians and populations living in Europe, South/Central Asia, and the Indian subcontinent. The analyses demonstrate that the genetic impacts of Mongolians likely resulted from the expansion of the Mongolian Empire in the 13th century. The genome will be of great help in further explorations of modern human evolution and genetic causes of diseases/traits specific to Mongolians.