Human nonsense-mediated RNA decay initiates widely by endonucleolysis and targets snoRNA host genes.

Eukaryotic RNAs with premature termination codons (PTCs) are eliminated by nonsense-mediated decay (NMD). While human nonsense RNA degradation can be initiated either by an endonucleolytic cleavage event near the PTC or through decapping, the individual contribution of these activities on endogenous substrates has remained unresolved. Here we used concurrent transcriptome-wide identification of NMD substrates and their 5'-3' decay intermediates to establish that SMG6-catalyzed endonucleolysis widely initiates the degradation of human nonsense RNAs, whereas decapping is used to a lesser extent. We also show that a large proportion of genes hosting snoRNAs in their introns produce considerable amounts of NMD-sensitive splice variants, indicating that these RNAs are merely by-products of a primary snoRNA production process. Additionally, transcripts from genes encoding multiple snoRNAs often yield alternative transcript isoforms that allow for differential expression of individual coencoded snoRNAs. Based on our findings, we hypothesize that snoRNA host genes need to be highly transcribed to accommodate high levels of snoRNA production and that the expression of individual snoRNAs and their cognate spliced RNA can be uncoupled via alternative splicing and NMD.

Consortium-based genome-wide meta-analysis for childhood dental caries traits.

Prior studies suggest dental caries traits in children and adolescents are partially heritable, but there has been no large-scale consortium genome-wide association study (GWAS) to date. We therefore performed GWAS for caries in participants aged 2.5-18.0 years from nine contributing centres. Phenotype definitions were created for the presence or absence of treated or untreated caries, stratified by primary and permanent dentition. All studies tested for association between caries and genotype dosage and the results were combined using fixed-effects meta-analysis. Analysis included up to 19 003 individuals (7530 affected) for primary teeth and 13 353 individuals (5875 affected) for permanent teeth. Evidence for association with caries status was observed at rs1594318-C for primary teeth [intronic within ALLC, odds ratio (OR) 0.85, effect allele frequency (EAF) 0.60, P 4.13e-8] and rs7738851-A (intronic within NEDD9, OR 1.28, EAF 0.85, P 1.63e-8) for permanent teeth. Consortium-wide estimated heritability of caries was low [h2 of 1% (95% CI: 0%: 7%) and 6% (95% CI 0%: 13%) for primary and permanent dentitions, respectively] compared with corresponding within-study estimates [h2 of 28% (95% CI: 9%: 48%) and 17% (95% CI: 2%: 31%)] or previously published estimates. This study was designed to identify common genetic variants with modest effects which are consistent across different populations. We found few single variants associated with caries status under these assumptions. Phenotypic heterogeneity between cohorts and limited statistical power will have contributed; these findings could also reflect complexity not captured by our study design, such as genetic effects which are conditional on environmental exposure.

A rare IL33 loss-of-function mutation reduces blood eosinophil counts and protects from asthma.

IL-33 is a tissue-derived cytokine that induces and amplifies eosinophilic inflammation and has emerged as a promising new drug target for asthma and allergic disease. Common variants at IL33 and IL1RL1, encoding the IL-33 receptor ST2, associate with eosinophil counts and asthma. Through whole-genome sequencing and imputation into the Icelandic population, we found a rare variant in IL33 (NM_001199640:exon7:c.487-1G>C (rs146597587-C), allele frequency = 0.65%) that disrupts a canonical splice acceptor site before the last coding exon. It is also found at low frequency in European populations. rs146597587-C associates with lower eosinophil counts (ß = -0.21 SD, P = 2.5×10-16, N = 103,104), and reduced risk of asthma in Europeans (OR = 0.47; 95%CI: 0.32, 0.70, P = 1.8×10-4, N cases = 6,465, N controls = 302,977). Heterozygotes have about 40% lower total IL33 mRNA expression than non-carriers and allele-specific analysis based on RNA sequencing and phased genotypes shows that only 20% of the total expression is from the mutated chromosome. In half of those transcripts the mutation causes retention of the last intron, predicted to result in a premature stop codon that leads to truncation of 66 amino acids. The truncated IL-33 has normal intracellular localization but neither binds IL-33R/ST2 nor activates ST2-expressing cells. Together these data demonstrate that rs146597587-C is a loss of function mutation and support the hypothesis that IL-33 haploinsufficiency protects against asthma.

A novel rare CUBN variant and three additional genes identified in Europeans with and without diabetes: results from an exome-wide association study of albuminuria.

AIMS/HYPOTHESIS: Identifying rare coding variants associated with albuminuria may open new avenues for preventing chronic kidney disease and end-stage renal disease, which are highly prevalent in individuals with diabetes. Efforts to identify genetic susceptibility variants for albuminuria have so far been limited, with the majority of studies focusing on common variants. METHODS: We performed an exome-wide association study to identify coding variants in a two-stage (discovery and replication) approach. Data from 33,985 individuals of European ancestry (15,872 with and 18,113 without diabetes) and 2605 Greenlanders were included. RESULTS: We identified a rare (minor allele frequency [MAF]: 0.8%) missense (A1690V) variant in CUBN (rs141640975, ß = 0.27, p = 1.3 × 10-11) associated with albuminuria as a continuous measure in the combined European meta-analysis. The presence of each rare allele of the variant was associated with a 6.4% increase in albuminuria. The rare CUBN variant had an effect that was three times stronger in individuals with type 2 diabetes compared with those without (pinteraction = 7.0 × 10-4, ß with diabetes = 0.69, ß without diabetes = 0.20) in the discovery meta-analysis. Gene-aggregate tests based on rare and common variants identified three additional genes associated with albuminuria (HES1, CDC73 and GRM5) after multiple testing correction (pBonferroni < 2.7 × 10-6). CONCLUSIONS/INTERPRETATION: The current study identifies a rare coding variant in the CUBN locus and other potential genes associated with albuminuria in individuals with and without diabetes. These genes have been implicated in renal and cardiovascular dysfunction. The findings provide new insights into the genetic architecture of albuminuria and highlight target genes and pathways for the prevention of diabetes-related kidney disease.

Fish Oil-Derived Fatty Acids in Pregnancy and Wheeze and Asthma in Offspring.

BACKGROUND: Reduced intake of n-3 long-chain polyunsaturated fatty acids (LCPUFAs) may be a contributing factor to the increasing prevalence of wheezing disorders. We assessed the effect of supplementation with n-3 LCPUFAs in pregnant women on the risk of persistent wheeze and asthma in their offspring. METHODS: We randomly assigned 736 pregnant women at 24 weeks of gestation to receive 2.4 g of n-3 LCPUFA (fish oil) or placebo (olive oil) per day. Their children formed the Copenhagen Prospective Studies on Asthma in Childhood2010 (COPSAC2010) cohort and were followed prospectively with extensive clinical phenotyping. Neither the investigators nor the participants were aware of group assignments during follow-up for the first 3 years of the children's lives, after which there was a 2-year follow-up period during which only the investigators were unaware of group assignments. The primary end point was persistent wheeze or asthma, and the secondary end points included lower respiratory tract infections, asthma exacerbations, eczema, and allergic sensitization. RESULTS: A total of 695 children were included in the trial, and 95.5% completed the 3-year, double-blind follow-up period. The risk of persistent wheeze or asthma in the treatment group was 16.9%, versus 23.7% in the control group (hazard ratio, 0.69; 95% confidence interval [CI], 0.49 to 0.97; P=0.035), corresponding to a relative reduction of 30.7%. Prespecified subgroup analyses suggested that the effect was strongest in the children of women whose blood levels of eicosapentaenoic acid and docosahexaenoic acid were in the lowest third of the trial population at randomization: 17.5% versus 34.1% (hazard ratio, 0.46; 95% CI, 0.25 to 0.83; P=0.011). Analyses of secondary end points showed that supplementation with n-3 LCPUFA was associated with a reduced risk of infections of the lower respiratory tract (31.7% vs. 39.1%; hazard ratio, 0.75; 95% CI, 0.58 to 0.98; P=0.033), but there was no statistically significant association between supplementation and asthma exacerbations, eczema, or allergic sensitization. CONCLUSIONS: Supplementation with n-3 LCPUFA in the third trimester of pregnancy reduced the absolute risk of persistent wheeze or asthma and infections of the lower respiratory tract in offspring by approximately 7 percentage points, or one third. (Funded by the Lundbeck Foundation and others; ClinicalTrials.gov number, NCT00798226 .).

Cadherin-related Family Member 3 Genetics and Rhinovirus C Respiratory Illnesses.

RATIONALE: Experimental evidence suggests that CDHR3 (cadherin-related family member 3) is a receptor for rhinovirus (RV)-C, and a missense variant in this gene (rs6967330) is associated with childhood asthma with severe exacerbations. OBJECTIVES: To determine whether rs6967330 influences RV-C infections and illnesses in early childhood. METHODS: We studied associations between rs6967330 and respiratory infections and illnesses in the COPSAC2010 (Copenhagen Prospective Studies on Asthma in Childhood 2010) and COAST (Childhood Origins of Asthma Birth Cohort Study) birth cohorts, where respiratory infections were monitored prospectively for the first 3 years of life. Nasal samples were collected during acute infections in both cohorts and during asymptomatic periods in COAST and analyzed for RV-A, RV-B, and RV-C, and other common respiratory viruses. MEASUREMENTS AND MAIN RESULTS: The CDHR3 asthma risk allele (rs6967330-A) was associated with increased risk of respiratory tract illnesses (incidence risk ratio [IRR] = 1.14 [95% confidence interval, 1.05-1.23]; P = 0.003). In particular, this variant was associated with risk of respiratory episodes with detection of RV-C in COPSAC2010 (IRR = 1.89 [1.14-3.05]; P = 0.01) and in COAST (IRR = 1.37 [1.02-1.82]; P = 0.03) children, and in a combined meta-analysis (IRR = 1.51 [1.13-2.02]; P = 0.006). In contrast, the variant was not associated with illnesses related to other viruses (IRR = 1.07 [0.92-1.25]; P = 0.37). Consistent with these observations, the CDHR3 variant was associated with increased detection of RV-C, but not of other viruses during scheduled visits at specific ages. CONCLUSIONS: The CDHR3 asthma risk allele is associated specifically with RV-C illnesses in two birth cohorts. This clinical evidence supports earlier molecular evidence indicating that CDHR3 functions as an RV-C receptor, and raises the possibility of preventing RV-C infections by targeting CDHR3.

UPF2-Dependent Nonsense-Mediated mRNA Decay Pathway Is Essential for Spermatogenesis by Selectively Eliminating Longer 3'UTR Transcripts.

During transcription, most eukaryotic genes generate multiple alternative cleavage and polyadenylation (APA) sites, leading to the production of transcript isoforms with variable lengths in the 3' untranslated region (3'UTR). In contrast to somatic cells, male germ cells, especially pachytene spermatocytes and round spermatids, express a distinct reservoir of mRNAs with shorter 3'UTRs that are essential for spermatogenesis and male fertility. However, the mechanisms underlying the enrichment of shorter 3'UTR transcripts in the developing male germ cells remain unknown. Here, we report that UPF2-mediated nonsense-mediated mRNA decay (NMD) plays an essential role in male germ cells by eliminating ubiquitous genes-derived, longer 3'UTR transcripts, and that this role is independent of its canonical role in degrading "premature termination codon" (PTC)-containing transcripts in somatic cell lineages. This report provides physiological evidence supporting a noncanonical role of the NMD pathway in achieving global 3'UTR shortening in the male germ cells during spermatogenesis.

A genome-wide association meta-analysis of diarrhoeal disease in young children identifies FUT2 locus and provides plausible biological pathways.

More than a million childhood diarrhoeal episodes occur worldwide each year, and in developed countries a considerable part of them are caused by viral infections. In this study, we aimed to search for genetic variants associated with diarrhoeal disease in young children by meta-analyzing genome-wide association studies, and to elucidate plausible biological mechanisms. The study was conducted in the context of the Early Genetics and Lifecourse Epidemiology (EAGLE) consortium. Data about diarrhoeal disease in two time windows (around 1 year of age and around 2 years of age) was obtained via parental questionnaires, doctor interviews or medical records. Standard quality control and statistical tests were applied to the 1000 Genomes imputed genotypic data. The meta-analysis (N = 5758) followed by replication (N = 3784) identified a genome-wide significant association between rs8111874 and diarrhoea at age 1 year. Conditional analysis suggested that the causal variant could be rs601338 (W154X) in the FUT2 gene. Children with the A allele, which results in a truncated FUT2 protein, had lower risk of diarrhoea. FUT2 participates in the production of histo-blood group antigens and has previously been implicated in the susceptibility to infections, including Rotavirus and Norovirus Gene-set enrichment analysis suggested pathways related to the histo-blood group antigen production, and the regulation of ion transport and blood pressure. Among others, the gastrointestinal tract, and the immune and neuro-secretory systems were detected as relevant organs. In summary, this genome-wide association meta-analysis suggests the implication of the FUT2 gene in diarrhoeal disease in young children from the general population.

NKG2D gene variation and susceptibility to viral bronchiolitis in childhood.

BACKGROUND: Genetic factors associated with bronchiolitis are inadequately characterized. We therefore inspected a selected subpopulation of our previous genome-wide association study (GWAS) of bronchiolitis for overlap with known quantitative trait loci (QTLs) to identify susceptibility loci that potentially affect mRNA and protein levels. METHODS: GWAS included a Finnish-Swedish case-control population (n = 187), matched for age and site. We integrated GWAS variants (p < 10-4) with QTL data. We subsequently verified allele-specific expression of identified QTLs by flow cytometry. Association of the resulting candidate loci with bronchiolitis was tested in three additional cohorts from Finland and Denmark (n = 1201). RESULTS: Bronchiolitis-susceptibility variant rs10772271 resided within QTLs previously associated with NKG2D (NK group 2, member D) mRNA and protein levels. Flow cytometric analysis confirmed the association with protein level in NK cells. The GWAS susceptibility allele (A) of rs10772271 (odds ratio [OR] = 2.34) corresponded with decreased NKG2D expression. The allele was nominally associated with bronchiolitis in one Finnish replicate (OR = 1.50), and the other showed directional consistency (OR = 1.43). No association was detected in Danish population CONCLUSIONS: The bronchiolitis GWAS susceptibility allele was linked to decreased NKG2D expression in the QTL data and in our expression analysis. We propose that reduced NKG2D expression predisposes infants to severe bronchiolitis.

Shared genetic variants suggest common pathways in allergy and autoimmune diseases.

BACKGROUND: The relationship between allergy and autoimmune disorders is complex and poorly understood. OBJECTIVE: We sought to investigate commonalities in genetic loci and pathways between allergy and autoimmune diseases to elucidate shared disease mechanisms. METHODS: We meta-analyzed 2 genome-wide association studies on self-reported allergy and sensitization comprising a total of 62,330 subjects. These results were used to calculate enrichment for single nucleotide polymorphisms (SNPs) previously associated with autoimmune diseases. Furthermore, we probed for enrichment within genetic pathways and of transcription factor binding sites and characterized commonalities in variant burden on tissue-specific regulatory sites by calculating the enrichment of allergy SNPs falling in gene regulatory regions in various cells using Encode Roadmap DNase-hypersensitive site data. Finally, we compared the allergy data with those of all known diseases. RESULTS: Among 290 loci previously associated with 16 autoimmune diseases, we found a significant enrichment of loci also associated with allergy (P = 1.4e-17) encompassing 29 loci at a false discovery rate of less than 0.05. Such enrichment seemed to be a general characteristic for autoimmune diseases. Among the common loci, 48% had the same direction of effect for allergy and autoimmune diseases. Additionally, we observed an enrichment of allergy SNPs falling within immune pathways and regions of chromatin accessible in immune cells that was also represented in patients with autoimmune diseases but not those with other diseases. CONCLUSION: We identified shared susceptibility loci and commonalities in pathways between allergy and autoimmune diseases, suggesting shared disease mechanisms. Further studies of these shared genetic mechanisms might help in understanding the complex relationship between these diseases, including the parallel increase in disease prevalence.

A novel common variant in DCST2 is associated with length in early life and height in adulthood.

Common genetic variants have been identified for adult height, but not much is known about the genetics of skeletal growth in early life. To identify common genetic variants that influence fetal skeletal growth, we meta-analyzed 22 genome-wide association studies (Stage 1; N = 28 459). We identified seven independent top single nucleotide polymorphisms (SNPs) (P < 1 × 10(-6)) for birth length, of which three were novel and four were in or near loci known to be associated with adult height (LCORL, PTCH1, GPR126 and HMGA2). The three novel SNPs were followed-up in nine replication studies (Stage 2; N = 11 995), with rs905938 in DC-STAMP domain containing 2 (DCST2) genome-wide significantly associated with birth length in a joint analysis (Stages 1 + 2; ß = 0.046, SE = 0.008, P = 2.46 × 10(-8), explained variance = 0.05%). Rs905938 was also associated with infant length (N = 28 228; P = 5.54 × 10(-4)) and adult height (N = 127 513; P = 1.45 × 10(-5)). DCST2 is a DC-STAMP-like protein family member and DC-STAMP is an osteoclast cell-fusion regulator. Polygenic scores based on 180 SNPs previously associated with human adult stature explained 0.13% of variance in birth length. The same SNPs explained 2.95% of the variance of infant length. Of the 180 known adult height loci, 11 were genome-wide significantly associated with infant length (SF3B4, LCORL, SPAG17, C6orf173, PTCH1, GDF5, ZNFX1, HHIP, ACAN, HLA locus and HMGA2). This study highlights that common variation in DCST2 influences variation in early growth and adult height.

C/EBPα is required for long-term self-renewal and lineage priming of hematopoietic stem cells and for the maintenance of epigenetic configurations in multipotent progenitors.

Transcription factors are key regulators of hematopoietic stem cells (HSCs) and act through their ability to bind DNA and impact on gene transcription. Their functions are interpreted in the complex landscape of chromatin, but current knowledge on how this is achieved is very limited. C/EBPα is an important transcriptional regulator of hematopoiesis, but its potential functions in HSCs have remained elusive. Here we report that C/EBPα serves to protect adult HSCs from apoptosis and to maintain their quiescent state. Consequently, deletion of Cebpa is associated with loss of self-renewal and HSC exhaustion. By combining gene expression analysis with genome-wide assessment of C/EBPα binding and epigenetic configurations, we show that C/EBPα acts to modulate the epigenetic states of genes belonging to molecular pathways important for HSC function. Moreover, our data suggest that C/EBPα acts as a priming factor at the HSC level where it actively promotes myeloid differentiation and counteracts lymphoid lineage choice. Taken together, our results show that C/EBPα is a key regulator of HSC biology, which influences the epigenetic landscape of HSCs in order to balance different cell fate options.

Temporal mapping of CEBPA and CEBPB binding during liver regeneration reveals dynamic occupancy and specific regulatory codes for homeostatic and cell cycle gene batteries.

Dynamic shifts in transcription factor binding are central to the regulation of biological processes by allowing rapid changes in gene transcription. However, very few genome-wide studies have examined how transcription factor occupancy is coordinated temporally in vivo in higher animals. Here, we quantified the genome-wide binding patterns of two key hepatocyte transcription factors, CEBPA and CEBPB (also known as C/EBPalpha and C/EBPbeta), at multiple time points during the highly dynamic process of liver regeneration elicited by partial hepatectomy in mouse. Combining these profiles with RNA polymerase II binding data, we find three temporal classes of transcription factor binding to be associated with distinct sets of regulated genes involved in the acute phase response, metabolic/homeostatic functions, or cell cycle progression. Moreover, we demonstrate a previously unrecognized early phase of homeostatic gene expression prior to S-phase entry. By analyzing the three classes of CEBP bound regions, we uncovered mutually exclusive sets of sequence motifs, suggesting temporal codes of CEBP recruitment by differential cobinding with other factors. These findings were validated by sequential ChIP experiments involving a panel of central transcription factors and/or by comparison to external ChIP-seq data. Our quantitative investigation not only provides in vivo evidence for the involvement of many new factors in liver regeneration but also points to similarities in the circuitries regulating self-renewal of differentiated cells. Taken together, our work emphasizes the power of global temporal analyses of transcription factor occupancy to elucidate mechanisms regulating dynamic biological processes in complex higher organisms.

UBL5 is essential for pre-mRNA splicing and sister chromatid cohesion in human cells.

UBL5 is an atypical ubiquitin-like protein, whose function in metazoans remains largely unexplored. We show that UBL5 is required for sister chromatid cohesion maintenance in human cells. UBL5 primarily associates with spliceosomal proteins, and UBL5 depletion decreases pre-mRNA splicing efficiency, leading to globally enhanced intron retention. Defective sister chromatid cohesion is a general consequence of dysfunctional pre-mRNA splicing, resulting from the selective downregulation of the cohesion protection factor Sororin. As the UBL5 yeast orthologue, Hub1, also promotes spliceosome functions, our results show that UBL5 plays an evolutionary conserved role in pre-mRNA splicing, the integrity of which is essential for the fidelity of chromosome segregation.

Amplification of pico-scale DNA mediated by bacterial carrier DNA for small-cell-number transcription factor ChIP-seq.

BACKGROUND: Chromatin-Immunoprecipitation coupled with deep sequencing (ChIP-seq) is used to map transcription factor occupancy and generate epigenetic profiles genome-wide. The requirement of nano-scale ChIP DNA for generation of sequencing libraries has impeded ChIP-seq on in vivo tissues of low cell numbers. RESULTS: We describe a robust, simple and scalable methodology for ChIP-seq of low-abundant cell populations, verified down to 10,000 cells. By employing non-mammalian genome mapping bacterial carrier DNA during amplification, we reliably amplify down to 50 pg of ChIP DNA from transcription factor (CEBPA) and histone mark (H3K4me3) ChIP. We further demonstrate that genomic profiles are highly resilient to changes in carrier DNA to ChIP DNA ratios. CONCLUSIONS: This represents a significant advance compared to existing technologies, which involve either complex steps of pre-selection for nucleosome-containing chromatin or pre-amplification of precipitated DNA, making them prone to introduce experimental biases.

spliceR: an R package for classification of alternative splicing and prediction of coding potential from RNA-seq data.

BACKGROUND: RNA-seq data is currently underutilized, in part because it is difficult to predict the functional impact of alternate transcription events. Recent software improvements in full-length transcript deconvolution prompted us to develop spliceR, an R package for classification of alternative splicing and prediction of coding potential. RESULTS: spliceR uses the full-length transcript output from RNA-seq assemblers to detect single or multiple exon skipping, alternative donor and acceptor sites, intron retention, alternative first or last exon usage, and mutually exclusive exon events. For each of these events spliceR also annotates the genomic coordinates of the differentially spliced elements, facilitating downstream sequence analysis. For each transcript isoform fraction values are calculated to identify transcript switching between conditions. Lastly, spliceR predicts the coding potential, as well as the potential nonsense mediated decay (NMD) sensitivity of each transcript. CONCLUSIONS: spliceR is an easy-to-use tool that extends the usability of RNA-seq and assembly technologies by allowing greater depth of annotation of RNA-seq data. spliceR is implemented as an R package and is freely available from the Bioconductor repository ( http://www.bioconductor.org/packages/2.13/bioc/html/spliceR.html).

The trans-ancestral genomic architecture of glycemic traits.

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Distinct immune phenotypes in infants developing asthma during childhood.

Early exposure to environmental triggers may elicit trajectories to chronic inflammatory disease through deregulated immune responses. To address relations between early immune competence and development of childhood asthma, we performed functional immune profiling of 186 parameters in blood of 541 18-month-old infants and examined links between their response phenotype and development of transient or persistent disease at 6 years of age. An abnormal neutrophil-linked antiviral response was associated with increased risk of transient asthma. Children who exhibited persistent asthma at year 6 showed enhanced interleukin-5 (IL-5) and IL-13 production in stimulated T cells at 18 months of age, which was associated with early life bacterial colonization of the airways. These findings highlight the early appearance of distinct immune characteristics in infants developing different asthma endotypes during childhood.